t;vcf/1000g.phase1.20110521.CFH.anno.vcf.gz&quot;, package = &quot;seqminer&quot;)
geneFile = system.file(&quot;vcf/refFlat_hg19_6col.txt.gz&quot;, package = &quot;seqminer&quot;)
geneName &lt;- &quot;CFH&quot;
cfh &lt;- readVCFToMatrixByGene(vcf.1000g, geneFile, geneName, &quot;Nonsynonymous&quot;)
</code></pre>

<pre><code>## Regcomp OK
</code></pre>

<pre><code class="r">dim(cfh)
</code></pre>

<pre><code>## [1] 47 85
</code></pre>

<pre><code class="r">
# Simulation phenotype under the null model
pheno 