从ucm类对象中提取p值

时间:2018-01-15 13:31:47

标签: r p-value

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

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

1 个答案:

答案 0 :(得分:1)

不幸的是,包不会暴露对象的 p - 值(这有点疏忽)。 print.ucm方法每次都会重新计算该值。如果你想得到对象的 p - 值,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"))

(代码稍作清理。)