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:[email protected]
[email protected]
======================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
