Running ‘DAS-scale.R’ [0s/0s]
  Running ‘DASvar-fallback.R’ [3s/3s]
  Running ‘PsiFunction.R’ [4s/4s]
  Running ‘asymptoticEfficiency.R’ [2s/2s]
  Running ‘compare-methods.R’ [5s/5s]
  Comparing ‘compare-methods.Rout’ to ‘compare-methods.Rout.save’ ...93,94c93,127
< Timing stopped at: 0.172 0.015 0.186
< Error : invalid assignment for reference class field 'b.s', should be from class "numeric" or a subclass (was class "dgeMatrix")
---
> Robust linear mixed model fit by DASvar 
> Formula: formula 
>    Data: data 
> 
> Scaled residuals: 
>     Min      1Q  Median      3Q     Max 
> -1.2872 -0.6775  0.1308  0.7136  1.7416 
> 
> Random effects:
>  Groups   Name        Variance Std.Dev.
>  Batch    (Intercept) 2184     46.73   
>  Residual             2933     54.16   
> Number of obs: 30, groups: Batch, 6
> 
> Fixed effects:
>             Estimate Std. Error t value
> (Intercept)  1526.32      22.04   69.25
> 
> Robustness weights for the residuals: 
>  25 weights are ~= 1. The remaining 5 ones are
>     2     3     5    18    23 
> 0.988 0.988 0.984 0.770 0.988 
> 
> Robustness weights for the random effects: 
> [1] 1.000 1.000 1.000 1.000 0.999 1.000
> 
> Rho functions used for fitting:
>   Residuals:
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     sig: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>   Random Effects, variance component 1 (Batch):
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     vcp: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>       A       B       C       D       E       F 
> -19.822   1.081  33.966 -30.283  67.013 -51.900 
97,98c130,164
< Timing stopped at: 0.057 0.001 0.058
< Error : invalid assignment for reference class field 'b.s', should be from class "numeric" or a subclass (was class "dgeMatrix")
---
> Robust linear mixed model fit by DAStau 
> Formula: formula 
>    Data: data 
> 
> Scaled residuals: 
>     Min      1Q  Median      3Q     Max 
> -1.2656 -0.6593  0.1234  0.6972  1.7041 
> 
> Random effects:
>  Groups   Name        Variance Std.Dev.
>  Batch    (Intercept) 2359     48.57   
>  Residual             3059     55.31   
> Number of obs: 30, groups: Batch, 6
> 
> Fixed effects:
>             Estimate Std. Error t value
> (Intercept)  1526.34      22.82   66.88
> 
> Robustness weights for the residuals: 
>  25 weights are ~= 1. The remaining 5 ones are
>     2     3     5    18    23 
> 0.996 0.996 0.992 0.787 0.994 
> 
> Robustness weights for the random effects: 
> [1] 1 1 1 1 1 1
> 
> Rho functions used for fitting:
>   Residuals:
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     sig: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>   Random Effects, variance component 1 (Batch):
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     vcp: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>       A       B       C       D       E       F 
> -19.761   1.048  33.601 -29.717  66.229 -51.399 
105c171
<        
---
>                     IRWLS    
107c173
<        
---
> (Intercept)         1526 (22)
109a176
> (Intercept) | Batch 46.7     
110a178
> sigma               54.2     
112d179
< sigma  
115c182
<        
---
>                     IRWLS      
117c184
<        
---
> (Intercept)         1526 (22.8)
119a187
> (Intercept) | Batch 48.6       
120a189
> sigma               55.3       
122d190
< sigma  
231,232c299,347
< Timing stopped at: 0.075 0 0.074
< Error : invalid assignment for reference class field 'b.s', should be from class "numeric" or a subclass (was class "dgeMatrix")
---
> Robust linear mixed model fit by DASvar 
> Formula: formula 
>    Data: data 
> 
> Scaled residuals: 
>     Min      1Q  Median      3Q     Max 
> -2.2164 -0.6474  0.0528  0.5885  3.4337 
> 
> Random effects:
>  Groups   Name        Variance Std.Dev.
>  plate    (Intercept) 0.8798   0.9380  
>  sample   (Intercept) 4.3155   2.0774  
>  Residual             0.2954   0.5435  
> Number of obs: 144, groups: plate, 24; sample, 6
> 
> Fixed effects:
>             Estimate Std. Error t value
> (Intercept)  23.0116     0.8929   25.77
> 
> Robustness weights for the residuals: 
>  122 weights are ~= 1. The remaining 22 ones are summarized as
>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
>   0.392   0.669   0.843   0.822   0.978   0.999 
> 
> Robustness weights for the random effects: 
>  25 weights are ~= 1. The remaining 5 ones are
>     7    13    19    24    30 
> 0.893 0.853 0.893 0.982 0.905 
> 
> Rho functions used for fitting:
>   Residuals:
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     sig: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>   Random Effects, variance component 1 (plate):
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     vcp: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>   Random Effects, variance component 2 (sample):
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     vcp: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>          a          b          c          d          e          f          g 
>  4.436e-01  4.436e-01  1.241e-01  2.030e-01 -1.094e-06 -2.536e-01 -8.075e-01 
>          h          i          j          k          l          m          n 
>  4.872e-01 -4.441e-01 -4.541e-01  5.688e-01  2.945e-01  8.498e-01  2.945e-01 
>          o          p          q          r          s          t          u 
>  5.481e-01  2.012e-02 -1.615e-01 -1.623e-01 -8.075e-01  5.688e-01 -4.951e-01 
>          v          w          x          A          B          C          D 
> -1.628e-01 -4.495e-01 -7.084e-01  5.551e-01 -2.672e-01  4.955e-01 -3.700e-02 
>          E          F 
> -2.667e-02 -7.956e-01 
235,236c350,396
< Timing stopped at: 0.074 0 0.074
< Error : invalid assignment for reference class field 'b.s', should be from class "numeric" or a subclass (was class "dgeMatrix")
---
> Robust linear mixed model fit by DAStau 
> Formula: formula 
>    Data: data 
> 
> Scaled residuals: 
>    Min     1Q Median     3Q    Max 
> -2.161 -0.642  0.057  0.574  3.333 
> 
> Random effects:
>  Groups   Name        Variance Std.Dev.
>  plate    (Intercept) 0.9054   0.9515  
>  sample   (Intercept) 4.5338   2.1293  
>  Residual             0.3114   0.5581  
> Number of obs: 144, groups: plate, 24; sample, 6
> 
> Fixed effects:
>             Estimate Std. Error t value
> (Intercept)  23.0002     0.9148   25.14
> 
> Robustness weights for the residuals: 
>  126 weights are ~= 1. The remaining 18 ones are summarized as
>    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
>   0.404   0.684   0.838   0.802   0.920   0.997 
> 
> Robustness weights for the random effects: 
>  25 weights are ~= 1. The remaining 5 ones are
>     7    13    19    24    30 
> 0.906 0.866 0.906 0.989 0.926 
> 
> Rho functions used for fitting:
>   Residuals:
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     sig: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>   Random Effects, variance component 1 (plate):
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     vcp: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>   Random Effects, variance component 2 (sample):
>     eff: smoothed Huber (k = 1.345, s = 10) 
>     vcp: smoothed Huber, Proposal 2 (k = 1.345, s = 10) 
>         a         b         c         d         e         f         g         h 
>  0.450838  0.450838  0.128411  0.205545  0.000502 -0.256535 -0.815539  0.490628 
>         i         j         k         l         m         n         o         p 
> -0.446899 -0.456519  0.575339  0.297994  0.858936  0.297994  0.554159  0.020649 
>         q         r         s         t         u         v         w         x 
> -0.163812 -0.164206 -0.815539  0.575339 -0.501182 -0.164248 -0.451685 -0.717135 
>         A         B         C         D         E         F 
>  0.559397 -0.265494  0.499107 -0.034026 -0.023533 -0.794300 
243c403
<        
---
>                      IRWLS     
245c405
<        
---
> (Intercept)          23 (0.893)
247a408,409
> (Intercept) | plate  0.938     
> (Intercept) | sample 2.077     
248a411
> sigma                0.543     
250d412
< sigma  
253c415
<        
---
>                      IRWLS     
255c417
<        
---
> (Intercept)          23 (0.915)
257a420,421
> (Intercept) | plate  0.952     
> (Intercept) | sample 2.129     
258a423
> sigma                0.558     
260d424
< sigma  
  Running ‘expectations.R’ [0s/0s]
  Running ‘fitDatasets.R’ [3s/3s]
Running the tests in ‘tests/fitDatasets.R’ failed.
Last 13 lines of output:
  Fixed Effects:
  (Intercept)  
        1.899  
  > fitDatasets_rlmer_custom <- function(datasets) {
  +     return(fitDatasets_rlmer(datasets,
  +                              method = "DASvar",
  +                              tuningParameter = c(1.345, 2.28, 1.345, 2.28, 5.14, 5.14),
  +                              label = "fitDatasets_rlmer_custom"))
  + }
  > test(fitDatasets_rlmer_custom(oneWay))
  Error : invalid assignment for reference class field 'b.s', should be from class "numeric" or a subclass (was class "dgeMatrix")
  Error in test(fitDatasets_rlmer_custom(oneWay)) : 
    !is(result[[1]], "try-error") is not TRUE
  Calls: test -> stopifnot
  Execution halted