tests/testthat/test-fastZI.R

In linogaliana/gravity: Gravity model with large scale data

# # Imagine 10 variables: 6 in selection, 4 in count
# params <- rnorm(10)
# X <- replicate(4, rnorm(200))
# Z <- cbind(
#   X,
#   replicate(2, rnorm(200))
# )
# Y <- rpois(200, lambda = 3)
# 
# 
# weights = rep(1, nrow(X))
# offsetx = rep(0, nrow(X))
# offsetz = rep(0, nrow(X))
# 
# # gravity:::fastZI(X, Z, Y)
# # gravity:::fastZI(X, Z, Y, link = "logit")
# 
# 
# 
# 
# 
#