"グループ化変数に純粋な数値の項目を含めることはできません"

Question

グループで分析しているので、私は性別（0 =男性、1 =女性）を使用するグループ変数を作成する必要があります。私が最初にしたことは、この変数のベクトルを作成することでした（手動でこれを行うように指示しました）が、「変数をグループ化するには純粋な数値の項目を含めてはいけません」というエラーがありました。その後、私はベクトルを論理的（TRUE/FALSE）に変換しましたが、どういうわけか私はまだこのエラーを受け取ります。

私はこのエラーが発生したときに何が問題になるのでしょうか。以下に添付

は私のデータセットの頭にコードです：

structure(c(7, 8, 7, 5, 6, 6, 4.9, NA, 6.9, 5.1, 5.8, NA, NA, 
     NA, 7, 3, 7, NA, NA, NA, 6.7, 4.1, 5.9, NA, NA, NA, 5, 6, 7, 
     NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 8, NA, NA, NA, 6.2, 
     4.3, 6.3, NA, NA, NA, 7, 5, 7, NA, NA, NA, 6.5, NA, NA, NA, NA, 
     NA, 6, NA, 7, NA, NA, NA, NA, NA, 5, NA, NA, NA, NA, NA, 7, NA, 
     NA, NA, NA, NA, 6.1, NA, NA, NA, NA, NA, 7, NA, NA, NA, NA, NA, 
     NA, NA, 16, 0.001, 12, 11, 11, 0.001, 0.001, 0.001, 12, 12, 12, 
     0.001, 0.001, 0.001, 12, 12, 12, 0.001, 0.001, 0.001, 15, 12, 
     12, 0.001, 0.001, 0.001, 16, 0.001, 12, 0.001, 0.001, 0.001, 
     0.001, 0.001, 15, 0.001, 0.001, 0.001, 0.001, 0.001, 16, 0.001, 
     0, 1, 0, 0, 1, 0), .Dim = c(6L, 24L), .Dimnames = list(c("800009", 
                   "800012", "800015", "800033", "800042", "800045"), c("gener_sat_T0", 
                                 "sel_T0", "gener_sat_T1", "sel_T1", "gener_sat_T2", "sel_T2", 
                                 "gener_sat_T3", "sel_T3", "gener_sat_T4", "sel_T4", "gener_sat_T5", 
                                 "sel_T5", "gener_sat_T6", "sel_T6", "gener_sat_T7", "sel_T7", 
                                 "dT1", "dT2", "dT3", "dT4", "dT5", "dT6", "dT7", "female"))) 

その後、私は何をしようとしていますがCTモデルに適合している（非グループデータに前にそれを使用し、それがうまく働いています） 。

ありがとうございました！あなたのコードが直接再現性がなかったとはいえ

Answer 1

[OK]を、私は、あなたの計算をやり直し、私はいくつかの変更を加えて、今では動作します：

だから最後にそれがグループ化変数の文字ベクトルを作るだけの問題だ

# create the structure object (data_wide), and change it to remove the 
# grouping: 

data_wide = structure(c(7, 8, 7, 5, 6, 6, 4.9, NA, 6.9, 5.1, 5.8, NA, NA, 
    NA, 7, 3, 7, NA, NA, NA, 6.7, 4.1, 5.9, NA, NA, NA, 5, 6, 7, 
    NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 8, NA, NA, NA, 6.2, 
    4.3, 6.3, NA, NA, NA, 7, 5, 7, NA, NA, NA, 6.5, NA, NA, NA, NA, 
    NA, 6, NA, 7, NA, NA, NA, NA, NA, 5, NA, NA, NA, NA, NA, 7, NA, 
    NA, NA, NA, NA, 6.1, NA, NA, NA, NA, NA, 7, NA, NA, NA, NA, NA, 
    NA, NA, 16, 0.001, 12, 11, 11, 0.001, 0.001, 0.001, 12, 12, 12, 
    0.001, 0.001, 0.001, 12, 12, 12, 0.001, 0.001, 0.001, 15, 12, 
    12, 0.001, 0.001, 0.001, 16, 0.001, 12, 0.001, 0.001, 0.001, 
    0.001, 0.001, 15, 0.001, 0.001, 0.001, 0.001, 0.001, 16, 0.001), .Dim = 
c(6L, 23L), 
    .Dimnames = list(c("800009", "800012", "800015", "800033", "800042", 
"800045"), 
        c("gener_sat_T0", "sel_T0", "gener_sat_T1", "sel_T1", 
"gener_sat_T2", 
         "sel_T2", "gener_sat_T3", "sel_T3", "gener_sat_T4", 
"sel_T4", "gener_sat_T5", 
         "sel_T5", "gener_sat_T6", "sel_T6", "gener_sat_T7", 
"sel_T7", 
         "dT1", "dT2", "dT3", "dT4", "dT5", "dT6", "dT7"))) 


CTMODEL <- ctModel(n.latent = 2, n.manifest = 2, Tpoints = 8, 
       manifestNames = c("gener_sat", "sel"), 
       latentNames = c("gener_sat", "sel"), LAMBDA = diag(2)) 
fem = c("f", "m", "f", "f", "m", "f") # grouping, which needs to be a 
# character vector 

fit_CTMODEL <- ctMultigroupFit(dat = data_wide, groupings=fem, ctmodelobj = 
CTMODEL) # dat instead of datawide 

。

追加：コードが実行されますが、じんましん様々なエラー：

Not all eigenvalues of Hessian are greater than 0 
Fit attempt generated errors 
Retry limit reached 

私はあなたが右であることが原因で、モデルのだと思うし、あなたに解決策を残す:)