exppar<-2 #指数分布的参数,指数分布的期望为exppar，方差为exppar^2
randexp<-rexp(100,1/exppar) #随机生成100个服从指数分布的样本
sortrandexp<-sort(randexp) #将样本按从小到大排列
dfyaxis<- 1-exp(-sortrandexp/exppar) #求每个样本对应的分布函数值(0,1)
sfyaxis<- 1-dfyaxis #求每个样本对应的生存函数值S(x)=1-F(x)
#服从指数分布的原数据的对数经验生存函数
plot(log(sortrandexp),sfyaxis,type="p",xlab="随机样本值的自然对数",ylab="经验生存函数值")
只有100个点，最后的图有点散
扫码加好友，拉您进群
全部版块
我的主页
收藏
