使用 postfile 和 post 命令

Question

使用post 命令时，出现以下错误：

post command requires expressions be bound in parenthesis

我的程序生成一个矩阵来存储每个模拟的回归系数，然后使用post 命令声明为float 并将矩阵的输出放在括号中（测试版）。

代码示例：

*Priors
set more off
global nmc=10
global l = 4  /* number of lags */ 
global cnt=150 /* number of countries */ 
set seed 10101

* Gen empty beta matrix 
matrix betas = J(153,$nmc+1,.) 

*** THIS IS WHERE MONTECARLO STARTS*** 
program bootStrapCH5, rclass 
tempname sim
postfile `sim' betas using results, replace  /* As trial I'll create only the betas matrix for now.  */ 
*postfile `sim' betas alpha_mean b1_mean b2_mean b3_mean b4_mean se_alpha se1 se2 se3 se4 using results, replace 

quietly {
forvalues i = 1/$nmc {

    * Fixed effects regression. 
    reg gdp_growth_wb L(1/4).gdp_growth_wb i.id
    matrix B1= e(b)
    mat li B1
    predict g_hat,xb
    gen e_hat= gdp_growth_wb - g_hat
    *gen flag=e(sample)


    * Generate the "wild" errors for the forecasts 
    gen eta=rnormal()
    gen e_star=e_hat*eta

        **RECURSION
    levelsof id, local(codes)

    capture noisily replace y_star= _b[_cons] + _b[L.gdp_growth_wb]*L.y_star + ///
        _b[L2.gdp_growth_wb]*L2.y_star + _b[L3.gdp_growth_wb]*L3.y_star + ///
        _b[L4.gdp_growth_wb]*L4.y_star + e_star if (id==1 & Dini4forward==1)

    forvalues cc= 2(1)150 {
        capture noisily replace y_star= _b[_cons] + _b[`cc'.id] + _b[L.gdp_growth_wb]*L.y_star + ///
        _b[L2.gdp_growth_wb]*L2.y_star + _b[L3.gdp_growth_wb]*L3.y_star + ///
        _b[L4.gdp_growth_wb]*L4.y_star + e_star if (id==`cc' & Dini4forward==1)

            }

    *Regression with new sample: y_star
    reg y_star L(1/4).y_star i.id
    matrix b= e(b)' 
    matrix betas= (betas , b)
    matrix list betas

    post `sim' float (betas)

    }
    }
    postclose `sim'
    end

    *Execute program 
    bootStrapCH5
    use results, clear
    summarize 

我也尝试了另一种方法：

post `sim'  (betas)

得到了错误：

> type mismatch
post:  above message corresponds to expression 1, variable betas

非常感谢有关如何解决此问题的任何想法。

Answer 1

我对@9​​87654321@ 不是很熟悉，但我认为一个问题可能是您尝试使用post 将kx2 矩阵插入到循环内的单个变量中。

当您使用以下命令发起 postfile 时：

postfile `sim' betas using results

您已经声明了一个包含单个变量 betas 的 Stata 数据集。

所以，不要使用

post `sim' float (betas)

你可以试试：

tempname sim
postfile `sim' float (betas1 betas2) using results, replace
forvalues i = 1/$nmc {

    * Some code. . . 
    local rows = rowsof(betas)
    forvalues i = 1/`r' {
        post `sim' (betas[`i',1]) (betas[`i',2])
    }
    * some other code. . .
}

或类似的东西来声明一个具有适当数量的变量的文件，您打算在 posting 上向数据集发送这些变量。

此外，我不确定您是否可以直接post 一个矩阵（我可能错了）。如果你不能，那么你可以在循环内嵌套一个forvalues 循环，你当前必须遍历betas 的元素并单独发布它们 - 正如我在上面的示例中所做的那样。

最后，您尝试在您的post 命令中将betas 的值转换为数据类型float。我相信存储类型需要在postfile 命令中声明（但同样，我可能错了）。您引用的第一个错误（括号中的表达式）是在post 命令中包含float 的直接结果。

底线 - 我怀疑第一个错误是由于在您尝试 post 数据时声明数据类型，第二个错误（类型不匹配）是尝试将 kx2 矩阵插入变量的结果。尝试（错误地）从矩阵创建数据时，请参阅下面的类型不匹配示例：

clear *
mat a = (1\2)
set obs 2
gen x = a

虽然我承认我会期望错误更类似于此：

mat a = (1\2)
set obs 2
gen x = a*2

返回此上下文中不允许的矩阵的矩阵运算符

还可以查看 svmat 从矩阵创建数据。