pscl_analyze.Rd
Run pattern scaling on the input temperature fields. This amounts to doing a linear regression on each grid cell, with the global temperature as the sole regressor.
pscl_analyze(tin, tgavin)
tin | Matrix of temperature fields (ntime x ngrid) |
---|---|
tgavin | Matrix of global mean temperature over time (ntime x 1) |
List of r, w, and b matrices (see details)
The return value is a list with the residuals and pattern coefficients as r (resid), w (slope), and b (intercept).
dim(r) == dim(tin) == ntime x nlat*nlon dim(w) == dim(b) == nlat*nlon