[adegenet-forum] DAPC

Jombart, Thibaut t.jombart at imperial.ac.uk
Thu May 19 10:51:41 CEST 2011


Hello,

you may have forgotten to provide group memberships. What does "pop(obj)" return? If it is missing, you need to provide something for the argument "pop". See ?find.clusters to identify clusters without priors.

Cheers

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 Martin Llewellyn [llewellynmartin at hotmail.com]
Sent: 18 May 2011 16:46
To: adegenet-forum at r-forge.wu-wien.ac.at
Subject: [adegenet-forum] DAPC

Hi Thimbaut and adegenet users

Thanks for your help on getting the package running. All working now.

I am attempting a DAPC on a genind object containing no missing data with 264 isolates and 24 microsat loci with the arguement:

dapc <- dapc(obj, pop=NULL, n.pca=NULL, n.da=NULL, scale=TRUE,scale.method=c("sigma"), truenames=TRUE, all.contrib=FALSE,pca.select=c("nbEig"), perc.pca=NULL)

Choose the number PCs to retain (>=1): (I put 2 - 15)

Message returned is:-

Error in svd(X, nu = 0) : infinite or missing values in 'x'
In addition: Warning message:
In .local(x, ...) : Some scaling values are null.
Corresponding alleles are removed.

Any advice most appreciated

Thanks

M
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20110519/e5c29cb4/attachment.htm>


More information about the adegenet-forum mailing list