x=30;h=10;p=0
> for(p in 0:5){
+ n=h+h*p
+ f=function(x,n,a,mean,sd){
+ i=1;j=0;k=0;
+ while(i<=n){
+ w=shapiro.test(rnorm(x,mean,sd))
+ u=w$p.value
+ if(u>a) j=j+1
+ else if(u<a) k=k+1
+ i=i+1
+ }
+ tp=j/n
+ Bp=c(n,tp)
+ print(Bp)
+ }
+ f(x,n,0.05,0,1)
+ p=p+1
+ }
[1] 10 1
[1] 20 1
[1] 30.0000000 0.9333333
[1] 40.000 0.975
[1] 50.00 0.96
[1] 60.0000000 0.9666667
结果以出来后,想在一张图上画散点图,怎么编程


雷达卡






京公网安备 11010802022788号







