誰かがここで私をまっすぐにしてくれることを願っています。
ビネットからhttp://www.jstatsoft.org/v59/i12/paper、9ページ...
library("mixAK")
data("PBC910", package = "mixAK")
tail(PBC910)[, c(1, 2, 3, 6:9)]
...ドキュメント化された出力 (9 ページ) が得られるので、これまでのところすべて順調です。次に、被験者ごとに1つの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)")
...明確な反転を示しています。
どこかで明らかな何かを見逃したに違いない....