# NOT RUN {
set.seed(27)
library(systemfit)
df <- data.frame(
X = rnorm(100),
Y = rnorm(100),
Z = rnorm(100),
W = rnorm(100)
)
fit <- systemfit(formula = list(Y ~ Z, W ~ X), data = df, method = "SUR")
tidy(fit)
tidy(fit, conf.int = TRUE)
# }
Run the code above in your browser using DataLab