<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><span style="font-family: 'Segoe UI', Helvetica, Arial, sans-serif; ">Dear,</span></div>
<div><span style="font-family: 'Segoe UI', Helvetica, Arial, sans-serif; "><br>
</span></div>
<div><span style="font-family: 'Segoe UI', Helvetica, Arial, sans-serif; ">I'm using Adegenet to do some population genomic analyses on a dataset of around 30.000 SNPs. I've read the data into a genind object from a genepop file. For get a first idea of the
data, I used the function summary and compared Hobs with Hexp</span></div>
<div><span style="font-family: 'Segoe UI', Helvetica, Arial, sans-serif; "><br>
</span></div>
<div>
<p class="p1"><span class="s1">d.genind </span><span class="s2"><-</span><span class="s1">
</span><span class="s2">read.genepop(</span>"filt.popgen.gen"<span class="s2">)</span></p>
</div>
<div>
<p class="p1">summary <span class="s1"><-</span> <span class="s1">summary(</span>d.genind<span class="s1">)</span></p>
<p class="p1"><span class="s1"></p>
<p class="p1"><span class="s1">barplot(</span><span class="s2">summary</span><span class="s1">$</span><span class="s2">Hexp </span><span class="s1">- </span><span class="s2">summary</span><span class="s1">$</span><span class="s2">Hobs</span><span class="s1">)</span></p>
<p class="p1"><span class="s1"><br>
</span></p>
<p class="p1"><span class="s1">Here I see that some loci have much higher Hobs than Hexp; so I plotted a qqplot to get a better picture:</span></p>
<p class="p1"><span class="s1"><br>
</span></p>
<p class="p1"><span class="s1">qqplot(</span>summary<span class="s1">$</span>Hobs<span class="s1">,</span>summary<span class="s1">$</span>Hexp)</p>
<p class="p1"><br>
</p>
<p class="p1">Here I see that I have loci for which Hobs is larger than 0.5, which is not possible. In the next step, I would therefore like to remove loci that have a Hobs larger than 0.5, but I don't know how. The most easy would be to state a logical, something
that states 'remove all loci that have Hobs larger then 0.5'. But I don't know how to do this. Is this possible for genind objects? And if so, how?</p>
<div><br>
</div>
<p class="p1">I know that keeping loci into the genind object can be done by (if you have a list of loci names):</p>
<p class="p1"><br>
</p>
<p class="p1">d.genind[loc=c("x1","x2",etc.)]</p>
<p class="p1"><br>
</p>
<p class="p1">But how can one easily remove loci from a genind object? Putting a minus '-' before c() returns an error. </p>
<p class="p1"><br>
</p>
<p class="p1">Finally I had a question about DAPC. It is possible to state the number of axis for PCA (n.pca). However, adding more and more axis increases the power to discriminate the populations. How can one know how many axis to add? If you add not enough
axis, than you will see no structure. If you add to much, then it looks like you have been overseperating the populations. </p>
<p class="p1"><br>
</p>
<p class="p1">Thank you very much for your answers.</p>
<p class="p1"><br>
</p>
<p class="p1">Sincerely,</p>
<p class="p1">Frederik</p>
</span>
<p></p>
</div>
</div>
</body>
</html>