Pycnophylactic disaggregation with spatial smoothness on rates
disaggregate_positive(
W,
y_r,
pop,
adj,
lambda_smooth = 1,
lambda_entropy = 0,
weights = NULL
)A list
m x n aggregation matrix
length m vector (observed regional counts)
length n vector (population per cell)
n x n adjacency matrix (1 if neighbours, 0 otherwise) or a sparse Matrix from the Matrix package
weight on spatial smoothness
weight on entropy (0 = pure QP, >0 = nonlinear)
optional n x n matrix of edge weights (e.g., inverse distance)
Solves: min ||W y_s - y_r||^2 + lambda_smooth * sum_(i~j) w_ij (r_i - r_j)^2 + lambda_entropy * sum(y_s * log(y_s)) s.t. y_s >= 0