Running examples in ‘monomvn-Ex.R’ failed
The error most likely occurred in:

> ### Name: monomvn
> ### Title: Maximum Likelihood Estimation for Multivariate Normal Data with
> ###   Monotone Missingness
> ### Aliases: monomvn
> ### Keywords: multivariate regression
> 
> ### ** Examples
> 
> ## standard usage, duplicating the results in
> ## Little and Rubin, section 7.4.3 -- try adding 
> ## verb=3 argument for a step-by-step breakdown
> data(cement.miss)
> out <- monomvn(cement.miss)
> out

Call:
monomvn(y = cement.miss)

Methods used (p=0.9):
2	complete
3	lsr

> out$mu
        [,1]
x1  6.655166
x2 49.965259
x3 11.769231
x4 27.047089
y  95.423077
> out$S
          x1         x2        x3         x4          y
x1  23.80287   22.26922 -26.97542  -12.28292   50.86578
x2  22.26922  258.86275 -17.13549 -273.80147  211.90393
x3 -26.97542  -17.13549  41.02564  -10.39915  -51.51923
x4 -12.28292 -273.80147 -10.39915  319.38777 -206.48170
y   50.86578  211.90393 -51.51923 -206.48170  226.31359
> 
> ##
> ## A bigger example, comparing the various methods
> ##
> 
> ## generate N=100 samples from a 10-d random MVN
> xmuS <- randmvn(100, 20)
> 
> ## randomly impose monotone missingness
> xmiss <- rmono(xmuS$x)
> 
> ## plsr
> oplsr <- monomvn(xmiss, obs=TRUE)
> oplsr

Call:
monomvn(y = xmiss, obs = TRUE)

Methods used (p=0.9):
1	complete
17	lsr
2	plsr-CV, ncomp range: [3,4]

> Ellik.norm(oplsr$mu, oplsr$S, xmuS$mu, xmuS$S)
[1] -35824.07
> 
> ## calculate the complete and observed RMSEs
> n <- nrow(xmiss) - max(oplsr$na)
> x.c <- xmiss[1:n,]
> mu.c <- apply(x.c, 2, mean)
> S.c <- cov(x.c)*(n-1)/n
> Ellik.norm(mu.c, S.c, xmuS$mu, xmuS$S)
[1] NA
> Ellik.norm(oplsr$mu.obs, oplsr$S.obs, xmuS$mu, xmuS$S)
[1] NA
> 
> ## plcr
> opcr <- monomvn(xmiss, method="pcr")
> Ellik.norm(opcr$mu, opcr$S, xmuS$mu, xmuS$S)
[1] -35819.55
> 
> ## ridge regression
> oridge <- monomvn(xmiss, method="ridge")
> Ellik.norm(oridge$mu, oridge$S, xmuS$mu, xmuS$S)
[1] -35831.75
> 
> ## lasso
> olasso <- monomvn(xmiss, method="lasso")
> Ellik.norm(olasso$mu, olasso$S, xmuS$mu, xmuS$S)
[1] -35876.26
> 
> ## lar
> olar <- monomvn(xmiss, method="lar")
Error in predict.lars(object, type = "coefficient", ...) : 
  Argument s out of range
Calls: monomvn ... regress.lars -> coef -> coef.lars -> predict -> predict.lars
Execution halted