用nlsur做aids但是一直报错r198，检查不出错误，求大佬们帮忙看看

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. program nlsuraids
   
2.            version 15
   
3.                    syntax varlist(min=16 max=16) if, at(name)
   
4.            tokenize `varlist'
   
5.            args w1 w2 w3 w4 w5 w6 w7 lnp1 lnp2 lnp3 lnp4 lnp5 lnp6 lnp7 lnp8 lnm
   
6.            tempname a1 a2 a3 a4 a5 a6 a7 a8
   
7.                    scalar `a1' = `at'[1,1]
   
8.            scalar `a2' = `at'[1,2]
   
9.            scalar `a3' = `at'[1,3]
   
10.                    scalar `a4' = `at'[1,4]
    
11.                    scalar `a5' = `at'[1,5]
    
12.                    scalar `a6' = `at'[1,6]
    
13.                    scalar `a7' = `at'[1,7]  
    
14.                    scalar `a8' =  1 - `a1' - `a2' - `a3'- `a4'- `a5'- `a6'- `a7'
    
15.                    tempname b1 b2 b3 b4 b5 b6 b7
    
16.                    scalar `b1' = `at'[1,8]
    
17.            scalar `b2' = `at'[1,9]
    
18.            scalar `b3' = `at'[1,10]
    
19.                    scalar `b4' = `at'[1,11]
    
20.                    scalar `b5' = `at'[1,12]
    
21.                    scalar `b6' = `at'[1,13]
    
22.                    scalar `b7' = `at'[1,14]
    
23.                    tempname g11 g12 g13 g14 g15 g16 g17 g18
    
24.            tempname g21 g22 g23 g24 g25 g26 g27 g28
    
25.            tempname g31 g32 g33 g34 g35 g36 g37 g38
    
26.            tempname g41 g42 g43 g44 g45 g46 g47 g48
    
27.                    tempname g51 g52 g53 g54 g55 g56 g57 g58
    
28.                    tempname g61 g62 g63 g64 g65 g66 g67 g68
    
29.                    tempname g71 g72 g73 g74 g75 g76 g77 g78
    
30.                    tempname g81 g82 g83 g84 g85 g86 g87 g88
    
31.                    scalar `g11' = `at'[1,15]
    
32.            scalar `g12' = `at'[1,16]
    
33.            scalar `g13' = `at'[1,17]
    
34.            scalar `g14' = `at'[1,18]
    
35.                    scalar `g15' = `at'[1,19]
    
36.                    scalar `g16' = `at'[1,20]
    
37.                    scalar `g17' = `at'[1,21]
    
38.                    scalar `g18' = -`g11'-`g12'-`g13'-`g14'-`g15'-`g16'-`g17'
    
39.            scalar `g21' = `g12'
    
40.            scalar `g22' = `at'[1,22]
    
41.            scalar `g23' = `at'[1,23]
    
42.            scalar `g24' = `at'[1,24]
    
43.                    scalar `g25' = `at'[1,25]
    
44.                    scalar `g26' = `at'[1,26]
    
45.                    scalar `g27' = `at'[1,27]
    
46.                    scalar `g28' = -`g21'-`g22'-`g23'-`g24'-`g25'-`g26'-`g27'
    
47.            scalar `g31' = `g13'
    
48.            scalar `g32' = `g23'
    
49.            scalar `g33' = `at'[1,28]
    
50.            scalar `g34' = `at'[1,29]
    
51.                    scalar `g35' = `at'[1,30]
    
52.                    scalar `g36' = `at'[1,31]
    
53.                    scalar `g37' = `at'[1,32]
    
54.                    scalar `g38' = -`g31'-`g32'-`g33'-`g34'-`g35'-`g36'-`g37'
    
55.            scalar `g41' = `g14'
    
56.            scalar `g42' = `g24'
    
57.            scalar `g43' = `g34'
    
58.            scalar `g44' = `at'[1,33]
    
59.                    scalar `g45' = `at'[1,34]
    
60.                    scalar `g46' = `at'[1,35]
    
61.                    scalar `g47' = `at'[1,36]
    
62.                    scalar `g48' = -`g41'-`g42'-`g43'-`g44'-`g45'-`g46'-`g47'
    
63.                    scalar `g51' = `g15'
    
64.            scalar `g52' = `g25'
    
65.            scalar `g53' = `g35'
    
66.            scalar `g54' = `g45'
    
67.                    scalar `g55' = `at'[1,37]
    
68.                    scalar `g56' = `at'[1,38]
    
69.                    scalar `g57' = `at'[1,39]
    
70.                    scalar `g58' = -`g51'-`g52'-`g53'-`g54'-`g55'-`g56'-`g57'   
    
71.                    scalar `g61' = `g16'
    
72.            scalar `g62' = `g26'
    
73.            scalar `g63' = `g36'
    
74.            scalar `g64' = `g46'
    
75.                    scalar `g65' = `g56'
    
