作者xung86 (PUSH!)
看板Statistics
标题[问题] 统软问题
时间Thu Nov 2 01:14:58 2006
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
1F:推 WANG3213:学弟,抓到你了,作业要自己作喔~~~ 11/02 02:31
2F:推 WANG3213:来去跟历史老师告状XD 11/02 02:37