&lt;- rnorm(ncol(cfh))

# Apply SKAT method
library(SKAT)
obj &lt;- SKAT_Null_Model(pheno ~ 1, out_type = &quot;C&quot;)
SKAT(t(cfh), obj)$p.value
</code></pre>

<pre><code>## Warning: 37 SNPs with either high missing rates or no-variation are
## excluded!
</code></pre>

<pre><code>## Warning: Genotypes of some variants are not the number of minor alleles!
</code></pre>

<pre><code>## [1] 0.2591
</code></pre>

<h2 id="toc_4">Contact</h2>

<p>SeqMiner is developed with joint effort from <a href="mailto:zhanx