source('../../Data/get_lunatic_data.R') DF = get_lunatic_data() DF$logPOP = log(DF$POP) ## POP ## m = glm(NBR ~ POP, family=poisson, data=DF) print(summary(m)) ## logPOP ## m = glm(NBR ~ logPOP, family=poisson, data=DF) print(summary(m)) ## PDEN (not as good a predictor as logPOP as the AIC is much larger for this model): ## ##m = glm(NBR ~ PDEN, family=poisson, data=DF) ##print(summary(m)) m = glm(NBR ~ logPOP + DIST, family=poisson, data=DF) print(summary(m)) m = glm(NBR ~ logPOP + I(1/DIST), family=poisson, data=DF) print(summary(m)) ## ## m1 = lm(PHOME ~ PDEN, data=DF) print(summary(m1)) m2 = lm(PHOME ~ DIST, data=DF) print(summary(m2))