benz1985 发表于 2012-10-4 23:13 
交互项控制约束,
* add constraints for pz interactions
. forvalues h=1(1)$ndem {
你试下显示出所有的constraint #,看看有没有超出范围
 forvalues h=1  (1) $ndem {
        di "h=`h'"
          forvalues j= 1 (1) $neq {
                  di "j=`j'"
                   local jplus1 = `j'+1
                   forvalues k=`jplus1' (1) $neq {
                                   di "k=`k'"
                           constraint `j'`k'`h' [s`j']np`k'z`h'=[s`k']np`j'z`h'  
                                                   di "constraint # = `j'`k'`h' [s`j']np`k'z`h'"
                          global conlist "$conlist `j'`k'`h'"
                                                  di
                                                  di
                   }
           }
   }
这个能过
 global ndem "5" 
 global neq "5"
这个就不行
 global ndem "10"
 global neq "5"