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:
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...
More information about the adegenet-forum