1

実行したい分析があり、対比を計画しました (事後比較ではありません!) 治療群間で行いたいと考えています。治療グループ変数には (k =) 4 つのレベルがあります。合計 3 つの異なる比較を行う予定です。そのため、私の理解が正しければ、比較は k-1 であるため、計算される p 値を調整する必要はありません。

これを行うには、Rのmultcomporパッケージを使用したいと思います。lsmeans私の質問は次のとおりです。信頼区間(およびp値)を調整せずに、この計画された比較を行うことが可能かどうかを誰かが知っていますか? 私が見たビネットと私が見た例からわかる限り、summary.glht()関数はデフォルトとして調整を行い、どのオプションがこれを元に戻すかは明確ではありません.

再現可能な例が必要な場合は、 http : //www.ats.ucla.edu/stat/r/faq/testing_contrasts.htmで見つけた次の例を使用できます。

    library(multcomp)

    hsb <- read.csv("http://www.ats.ucla.edu/stat/data/hsb2.csv")

    m1 <- lm(read ~ socst + factor(ses) * factor(female), data = hsb)
    summary(m1)

    K <- matrix(c(0, 0, 1, -1, 0, 0, 0), 1)
    t <- glht(m1, linfct = K)
    summary(t)
4

1 に答える 1