関数を使用して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")