2015年10月15日

#生成随机数编号为1
set.seed(1)
x <- seq(-5,5,length.out = 1000)     #生成序列,序列长度1000,可以理解为序列的精细度       
y1 <- dnorm(x,0,1)                   #生成N(O,1)的PDF

#?xaxs"i";lwd是线宽,递增
plot(x,y1,col="red",xlim=c(-5,5),ylim=c(0,0.45),type='l',
     xaxs="i", yaxs="i",ylab='Density',xlab='x',
     main="The Normal&t Density Distribution",lwd=5) 

#生成t分布PDF;lty是虚线密度,1默认为实线
lines(x,dt(x,5,0),col="green",lty=5,lwd=2)        
lines(x,dt(x,10,0),col="blue",lty=5,lwd=2)
lines(x,dt(x,20,0),col="grey",lty=5,lwd=2)
lines(x,dt(x,100,0),col="black",lty=5,lwd=2)
lines(x,dt(x,2,0),col="orange",lty=5,lwd=2)

#图例
legend("topright",legend=paste("",c("N(0,1)","t(2)","t(5)","t(10)","t(20)","t(100)")),
       col=c("red","orange","green","blue","grey","black"),
       lty = c(1,5,5,5,5,5),lwd = c(5,2,2,2,2,2),cex=0.48,bty = "0") 2015年10月15日 - paul_sweet@126 - MOMO


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值