<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Hi, <br><br>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:<br><br>setwd ("C:\\Users\\......") <br>library(ape)<br>library(ade4)<br>library(seqinr)<br>library(graphComp)<br>library (MASS)<br>library(adegenet)<br>read.dna ("xxx.fasta", format="fasta")->x<br>DNAbin2genind(x, polyThres = 0.01)->xSNP<br>find.clusters(xSNP, choose.n.clust=FALSE,criterion="diffNgroup", clust=NULL)->cls_x<br><br>In this point, I should choose the number PC's to retain (>=1) but any chosen number I always obtain the following response:<br><br><b>Error en kmeans(XU, centers = nbClust[i], iter.max = n.iter, nstart = n.start) : <br> more cluster centers than distinct data points.</b><br><br>I know my data set show little variation, I don't know if it could produce such a result...<br><br>Please, someone can help me?<br><br><br>Thanks in advance!!<br><br><br><br> </div></body>
</html>