Prezad@s

Alguma ideia de solução na seguinte dificuldade:

 

library(survival)

fit <- survreg(Surv(futime, fustat) ~ rx + age, data = ovarian, dist = "weibull")

 

Error in coxph.wtest(t(x) %*% (wt * x), c((wt * eta + weights * deriv$dg) %*%  : object 'Ccoxph_wtest' not found

 

fit <- survreg(Surv(futime, fustat) ~ 1, data = ovarian,  dist = "weibull")

summary(fit)

 
Call:
survreg(formula = Surv(futime, fustat) ~ 1, data = ovarian, dist = "weibull")
             Value Std. Error      z         p
(Intercept)  7.111      0.293 24.292 2.36e-130
Log(scale)  -0.103      0.254 -0.405  6.86e-01
 
Scale= 0.902 
 
Weibull distribution
Loglik(model)= -98   Loglik(intercept only)= -98
Number of Newton-Raphson Iterations: 5 
n= 26 

 

Ou seja, a função só funciona no ajuste do modelo sem variáveis explanatórias. Observação: o pacote survival tem a função coxph.wtest(.).

Desde já agradeço



Gilênio Borges Fernandes
Professor Associado IV (Aposentado)
Professor Adjunto A (Substituto)
Universidade Federal da Bahia
Instituto de Matemática e Estatística
Departamento de Estatística
Av. Adhemar de Barros, s/n – Ondina.
40.170-110 - Salvador - BA, Brasil
Tel.: (071)3283-6340/6341/6337  Fax:  (071)3283-6336
Skype: gilenio.fernandes
Lattes: http://lattes.cnpq.br/6764860618464860