[ferret_users] Spatial Correlation with regressxy

[Molulaqhooa Maoyi <molulaqhooa.maoyi@xxxxxxxxx>]

Dear Ferret Users.

I am struggling to find the spatial correlation between two datasets. I keep on getting the error  SET GRID Q" after defining Q. **ERROR: unknown command: missing
missing. I am not sure what is missing.

Initially, I defined Q without specifying gt=precip[d=1] because both datasets are on the same grid. However, that did not work and tried regridding again by adding gt=precip[d=1]

but that did not work still. Any suggestions?

use final.gpcp.nc
use final.20th.nc
 
let P = precip[d=1]
let Q = prate[d=2, gt=precip[d=1]]
 
SET GRID Q

 exit

go regressxy

... Linear Regression In the XY Plane
... Instructions:
Use the LET command to define new variables
Define the variable P as your independent (X) variable
Define the variable Q as your   dependent (Y) variable
Results will be variables "SLOPE", "INTERCEP" and "RSQUARE"
QHAT will be the regression estimate
Note: If "Y" is your independent variable then
...   "SET GRID Q" after defining Q.
...
 **ERROR: unknown command: missing
missing
Command file, command group, or REPEAT execution aborted
yes? sh d
     currently SET data sets:
    1> ./final.gpcp.nc
 name     title                             I         J         K         L         M         N
 PRECIP   Average Monthly Rate of Precipi  1:720     1:360     ...       1:1       ...       ...
 
    2> ./final.20th.nc
 name     title                             I         J         K         L         M         N
 TIME_BNDS_2
          Time Boundaries                  1:2       ...       ...       1:1       ...       ...
 PRATE    Monthly Precipitation Rate at s  1:720     1:360     ...       1:1       ...       ...

Kind Regards

--

Linda

"Try and Fail than to Fail to Try" - Israel Maoyi