D_out <- D_regularized(
data = iris[iris$Species == "versicolor" | iris$Species == "virginica", ],
mv.vars = c("Sepal.Length", "Sepal.Width", "Petal.Length", "Petal.Width"),
group.var = "Species", group.values = c("versicolor", "virginica"),
out = TRUE,
size = 15
)
pcc(
data = D_out$pred.dat,
pred.var = "pred",
group.var = "group",
group.values = c("versicolor", "virginica")
)
Run the code above in your browser using DataLab