DESCRIPTION

Module r.futures.potential implements POTENTIAL submodel as a part of FUTURES land change model. POTENTIAL is implemented using a set of coefficients that relate a selection of site suitability factors to the probability of a place becoming developed. This is implemented using the parameter table in combination with maps of those site suitability factors (mapped predictors). The coefficients are obtained by conducting multilevel logistic regression in R with package lme4 where the coefficients may vary by subregions. The best model is selected automatically using dredge function from package MuMIn (which has numerous caveats).

Module r.futures.potential can run it two modes. Without the -d flag, it uses all the given predictors to construct the model. With -d flag, it evaluates all the different combinations of predictors and picks the best one based on AIC.

Format

The format of the output file is a CSV file (but with tabs as delimiters). The header contains the names of the predictor maps and the first column contains the identifiers of the subregions. The order of columns is important, the second column represents intercept, the third development pressure and then the predictors. Therefore the development pressure column must be specified as the first column in option columns.
ID	Intercept	devpressure_0_5	slope	road_dens_perc	forest_smooth_perc ...
37037	-1.873	12.595	-0.0758	0.0907	-0.0223	...
37063	-2.039	12.595	-0.0758	0.0907	-0.0223	...
37069	-1.795	12.595	-0.0758	0.0907	-0.0223	...
37077	-1.264	12.595	-0.0758	0.0907	-0.0223	...
37085	-1.925	12.595	-0.0758	0.0907	-0.0223	...
...

NOTES

Note that this module is designed to automate the FUTURES workflow by brute-force selection of model, which has numerous caveats.

In case there is only one subregion, R function glm is used instead of glmer.

EXAMPLES

SEE ALSO

FUTURES, r.futures.pga, r.futures.devpressure, r.futures.potsurface, r.futures.demand, r.futures.calib, r.sample.category

REFERENCES

AUTHOR

Anna Petrasova, NCSU GeoForAll

Last changed: $Date$