我正在使用该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")