[問題] 統軟問題
pois<- function(n,rep)
{
i<- 1
z<- vector(mode="numeric",length=rep)
y<- vector(mode="numeric",length=rep)
while(i <=rep)
{lambda<-20
u<-rpois(n,lambda)
mu<-mean(u)
z[i]<-sqrt(n/lambda)*(mu-lambda)
i<-i+1
}
x<-seq(-3,3,length=rep)
plot(x,pnorm(x),type="l",xlab="x",ylab="Probability",col=4,lty=1)
for (j in (1:rep))
{
y[j]<-sum(z<=x[j])/rep
}
lines(x,y,lty=2,col=2)
ch<-c("Normal CDF","Approx. CDF")
legend(-3,0.9,ch,col=c(4,2),lty=1:2)
print(shapiro.test(z))
}
請問一下分配要改成bonomial要怎麼改呢?
改了好幾次圖都不是正確的!!
有大大可以幫我解決嗎?
--
http://www.wretch.cc/album/xung86
我要廣末住我家
因為妳就是她
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.171.53.152
推
11/02 02:31, , 1F
11/02 02:31, 1F
推
11/02 02:37, , 2F
11/02 02:37, 2F