tests/testthat/prediction.R

In semtree: Recursive Partitioning for Structural Equation Models

skip_on_cran()

N<-100
x<-rnorm(N)
y<-rnorm(N)
pred1 <- ordered( sample(c(0,1,2),N,replace=TRUE) )
pred2 <- sample(0:10, N, replace=TRUE)

y2 <- ifelse(pred2>5,0.2*y+0.8*x,x)              

model <- "x~~y; x~~x;y~~y"

df <- data.frame(x,y=y2,pred1,pred2)

fitted_model <- lavaan(model,df)
tree <- semtree(fitted_model, df)
plot(tree)

forest <- semforest(fitted_model, df)

pars_pred <- predict(forest, df, "pars")


#pars_pred <- predict(strip(forest), df, "pars")