<div dir="ltr">Hi Carol, <div><br></div><div>Roman's solution should have worked, but adding random noise to locations may not be ideal anyway. You may prefer using a connection network which can accommodate duplicated locations, such as neighbourhood defined by distance interval (see ?chooseCN).</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"><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<br><a href="http://sites.google.com/site/thibautjombart/" target="_blank">sites.google.com/site/thibautjombart/</a><br>Twitter: @TeebzR<br>+44(0)20 7594 3658</div></div>
<br><div class="gmail_quote">On 3 August 2017 at 23:53, Roman Luštrik <span dir="ltr"><<a href="mailto:roman.lustrik@biolitika.si" target="_blank">roman.lustrik@biolitika.si</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-family:trebuchet ms,sans-serif;font-size:12pt;color:#000000"><div>As the error is reporting, you should use jitter() on your coordinates. Either specify new jittered coordinates directly to xy argument of spca(), or, jitter your genind object which you can access using `other(x)$xy`.</div><div><br></div><div>Here's an example of the former approach:</div><div><br></div><div><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> library(adegenet)</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> data(spcaIllus)</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> x <- spcaIllus[[1]]</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> out <- spca(x, ask = FALSE, type = 1, scannf = FALSE)</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> xy <- other(x)$xy</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> newxy <- xy[rep(1:(nrow(xy)/10), length.out = nrow(xy)), ]</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> out.clumped <- spca(x, ask = FALSE, type = 1, scannf = FALSE, xy = newxy)</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">Error in chooseCN(xy = xy, ask = ask, type = type, plot.nb = plot.nb, : </span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif"> Duplicate locations detected and incompatible with graph type 1-4.</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">Please choose another graph (5-7) or add random noise to locations (see ?jitter).</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> out.jittered <- spca(x, ask = FALSE, type = 1, scannf = FALSE, xy = jitter(newxy))</span><br><span style="font-family:"courier new",courier,monaco,monospace,sans-serif">> </span></div><div><br></div><div>Cheers,</div><div>Roman</div><div><br></div><div>----<br>In god we trust, all others bring data.<br>> Zahtevaj IJZ na <a href="https://kurc.biolitika.si" target="_blank">https://kurc.biolitika.si</a></div><br><hr id="m_421054250945395128zwchr"><div><b>From: </b>"Carol Ann Reeb" <<a href="mailto:creeb@stanford.edu" target="_blank">creeb@stanford.edu</a>><br><b>To: </b><a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">adegenet-forum@lists.r-forge.<wbr>r-project.org</a><br><b>Sent: </b>Friday, August 4, 2017 12:14:01 AM<br><b>Subject: </b>[adegenet-forum] Jitter & sPCA<br></div><br><div><span class="">


<div id="m_421054250945395128divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>I'm having trouble  jittering x-coordinates in order to run sPCA.   I'm still learning R so I don't know what code would work in adegent.  I did find R usage:  Jitter(x, factor=1, amount=NULL).  But that does not tell me how to call the proper files.  I
 did find code from a supplement from a manuscript comparing MENE to sPCA comparison.  And although that routine ran, the program would still not allow me to choose a network to plot sPCA coordinates because it said they were still identical in the xy file.
  Any help, suggestions would be greatly appreciated.</p>
<p><br>
</p>
<p>Carol</p>
</div>


<br></span>______________________________<wbr>_________________<br>adegenet-forum mailing list<br><a href="mailto:adegenet-forum@lists.r-forge.r-project.org" target="_blank">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" target="_blank">https://lists.r-forge.r-<wbr>project.org/cgi-bin/mailman/<wbr>listinfo/adegenet-forum</a><br></div></div></div><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>