The following object is masked from 'package:base':
backsolve
par(mfrow =c(1,2))data(engel);plot(engel)plot(log10(foodexp)~log10(income),data = engel,main="'engel' data ")taus <-c(0.15,0.25,0.5,0.75,0.95,0.99)rqs<-as.list(taus)for (i inseq(along=taus)){ rqs[[i]]=rq(log10(foodexp)~log10(income),tau=taus[i],data=engel)lines(log10(engel$income),fitted(rqs[[i]]),col=i+1)}legend("bottomright",paste("tau=",taus),inset=0.02,col=2:(length(taus)+1),lty=1)
13.2 生存分析的Cox模型
T为随机变量,在活过时间t的概率为
\[
S(t)=P(T>t)=1-P(T\leq t),t>0
\]
其中\(S(t)\) 为生存函数
library(survival)
Attaching package: 'survival'
The following object is masked from 'package:quantreg':
untangle.specials