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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================