Rglimclim-package {Rglimclim}R Documentation

Generalized Linear Modelling for Daily Climate Time Series

Description

Modelling and simulation of daily multi-site time series, based on generalised linear models and designed in particular for climatological downscaling applications. The references below give an overview of the theory upon which much of this work is based.

The mathematical specification of the models can be found in Section 3 and in the Appendix of the PDF package manual.

Details

This package is designed to allow flexibility in fitting and simulating fairly complicated models for daily climate/weather data, and to avoid some of the tedious manipulation of data files that would be necessary if standard software packages were used to perform this kind of analysis. The package is written to take advantage of the fact that all exercises involving fitting, and simulating, GLMs for daily weather series must inevitably share common features, as follows:

Climate datasets often have other unusual features relating to measurement methods. For example, in daily raingauge data any non-zero amount that is less than some small threshold may be recorded as a ‘trace’ amount, because it is too small to be measured accurately. Such features pose potential problems for statistical analysis. The software aims to provide methods for dealing with them.

The package is essentially an R port of the original GLIMCLIM package (Chandler, 2002) which was written in Fortran 77. The Fortran code for model fitting and simulation has been retained since it is considerably faster than, for example, the glm command in R itself. However, the R interface makes it much easier to save and update models, and to transfer and visualise information. In particular, GLIMCLIM made extensive use of definition files to define sites and models. All of these had to be prepared manually, with a consequent risk of error. Rglimclim contains routines that will read old GLIMCLIM definition files; however, it is intended that the primary means of defining sites, models etc. will be via manipulation of R objects, and hence that definition files will eventually become defunct.

Author(s)

Richard E. Chandler

Maintainer: Richard E. Chandler <richard@stats.ucl.ac.uk>

References

Chandler, R.E. (2002). GLIMCLIM: Generalised linear modelling for daily climate time series (software and user guide). Research Report No. 227, Department of Statistical Science, University College London (available from http://www.ucl.ac.uk/statistics/research/reports).

Chandler, R.E. (2005). On the use of generalized linear models for interpreting climate variability. Environmetrics 16(7): 699-715.

Chandler, R.E. and Bate, S. (2007). Inference for clustered data using the independence log-likelihood. Biometrika 94: 167-183.

Yan, Z., Bate, S., Chandler, R.E., Isham, V. and Wheater, H. (2002). An analysis of daily maximum windspeed in northwestern Europe using generalized linear models. J. Climate 15(15): 2073-2088.

Yang, C., Chandler, R.E., Isham, V. and Wheater, H.S. (2005). Spatial-temporal rainfall simulation using Generalized Linear Models. Water Resources Research 41, doi:10.1029/2004WR003739.

See Also

See the help for GLCfit (the main model fitting routine) to get started. A full description of the package can be found in the PDF package manual.


[Package Rglimclim version 1.2-2 Index]