This function both fits distributions and projects those distributions to new covariates..
pnp_gaussian(data, method, type = "regularized", object = NULL)
dataframe of covariates
one of either "fit" or "predict"
one of either "classical", "robust", or "regularized" (the default)
fitted object returned by a pnp_... function. Only needed when method = "predict"