作者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