作者vioan223 (vioan223)
看板R_Language
标题[问题] silhouette 使用方法
时间Sun Mar 15 13:38:18 2015
[问题类型]:
程式谘询
[软体熟悉度]:
新手(没写过程式，R 是我的第一次)
[问题叙述]:
我想使用silhouette这个package
我把资料改成我的(其他都没改)
但结果怪怪的
我想可能是有些变数必需修改
但因为我看不懂程式码
不知道从何改
想请问有人成功使用过这个package吗?
可以教我吗?
我会非常感激
谢谢
package的pdf档网址
http://ftp.yzu.edu.tw/CRAN/web/packages/cluster/cluster.pdf
在第66页
--------------------------------------------------------------------
我复制example的范例
把原先"ruspini"的资料
改成我自己的资料"abbc"
然後将原范例算距离的方法"canberra"
改成另一个算距离的方法"euclidean"
其他都没改 只改了这两部分
我最终的目的就是希望透过这个方法
得到我的资料应该分K群
不知道程式这样改之後是不是对的
再麻烦您帮我看看
谢谢
程式码如下:
pr4 <- pam(abbc, 4)
str(si <- silhouette(pr4))
(ssi <- summary(si))
plot(si) # silhouette plot
plot(si, col = c("red", "green", "blue", "purple"))# with cluster-wise
coloring
si2 <- silhouette(pr4$clustering, dist(abbc, "euclidean"))
summary(si2) # has small values: "euclidean"'s fault
plot(si2, nmax= 80, cex.names=0.6)
op <- par(mfrow= c(3,2), oma= c(0,0, 3, 0),
mgp= c(1.6,.8,0), mar= .1+c(4,2,2,2))
for(k in 2:6)
plot(silhouette(pam(abbc, k=k)), main = paste("k = ",k), do.n.k=FALSE)
mtext("PAM(abbc) as in Kaufman & Rousseeuw, p.101",
outer = TRUE, font = par("font.main"), cex = par("cex.main")); frame()
## the same with cluster-wise colours:
c6 <- c("tomato", "forest green", "dark blue", "purple2", "goldenrod4",
"gray20")
for(k in 2:6)
plot(silhouette(pam(abbc, k=k)), main = paste("k = ",k), do.n.k=FALSE,
col = c6[1:k])
par(op)
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 115.43.79.63
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1426397900.A.CA9.html
1F：→ andrew43: 该package的网址可否提供一下？ 03/15 15:15
你好 我附上pdf档的网址了 谢谢你
※ 编辑: vioan223 (115.43.79.63), 03/15/2015 18:50:52
2F：推 andrew43: 谢谢，提醒你那是cluster package。 03/15 19:57
3F：推 andrew43: 另外建议你提出具体问题，不然根本无从帮起。 03/15 20:01
※ 编辑: vioan223 (115.43.79.63), 03/15/2015 22:06:45
※ 编辑: vioan223 (115.43.79.63), 03/15/2015 22:07:20
4F：→ andrew43: 没有说到重点：abbc的内容是什麽？ 03/15 23:37
5F：→ andrew43: 请考虑让其它人能重现你卡关的方式来问。 03/15 23:38
abbc资料的内容是有九个数值属性 像是年龄、身高这种 共有1千多笔资料 谢谢你~
※ 编辑: vioan223 (115.43.79.63), 03/16/2015 08:22:42
6F：→ celestialgod: 发文前请看一下注释，"提供可重现错误的code" 03/16 12:26
7F：→ celestialgod: 建议提供partial data 版友才能比较快的提供协助 03/16 12:26