vendredi 13 mai 2016

Get the same pam average Silhouettes using different seeds in R

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