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

> ### Name: reSamGF
> ### Title: Gross flows variance estimation.
> ### Aliases: reSamGF
> 
> ### ** Examples
> 
> library(TeachingSampling)
Loading required package: dplyr

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

Loading required package: magrittr
> library(data.table)

Attaching package: ‘data.table’

The following objects are masked from ‘package:dplyr’:

    between, first, last

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

    %notin%

> # Colombia's electoral candidates in 2014
> candidates_t0 <- c("Clara","Enrique","Santos","Martha","Zuluaga","Blanco", "NoVoto")
> candidates_t1 <- c("Santos","Zuluaga","Blanco", "NoVoto")
> 
> N 	   <- 100000
> nCanT0  <- length(candidates_t0)
> nCanT1  <- length(candidates_t1)
> 
> # Initial probabilities 
> eta <- matrix(c(0.10, 0.10, 0.20, 0.17, 0.28, 0.1, 0.05),
+ 			   byrow = TRUE, nrow = nCanT0)
> # Transition probabilities
> P <- matrix(c(0.10, 0.60, 0.15, 0.15,
+ 			  0.30, 0.10, 0.25,	0.35,
+ 			  0.34, 0.25, 0.16, 0.25,
+ 			  0.25,	0.05, 0.35,	0.35,
+ 			  0.10, 0.25, 0.45,	0.20,
+ 			  0.12, 0.36, 0.22, 0.30,
+ 			  0.10,	0.15, 0.30,	0.45),
+ 	 byrow = TRUE, nrow = nCanT0)
> 
> citaMod <- matrix(, ncol = nCanT1, nrow = nCanT0)
> row.names(citaMod) <- candidates_t0
> colnames(citaMod) <- candidates_t1
> 
> for(ii in 1:nCanT0){ 
+  citaMod[ii,] <- c(rmultinom(1, size = N * eta[ii,], prob = P[ii,]))
+ }
> 
> # # Model I
> psiI   <- 0.9
> rhoRRI <- 0.9
> rhoMMI <- 0.5
> 
> citaModI <- matrix(nrow = nCanT0 + 1, ncol = nCanT1 + 1)
> rownames(citaModI) <- c(candidates_t0, "Non_Resp")
> colnames(citaModI) <- c(candidates_t1, "Non_Resp")
> 
> citaModI[1:nCanT0, 1:nCanT1] 		 <- P * c(eta) * rhoRRI * psiI  
> citaModI[(nCanT0 + 1), (nCanT1 + 1)]  <- rhoMMI * (1-psiI) 
> citaModI[1:nCanT0, (nCanT1 + 1)]  	 <-  (1-rhoRRI) * psiI * rowSums(P * c(eta))
> citaModI[(nCanT0 + 1), 1:nCanT1 ] 	 <-  (1-rhoMMI) * (1-psiI) * colSums(P * c(eta))
> citaModI   <- round_preserve_sum(citaModI * N)
> DBcitaModI <- createBase(citaModI)
> 
> # Creating auxiliary information
> DBcitaModI[,AuxVar := rnorm(nrow(DBcitaModI), mean = 45, sd = 10)]
> # Selects a sample with unequal probabilities
> res <- S.piPS(n = 1200, as.data.frame(DBcitaModI)[,"AuxVar"])
> sam <- res[,1]
> pik <- res[,2]
> DBcitaModISam <- copy(DBcitaModI[sam,])
> DBcitaModISam[,Pik := pik]
> 
> # Gross flows estimation
> estima <- estGF(sampleBase = DBcitaModISam, niter = 500, model = "II", colWeights = "Pik")
> # gross flows variance estimation
> varEstima <- reSamGF(sampleBase = DBcitaModISam, type = "Bootstrap", nRepBoot = 100,
+ 						model = "II", niter = 101,  colWeights = "Pik")
Error:
! `group_by_()` was deprecated in dplyr 0.7.0 and is now defunct.
ℹ Please use `group_by()` instead.
ℹ See vignette('programming') for more help
Backtrace:
     ▆
  1. ├─GFE::reSamGF(...)
  2. │ └─GFE (local) varParam(EstGF)
  3. │   └─... %>% as.data.frame()
  4. ├─base::as.data.frame(.)
  5. ├─dplyr::summarize_all(...)
  6. │ └─dplyr:::manip_all(...)
  7. │   ├─rlang::syms(tbl_nongroup_vars(.tbl))
  8. │   │ └─rlang:::map(x, sym)
  9. │   │   └─base::lapply(.x, .f, ...)
 10. │   └─dplyr::tbl_nongroup_vars(.tbl)
 11. │     ├─generics::setdiff(tbl_vars(x), group_vars(x))
 12. │     └─dplyr::tbl_vars(x)
 13. │       ├─dplyr:::new_sel_vars(tbl_vars_dispatch(x), group_vars(x))
 14. │       │ └─base::structure(...)
 15. │       └─dplyr:::tbl_vars_dispatch(x)
 16. └─dplyr::group_by_(., var2)
 17.   └─dplyr:::lazy_defunct("group_by")
 18.     └─lifecycle::deprecate_stop(...)
 19.       └─lifecycle:::deprecate_stop0(msg)
 20.         └─rlang::cnd_signal(...)
Execution halted