<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<!--[if !mso]><style id=owaParaStyle>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:PMingLiU;
panose-1:2 2 3 0 0 0 0 0 0 0;}
@font-face
{font-family:PMingLiU;
panose-1:2 2 3 0 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:PMingLiU;
panose-1:2 2 3 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin-top:0in;
margin-right:0in;
margin-bottom:10.0pt;
margin-left:0in;
line-height:115%;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.emailstyle17
{mso-style-name:emailstyle17;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.25in 1.0in 1.25in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:12.0pt;font-family:"Arial","sans-serif";color:black">Dear Thibaut,<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:12.0pt;font-family:"Arial","sans-serif";color:black"><o:p> </o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:12.0pt;font-family:"Arial","sans-serif";color:black">Thank you for your help. I used the way you advised me to convert the Strucutre file of my haploid data set into a genind object as follow: <o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:12.0pt;font-family:"Arial","sans-serif";color:black"><o:p> </o:p></span></p>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:12.0pt;font-family:"Arial","sans-serif";color:black">dat <- read.table("structure.str", sep="\t",head=TRUE)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">x <- df2genind(dat[,-(1:2)], ploidy=1, pop=dat[,2], ind.names=dat[,1])</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">and then<br>
<br>
###<br>
find.cluster(x)<br>
dapc(x)</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">…</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black"><br>
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">However,
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">the membership plot derived from th</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">is Structure file</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">(showing admixture individuals)
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">was different
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">those obtained from using
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif"">import2genind
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif"">to convert
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">Genepop and Fstat files
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">of the same data set where data were coded as
</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">a pseudo-diploid</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black"> (no admixture individuals). Did I convert
my structure data set wrong? Or, I should trust results from </span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">Genepop and Fstat</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">
data sets? Thank you.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">Regards,
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif";color:black">Arhau</span><span style="font-size:12.0pt;line-height:115%;font-family:"Arial","sans-serif""><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Jombart, Thibaut [mailto:t.jombart@imperial.ac.uk]
<br>
<b>Sent:</b> Monday, July 04, 2011 9:30 AM<br>
<b>To:</b> Huang,Cheng-Hua; adegenet-forum@r-forge.wu-wien.ac.at<br>
<b>Cc:</b> arhuahuang@gmail.com<br>
<b>Subject:</b> RE: [adegenet-forum] Structure-like membership plot<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">Dear Arhua,
<br>
<br>
read.structure is only meant for diploid data - will add a warning on this point in the doc for 1.3-1.
<br>
<br>
So your analysis based on read.structure here is wrong. The other analyses are surely OK. Otherwise use read.table/read.csv to read your structure data into R as a data.frame, and then df2genind with the proper arguments to convert the data into a genind object.<br>
<br>
Best<br>
<br>
Thibaut<o:p></o:p></span></p>
<div>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:0in;margin-bottom:.0001pt;line-height:normal">
<span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">-- <br>
######################################<br>
Dr Thibaut JOMBART<br>
MRC Centre for Outbreak Analysis and Modelling<br>
Department of Infectious Disease Epidemiology<br>
Imperial College - Faculty of Medicine<br>
St Mary’s Campus<br>
Norfolk Place<br>
London W2 1PG<br>
United Kingdom<br>
Tel. : 0044 (0)20 7594 3658<br>
t.jombart@imperial.ac.uk<br>
http://sites.google.com/site/thibautjombart/<br>
http://adegenet.r-forge.r-project.org/<o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div class="MsoNormal" align="center" style="margin-bottom:0in;margin-bottom:.0001pt;text-align:center;line-height:normal">
<span style="font-size:12.0pt;font-family:"Times New Roman","serif";color:black">
<hr size="2" width="100%" align="center">
</span></div>
<div id="divRpF735054">
<p class="MsoNormal" style="margin-bottom:12.0pt;line-height:normal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black"> adegenet-forum-bounces@r-forge.wu-wien.ac.at
[adegenet-forum-bounces@r-forge.wu-wien.ac.at] on behalf of Cheng-Hua Huang [arhua@ufl.edu]<br>
<b>Sent:</b> 01 July 2011 17:10<br>
<b>To:</b> adegenet-forum@r-forge.wu-wien.ac.at<br>
<b>Cc:</b> arhuahuang@gmail.com<br>
<b>Subject:</b> [adegenet-forum] Structure-like membership plot</span><span style="font-size:12.0pt;font-family:"Times New Roman","serif";color:black"><o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;color:black">Dear adegent users and Thibaut,</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;color:black">I would like to ask questions about using DAPC. I have used
</span><span style="color:black">import2genind</span><span style="font-size:12.0pt;line-height:115%;color:black"> and
</span><span style="font-size:10.0pt;line-height:115%;font-family:"Tahoma","sans-serif";color:black">df2genind</span><span style="font-size:12.0pt;line-height:115%;color:black"> for importing Structure, Genepop, and Fstat files of the same data set to find
clusters and to plot Struture-like membership information. Since I was working on a haploid, asexual fungus collected from three regions, I coded data sets as a pseudo-diploid for Genepop and Fstat files but used the original haploid data set for the Structure
format. Even though I got K=3 for biological interpretation based on BIC, the membership plot derived from the Structure input file was different from those obtained from the other two file formats. The result using the Structure input format showed many admixture
individuals, but the other two input files did not exhibit any admixture individuals (all with membership probability = 1). The latter makes sense to me because the fungus is asexual and genotype flow plays a main role in shaping population structure. Moreover,
recombination is not possible for this fungus. Thus, admixture individuals may not occur. What caused this difference between these input files? How do I interpret these results? Thank you for your help.</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;color:black">Regards,</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;color:black">Arhua</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;color:black"> </span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:black"> <o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
</body>
</html>