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")