[adegenet-forum] DAPC

[josé manuel lucas cánovas]

Hi, 

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(ape)
library(ade4)
library(seqinr)
library(graphComp)
library (MASS)
library(adegenet)
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>