作者littlehana (我應該又做錯了~)
看板NCTU-STAT94G
標題倖存第一題第二小題-R的code
時間Tue Mar 14 19:37:59 2006
A<-matrix(c(11,12,15,45,16,17,19,28,30,33,1,1,1,1,0,0,0,0,0,0),ncol=2)
time<-A[,1]
status<-A[,2]
AA<-survreg(Surv(time,status)~1,dist="weib")
AA
summary(AA)
muhat<-AA$coeff
sigma<-AA$scale
mletheta<-exp(-1*muhat)
mletheta
mlealpha<-(1/sigma)
mlealpha
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.114.149
1F：→ littlehana:這是使用survival regression的方法找mle 有疑問可以問 03/14 19:46
2F：→ littlehana:我 不過之後 老師應該會教到這個部份 03/14 19:47
3F：→ littlehana: 記得要先下載survival pkg喔~~ 03/14 19:47
4F：推 roussas:找不到package的可以來找我...=.=|| 03/15 00:08