1つのレベルと他のレベルの平均との間の差のコントラストマトリックス（R）を指定する方法は？

私はこのように見えることが回帰モデルを有する：

...またはR表記： y ~ x1 + x2 + x3 + x1:x2 + x1:x3 + x1:x2:x3

とX 2がカテゴリー変数で、X 3が数値であるとしましょう。複雑なのは、X 1に3つのレベルX 1 a、X 1 b、X 1 cがあり、標準的なコントラストの代わりに、テストする必要があることです。$X_1$$X_2$$X_3$$X_1$$X_{1a}, X_{1b}, X_{1c}$

• レベルについてインターセプトか有意レベルの平均切片は異なるX 1 、BおよびX 1個のC。$X_{1a}$$X_{1b}$$X_{1c}$
• 応答かどうかレベル間で有意に異なっているX 1 、A及びレベルの平均値X 1 、BおよびX 1個のC。$X_2$$X_{1a}$$X_{1b}$$X_{1c}$
• 傾きか否有意レベル間で異なるX 1 、A及びレベルの平均値X 1 、BおよびX 1個のC。$X_3$$X_{1a}$$X_{1b}$$X_{1c}$

この投稿に基づいて、私が欲しいマトリックスは...

 2
-1
-1

contrasts(mydata$x1)<-t(ginv(cbind(2,-1,-1)))$\beta_1$$beta_1$$X_1b$$X_1c$$X_3=0$$X_2$$X_1a$

誰かが細胞の平均とコントラストの関係に頭を抱える方法について何かアドバイスはありますか？ありがとう。このタイプのコントラストの標準的な名前はありますか？

ああ！Glen_bの回答に掲載されているリンクによると、一番下の行は、次のように、グループ平均の任意の比較をRスタイルのコントラスト属性に変換できることです。

1. 正方行列を作成します。行は因子のレベルを表し、列は対比を表します。最初のものを除いて、切片が表すものをモデルに伝えます。
2. 切片を総平均にしたい場合は、最初の列にゼロ以外の同じ値をすべて入力します。これは問題ではありません。切片をレベル平均の1つにしたい場合は、その行に数値を入力し、残りをゼロで埋めます。切片をいくつかのレベルの平均にしたい場合は、それらの行に数値を入れ、残りにゼロを入れます。加重平均にする場合は、異なる数値を使用します。それ以外の場合は、同じ数値を使用します。切片の列に負の値を入力することもでき、それはおそらく何かを意味しますが、それは他のコントラストを完全に変更します。
3. 他の列と比較して必要なレベルを示す正と負の値を残りの列に入力します。ゼロへの合計が重要な理由を忘れていますが、列が合計してゼロになるように値を調整します。
4. t()関数を使用して行列を転置します。
5. パッケージginv()から、MASSまたはsolve()転置行列の逆を取得するために使用します。
6. たとえば、最初の列をドロップしますmycontrast<-mycontrast[,-1]。これでapx p-1マトリックスが作成されましたが、インターセプトに入力した情報は、手順5でマトリックス全体にエンコードされました。
7. 要約出力のラベルをlm()他のデフォルト出力よりも読みやすくしたい場合は、それに応じてマトリックスの列に名前を付けます。(Intercept)ただし、切片は常に自動的に名前が付けられます。
8. マトリックスを問題の因子の新しいコントラストにします。例： contrasts(mydata$myfactor)<-mymatrix
9. ファイル名を指定して実行lm()（および数式を使用し、おそらく他の多くの機能）負荷にせずに、標準Rで通常通りglht、doByまたはcontrasts。

Glen_b、ありがとう、そしてUCLA Statistical Consulting Groupに感謝します。私の応用統計学教授はこのトピックに数日間手を振りましたが、実際に自分のコントラストマトリックスを実際に作成する方法がわかりませんでした。そして今、1時間のRを読んで遊んでいましたが、ようやく理解できたと思います。代わりにUCLAに応募すべきだったと思います。またはStackExchange大学。

$(0, 1, -1)^\top$

Rがヘルマートコーディングと呼ぶもので、これは「逆ヘルマート」と呼ばれます。変数の順序を変更するまでは同等です。