mixAK getProfiles 错误？

Question

我希望有人能在这里理顺我。

来自插图 http://www.jstatsoft.org/v59/i12/paper，第 9 页...

library("mixAK")
data("PBC910", package = "mixAK")
tail(PBC910)[, c(1, 2, 3, 6:9)]

...到目前为止一切都很好，因为这给出了记录的输出（第 9 页）。然后下面应该给每个主题一个data.frame（从第11页开始）。

ip <- getProfiles(t = "month", y = c("lbili", "platelet", "spiders",    
     "jspiders"), id = "id", data = PBC910)

所以要测试...

 ip[[1]]

正确地给出了第一个子集（尽管 id=2，而不是一个更有用的索引）。

然而，深入挖掘，子设置似乎不同步，至少对我而言。例如，通过检查结果...

ip[[11]];ip[[12]];ip[[13]];ip[[14]];ip[[15]]

这里显示...

ip[[12]]

我得到以下内容，它从前一个子集的最后一个条目开始。子集的末尾溢出了，有点乱。

              month      lbili       platelet  spiders    jspiders
           44 23.425051 -0.3566749   133       0          0.207410112
           45  0.000000 -0.2231436   295       0          0.003706764
           46  6.439425 -0.9162907   235       1          0.709859749
           47 12.024641 -0.2231436   268       0          0.146529978

绘制结果显示相同的效果...

 plotProfiles(ip = ip, data = PBC910, var = "lbili", tvar = "month", 
     main = "Log(bilirubin)", highlight = c(12), 
     xlab = "Time (months)", ylab = "Log(bilirubin)")

...显示出明显的逆转。

我一定是漏掉了一些明显的东西....

Answer 1

这是在 4.1 版软件包中错误创建的错误。它应该在 4.2 版中得到纠正（希望很快在 CRAN 上可用）。

阿诺斯特

Answer 2

同意 - 正如 Arnost 友善地指出的那样，4.2 版修复了这个错误。

对于任何其他关注有间隙纵向数据聚类的人来说，这个包就是蜜蜂的膝盖！

mixAK getProfiles 错误？

mixAK getProfiles bug?

r

cluster-analysis

mixture-model