如何将所作出的五次的数组形成的五张散点图做到同一个坐标系中

1. x=30;h=10;p=1;
   
2. for(s in 1:5){
   
3. n=h*s
   
4. f=function(x,n,a,mean,sd){
   
5. i=1;j=0;
   
6. while(i<=n){
   
7.   w=shapiro.test(rnorm(x,mean,sd))
   
8.   u=w$p.value
   
9.   if(u>a) j=j+1
   
10.    i=i+1
    
11. }
    
12. y=j/n
    
13. dev.new()
    
14. plot(y,n)
    
15. }
    
16. f(x,n,0.05,0,1)
    
17. }

复制代码

作出的五次的数组形成的五张散点图做到同一个坐标系中，或者是将每次得到的数组存储到一个矩阵或是数组中，然后做出整合后的散点图。

1. x=30;h=10;p=1;
   
2. op = par(mfrow = c(3, 2))
   
3. for(s in 1:5){
   
4.   n=h*s
   
5.   f=function(x,n,a,mean,sd){
   
6.     i=1;j=0;
   
7.     while(i<=n){
   
8.       w=shapiro.test(rnorm(x,mean,sd))
   
9.       u=w$p.value
   
10.       if(u>a)
    
11.         j=j+1
    
12.       i=i+1
    
13.     }
    
14.     y=j/n
    
15.     plot(y,n)
    
16.   }
    
17.   f(x,n,0.05,0,1)
    
18. }
    
19. par(op)

复制代码

1. x=30;h=10;p=1;mean=0;sd=1;a=0.01
   
2. y=rep(0,5)
   
3. n=rep(0,5)
   
4. for(s in 1:5){
   
5.   n[s]=h*s
   
6.   i=1;j=0;
   
7.   while(i<=n[s]){
   
8.     w=shapiro.test(rnorm(x,mean,sd))
   
9.     u=w$p.value
   
10.     if(u>a)
    
11.       j=j+1
    
12.     i=i+1
    
13.   }
    
14.   y[s]=j/n[s]
    
15. }
    
16.   plot(y,n)

复制代码

zhou1_20 发表于 2016-10-19 10:35

谢谢，初学，卡了几天了