medflex

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

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: neImpute.default
> ### Title: Expand the dataset and impute nested counterfactual outcomes
> ### Aliases: neImpute.default
> 
> ### ** Examples
> 
> data(UPBdata)
> 
> ## example using glm imputation model with binary exposure
> fit.glm <- glm(UPB ~ factor(attbin) + negaff + gender + educ + age, 
+                family = binomial, data = UPBdata)
> impData <- neImpute(fit.glm)
> head(impData)
  id attbin0 attbin1        att attcat    negaff initiator gender educ age
1  1       1       1  1.0005617      M  0.840461    myself      F    M  41
2  1       0       1  1.0005617      M  0.840461    myself      F    M  41
3  2       0       0 -0.7085889      L -1.257465      both      M    M  42
4  2       1       0 -0.7085889      L -1.257465      both      M    M  42
5  3       0       0 -0.7085889      L -1.202256      both      F    H  43
6  3 
...[truncated]...
 1      M  0.840461    myself      F    M  41 0.5106796    1
2  1  low high      1      M  0.840461    myself      F    M  41 0.4059912    1
3  2  low  low      0      L -1.257465      both      M    M  42 0.1999050    0
4  2 high  low      0      L -1.257465      both      M    M  42 0.2761580    0
5  3  low  low      0      L -1.202256      both      F    H  43 0.1668007    0
6  3 high  low      0      L -1.202256      both      F    H  43 0.2341212    0
> 
> # test with gam?
> if (requireNamespace("gam", quietly = TRUE)) {
+ library(gam)
+ impFit4 <- gam(UPB ~ att + negaff + gender + educ + age, family = binomial, data = UPBdata)
+ impData4 <- neImpute(impFit4)
+ head(impData2); head(impData4)
+ # check!
+ }
Loading required package: foreach
Loaded gam 1.22-7


Attaching package: ‘gam’

The following object is masked from ‘package:VGAM’:

    s

Error in a$specials[[i]] : subscript out of bounds
Calls: neImpute ... predict.glm -> predict.lm -> delete.response -> off
Execution halted