missoNet-package | R Documentation |
missoNet fits a joint multivariate regression and conditional
dependency (precision–matrix) model when some response entries are missing.
The method estimates a sparse coefficient matrix B
linking
predictors X
to multivariate responses Y
, together with a sparse
inverse covariance \Theta
for the residuals in
Y = \mathbf{1}\mu^{\mathsf{T}} + XB + E
, E \sim \mathcal{N}(0, \Theta^{-1})
.
Responses may contain missing values (e.g., MCAR/MAR); predictors must be
finite. The package provides cross-validation, prediction, publication-ready
plotting, and simple simulation utilities.
Key features
Joint estimation of B
(regression) and \Theta
(conditional network).
\ell_1
-regularization on both B
and \Theta
with user-controlled grids.
K-fold cross-validation with optional 1-SE model selections.
Heatmap and 3D surface visualizations for CV error or GoF across
(\lambda_B,\lambda_\Theta)
.
Fast prediction for new data using stored solutions.
Lightweight data generator for simulation studies.
Workflow
Fit a model across a grid of penalties with missoNet
or
select penalties via cv.missoNet
.
Visualize the CV error/GoF surface with plot.missoNet
.
Predict responses for new observations with predict.missoNet
.
missoNet
Fit models over user-specified penalty grids for
\lambda_B
and \lambda_\Theta
; returns estimated
\mu
, B
, \Theta
, and metadata (grids, GoF).
cv.missoNet
Perform k-fold cross-validation over a penalty grid;
stores est.min
and (optionally) est.1se.beta
,
est.1se.theta
.
plot.missoNet
S3 plotting method; heatmap or 3D scatter of CV error or GoF.
predict.missoNet
S3 prediction method; returns
\hat{Y} = \mathbf{1}\hat{\mu}^{\mathsf{T}} + X_\mathrm{new}\hat{B}
for a chosen solution.
generateData
Generate synthetic datasets with controllable dimensions, signal, and missingness mechanisms for benchmarking.
GPL-2.
Maintainer: Yixiao Zeng yixiao.zeng@mail.mcgill.ca [copyright holder]
Authors:
Celia Greenwood celia.greenwood@mcgill.ca [thesis advisor]
missoNet
, cv.missoNet
, plot.missoNet
,
predict.missoNet
, generateData
,
and browseVignettes("missoNet")
for tutorials.
sim <- generateData(n = 100, p = 8, q = 5, rho = 0.1, missing.type = "MCAR")
fit <- missoNet(X = sim$X, Y = sim$Z) # fit over a grid
plot(fit) # GoF heatmap
cvfit <- cv.missoNet(X = sim$X, Y = sim$Z, kfold = 5, compute.1se = TRUE)
plot(cvfit, type = "scatter", plt.surf = TRUE) # CV error surface
yhat <- predict(cvfit, newx = sim$X, s = "lambda.min")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.