从 ucm 类对象中提取 p 值答案

【问题标题】：Extract p-value from ucm class object从 ucm 类对象中提取 p 值
【发布时间】：2018-01-15 13:31:46
【问题描述】：

我正在使用rucm 包，我有一个名为multivar_ucm 的ucm 类对象。当我在 r 控制台中输入 multivar_ucm 并按回车键时，我会得到所有外生变量的 beta 估计值，以及在 r 控制台中打印的标准误差和 p 值。

我可以通过使用as.data.frame(multivar_ucm$est) 获取 beta 系数并将其写入 csv 文件。也在寻找一种类似的方法来导出相应的 p 值。

【问题讨论】：

欢迎来到 StackOverflow！请阅读有关how to ask a good question 的信息以及如何提供reproducible example。这将使其他人更容易帮助您。

【解决方案1】：

不幸的是，包没有暴露对象的 p 值（这有点疏忽）。 print.ucm 方法每次都会重新计算该值。如果要获取对象的p-value，you’ll need to do the same:

t_value = x$est / attr(x, "se")
p_value = 2 * pt(-abs(t_value), df = attr(x, "df"))

（代码稍微清理一下。）

【讨论】：