Displaying 2 results from an estimated 2 matches for "kseq".
Did you mean:
cseq
2006 Apr 19
1
determining optimal # of clusters for a given dataset (e.g. between 2 and K)
...ering the algorithms in "cluster" were developed decades ago. By the way, I have alot of missing values in my dataset, coded as "NA", so some software packages don't work.
Here is the code I've been using:
library(cluster)
avgsil <- c()
for (k in kseq){
clarares <- clara(data, k, rngR = TRUE)
savg <- clarares$silinfo$avg.width
print(c(k,savg))
avgsil[k] <- savg
}
k<-kseq
plot(k,avgsil[k])
lines(k,avgsil[k])
Sincerely,
Andrew McSweeny
grad student
Medical University of Ohio
[[alternative HTML version delete...
2006 Apr 15
0
clustering genes / automatically determining # of clusters
...ering the algorithms in "cluster" were developed decades ago. By the way, I have alot of missing values in my dataset, coded as "NA", so some software packages don't work.
Here is the code I've been using:
library(cluster)
avgsil <- c()
for (k in kseq){
clarares <- clara(data, k, rngR = TRUE)
savg <- clarares$silinfo$avg.width
print(c(k,savg))
avgsil[k] <- savg
}
k<-kseq
plot(k,avgsil[k])
lines(k,avgsil[k])
Sincerely,
Andrew McSweeny
grad student
Medical University of Ohio
[[alternative HTML version delete...