<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 14 (filtered medium)"><style><!--
/* Font Definitions */
@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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
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;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Tahoma","sans-serif";
        color:windowtext;
        font-weight:normal;
        font-style:normal;
        text-decoration:none none;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
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-GB link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>Hi,<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>I’m very new to adegenet, and trying to determine if it is sensible to create a distance matrix from a neighbour joining tree (njt) for correspondence analysis with a genetic matrix.  I couldn’t find any information on this from a search of the archives.  <o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>I have created a neighbour joining tree and then used it in an sPCA as follows:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>njt <- chooseCN(myind@other$xy,ask=FALSE,type=4) #create njt from xy<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>njt_ed<-edit.nb(njt,myind@other$xy,polys=rb_polys) #edit njt to insert “barriers”, creates a “forest”<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>myspca<-spca(myind, cn=njt_ed, scale=TRUE, type=1, plot.nb=TRUE, nfposi=40, nfnega=40, ask=FALSE, scannf=FALSE)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>I can now create a genetic distance matrix from the microsatellite data:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>dg<-dist(myind$tab)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>However, I don’t know how to create a distance matrix for the individuals based on the edited neighbour joining tree, rather than the xy coordinates (as in the tutorial) or the results of the sPCA which also represent genetic variation.  I want to compare distances for individuals based on the njt “forest”, so that I can compare them, e.g.:<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>plot(dg,?? geographic distance matrix using njt ??)<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>Perhaps it would be more sensible to use a different type of connection network?<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>If anyone can help I’d be very grateful, thanks.<o:p></o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'>Peter<o:p></o:p></span></p></div></body></html>