R Under development (unstable) (2024-02-07 r85873 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > library('diseasemapping') > data('casedata') > data('popdata') > popdata = terra::unwrap(popdata) > > therates = getRates(casedata, popdata, ~age*sex,breaks=c(seq(0, 80, by=10), Inf)) > therates Call: stats::glm(formula = formula1, family = family, data = newdata) Coefficients: (Intercept) age[70,80) age[80,Inf) age[50,60) -5.4245 0.2084 -0.1015 -1.9418 age[40,50) age[30,40) age[20,30) age[10,20) -1.7551 -3.0343 -5.1459 -8.1712 age[0,10) sexM age[70,80):sexM age[80,Inf):sexM -7.7040 0.3456 0.0736 0.6850 age[50,60):sexM age[40,50):sexM age[30,40):sexM age[20,30):sexM 0.8576 -0.3938 -0.9722 -0.3720 age[10,20):sexM age[0,10):sexM 0.6994 NA Degrees of Freedom: 16 Total (i.e. Null); 0 Residual Null Deviance: 34560 Residual Deviance: -9.313e-13 AIC: 155.9 > > proc.time() user system elapsed 5.28 0.12 5.39