2015年10月21日

  ＃改天写注释

set.seed(1)

x1<-rt(10,4,ncp = 0)

x2<-rt(20,4,ncp = 0)

x3<-rt(50,4,ncp = 0)

x4<-rt(100,4,ncp = 0)

x5<- seq(-8,8,length.out = 1000)

opar <- par(no.readonly = T)

par(mfrow=c(2,2))

hist(x1,freq = F,breaks=10,main="n=10,t(4)",xlab = "x")

lines(density(x1),col="red",lty=5,lwd=2)

hist(x2,freq = F,breaks=20,main="n=20,t(4)",xlab = "x")

lines(density(x2),col="red",lty=5,lwd=2)

hist(x3,freq = F,breaks=50,main="n=50,t(4)",xlab = "x")

lines(density(x3),col="red",lty=5,lwd=2)

hist(x4,freq = F,breaks=100,main="n=100,t(4)",xlab = "x")

lines(density(x4),col="red",lty=5,lwd=2)