
<HTML><HEAD>

<META content="text/html; charset=utf-8" http-equiv=Content-Type>

<META name=GENERATOR content="MSHTML 8.00.6001.23515"></HEAD>

<BODY style="MARGIN: 4px 4px 1px; FONT: 10pt Tahoma">

<DIV>

<P style="LINE-HEIGHT: normal; MARGIN: 0cm 0cm 0pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Tahoma','sans-serif'; FONT-SIZE: 10pt; mso-fareast-font-family: 'Times New Roman'; mso-ansi-language: EN-US; mso-fareast-language: DE" lang=EN-US>Dear Valeria,<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" /><o:p></o:p></SPAN></P>

<P style="LINE-HEIGHT: normal; MARGIN: 0cm 0cm 0pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Tahoma','sans-serif'; FONT-SIZE: 10pt; mso-fareast-font-family: 'Times New Roman'; mso-ansi-language: EN-US; mso-fareast-language: DE" lang=EN-US> <o:p></o:p></SPAN></P>

<P style="LINE-HEIGHT: normal; MARGIN: 0cm 0cm 0.75pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Tahoma','sans-serif'; FONT-SIZE: 10pt; mso-fareast-font-family: 'Times New Roman'; mso-ansi-language: EN-US; mso-fareast-language: DE" lang=EN-US>thank you very much for your quick answer. I’m aware of the problems STUCTURE has to analyze genetic data of continuous populations (see also <A href="http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2664.2008.01606.x/pdf">http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2664.2008.01606.x/pdf</A>). That is one reason I don’t want to use STUCTURE as the only cluster analysis. <SPAN style="mso-spacerun: yes"> </SPAN>I haven’t attempted to use BAPS yet, but I gave GENELAND a trial to include spatial information. Besides testing for IBD with a Mantel test, I also modified the geographic distances by resistance values etc. I inferred from a SDM. A spatial autocorrelations didn’t show a clear pattern of spatial relation (also in different distance classes). <SPAN style="mso-spacerun: yes"> </SPAN>A PCA indicates a big cloud around the center point. Each of the first two axes explained about 19 % of the variance.<o:p></o:p></SPAN></P>

<P style="LINE-HEIGHT: normal; MARGIN: 0cm 0cm 0.75pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Tahoma','sans-serif'; FONT-SIZE: 10pt; mso-fareast-font-family: 'Times New Roman'; mso-ansi-language: EN-US; mso-fareast-language: DE" lang=EN-US>Thanks to assure the correctness of my DAPC script. I set the maximum number of clusters to 50 to exclude a missing of structural shifts.<o:p></o:p></SPAN></P>

<P style="LINE-HEIGHT: normal; MARGIN: 0cm 0cm 0.75pt" class=MsoNormal><SPAN style="FONT-FAMILY: 'Tahoma','sans-serif'; FONT-SIZE: 10pt; mso-fareast-font-family: 'Times New Roman'; mso-ansi-language: EN-US; mso-fareast-language: DE" lang=EN-US>Nonetheless, I cannot explain the contrary results of structure indicating a panmictic population (4 parallel stripes) and DAPC assigning most individuals to one specific cluster. <o:p></o:p></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="mso-ansi-language: EN-US" lang=EN-US>Thanks again for your comments. I will have a look at BAPS.<o:p></o:p></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="mso-ansi-language: EN-US" lang=EN-US>Best wishes, <o:p></o:p></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN style="mso-ansi-language: EN-US" lang=EN-US>Jutta<o:p></o:p></SPAN></P></DIV>

<DIV>>>> Valeria Montano <mirainoshojo@gmail.com> 9/5/2013 10:59 >>><BR></DIV>

<DIV dir=ltr>

<DIV>

<DIV>Dear Jutta,<BR><BR></DIV>cluster analysis can be tricky when the samples analysed are distributed along a gradient and if there is no clear-cut subdivision, this can lead to contradictory results (have a look at this paper <A href="http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.192.3029&rep=rep1&type=pdf">http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.192.3029&rep=rep1&type=pdf</A>). You may want to consider using TESS or BAPS with the admixture model option. These two software allow including the geographic coordinates as a prior information and the admixture model is a way to model spatial gradients. If you tested the IBD with a Mantel test, just be careful that a significant mantel test is not directly due to IBD, geo to gen correlation can be significant for different spatial/migratory schemes. I think your DAPC is ok, a part from the fact that there is no need to use the find.clusters with the number of PCs indicated by the optim.a.score. This procedure is used to optimize the discriminant space among clusters in the DAPC. To assign individuals to clusters you can simply retrieve all the variance (even though in your case is almost the same given that you have 98%). Only thing, I would try with max number of clusters around 20, more than your sampling locations. You can also give sPCA a try.<BR><BR></DIV>

