<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)">
<!--[if !mso]><style>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: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:"Lucida Console";
        panose-1:2 11 6 9 4 5 4 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        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;}
pre
        {mso-style-priority:99;
        mso-style-link:"HTML Preformatted Char";
        margin:0in;
        margin-bottom:.0001pt;
        font-size:10.0pt;
        font-family:"Courier New";}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
        {mso-style-priority:99;
        mso-style-link:"Balloon Text Char";
        margin:0in;
        margin-bottom:.0001pt;
        font-size:8.0pt;
        font-family:"Tahoma","sans-serif";}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
span.HTMLPreformattedChar
        {mso-style-name:"HTML Preformatted Char";
        mso-style-priority:99;
        mso-style-link:"HTML Preformatted";
        font-family:"Courier New";}
span.gcg2ujhdeab
        {mso-style-name:gcg2ujhdeab;}
span.BalloonTextChar
        {mso-style-name:"Balloon Text Char";
        mso-style-priority:99;
        mso-style-link:"Balloon Text";
        font-family:"Tahoma","sans-serif";}
span.gcg2ujhdiab
        {mso-style-name:gcg2ujhdiab;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
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">Hello all,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I’m attempting to recreate a similar figure to the one seen in Sarno et al. (<a href="http://www.plosone.org/article/fetchObject.action?uri=info:doi/10.1371/journal.pone.0096074&representation=PDF">http://www.plosone.org/article/fetchObject.action?uri=info:doi/10.1371/journal.pone.0096074&representation=PDF</a>).
 To do so, I create the attached matrix (HGF.txt) that contains populations (rows) and Y-STR haplogroups (columns). Additionally, due to my inexperience with adegenet, I created the attached coordinates list in a separate file (XY.txt). I then did the following.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Function:<o:p></o:p></p>
<pre style="word-break:break-all">mydata <- <span style="font-family:"Lucida Console";color:blue;background:#E1E2E5">genpop(HGF,ploidy=1)</span><span style="font-family:"Lucida Console";color:black;background:#E1E2E5"><o:p></o:p></span></pre>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Result:<o:p></o:p></p>
<p class="MsoNormal">Object created<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Function:<o:p></o:p></p>
<pre style="word-break:break-all"><span class="gcg2ujhdeab"><span style="font-family:"Lucida Console";color:blue;background:#E1E2E5">spca(mydata,xy=xy,d1=0,d2=5)</span></span><span style="font-family:"Lucida Console";color:black;background:#E1E2E5"><o:p></o:p></span></pre>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Result:<o:p></o:p></p>
<p class="MsoNormal">Plot created<o:p></o:p></p>
<p class="MsoNormal"><img border="0" width="306" height="130" id="Picture_x0020_1" src="cid:image003.jpg@01D0B4E0.1EA1AB50"><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Now comes the problem. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Function:<o:p></o:p></p>
<pre style="word-break:break-all"><span class="gcg2ujhdeab"><span style="font-family:"Lucida Console";color:blue;background:#E1E2E5">myspca<-spca(mydata,xy=xy,type=5,d1=0,d2=5)</span></span><span style="font-family:"Lucida Console";color:black;background:#E1E2E5"><o:p></o:p></span></pre>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Result:<o:p></o:p></p>
<p class="MsoNormal">No “myspca” list is created and I get the following error<o:p></o:p></p>
<pre style="word-break:break-all"><span class="gcg2ujhdiab"><span style="font-family:"Lucida Console";color:#C5060B;background:#E1E2E5">Error in if (nf > rank) nf <- rank : <o:p></o:p></span></span></pre>
<pre style="word-break:break-all"><span class="gcg2ujhdiab"><span style="font-family:"Lucida Console";color:#C5060B;background:#E1E2E5">  missing value where TRUE/FALSE needed</span></span><span style="font-family:"Lucida Console";color:black;background:#E1E2E5"><o:p></o:p></span></pre>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I feel that if I can obtain the eig values for the dataset the rest will be smooth sailing; however, I am unable to find any similar instances on the forum (some people have mentioned this general error in regard to clustering but not spca).
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Any help is greatly appreciated. Thanks in advance <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Jonathan King</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Research and Development Lab Manager</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Institute of Applied Genetics</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">UNT Health Science Center</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">3500 Camp Bowie Blvd</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Fort Worth, TX 76107</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Office: 817-735-2773<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Phone: 817-735-2940</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><a href="https://www.unthsc.edu/graduate-school-of-biomedical-sciences/molecular-and-medical-genetics/laboratory-faculty-and-staff/">Research and Development Lab</a><u><span style="color:blue"><o:p></o:p></span></u></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>