[R-br] Problema com TL 50 em análise de sobrevivência

ASANTOS alexandresantosbr em yahoo.com.br
Terça Janeiro 22 13:12:08 BRST 2013


Boa tarde pessoal,


         Estou com um problema para calcular a TL50 (tempo de morte de 
50% do indivíduos de uma população), usando o modelo de sobrevivência 
paramétrico de Weibull, pois o valor extraído do summary do modelo com 
exp(Intercept) não esta batendo com o curve do modelo, para explicar 
melhor segue CRM:

require(survival)

##Dados obtidos na web
A <- read.table("http://www.mas.ncl.ac.uk/~nmf16/teaching/mas3311/A.dat")

attach(A)


## Ajuste do modelo

mod = survreg(Surv(time, status)~1,dist="weibull")
summary(mod)

###########################   Calculando os parâmetros alpha e tempo médio


##########   Valor de Alpha (1/Scale)


1/0.891 # =1.122334


##TL 50 exp do (Intercept)  do modelo


exp(3.887)  #= 48.76437

##Plotando


par(cex.lab=1.8,bty="n",yaxs="r",xaxs="r",lwd=2,mar=c(4.5,4.5,1,1),mex=1.3,family="")

curve(exp((-(48.76437)^(-1.122334))*(x^1.122334)),xaxt="n",yaxt="n",from=0,to=170,ylim=c(0,1),xlab="",ylab="")

x=seq(0,170,10)
y=seq(0,1,0.1)
axis(1,x,lwd=2,cex.axis=1.5)
axis(2,y,lwd=2,cex.axis=1.8)

#

abline(h=0.5, col="red")
abline(v=48.76437, col="red")

## O que não bate o valor encontrado com a curva


Alguém poderia de dar uma ajuda,



Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Avenida dos Ramires, s/n
Caixa Postal 744
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
======================================================================



Mais detalhes sobre a lista de discussão R-br