Closed BERENZ closed 2 years ago
Imagine we have an aggregated dataset where each row corresponds to two observations in the whole dataset. This is denoted by rep.int(2, nrow(farmsubmission))
rep.int(2, nrow(farmsubmission))
ModelPo <- estimate_popsize(formula = TOTAL_SUB ~ ., data = farmsubmission, pop.var = "analytic", model = "ztpoisson", method = "mle", weights = rep.int(2, nrow(farmsubmission)))
so the observed population should be sum(nrow(farmsubmission)* weights) = 24072 but it is not included included in the output
sum(nrow(farmsubmission)* weights) = 24072
Observed proportion: 65.6% (N obs = 12036)
Added to 0.3.0 roadmap.
Imagine we have an aggregated dataset where each row corresponds to two observations in the whole dataset. This is denoted by
rep.int(2, nrow(farmsubmission))
so the observed population should be
sum(nrow(farmsubmission)* weights) = 24072
but it is not included included in the output