2000论坛币求两个不同总体检验的SAS程序

我想请问一下高手，请问您是否可以编一个sas程序，满足以下要求：1.生成两个正态分布的总体A与B；2.从正态分布的总体B中，随机抽出N个含有相同样本量的样本，然后一一检验每个样本与正态分布总体A之间的差别，看每个样本是否来自总体A，这样可以得到N个P值，再把这N个P值一一列出来。谢谢高手！

1. /*产生总体A*/
   
2. data a;
   
3.         do i=1 to 10000;
   
4.                 x=0+1*rannor(0);
   
5.                 group=0;
   
6.                 output;
   
7.         end;
   
8.         drop i;
   
9. run;
   
10. 
    
11. /*产生总体B*/
    
12. data b;
    
13.         do i=1 to 10000;
    
14.                 x=0+1*rannor(0);
    
15.                 output;
    
16.         end;
    
17.         drop i;
    
18. run;
    
19. 
    
20. %macro test(n1);
    
21. %do i=1 %to &n1.;
    
22. /*从总体B中随机抽取n1个样本*/
    
23.         proc surveyselect data=b out=b&i.
    
24.         method = srs
    
25.                n=53
    
26.         seed =0;
    
27. run;
    
28. 
    
29. data b&i;
    
30.         set b&i;
    
31.         group=&i.;
    
32. run;
    
33. 
    
34. data test&i.;
    
35.         set a b&i.;
    
36. run;
    
37. 
    
38. /*将每个样本与总体A进行t检验*/
    
39. ods output ttests=p&i.;
    
40.         proc ttest data=test&i.;
    
41.                 class group;
    
42.                 var  x;
    
43.         quit;
    
44. ods output close;
    
45. 
    
46. /*将每个样本与总体A进行t检验的p值取出来*/
    
47. data p&i.;
    
48.         set p&i.;
    
49.         if Variances="Equal" then call symput("p&i.",Probt);
    
50. run;
    
51. %end;
    
52. 
    
53. /*将每个样本与总体t检验的p值输出到同一个数据集中*/
    
54. data p;
    
55.                 retain count;
    
56.         %do j=1 %to &n1.;
    
57.                  p_valute=&&p&j.;
    
58.                                  if p_valute>0.05 then count+1;/*count表示p值大于0.05的样本个数*/
    
59.                                  output;
    
60.         %end;
    
61. run;
    
62. 
    
63. %mend;
    
64. %test(10);

复制代码

不好意思，帮不上忙！

谢谢关注！！

你上个悬赏还没有给人家奖励，怎么还能够期待人家会继续帮你呢？

同楼上

谢谢楼上两位提醒，请大家看我原来悬赏贴的解释及结果！！我原来的悬赏贴是“200论坛币请教正态分布问题
本文来自: 人大经济论坛 SAS专版 版，详细出处参考： https://bbs.pinggu.org/forum.php?mod=viewthread&tid=1295597&page=1&from^^uid=2090664”

感觉楼主的描述有不妥之处，样本与A总体的差异是指什么样是什么？1.如果是指均数之间的差异，那么计算抽样误差时是利用A总体标准差还是每个样本的标准差啊？2。如果是想检验每个样本是否来自A总体，那么频数表组段的取值又因人而异了。

谢谢kutuomonk！我说得差异是指均数之间的差异；2.具体说是B1均数、B2均数、B3均数、B4均数（B1、B2、B3、B4······分别代表从总体B中抽出的各个样本）······分别与总体A均数的差别，这样可以得出P1、P2、P3、P4等（可以用样本均数与总体均数比较的t检验），谢谢高手指点！

请问kutuomonk，在“如果是想检验每个样本是否来自A总体，那么频数表组段的取值又因人而异了。”“频数表组段的取值”是什么意思！谢谢！！