Rのパッケージを使用してlsmeans/emmeans、treatA (バイナリ/因子変数) のレベル間の応答におけるペアワイズ比較のプロットを作成しています。lsmeans (コントラスト) を使用して差の推定値を取得できますが、信頼限界ではなく、推定値の SE のみが提供されます。信頼限界は、対比ではなく、個々の効果に対してのみ提供されます。平均差 (コントラスト) の CL を生成するのを手伝ってくれる人はいますか?

最終的にはこのようなプロットを作成したいと考えていますが、SE ではなく CL を使用します。

最初にテーブルをアップロードします。表には 9 行が含まれており、そのうちの 6 行は要因であり、残りの 3 行は 152 人の個人 (n01、n02、n03) の成長率の個別の測定値です。次に、要因を指定します。

次に、興味のある要素を含む新しいテーブル "r2" にデータフレームを溶かし、na.omit 関数で NA 値を削除します。

r2 は次のようになります。

その後、固定分散を設定し、次のように 2 つの gls モデルを適用して実行します。

モデルの要約と anova を取得できるため、モデルは正常に機能しているようです。次に、次のように、パッケージ emmeans の lsmeans 関数を使用して事後ペアワイズ比較を実行しようとします。

lsmeans は、2 因子モデル mix1 で正常に機能するようです。ただし、モデル mix3 で lsmeans を実行すると、次のエラーが表示されます。

crossprod(x, y) のエラー: 数値/複素行列/ベクトルの引数が必要です

モデルを行列に変換しようとしましたが、lsmeans関数に対して正しいオブジェクトではありません。また、係数を設定せずに列を数値のままにしてみましたが、同じエラーが表示されます。lsmeans 関数について読んでいるときに、それに関連する crossprod 関数が見つかりません。

例として "mpg" データを使用して、非関数形式で lsmeans 関数を呼び出すコードを書きました。出力は適切です (以下に示すように)。

コードを関数形式に変更して同じ出力を生成しようとしましたが、データの列は私が書いた関数で認識できませんでした。

私の関数から報告されたエラーは次のとおりです。

私のコードは正常に機能する関数形式ではありません:

私が書いた関数コードが機能しませんでした:

「mpg」データを使用して関数を呼び出します。

非機能形式コードからの出力:

私の関数フォーマットコードからのエラー: