[adegenet-forum] question scatter.dapc
Jombart, Thibaut
t.jombart at imperial.ac.uk
Thu May 19 10:40:45 CEST 2011
Dear Rudy,
the "d=..." simply indicates the mesh of the grid, which is not necessarily useful for a single graph, but may be used to compare different graphs (e.g. plane 1-2 and 3-4, to see how the distances in one plot relate to the other one). You can remove this indication using "cgrid=0" or remove the grid using "grid=FALSE".
About the inset, this is generated by the function add.scatter.eig. It scales so that the largest value uses most of the inset window. For instance:
plot(1:10)
add.scatter.eig(10:1,2,1,2)
add.scatter.eig(1:10,0,9,10,posi="bottomright")
About the inertia ellipses, they are graphical summaries of a cloud of points. This topic is actually being discussed in the adelist forum (forum for ade4):
http://listes.univ-lyon1.fr/wws/arc/adelist/2011-05/msg00003.html
They're not generally confidence ellipses. They are whenever the cloud of point is a sample from a bivariate normal distribution. In this case, "cellipse" (determines the size of the ellipse) indirectly controls the alpha threshold. 1.5 (default) corresponds to 67%; 2.5 corresponds to the magical 95%.
The exact formula is:
p = 1 − exp( − cellipse^2 /2 )
All the best
Thibaut
________________________________
From: adegenet-forum-bounces at r-forge.wu-wien.ac.at [adegenet-forum-bounces at r-forge.wu-wien.ac.at] on behalf of Rudy Jonker [mrjonker at gmail.com]
Sent: 18 May 2011 09:39
To: adegenet-forum at r-forge.wu-wien.ac.at
Subject: [adegenet-forum] question scatter.dapc
Dear Thibaut and group,
I have a question with respect to the graphs made with scatter.dapc. I noticed a d=n in the top right corner of the graph and was wondering what it stands for. I could not find it in the manual or in the article. I also have a question with respect to the eigenvalue inset in the graph. Are the bars scaled such that the first eigenvalue is always the full length of the graph and that the rest is then plotted relative to the first?
And then the inertia ellipse: is this a sort of 95% confidence interval on the position of the centroid of a cluster?
Thanks in advance,
Rudy
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20110519/80cbfc1b/attachment.htm>
More information about the adegenet-forum
mailing list