我在R中使用纵向医学数据进行了KML3d，在3个时间点测量了2个结果测量值(Eq5d评分+ oxford评分)。
预处理数据后，我使用以下代码构建模型:
trajectory_knee <-clusterLongData3d(data, timeInData = list(oxford = c(17, 19, 21), eq5d = c(18, 20, 22)), varNames = c("Oxford score", "Eq5d score"))
kml3d(trajectory_knee, nbClusters = 2:5, nbRedrawing = 4, toPlot = "both")
我的KML3d的目标是预测包含观测值的不同簇。然而，运行KML3D给了我calinski harabatz图，而我的目标是:
- 获取模型生成的聚类标签
- 绘制集群
- 使用BIC + plot轨迹找到最优簇数
然而，我不知道如何达到上述目标。有人能给我指个方向吗?
谢谢!
我试过使用
- NbClust
- Plotallcriterion
- 选择(轨迹)
但是这些都没有给我以前制定的目标的信息…
getClusters()是合适的函数。见https://www.jstatsoft.org/index.php/jss/article/view/v065i04