次の長いデータセットがあります (「dat」と呼ばれる R データフレーム)。1200 人の子供の約 10,000 の観測があり、子供ごとに最大 10 の観測があります。これらの観察結果は、出生から 5 歳までの幼児期のさまざまな時点で収集された子供の体重です。データセットの年齢は日単位 (変数 = agedays) で、体重は変数 wtkg のキログラム単位で与えられます。
structure(list(subjid = c(1001L, 1001L, 1001L, 1001L), sex = structure(c(2L, 2L, 2L, 2L), .Label = c("Female", "Male"), class = "factor"),
agedays = c(0L, 2L, 30L, 107L), wtkg = c(3.78, 3.64, 4.71,
6.5), Uobservations = c(10L, 10L, 10L, 10L), BMI_group = structure(c(2L,
2L, 2L, 2L), .Label = c("normal", "obese", "overweight",
"true_NA", "underweight"), class = "factor"), GWG_T3_cat2 = structure(c(2L,
2L, 2L, 2L), .Label = c("T3notexcessive", "T3excessive"), class = "factor")), .Names = c("subjid", "sex", "agedays", "wtkg", "Uobservations", "BMI_group", "GWG_T3_cat2"), row.names = c(NA, 4L), class = "data.frame")
SITAR パッケージを使用して、肥満の母親 (変数 = BMI_group) の間で (変数 = GWG_T3_cat2) (変数 = GWG_T3_cat2) 妊娠 3 期の過度の妊娠中の体重増加に応じて、子供の速度、テンポ、およびサイズを調査したいと考えています。 .
このリンクにあるモデルを実行してみました: https://www.rdocumentation.org/packages/sitar/versions/1.0.9
m1 <- sitar(x=agedays, y=wtkg, id=subjid, data=dat, df=2)
しかし、私はエラーが発生します:
Error in nlme.formula(y ~ fitnlme(x, s1, s2, a, b, c), fixed = s1 + s2 +:
Singularity in backsolve at level 0, block 1
誰かがこの問題の解決に役立つことができれば、本当に感謝しています。