Running ‘DAS-scale.R’
  Running ‘DASvar-fallback.R’ [8s/21s]
  Running ‘PsiFunction.R’ [13s/31s]
  Running ‘asymptoticEfficiency.R’ [6s/16s]
  Running ‘compare-methods.R’ [17s/43s]
  Comparing ‘compare-methods.Rout’ to ‘compare-methods.Rout.save’ ...191,199d190
< Warning messages:
< 1: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
<   convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
< 2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
<   Model failed to converge with max|grad| = 0.494878 (tol = 0.002, component 1)
< 3: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
<   convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
< 4: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
<   Model failed to converge with max|grad| = 0.494878 (tol = 0.002, component 1)
203a195
> boundary (singular) fit: see help('isSingular')
240a233
> boundary (singular) fit: see help('isSingular')
302,310d294
< Warning messages:
< 1: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
<   convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
< 2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
<   Model failed to converge with max|grad| = 1.75912 (tol = 0.002, component 1)
< 3: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
<   convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
< 4: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
<   Model failed to converge with max|grad| = 1.75912 (tol = 0.002, component 1)
355c339
<  4.436e-01  4.436e-01  1.241e-01  2.030e-01 -1.084e-06 -2.536e-01 -8.075e-01 
---
>  4.436e-01  4.436e-01  1.241e-01  2.030e-01 -1.094e-06 -2.536e-01 -8.075e-01 
406c390
<  0.450838  0.450838  0.128411  0.205546  0.000502 -0.256535 -0.815539  0.490628 
---
>  0.450838  0.450838  0.128411  0.205545  0.000502 -0.256535 -0.815539  0.490628 
412c396
<  0.559398 -0.265494  0.499107 -0.034026 -0.023533 -0.794300 
---
>  0.559397 -0.265494  0.499107 -0.034026 -0.023533 -0.794300 
441,449d424
< Warning messages:
< 1: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
<   convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
< 2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
<   Model failed to converge with max|grad| = 1.09789 (tol = 0.002, component 1)
< 3: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
<   convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
< 4: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
<   Model failed to converge with max|grad| = 1.09789 (tol = 0.002, component 1)
  Running ‘expectations.R’
  Running ‘fitDatasets.R’ [54s/145s]
  Running ‘genericFunctions.R’ [8s/20s]
  Running ‘getME.R’ [8s/22s]
  Comparing ‘getME.Rout’ to ‘getME.Rout.save’ ... OK
  Running ‘multipleGroupingFactorsTestData.R’
  Running ‘offset.R’ [8s/23s]
  Running ‘psi-rho-funs.R’ [7s/19s]
  Running ‘rlmerMod.R’
  Running ‘simulationStudies.R’ [7s/18s]
  Running ‘subset.R’
  Running ‘tau.R’
  Running ‘testMatrices.R’
  Running ‘weights.R’ [8s/20s]
Running the tests in ‘tests/weights.R’ failed.
Complete output:
  > require(robustlmm)
  Loading required package: robustlmm
  Loading required package: lme4
  Loading required package: Matrix
  > 
  > checkEquality <- function(cW, rW, tolerance = 1e-8) {
  +     stopifnot(
  +         all.equal(
  +             coef(cW),
  +             coef(rW),
  +             tolerance = tolerance,
  +             check.attributes = FALSE
  +         ),
  +         all.equal(
  +             fixef(cW),
  +             fixef(rW),
  +             tolerance = tolerance,
  +             check.attributes = FALSE
  +         ),
  +         all.equal(
  +             ranef(cW) ,
  +             ranef(rW),
  +             tolerance = 100 * tolerance,
  +             check.attributes = FALSE
  +         ),
  +         all.equal(
  +             fitted(cW) ,
  +             fitted(rW),
  +             tolerance = tolerance,
  +             check.attributes = FALSE
  +         ),
  +         all.equal(
  +             predict(cW) ,
  +             predict(rW),
  +             tolerance = tolerance,
  +             check.attributes = FALSE
  +         ),
  +         all.equal(
  +             coef(summary(cW)) ,
  +             coef(summary(rW)),
  +             tolerance = 100 * tolerance,
  +             check.attributes = FALSE
  +         )
  +     )
  + }
  > 
  > testBattery <- function(formula, data, tolerance) {
  +     nobs <- nrow(data)
  + 
  +     test <- function(weights) {
  +         cW <- lmer(formula, data, weights = weights)
  +         rW <-
  +             rlmer(
  +                 formula,
  +                 data,
  +                 weights = weights,
  +                 rho.e = cPsi,
  +                 rho.b = cPsi,
  +                 init = lmerNoFit
  +             )
  +         checkEquality(cW, rW, tolerance)
  +     }
  + 
  +     test(rep(2, nobs))
  +     test(rep(0.5, nobs))
  + 
  +     set.seed(133)
  +     test(runif(nobs))
  +     test(rexp(nobs))
  + }
  > 
  > testBattery(Yield ~ 1 | Batch, Dyestuff, 1e-8)
  Error in checkEquality(cW, rW, tolerance) : 
    coef(cW) and coef(rW) are not equal:
    Component "Batch": Component "(Intercept)": Mean relative difference: 0.002937986
  Calls: testBattery -> test -> checkEquality -> stopifnot
  In addition: Warning messages:
  1: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
    convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
  2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
    Model failed to converge with max|grad| = 2.53309 (tol = 0.002, component 1)
  3: In optwrap(optimizer, devfun, getStart(start, rho$pp), lower = rho$lower,  :
    convergence code -2 from nloptwrap: NLOPT_INVALID_ARGS: Invalid arguments (e.g. lower bounds are bigger than upper bounds, an unknown algorithm was specified, etcetera).
  4: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
    Model failed to converge with max|grad| = 2.53309 (tol = 0.002, component 1)
  Execution halted