我是 R 新手。我目前正在尝试基于sysid
R 包中的工具变量实现回归。我选择了这个包,因为它可以预测我的乐器。
我找到了一种合适的方法(“iv”是这里的功能)来解决我的问题。但是 R 函数没有将“预测仪器”作为其返回参数之一返回。我对那个预测变量非常感兴趣。有没有办法让这个变量作为参数?
我已经尝试创建这个函数的克隆,但是它有很多来自sysid
包的依赖函数,所以它失败了。我还尝试使用“source”命令在我的 R 代码中链接这个修改后的函数,但其余的库与我当前的脚本脱链。请为我提供任何解决方案以获取预测的仪器。源代码如下:
https ://rdrr.io/cran/sysid/src/R/iv.R 。
iv4 <- function(z,order=c(0,1,0)){
na <- order[1]; nb <- order[2]
# Steps 1-2
mod_iv <- iv(z,order)
# Step 3
w <- resid(mod_iv)
mod_ar <- ar(w,aic = F,order.max =na+nb)
Lhat <- signal::Arma(b=c(1,-mod_ar$ar),a=1)
# Step 4
x2 <- matrix(sim(mod_iv$sys,inputData(z)))
ivcompute(z,x2,order,Lhat)
}
我想要返回预测的仪器- Lhat。我也欢迎使用任何其他可以做同样的包或回归方法的建议(预测工具)。