代替タイトル: モデル マトリックスと係数のセットは、変数の数が異なることを示しています

私はR用のマウスパッケージを使用していくつかの分析を行っています。を使用して ( miraオブジェクトに保持されている) 2 つのモデルを比較したかったのですが、次のエラーが発生し続けます。pool.compare()

Error in model.matrix(formula, data) %*% coefs : non-conformable arguments

二項演算子は、 R%*%での行列の乗算を示します。

この式model.matrix(formula, data)は、「指定された数式とデータを使用した回帰のようなモデルの設計行列」を生成します ( model.matrix {stats}のRドキュメントから)。

エラー メッセージでは、はmipoオブジェクトcoefsから描画され、要素は「完全なデータ推定値の平均。多重代入推定値」です。( mipo-class {mice}のドキュメントから)。 est1$qbarest1qbar

私の場合

• est1$qbar長さ 36 の数値ベクトル
• data82 個の変数の 918 個の観測値を含む data.frame です。
• formula私のモデルの式を含むクラス「式」です
• model.matrix(formula, data)次元が 918 x 48 の行列です。

このエラーを解決/防止するにはどうすればよいですか?

R のMICEパッケージを使用して、一連の変数に対して複数の代入を実行したいと考えています。

の最小値gcsは3で、最大値は15で、小数ではない可能性があります。MICEでこれらの制約を設定するにはどうすればよいですか? にも同じことが言えますが、これには下限が0hfしかありません。

R でマウス アルゴリズムを実行して得た mids オブジェクトがあります。plot.mids を使用して収束を評価しようとしています。指定された 2 つの変数の収束診断をプロットできるようにしたいと考えています。次に例を示します。

ここでは、変数 chl のプロットのみが表示されます。戻る矢印をクリックすると他のプロットも表示されますが、レポートの同じプロットに chl と bmi を一緒に表示したいと思います。y=式で変数を指定してみましたが、使い方がわかりません。

それは、chl の手段を bmi の手段に対してプロットしようとしているようでした (私が思うに)、これは私が求めているものではありません。

コンポーネントが欠落しているデータがあるため、マウスアルゴリズムを実行しました（パッケージからmice）。この関数は .mids オブジェクトを返します。これをトレーニング データセットとテスト データセットに分割して、モデルの適合性を評価します。トレーニング データとテスト データも .mids 形式にしてpool、ルービンの規則に従って標準誤差を調整するなど、他のさまざまな機能と組み合わせることができるようにしたいと考えています。

これは、トレーニングセットを取得するためにデータから行を削除するだけの私の試みです:

しかし、このデータでモデルを実行しようとすると:

9 行を 7 行に置き換えようとしているというエラーが表示されます (おそらく、他のものも調整せずに train$data の行数を減らしたためです)。

どんな助けでも大歓迎です。

連続ベクトルに欠損値があります

私は郵便番号を持っています

郵便番号を使用して y を代入したいと思います。私が試してみました

しかし、クロス積を実行できないことを示すエラーが表示されます (rstudio セッションがハングしているため、特定のメッセージを投稿できません)。

マウスまたは別のパッケージでこれを行う方法についてのアイデアはありますか?

編集:dputデータの一部

MICEパッケージで R に代入した後、分割表を生成したいと考えています。フィットはリスト内のテーブルを表示しますが、 pool()それらを使用すると、次のエラーがスローされます:Error in pool(fit) : Object has no coef() method.何が間違っていますか?

この基本的な例は、エラーを再現します。