<DIV>Hope this helps<BR><BR></DIV>

<DIV>Ciao<BR><BR></DIV>

<DIV>Valeria<BR></DIV></DIV>

<DIV class=gmail_extra><BR><BR>

<DIV class=gmail_quote>On 4 September 2013 15:03, Jutta Geismar <SPAN dir=ltr><<A href="mailto:Jutta.Geismar@senckenberg.de" target=_blank>Jutta.Geismar@senckenberg.de</A>></SPAN> wrote:<BR>

<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class=gmail_quote>

<DIV style="MARGIN: 4px 4px 1px; FONT: 10pt Tahoma">

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>Dear Mr Jombart and DAPC users,<U></U><U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><U></U><FONT size=3 face=Calibri></FONT><U></U></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>I used DAPC to analyze genetic structure in a small region with 20 microsatellite markers. I analyzed 330 individuals (14 sampling sites) and found little genetic differences (FST, D Jost), but a significant isolation by distance pattern. A cluster analysis in STRUCTURE resulted in four clusters (STRUCTURE Harvester) but all individuals had more or less equal posterior probability in all of the four inferred clusters. Therefore I assume a panmictic population structure. Since STRUCTURE is known for some problems analyzing datasets under IBD I analyzed the data with DAPC. DAPC resulted in 3 or 4 clusters (and tested up until K=7 to be sure), but in both cases these were randomly distributed among all individuals without a geographic context. Only 94 individuals were not assigned to one cluster with more than 90% and therefore would be counted as “admixed” (example in DAPC tutorial). For me the results of STRUCTURE and DAPC are in conflict to each other, but I don’t know how a panmictic population would look like in DAPC. Distances between sites are small and it is very likely that gene flow occurs among my sampling points, which might cause problems in genetic cluster analyses. I don’t know if I made any mistake in my thinking, that’s why I want to explain my procedure briefly:<U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 0pt 36pt"><SPAN lang=EN-US><SPAN><FONT size=3 face=Calibri>1.</FONT><SPAN style="FONT: 7pt 'Times New Roman'"> </SPAN></SPAN></SPAN><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>I used dapc and chose 1/3 of the sample size as PC (as suggested) and counted DAs in the plot (100% of the variability was included, 110 PC, 13 DA)<U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 0pt 36pt"><SPAN lang=EN-US><SPAN><FONT size=3 face=Calibri>2.</FONT><SPAN style="FONT: 7pt 'Times New Roman'"> </SPAN></SPAN></SPAN><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>To reduce variability I used optim.a.score (smart FALSE). The best a-score was around 0.2 (PC 61)<U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 0pt 36pt"><SPAN lang=EN-US><SPAN><FONT size=3 face=Calibri>3.</FONT><SPAN style="FONT: 7pt 'Times New Roman'"> </SPAN></SPAN></SPAN><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>After that I wanted to estimate the number of clusters by find.clusters and used the a-score as number of PCs and repeated the dapc (conserved variance was still 98%, 61 PCs, 2 DA) <U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt 36pt"><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>I chose k in the BIC values after which the decrease was less compared to the previous, but not the lowest k.<U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>If I have some mistakes in my procedure I would appreciate some advice. But also if the procedure is okay I cannot explain the contrariness of these two analyses. <U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>Thanks a lot in advance for some help.<U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>Jutta Geismar <U></U><U></U></FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>PhD student</FONT></FONT></SPAN></P>

<P style="MARGIN: 0cm 0cm 10pt" class=MsoNormal><SPAN lang=EN-US><FONT size=3><FONT face=Calibri>Germany<U></U><U></U></FONT></FONT></SPAN></P></DIV><BR>_______________________________________________<BR>adegenet-forum mailing list<BR><A href="mailto:adegenet-forum@lists.r-forge.r-project.org">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" target=_blank>https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/adegenet-forum</A><BR></BLOCKQUOTE></DIV><BR></DIV></BODY></HTML>