dat <- simdat(model = "mean")
form <- Y - f ~ 1
Y_l <- dat[dat$set_label == "labeled", all.vars(form)[1]] |> matrix(ncol = 1)
f_l <- dat[dat$set_label == "labeled", all.vars(form)[2]] |> matrix(ncol = 1)
f_u <- dat[dat$set_label == "unlabeled", all.vars(form)[2]] |> matrix(ncol = 1)
ppi_mean(Y_l, f_l, f_u)
Run the code above in your browser using DataLab