我正在使用rucm
包,我有一个名为multivar_ucm的ucm类对象。当我在r控制台中输入multivar_ucm并点击回车时,我得到所有外生变量的beta估计值以及r控制台中打印的标准误差和p值。
我可以使用as.data.frame(multivar_ucm$est)
获取beta系数并将其写入csv文件中。寻找类似的方法来输出相应的p值。
答案 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"))
(代码稍作清理。)