<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Coraline
<br>
<br>
The colorplot is basically translating up to 3 quantitative variables into colors using RGB coding. So each axis corresponds to a channel, but the problem is our eyes are not good at knowing what color comes out when combining values RBG channels. To have a
proper legend for 3 axes we would need some kind of color triangle. And a different system for 2 axes.. nothing impossible to code, but cumbersome.
<br>
<br>
As for your second question, your approach is the best way to go. But there's a trick.<br>
The function add.scatter.eig, which adds the barplots of eigenvalues to the plot, changes the coordinate system, so that you can't plot anything afterwards. You need to disable them,e g.:<br>
<br>
library(adegenet)<br>
example(dapc)<br>
<div>scatter(dapc1, col="transparent", scree.da=FALSE)<br>
points(dapc1$ind.coord[,1],dapc1$ind.coord[,2], col=funky(100))<br>
<div style="font-family:Tahoma; font-size:13px">
<div class="BodyFragment"><font size="2"><span style="font-size:10pt">
<div class="PlainText"><br>
Note that if some points are missing around the edges, you may need a par(xpd=TRUE) right before your 'points(..)'.<br>
<br>
Cheers<br>
Thibaut<br>
<br>
==============================<br>
Dr Thibaut Jombart<br>
MRC Centre for Outbreak Analysis and Modelling<br>
Department of Infectious Disease Epidemiology<br>
Imperial College - School of Public Health<br>
Norfolk Place, London W2 1PG, UK<br>
Tel. : 0044 (0)20 7594 3658<br>
http://sites.google.com/site/thibautjombart/<br>
http://sites.google.com/site/therepiproject/<br>
http://adegenet.r-forge.r-project.org/<br>
Twitter: @thibautjombart<br>
<br>
<br>
</div>
</span></font></div>
</div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF839377"><font face="Tahoma" color="#000000" size="2"><b>From:</b> adegenet-forum-bounces@lists.r-forge.r-project.org [adegenet-forum-bounces@lists.r-forge.r-project.org] on behalf of BICHET CORALINE [coraline.bichet@univ-lyon1.fr]<br>
<b>Sent:</b> 20 August 2015 10:29<br>
<b>To:</b> adegenet-forum@lists.r-forge.r-project.org<br>
<b>Subject:</b> [adegenet-forum] colorplot and scatter<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div>Dear forum</div>
<div><br>
</div>
<div>I have two questions about the adegenet functions "colorplot" and "scatter"</div>
<div>I precise that I am a beginner with the package adegenet, so excuse me in advance if my questions seem a little stupid...</div>
<div><br>
</div>
<div>1. Colorplot</div>
<div><br>
</div>
<div>I made a colorplot to represent the result of a sPCA. The colorplot plots the different populations according to their localisation and the colors correspond to the global sPCA score. I would like to add a legend in a box representing the correspondance
between colors and sPCA score (as the image below). I search... but I did not find anything to do this... </div>
<div><br>
</div>
<div>
<p><img src="cid:426d0c7f-6e00-4dd2-a001-eeab8f80089a"></p>
<div><br>
</div>
<div><br>
</div>
<div>2. Scatter</div>
<div><br>
</div>
<div>I used the function scatter to represent the results of a DAPC. The colored dots in the plot represent the cluster, but I would like to know if it is possible that the dots represent populations instead of clusters. The clusters would be still represented
with the colored bars and the ellipses. I only want to change the dot colors. I tried to manually add dots in the plot (with the function "points"), using the coordinates of each indviduals in the DAPC axes. But the dots did not appeared... maybe due to a
scale problem that I cannot solve... </div>
<div>Here the function that I use for one population. The data frame "ind" contains the DAPC LD1 and LD2, and the population names of the different individuals :</div>
</div>
<div>
<div><br>
</div>
<div><span style="font-size:13.3333330154419px">points(ind$LD1[ind$pop=="Aussois"]), scale(ind$LD2[ind$pop=="Aussois"]), col="black", pch=16, cex = 1)</span></div>
<div><span style="font-size:13.3333330154419px"><br>
</span></div>
<div><span style="font-size:13.3333330154419px"><br>
</span></div>
<div>
<div>I hope that my questions are comprehensible. Thank you very much in advance for your help. Don't hesitate to let me know if you need further information to more understand my problems.</div>
<div><br>
</div>
<div>Thanks a lot!</div>
</div>
<div><br>
<div style="font-family:Tahoma; font-size:13px"><i><font size="4">Coraline Bichet</font></i>
<div>
<div><font size="2">coraline.bichet@univ-lyon1.fr</font></div>
<div><font size="2">+33(0)472433584</font></div>
<div>UMR-CNRS 5558, Laboratoire de Biométrie et Biologie Evolutive (LBBE) <br>
<div>Université Claude Bernard Lyon 1, bâtiment Mendel</div>
<div>43 boulevard du 11 novembre 1918</div>
<div>69622 Villeurbanne Cedex, FRANCE</div>
<div><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>