[adegenet-forum] DAPC

josé manuel lucas cánovas josemanuel.lc at hotmail.com
Wed Mar 28 12:19:04 CEST 2012


I'm trying to do a DAPC with a sequence matrix of mtDNA. I have 89 sequences (they collapsed in 7 haplotypes) and this is my batch:

setwd ("C:\\Users\\......") 
library (MASS)
read.dna ("xxx.fasta", format="fasta")->x
DNAbin2genind(x, polyThres = 0.01)->xSNP
find.clusters(xSNP, choose.n.clust=FALSE,criterion="diffNgroup", clust=NULL)->cls_x

In this point, I should choose the number PC's to retain (>=1) but any chosen number I always obtain the following response:

Error en kmeans(XU, centers = nbClust[i], iter.max = n.iter, nstart = n.start) : 
  more cluster centers than distinct data points.

I know my data set show little variation, I don't know if it could produce such a result...

Please, someone can help me?

Thanks in advance!!

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20120328/1e529013/attachment.html>

More information about the adegenet-forum mailing list