如何进行分组

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有两列数a 和 b，我想根据b来从小到大分成5组，然后将a按照这5组来分析，比如a在每一组的mean一类的统计量。

现在问题是我想将b分成很多组，然后通过centile得出每组的数再按数进行keep太繁琐。求code

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：从小到大 cent code Tile mean 如何

1. egen var_group=cut(b), group(10)
   
2. egen meanv=mean(a),by(var_group)

复制代码

ywh19860616 发表于 2013-12-20 08:40

这儿mean给出的是估计平均数，默认的置信水平为95%，所以也就是说，一组数用mean求得的均值并不是其算数平均数。

瘠北幼农 发表于 2013-12-20 09:56
这儿mean给出的是估计平均数，默认的置信水平为95%，所以也就是说，一组数用mean求得的均值并不是其算数平 ...

谢谢提醒，我试了两组数据，结果是正确的。

1. . clear
   
2. 
   
3. . input num1 num2
   
4. 
   
5.            num1        num2
   
6.   1. 1 1
   
7.   2. 2 2
   
8.   3. 1 1
   
9.   4. 3 3
   
10.   5. 2 2
    
11.   6. 1 1
    
12.   7. 4 11
    
13.   8. end
    
14. 
    
15. . egen id1=mean(num1)
    
16. 
    
17. . egen id2=mean(num2)
    
18. 
    
19. . list id1 id2
    
20. 
    
21.      +-----------+
    
22.      | id1   id2 |
    
23.      |-----------|
    
24.   1. |   2     3 |
    
25.   2. |   2     3 |
    
26.   3. |   2     3 |
    
27.   4. |   2     3 |
    
28.   5. |   2     3 |
    
29.      |-----------|
    
30.   6. |   2     3 |
    
31.   7. |   2     3 |
    
32.      +-----------+
    
33. 
    
34. .

复制代码

1. . clear all
   
2. 
   
3. . set obs 200
   
4. obs was 0, now 200
   
5. 
   
6. . gen num1=uniform()
   
7. 
   
8. . gen num2=uniform()
   
9. 
   
10. . egen id1=mean(num1)
    
11. 
    
12. . egen id2=mean(num2)
    
13. 
    
14. . list id1 id2 in 1/5
    
15. 
    
16.      +----------------------+
    
17.      |       id1        id2 |
    
18.      |----------------------|
    
19.   1. | .49207754   .4918756 |
    
20.   2. | .49207754   .4918756 |
    
21.   3. | .49207754   .4918756 |
    
22.   4. | .49207754   .4918756 |
    
23.   5. | .49207754   .4918756 |
    
24.      +----------------------+
    
25. 
    
26. . summarize num1 num2
    
27. 
    
28.     Variable |       Obs        Mean    Std. Dev.       Min        Max
    
29. -------------+--------------------------------------------------------
    
30.         num1 |       200    .4920775    .2874816   .0057233    .984407
    
31.         num2 |       200    .4918756    .2802483   .0007097   .9931726
    
32. 
    
33. .

复制代码

ywh19860616 发表于 2013-12-20 08:40

非常感谢~

ywh19860616 发表于 2013-12-20 10:11
谢谢提醒，我试了两组数据，结果是正确的。

那为什么我上次算出来错的呢？导致我只能换方法了！

瘠北幼农 发表于 2013-12-20 15:20
那为什么我上次算出来错的呢？导致我只能换方法了！

嗯，谢谢提醒，以前没有注意过这个问题，一直默认他就是几何平均。

ywh19860616 发表于 2013-12-21 09:34
嗯，谢谢提醒，以前没有注意过这个问题，一直默认他就是几何平均。

给出的是估计的平均数吧？那算数平均数该怎么算呢？

瘠北幼农 发表于 2013-12-21 10:39
给出的是估计的平均数吧？那算数平均数该怎么算呢？

1. . set obs 200
   
2. obs was 0, now 200
   
3. 
   
4. . gen var1=uniform()
   
5. 
   
6. . ameans var1
   
7. 
   
8.     Variable |    Type        Obs        Mean       [95% Conf. Interval]
   
9. -------------+----------------------------------------------------------
   
10.         var1 | Arithmetic     200    .4920775        .4519916   .5321635
    
11.              |  Geometric     200    .3624659        .3162113   .4154865
    
12.              |   Harmonic     200    .1550841        .1098098   .2638821
    
13. ------------------------------------------------------------------------
    
14. 
    
15. . egen aa=mean(var1)
    
16. 
    
17. . list aa in 1/5
    
18. 
    
19.      +-----------+
    
20.      |        aa |
    
21.      |-----------|
    
22.   1. | .49207754 |
    
23.   2. | .49207754 |
    
24.   3. | .49207754 |
    
25.   4. | .49207754 |
    
26.   5. | .49207754 |
    
27.      +-----------+
    
28. 
    
29. . summarize aa
    
30. 
    
31.     Variable |       Obs        Mean    Std. Dev.       Min        Max
    
32. -------------+--------------------------------------------------------
    
33.           aa |       200    .4920775           0   .4920775   .4920775
    
34. 
    
35. .

复制代码