1

Survival Packageの条件付きロジスティック回帰の例を見てみましょう

そして、次のコマンドを使用して

library(survival)

data(logan)

resp <- levels(logan$occupation)
n <- nrow(logan)
indx <- rep(1:n, length(resp))
logan2 <- data.frame(logan[indx,],
                     id = indx,
                     tocc = factor(rep(resp, each=n)))
logan2$case <- (logan2$occupation == logan2$tocc)
B <- clogit(case ~ tocc + tocc:education + strata(id), logan2)

これで回帰パラメータを生成できますが、具体的toccfarmには -1.896 の値が必要だとしましょう。

これだけを出力したり、xsave this として保存したりするにはどうすればよいでしょうか。

使用する場合

B$coefficients

すべての回帰係数を取得します。

私は次のようなことを試しました

B$coefficients[1,]
B$coefficients(term=1)
B$coefficients("toccfarm")

しかし、どれも機能していません

4

1 に答える 1

-1

別のアプローチとして、summary関数を使用することもできます。要約すると、モデルの係数が行列として取られていることがわかります。

> is.matrix(summary(B)$coefficients)
[1] TRUE

この時点summary(B)$coefficientsで、オブジェクトに格納してから、必要に応じてサブセット化できます。

summary(B)$coefficients[1,1]
于 2015-05-29T13:26:21.027 に答える