<div dir="ltr">Hi again, <div><br></div><div>such plot typically indicates no clustering. Just to confirm: are we talking about 7 rows and 100,000 columns?</div><div><br></div><div>If so, your pools are technically your statistical individuals, and the method explore clustering solutions for 1-6 clusters for 7 individuals, which won't go far - not enough individuals to detect clustering really. Apologies if I misunderstood.</div><div><br></div><div>Best</div><div>Thibaut</div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><br>--<br>Dr Thibaut Jombart<br>Lecturer, Department of Infectious Disease Epidemiology, Imperial College London<br>Head of RECON: <a href="http://repidemicsconsortium.org" target="_blank">repidemicsconsortium.org</a><br>WHO Consultant - outbreak analysis</div><div><a href="https://thibautjombart.netlify.com" target="_blank">https://thibautjombart.netlify.com</a><br>Twitter: @TeebzR<br>+44(0)20 7594 3658</div></div></div></div>
<br><div class="gmail_quote">On 2 February 2018 at 09:07, Benjamin Dauphin <span dir="ltr"><<a href="mailto:benjamin.dauphin@wsl.ch" target="_blank">benjamin.dauphin@wsl.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Mark,<br>
<br>
Thanks for response. I’ve run find.clusters() with the matrix of allele frequencies as input file, and then run the DAPC using still the matrix (not the genind or genlight object) by assigning the group generated with kmeans (grp$grp). It works but I have a strange “inverted parabolic curve" for the kmean analysis.<br>
Is it a common picture for pooldseq data?<br>
<br>
Thanks,<br>
Ben<br>
<br>
<br><br>
<br>
> On 1 Feb 2018, at 18:01, Mark Coulson <<a href="mailto:Mark.Coulson.ic@uhi.ac.uk">Mark.Coulson.ic@uhi.ac.uk</a>> wrote:<br>
><br>
> Hi Ben,<br>
><br>
> I have used allelotype data with the input as a matrix of the frequency of the A allele in each group to run DAPC and it worked well. However, my groups were defined already but could the same type of input not be used to find.clusters?<br>
><br>
> Mark<br>
><br>
><br>
> -----Original Message-----<br>
> From: <a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org">adegenet-forum-bounces@lists.<wbr>r-forge.r-project.org</a> [mailto:<a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org">adegenet-forum-<wbr>bounces@lists.r-forge.r-<wbr>project.org</a>] On Behalf Of Benjamin Dauphin<br>
> Sent: 31 January 2018 09:18<br>
> To: <a href="mailto:adegenet-forum@lists.r-forge.r-project.org">adegenet-forum@lists.r-forge.<wbr>r-project.org</a><br>
> Subject: [adegenet-forum] Kmeans and DAPC on poolSeq data<br>
><br>
> Dear all,<br>
><br>
> I am newly working on pool sequencing data and I simply wonder if I can use kmeans (find.cluster) and DAPC to investigate population structure from poolseq data (allele frequencies)? How find.clusters can deal with allele frequencies?<br>
><br>
> Dataset: 7 pools and 100’000 SNPs<br>
><br>
> Any comment or help would be much appreciated.<br>
> Best regards<br>
> Ben<br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> adegenet-forum mailing list<br>
> <a href="mailto:adegenet-forum@lists.r-forge.r-project.org">adegenet-forum@lists.r-forge.<wbr>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-<wbr>project.org/cgi-bin/mailman/<wbr>listinfo/adegenet-forum</a><br>
> Inverness College UHI, a partner in the University of the Highlands and Islands <a href="http://www.inverness.uhi.ac.uk" rel="noreferrer" target="_blank">www.inverness.uhi.ac.uk</a> Board of Management of Inverness College (known as Inverness College UHI), Scottish Charity No SC021197.<br>
<br>
<br>______________________________<wbr>_________________<br>
adegenet-forum mailing list<br>
<a href="mailto:adegenet-forum@lists.r-forge.r-project.org">adegenet-forum@lists.r-forge.<wbr>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-<wbr>project.org/cgi-bin/mailman/<wbr>listinfo/adegenet-forum</a><br></blockquote></div><br></div>