求助stata大神：做广义倾向得分匹配出现valid syntax

ml19930812

1915

收藏 2017-03-12

我用doseresponse，在stata中输入：gen cut=0.01 if subsidy<=0.01
replace cut=0.02 if subsidy>0.01 & subsidy<=0.02
replace cut=0.5 if subsidy>0.02
tab cut
matrix define tp=(10\20\30\40\50\60\70\80\90\100)
doseresponse size prof drd age export gov,outcome(tech) t(subsidy) gpscore(pscore) predict(hat_treat) sigma(sd) cutpoints(cut) index(p50) nq_gps(5) dose_response(dose_response) t_transf(ln) tpoints(tp) delta(1) reg_type_t(quadratic) reg_type_gps(quadratic) interaction(1) bootstrap(yes) boot_reps(100)
filename("output") analysis(yes) graph("graph_output") detail
可是输入的命令时候呈现结果错误 invalid syntax：

请问是什么原因？谢谢各位大神！

扫码加我拉你入群

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

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

全部回复

ml19930812

2017-3-12 23:39:36

输出结果如下：

********************************************
ESTIMATE OF THE GENERALIZED PROPENSITY SCORE
********************************************

Generalized Propensity Score
By agreement we assume that the logarithm of 0 is 0

******************************************************
Algorithm to estimate the generalized propensity score
******************************************************

Estimation of the propensity score

The log transformation of the treatment variable subsidy is used

                           T
-------------------------------------------------------------
   Percentiles    Smallest
1% -8.316958    -8.917275
5% -6.945942    -8.721264
10%    -6.03387    -8.645988       Obs                672
25% -5.079773    -8.529114    Sum of Wgt.       672

50% -4.222722                         Mean       -4.167873
                                 Largest    Std. Dev.    1.618843
75% -3.385871             0
90% -2.359686             0    Variance    2.620653
95%          0                0    Skewness    .3090964
99%          0                0    Kurtosis    3.966695

invalid syntax
r(198);

.
. filename("output") analysis(yes) graph("graph_output") detail
command filename is unrecognized
r(199);