我使用systemfit()
函数进行看似无关的回归(SUR)并收到以下错误:
Error in solve(sigma, tol = solvetol) :
Lapack routine dsptrf returned error code 1
然而,我无法找到错误表明正在进行的有意义的解释。下面是一些模拟代码,用于显示我正在使用的功能(模拟代码不会产生错误)。感谢您对此错误的想法。
y <- sample(seq(1:4), 100, replace = TRUE)
x1 <- sample(seq(0:1), 100, replace = TRUE) -1
x2 <- sample(seq(0:1), 100, replace = TRUE) - 1
x3 <- sample(seq(1:4), 100, replace = TRUE)
frame <- as.data.frame(cbind(y,x1,x2, x3))
mod_1 <- y ~ x1 + x3 + x1:x3
mod_2 <- y ~ x2 + x3 + x2:x3
output <- systemfit(list(mod_1, mod_2), data = frame, method = "SUR")