R返回從GLM cbind

Question

統計和R小白內的p值知道是否有有沒有辦法從GLM添加p值到從下面的命令而產生的輸出的結束：

exp(cbind(OR = coef(mod1), confint(mod1))) 

也許像這樣：

summary(mod1)$coefficients[,4] 

我知道這是一個'化妝品'的問題，但它會很方便。

感謝

Answer 1

您可以保存彙總的結果（MOD1），然後使用coefficients訪問係數表。

您可以編寫一個函數，會爲你做的全過程......

OR.summary <- function(x){ 
# get the summary 
    xs <- summary(x) 
# and the confidence intervals for the coefficients 
    ci = confint(x) 
# the table from the summary object 
    coefTable <- coefficients(xs) 
# replace the Standard error/test statistic columns with the CI 
    coefTable[,2:3] <- ci 
# rename appropriatly 
    colnames(coefTable)[2:3] <- colnames(ci) 
# exponentiate the appropriate columns 
    coefTable[,1:3] <- exp(coefTable[,1:3]) 
# return the whole table.... 
    coefTable 

} 

更健壯的方法是使用包像rms ....