stata 出现weights not allowed

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我想求元素个数为10的列向量x的累计的算术平均值：set obs 10
gen m=x
forvalues i=1/10{ replace m[i+1]=mean(x[1]-x[i+1])
}
为什么总是显示weights not allowed？

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：weights allowed Weight Allow weigh

两个错误：
1. mean()是egen的函数，不能用replace.
2. m[i+1]不能用，如果确实要用应该为

1. replace m = ... in `i'

复制代码

楼主的问题，可以这样实现：

1. . clear
   
2. 
   
3. . input x
   
4. 
   
5.              x
   
6.   1. 1
   
7.   2. 3
   
8.   3. 4
   
9.   4. 2
   
10.   5. 7
    
11.   6. 9
    
12.   7. 10
    
13.   8. 3
    
14.   9. 6
    
15. 10. 8
    
16. 11. end
    
17. 
    
18. .
    
19. . list
    
20. 
    
21.      +----+
    
22.      |  x |
    
23.      |----|
    
24.   1. |  1 |
    
25.   2. |  3 |
    
26.   3. |  4 |
    
27.   4. |  2 |
    
28.   5. |  7 |
    
29.      |----|
    
30.   6. |  9 |
    
31.   7. | 10 |
    
32.   8. |  3 |
    
33.   9. |  6 |
    
34. 10. |  8 |
    
35.      +----+
    
36. 
    
37. .
    
38. . gen m = .
    
39. (10 missing values generated)
    
40. 
    
41. . forvalues i = 1/10 {
    
42.   2.         egen b = mean(x) if _n <= `i'
    
43.   3.         replace m = b in  `i'
    
44.   4.         drop b
    
45.   5. }
    
46. (9 missing values generated)
    
47. (1 real change made)
    
48. (8 missing values generated)
    
49. (1 real change made)
    
50. (7 missing values generated)
    
51. (1 real change made)
    
52. (6 missing values generated)
    
53. (1 real change made)
    
54. (5 missing values generated)
    
55. (1 real change made)
    
56. (4 missing values generated)
    
57. (1 real change made)
    
58. (3 missing values generated)
    
59. (1 real change made)
    
60. (2 missing values generated)
    
61. (1 real change made)
    
62. (1 missing value generated)
    
63. (1 real change made)
    
64. (1 real change made)
    
65. 
    
66. . list
    
67. 
    
68.      +---------------+
    
69.      |  x          m |
    
70.      |---------------|
    
71.   1. |  1          1 |
    
72.   2. |  3          2 |
    
73.   3. |  4   2.666667 |
    
74.   4. |  2        2.5 |
    
75.   5. |  7        3.4 |
    
76.      |---------------|
    
77.   6. |  9   4.333333 |
    
78.   7. | 10   5.142857 |
    
79.   8. |  3      4.875 |
    
80.   9. |  6          5 |
    
81. 10. |  8        5.3 |
    
82.      +---------------+

复制代码

andruw 发表于 2015-7-3 15:46
两个错误：
1. mean()是egen的函数，不能用replace.
2. m不能用，如果确实要用应该为楼主的问题，可以这样 ...

谢谢~~但是为什么我显示的是这个啊。。求解答~

andystata 发表于 2015-7-4 10:44
谢谢~~但是为什么我显示的是这个啊。。求解答~

1. if_n<=`i'

复制代码

应为

1. if _n <= `i'

复制代码

注意 if 与_n之间有空格。

andruw 发表于 2015-7-4 16:04
应为注意 if 与_n之间有空格。

嗯嗯 谢谢 问题是解决了  但是我的数据有一千万个 这样算好慢。。请问还有办法吗？

andruw 发表于 2015-7-4 16:04
应为注意 if 与_n之间有空格。

在你的启示下，我是这样写的
gen m=.
forvalues i=1/10000000{
  gen b=v/`i'
   replace m=b in `i'
   drop b
}
但是也还是好慢。。

andystata 发表于 2015-7-4 16:12
在你的启示下，我是这样写的
gen m=.
forvalues i=1/10000000{

按照你的代码：

1. gen m = .
   
2. forvalues i = 1/10000000 {
   
3.    scalar b = v/`i'
   
4.    replace m = b in `i'
   
5. }

复制代码

生成scalar比生成variable稍微快一点。

andystata 发表于 2015-7-4 16:12
在你的启示下，我是这样写的
gen m=.
forvalues i=1/10000000{

gen v=sum(x)

gen meanx=sum(x)/_n

andystata 发表于 2015-7-4 16:12
在你的启示下，我是这样写的
gen m=.
forvalues i=1/10000000{

还有一个加速的办法，就是用mata写个小程序。