2015-05-29 6 views
1

Давайте условный logisitc пример регрессии от Survival Packageполучить индивидуальное значение Коэффициент из clogit в R

И используя следующие команды

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.

Как мы будем выводить это или сохранить его как x Сохраните это.

Когда мы используем

B$coefficients 

мы получаем все коэффициенты регрессии.

Я пытался что-то вроде

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

Но никто из них не работал

+2

'B $ коэффициентов [1]' или 'коэффици- (B) [1]' или 'B $ коэффициенты [ 'toccfarm']' должен работать – akrun

+0

да я только что нашел пример коэф() функция – lukeg

ответ

-1

Другой подход может быть с использованием функции summary. Вы можете видеть, что с суммированием коэффициенты модели берутся как матрица.

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

На данный момент вы можете хранить summary(B)$coefficients в объекте, а затем подмножество, как вы хотите.

summary(B)$coefficients[1,1] 

 Смежные вопросы

  • Нет связанных вопросов^_^