if(requireNamespace("OpenMx", quietly = TRUE)){
library(OpenMx)
dat <- data.frame(x = iris$Petal.Length)
mixmod <- mx_profiles(dat,
classes = 2)
res <- BCH(mixmod, "y ~ 1", data = data.frame(y = iris$Sepal.Length))
}
Run the code above in your browser using DataLab