Running examples in 'pbkrtest-Ex.R' failed
The error most likely occurred in:

> ### Name: get_ddf_Lb
> ### Title: Adjusted denomintor degress freedom for linear estimate for
> ###   linear mixed model.
> ### Aliases: Lb_ddf ddf_Lb get_Lb_ddf get_Lb_ddf.lmerMod get_ddf_Lb
> ###   get_ddf_Lb.lmerMod
> ### Keywords: inference models
> 
> ### ** Examples
> 
> 
> (fmLarge <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
Linear mixed model fit by REML ['lmerMod']
Formula: Reaction ~ Days + (Days | Subject)
   Data: sleepstudy
REML criterion at convergence: 1743.628
Random effects:
 Groups   Name        Std.Dev. Corr
 Subject  (Intercept) 24.737       
          Days         5.923   0.07
 Residual             25.592       
Number of obs: 180, groups:  Subject, 18
Fixed Effects:
(Intercept)         Days  
     251.41        10.47  
> ## removing Days
> (fmSmall <- lmer(Reaction ~ 1 + (Days|Subject), sleepstudy))
Linear mixed model fit by REML ['lmerMod']
Formula: Reaction ~ 1 + (Days | Subject)
   Data: sleepstudy
REML criterion at convergence: 1769.845
Random effects:
 Groups   Name        Std.Dev. Corr 
 Subject  (Intercept) 25.53         
          Days        11.93    -0.18
 Residual             25.59         
Number of obs: 180, groups:  Subject, 18
Fixed Effects:
(Intercept)  
      257.8  
> anova(fmLarge,fmSmall)
refitting model(s) with ML (instead of REML)
Data: sleepstudy
Models:
fmSmall: Reaction ~ 1 + (Days | Subject)
fmLarge: Reaction ~ Days + (Days | Subject)
        Df    AIC    BIC  logLik deviance  Chisq Chi Df Pr(>Chisq)    
fmSmall  5 1785.5 1801.4 -887.74   1775.5                             
fmLarge  6 1763.9 1783.1 -875.97   1751.9 23.537      1  1.226e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> 
> KRmodcomp(fmLarge, fmSmall)  ## 17 denominator df's
F-test with Kenward-Roger approximation; computing time: 1.00 sec.
large : Reaction ~ Days + (Days | Subject)
small : Reaction ~ 1 + (Days | Subject)
        stat    ndf    ddf F.scaling   p.value    
Ftest 45.843  1.000 17.000         1 3.268e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> get_Lb_ddf(fmLarge, c(0,1)) ## 17 denominator df's
[1] 17
> 
> # Notice: The restriction matrix L corresponding to the test above
> # can be found with
> L <- model2restrictionMatrix(fmLarge, fmSmall)
> L
1 x 2 sparse Matrix of class "dgCMatrix"
Error in isFALSE(suppRows) : could not find function "isFALSE"
Calls: <Anonymous> -> <Anonymous> -> printSpMatrix2
Execution halted