能不能用R语言按下面编程形式将正态分布改为指数分布,画出指数分布概率密度和分布函数?
sigma=1
u=c(-2,0,2)
#计算和绘图
x=seq(-6,6,0.1)
t1=t2=list()
for(i in 1:3){
t1[[i]]=dnorm(x,u[i],sigma)
t2[[i]]=pnorm(x,u[i],sigma)
}
par(mar=c(2,2,2,1))
plot(x,t1[[1]],xlim=c(-6,6),type="l",lty=2,col=2) #概率分布图
lines(x,t1[[2]])
lines(x,t1[[3]],col=4,lty=3)
text(u,0.35,paste("u=",u,sep=""),col=c(2,1,4))
plot(x,t2[[1]],xlim=c(-6,6),type="l",lty=2,col=2) #累计分布图
lines(x,t2[[2]])
lines(x,t2[[3]],col=4,lty=3)
text(u,0.5,paste("u=",u,sep=""),col=c(2,1,4))