如何在R中处理KML3D结果？

Question

我使用纵向医学数据在 R 中完成了 KML3d，并在 3 个时间点测量了 2 个结果测量值（Eq5d 分数 + 牛津分数）。

对数据进行预处理后，我使用以下代码来构建模型：

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 habatz 图，而我的目标是：

• 获取模型生成的聚类标签
• 绘制集群
• 使用 BIC + 绘制轨迹来找到最佳簇数

但是，我不知道如何才能达到上述目标。有人可以帮助我/引导我走向正确的方向吗？

谢谢！

我尝试使用

• 铌簇
• 绘图所有标准
• 选择（轨迹）

但这些都没有给我有关先前制定的目标的信息......

Answer 1

getClusters() 是合适的函数。请参阅https://www.jstatsoft.org/index.php/jss/article/view/v065i04