For each k, I use different seeds and add each pam average Silhouette to sil array. the code is like bellow:
sil <- numeric(5)
for(k in 3:10){
for (i in 1:5){
set.seed(i)
PAM<-pam(RFML_table,k)
sil[[i]]<-PAM$silinfo$avg.width
}
minSil[[k]]<-min(sil)
maxSeed[[k]]<-which.max(sil)
}
kClust<-which.max(minSil)
kSeed<-maxSeed[kClust]
However, each time, my sil array's elements are equal. For example, sil elements for k=3 are 0.8695872 0.8695872 0.8695872 0.8695872 0.8695872 . My question is why I get same elements in each k. I expect different sil elements because pam results may be different for each seed and, as a consequence, average Silhoettes for pams may be different!
Aucun commentaire:
Enregistrer un commentaire