juanzi
一、参数回归模型条件下(误差项服从指数分布)存在异常点的情形:
参数回归模型为 status<-c(1:100)
time<-c(1:100)
x1<-runif(100,min=0.5,max=1.5)
x2<-rbinom(100,1,0.4)
e<-rexp(100,rate=1)
logt<-2*x1+0.6*x2+e
z<-runif(100,min=2,max=3)
logc<-1.6*z
for (i in 1:100){if (logt<logc) {time<-logt;status<-1;} else {time<-logc;status<-0;}}
异常点的比例假定为5%该如何设定?
二、满足PH假定的COX数据
生存时间还是服从指数分布,自变量同上。就是还想在不同截尾比例下考察看看截尾百分位数估计量的无偏性和有效性。截尾的比例(分别是20%、40%、80%)。不知道该如何控制截尾比例?