# \donttest{
if (requireNamespace("lavaan", quietly = TRUE)) {
set.seed(42)
n <- 100
nitrogen <- rnorm(n, 60, 10)
phosphorus <- rnorm(n, 30, 5)
yield <- 0.5 * nitrogen + 0.3 * phosphorus + rnorm(n, 0, 5)
df <- data.frame(yield = yield, nitrogen = nitrogen, phosphorus = phosphorus)
model <- 'yield ~ nitrogen + phosphorus'
result <- sem_analysis(df, model, plot = FALSE)
}
# }
Run the code above in your browser using DataLab