if(requireNamespace("OpenMx", quietly = TRUE)){
library(OpenMx)
iris_sample <- iris[c(1:5, 145:150), c("Sepal.Length", "Sepal.Width")]
names(iris_sample) <- c("x", "y")
res <- mx_profiles(iris_sample, classes = 2)
plot_bivariate(res, rawdata = FALSE)
}
Run the code above in your browser using DataLab