<html>
<body>
Dear Thibaut Jombart and List members,<br><br>
I have three questions concerning the mapping of DAPC results:<br><br>
1) after find.clusters, it appear that individuals from the 5 samples
sites are genetically grouped into 3 clusters. I would like to map this
(as in Fig 7, Jombart et al BMC Genetics), but given that I have 5
localities (xy) and only 3 clusters, isn't it normal that I do not have
the same number of row numbers in xy and X? (?? dem to post from<b>
</b>Bernhard Stoeckle??). How do I resolve this error (Erreur dans
colorplot.default(xy, -dapc.a$grp.coord, cex = 3, add = TRUE, :
<br>
xy and X have different row numbers) ?<br><br>
> dim(xy)<br>
[1] 5 2<br>
> dim(dapc.a$grp.coord)<br>
[1] 3 2<br><br>
> if(require(maps)){<br>
+ <br>
+ xy <- cbind(data$other$popInfo$Longitude,
data$other$popInfo$Latitude)<br>
+ <br>
+ par(mar=rep(.1,4))<br>
+ map(fill=TRUE, col="lightgrey")<br>
+ colorplot(xy, -dapc.a$grp.coord, cex=3, add=TRUE, trans=FALSE)<br>
+ }<br>
Erreur dans colorplot.default(xy, -dapc.a$grp.coord, cex = 3, add =
TRUE, : <br>
xy and X have different row numbers<br><br>
2) where can one download maps from countries other than those provided
in the map package?<br><br>
3) finally, when there are 2 BIC values that are very similar, should one
choose for the smaller number of clusters even if it has a higher BIC
value?<br><br>
Thank you for your help!<br><br>
Sincerely yours,<br><br>
Sofie <br>
<x-sigsep><p></x-sigsep>
Sofie Vandewoestijne<br>
Postdoctoral fellow<br>
Biodiversity Research Centre<br>
Catholic University Louvain<br>
Place Croix du Sud 5<br>
B-1348 Louvain-la-Neuve<br>
Belgium<br><br>
Tel: 32 (0)10 473454<br>
Fax: 32 (0)10 473490<br><br>
<font face="Webdings" size=5 color="#009999"><b>P</font>
<font color="#009999">
</font><font face="Verdana" size=1 color="#008080">Please consider your
environmental responsibility before printing this e-mail<br>
</font></b></body>
</html>