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"