Stata：提取p值并将其保存在列表中

Question

这可能是一个琐碎的问题，但是作为R用户来到Stata以来，我到目前为止未能找到正确的Google术语来找到答案。 我要执行以下步骤：

• 做一堆测试（例如，lrtest导致一个foreach循环）
• 从每个测试中提取p值并将其保存在某种列表中
• 列出我可以做的进一步操作（例如执行多重比较校正）

所以我想知道如何从命令结果中提取p值（或类似值），以及如何将它们保存到我可以使用的类似矢量的对象中。 这是一些执行类似操作的R代码：

myData <- data.frame(a=rnorm(10), b=rnorm(10), c=rnorm(10)) ## generate some data
pValue <- c()
for (variableName in c("b", "c")) {
    myModel <- lm(as.formula(paste("a ~", variableName)), data=myData) ## fit model
    pValue <- c(pValue, coef(summary(myModel))[2, "Pr(>|t|)"]) ## extract p-value and save in vector
}
pValue * 2 ## do amazing multiple comparison correction

在我看来，Stata似乎比R具有更少的“编程”心态。如果您对可以编程的R用户有任何一般的Stata文献建议，也将不胜感激。

Answer 1

这是一种将p值保存在矩阵中，然后可以操纵矩阵的方法，也许可以使用Mata或Stata中的标准矩阵操纵。

matrix storeMyP = J(2, 1, .)  //create empty matrix with 2 (as many variables as we are looping over) rows, 1 column
matrix list storeMyP  //look at the matrix

loc n = 0 //count the iterations

foreach variableName of varlist b c {

  loc n = `n' + 1  //each iteration, adjust the count

  reg a `variableName'
  test `variableName' //this does an F-test, but for one variable it's equivalent to a t-test (check: -help test- there is lots this can do
  matrix storeMyP[`n', 1] = `r(p)'  //save the p-value in the matrix 

}
matrix list storeMyP  //look at your p-values 
matrix storeMyP_2 = 2*storeMyP //replicating your example above 

发生这样的情况是，Stata在估算和测试命令之后会自动存储某些数量。 当帮助文件说此命令将以下值存储在r（）中时，您将其用单引号引起来。 使用svmat storeMyP将矩阵列转换为变量也可能很有趣，或者请参阅help svmat以获取更多信息。