76.                    scalar `g66' = `at'[1,40]
    
77.                    scalar `g67' = `at'[1,41]
    
78.                    scalar `g68' = -`g61'-`g62'-`g63'-`g64'-`g65'-`g66'-`g67'
    
79.                    scalar `g71' = `g17'
    
80.            scalar `g72' = `g27'
    
81.            scalar `g73' = `g37'
    
82.            scalar `g74' = `g47'
    
83.                    scalar `g75' = `g57'
    
84.                    scalar `g76' = `g67'
    
85.                    scalar `g77' = `at'[1,42]
    
86.                    scalar `g78' = -`g71'-`g72'-`g73'-`g74'-`g75'-`g76'-`g77'
    
87.                    scalar `g81' = `g18'
    
88.            scalar `g82' = `g28'
    
89.            scalar `g83' = `g38'
    
90.            scalar `g84' = `g48'
    
91.                    scalar `g85' = `g58'
    
92.                    scalar `g86' = `g68'
    
93.                    scalar `g87' = `g78'
    
94.                    scalar `g88' = -`g81'-`g82'-`g83'-`g84'-`g85'-`g86'-`g87'
    
95.         quietly {
    
96.                     tempvar lnpindex
    
97.                     gen double `lnpindex' = 5 + `a1'*`lnp1' + `a2'*`lnp2' + `a3'*`lnp3' + `a4'*`lnp4'+ `a5'*`lnp5' + `a6'*`lnp6' + `a7'*`lnp7' + `a8'*`lnp8'
    
98.            forvalues i = 1/8 {
    
99.                       forvalues j = 1/8 {
    
100.                       replace `lnpindex' = `lnpindex' + 0.5*`g`i'`j''*`lnp`i''*`lnp`j''
     
101.                                         }
     
102.                              }
     
103.            replace `w1' = `a1' + `g11'*`lnp1' + `g12'*`lnp2' + `g13'*`lnp3' + `g14'*`lnp4' + `g15'*`lnp5' + `g16'*`lnp6' + `g17'*`lnp7' + `g18'*`lnp8' + `b1'*(`lnm' - `lnpindex')
     
104.            replace `w2' = `a2' + `g21'*`lnp1' + `g22'*`lnp2' + `g23'*`lnp3' + `g24'*`lnp4' + `g25'*`lnp5' + `g26'*`lnp6' + `g27'*`lnp7' + `g28'*`lnp8' + `b2'*(`lnm' - `lnpindex')
     
105.            replace `w3' = `a3' + `g31'*`lnp1' + `g32'*`lnp2' + `g33'*`lnp3' + `g34'*`lnp4' + `g35'*`lnp5' + `g36'*`lnp6' + `g37'*`lnp7' + `g38'*`lnp8' + `b3'*(`lnm' - `lnpindex')
     
106.                    replace `w4' = `a4' + `g41'*`lnp1' + `g42'*`lnp2' + `g43'*`lnp3' + `g44'*`lnp4' + `g45'*`lnp5' + `g46'*`lnp6' + `g47'*`lnp7' + `g48'*`lnp8' + `b4'*(`lnm' - `lnpindex')
     
107.                    replace `w5' = `a5' + `g51'*`lnp1' + `g52'*`lnp2' + `g53'*`lnp3' + `g54'*`lnp4' + `g55'*`lnp5' + `g56'*`lnp6' + `g57'*`lnp7' + `g58'*`lnp8' + `b5'*(`lnm' - `lnpindex')
     
108.                    replace `w6' = `a6' + `g61'*`lnp1' + `g62'*`lnp2' + `g63'*`lnp3' + `g64'*`lnp4' + `g65'*`ln65' + `g66'*`lnp6' + `g67'*`lnp7' + `g68'*`lnp8' + `b6'*(`lnm' - `lnpindex')
     
109.                    replace `w7' = `a7' + `g71'*`lnp1' + `g72'*`lnp2' + `g73'*`lnp3' + `g74'*`lnp4' + `g75'*`lnp5' + `g76'*`lnp6' + `g77'*`lnp7' + `g78'*`lnp8' + `b7'*(`lnm' - `lnpindex')
     
110.                    }
     
111. end
     
112. nlsur aids @ w1 w2 w3 w4 w5 w6 w7 lnp1 lnp2 lnp3 lnp4 lnp5 lnp6 lnp7 lnp8 lnexp, parameters(a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 g11 g12 g13 g14 g15 g16 g17 g22 g23 g24 g25 g26 g27 g33 g34 g35 g36 g37 g44 g45 g46 g47 g55 g56 g57 g66 g67 g77)neq(7) ifgnls
     
113.           

复制代码

就是上面这些代码，用aids做8大类消费，一致报错r198，求各位帮忙看看，多谢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：NLSUR AIDS SUR nls IDS aids模型 建模 报错

顶一下

请问您解决了么？也一直报错198

我也一直报错198，反复检查了，不知道问题在哪。请问楼主解决了嘛？