<div dir="ltr"><div>Hi Rachel,</div><div><br></div><div>I believe that the statement for why DAPC can not be scripted is because the steps to determine the number of principal components and discriminant functions are interactive. If you know this information you can send the DAPC step to a server without a monitor. Not sure DAPC is going to help you find SNPs associated with phenotypes though.</div><div><br></div><div>Good luck!</div><div>Brian<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 6, 2019 at 11:14 AM Rachel Turba <<a href="mailto:rturba@outlook.com">rturba@outlook.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div class="gmail-m_-7965030034425534424WordSection1">
<p class="MsoNormal">Hi Brian,</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thank you for the help. I was avoiding to subset my data because I am interested in investigating SNPs that could be related to a change in morph between these populations.</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Cheers,</p>
<p class="MsoNormal">Rachel</p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-7965030034425534424divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> brian knaus <<a href="mailto:briank.lists@gmail.com" target="_blank">briank.lists@gmail.com</a>><br>
<b>Sent:</b> Tuesday, March 5, 2019 2:05:03 PM<br>
<b>To:</b> Rachel Turba<br>
<b>Cc:</b> <a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.r-project.org</a><br>
<b>Subject:</b> Re: [adegenet-forum] Scripting adegenet with whole-genome data</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>Hi Rachel,</div>
<div><br>
</div>
<div>Demographic patterns are generally thought to occur throughout the genome. In which case you shouldn't need the entire genome but can use a representative subset of your variants. We've demonstrated how to accomplish this at the below link.</div>
<div><br>
</div>
<div><a href="https://grunwaldlab.github.io/Population_Genetics_in_R/gbs_analysis.html#subsetting-a-vcfr-object-to-200-random-variants" target="_blank">https://grunwaldlab.github.io/Population_Genetics_in_R/gbs_analysis.html#subsetting-a-vcfr-object-to-200-random-variants</a><br>
</div>
<div><br>
</div>
<div>Good luck!</div>
<div>Brian<br>
</div>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Tue, Mar 5, 2019 at 11:24 AM Rachel Turba <<a href="mailto:rturba@outlook.com" target="_blank">rturba@outlook.com</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="EN-US">
<div class="gmail-m_-7965030034425534424gmail-m_-4436784124459730377WordSection1">
<p class="MsoNormal">Hi all,</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have researched this on the forum but was not successful in finding help, so I apologize in advance if this has been addressed already.</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">In this tutorial (<a href="https://grunwaldlab.github.io/Population_Genetics_in_R/DAPC.html" target="_blank">https://grunwaldlab.github.io/Population_Genetics_in_R/DAPC.html</a>) it says that adegenet cannot be scripted a priori, but my
data has 6Gb and takes a very long time to generate the DAPC object. So I have to submit it as a script to the university cluster. I get the plot of variance explained by PC but then everything halts because I do not provide the number of PCs to be retained
for continuing the analysis. Is there a way around this? How does anyone handle large datasets?</p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thank you in advance,</p>
<p class="MsoNormal">Rachel</p>
</div>
</div>
_______________________________________________<br>
adegenet-forum mailing list<br>
<a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum" rel="noreferrer" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum</a></blockquote>
</div>
</div>
</div>
_______________________________________________<br>
adegenet-forum mailing list<br>
<a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum" rel="noreferrer" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum</a></blockquote></div>