Re: [ferret_users] Spatial Correlation with regressxy

To: ferret_users@xxxxxxxx, Molulaqhooa Maoyi <molulaqhooa.maoyi@xxxxxxxxx>

Subject: Re: [ferret_users] Spatial Correlation with regressxy

From: "'Dr Kamoru Abiodun Lawal' via _OAR PMEL Ferret Users" <ferret_users@xxxxxxxx>

Date: Mon, 19 Aug 2019 08:40:12 +0000 (UTC)

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@yahoo.co.uk header.s=s2048 header.b=HWWjZc1J; arc=pass (i=2 dkim=pass dkdomain=yahoo.co.uk); spf=neutral (google.com: 52.222.5.176 is neither permitted nor denied by domain of lawal_ka@xxxxxxxxxxx) smtp.mailfrom=lawal_ka@xxxxxxxxxxx; dmarc=pass (p=REJECT sp=REJECT dis=NONE) header.from=yahoo.co.uk

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@yahoo.co.uk header.s=s2048 header.b=HWWjZc1J; arc=pass (i=2 dkim=pass dkdomain=yahoo.co.uk); spf=neutral (google.com: 52.222.5.176 is neither permitted nor denied by domain of lawal_ka@xxxxxxxxxxx) smtp.mailfrom=lawal_ka@xxxxxxxxxxx; dmarc=pass (p=REJECT sp=REJECT dis=NONE) header.from=yahoo.co.uk

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@yahoo.co.uk header.s=s2048 header.b=HWWjZc1J; arc=pass (i=1 spf=pass spfdomain=yahoo.co.uk dkim=pass dkdomain=yahoo.co.uk dmarc=pass fromdomain=yahoo.co.uk); spf=neutral (google.com: 65.153.203.168 is neither permitted nor denied by domain of lawal_ka@xxxxxxxxxxx) smtp.mailfrom=lawal_ka@xxxxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@yahoo.co.uk header.s=s2048 header.b=HWWjZc1J; spf=pass (google.com: domain of lawal_ka@xxxxxxxxxxx designates 77.238.179.82 as permitted sender) smtp.mailfrom=lawal_ka@xxxxxxxxxxx; dmarc=pass (p=REJECT sp=REJECT dis=NONE) header.from=yahoo.co.uk

Arc-message-signature: i=4; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=list-archive:list-help:list-post:list-id:mailing-list:precedence :message-id:mime-version:subject:references:in-reply-to:to:reply-to :from:date:dkim-signature; bh=cvRh/cAKDu7YbNnGPdu12v4k9wuSh0DsEcLOaa+ID8w=; b=hRpYZB2sGEbxCWJClwFCGMqdS5jpSmmquK+S6e5bN2rxeQuGH6LRJcVTQ6dYPQURmN 5pw045rxMffxPQhUrEORqPaH5re4hdppKwSLy1uLbIM8KbUwUwtmKuVLc+5OBLU5JPQm w2gOBPEElo0i96jmb+h9sKyCBZtmIeWwWRQm83BcW9ZAuCmRucEsuin49FVqtpYaGaGC Ze34JTmVO9LHObFNbRGeU3A2GJBF5WJl0DhpbZQnCHtw//HpJAiyj/EG6fCWgAJDsD+Z D731ZnZEi3F5rn76YrnxU+rxFEiNf8BXf0o7KUyldLS2d401ucD2f6JMl1sEs1cQguNd vnEg==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:mime-version:subject:references:in-reply-to:to:reply-to :from:date:dkim-signature; bh=cvRh/cAKDu7YbNnGPdu12v4k9wuSh0DsEcLOaa+ID8w=; b=eNL3EbdYgJNkb3eYk/plf50+GAgvr5Kas+uEbbvtO8pCVADr6UHLNUxWjkm0CiXKXF GnM9uw1vuuREUrWcXMd3dR54SZly3l2g9zzW3Taw7PJ+Xi9wR1zgqOPDb2AXviL1ZJUH cCzhWAr8eAFu+03O7Il6NnNcoKzY8OB+OvTNnpYaOV4FliRUhci45yToMq6rzTQQENrg hO50tyy1/f+TI94mRp0knDnm5slxTPX78hJqtirh+KSXMjmeO53m+aX2SxEXuiRUcUSi qLWFw6SECLhHYYZAUfxSMh3sfQDnOOVmM5YIEuDx36Y1FuLndTQJbwwsdU56NpdP4Jg4 2vfg==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:subject:references:in-reply-to:message-id:to:reply-to :from:date:dkim-signature; bh=cvRh/cAKDu7YbNnGPdu12v4k9wuSh0DsEcLOaa+ID8w=; b=s8PTGWq1Frp5iueXxlw5b6GKAwYbr1+5dYJT+HrVClUFsVaibEQeNssanpMzURVolV 0tnwQPKCDLLN2NaxDNxAeqkaeLkIW1g9+xYWETcTTQTnH0aXCuuI1ySs4b3NUJezTJ72 zEKxStkIsdd+u8vIlEcyCJhuy8tPpK0PsA8JzSxzVhFbSJgMkNSeDTAAYeaVzc5FIN5R fGQop4Nq17fB+gG+rwnzdUNa6EO/QZu8EmnJYcYgiL642rIMvg5lcIMaD+6m3i+TM4Hi XIjzZSj8kUdMYtxMBSyFuo9NU8pA41waWfeydNVYZHLG7F/w3xo9PXGNFEUPxb5r7Uyb 6uuA==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:subject:references:in-reply-to:message-id:to:reply-to :from:date:dkim-signature; bh=cvRh/cAKDu7YbNnGPdu12v4k9wuSh0DsEcLOaa+ID8w=; b=urSaR/nBFEMttJ9wr3SsliU1nI2q9fThHBvjYepoFrpokYAUhthp4LmT7OZlYlPiQr Anli0lLekElNjt2KsRqMPzT7A3CE4kOxgYDvv+tCQCYSwkMFEX+jEJgayN/7km5yrfCQ PYhCdsOLnUYyQb6f+zL9pgbu0Pm893TlGHTq9IGAJ1MY3crxE8a59tL8kqODrnLuWuyU JnFTF8H/vy2CQg7H/QqknO9CmlDXRIANHOFBt34efXhcC0FSYPjZgsMc0lvZT8oLZOrv RrVy8X1FhbSAFehCg4q3e4PT+qZkNo+6YOTXvDDb/W0Lg8hBTX2txNXjujHVG3SwdV5S 0I0w==

Arc-seal: i=4; a=rsa-sha256; t=1566204022; cv=pass; d=google.com; s=arc-20160816; b=hQ8XdEoez/F3vwFx+lMLRtKTCDEXuV/ZWb92N8szfMVaDNjb9w4p7yctOLFqqnb+gA ss5gtkduKEZzURcoXh1jsQPOjIDSOdFwERx3c4lWL4FzLbVseJCEh6FQhT4KmQtwU5xQ RvAKNgrR923Pq+8Jmsamy37U2IvLvXwuXwpzBxX+R/bGnXpKIYhtmFxmaKctpML/Y3GJ 5cMtz/O0qQ5zYllW1fNMuvwOmKXSAQuSHxy+HY8koJ6d7S4iPqzYv7UWMRCJ9DQg2C5U airvvG4OFv4YfSVR6T1pfSx4Npaalzh9fF19wbfXdFMo15XG7A4h6mwHwwZrussFM9UT dKZQ==

Arc-seal: i=3; a=rsa-sha256; t=1566204021; cv=pass; d=google.com; s=arc-20160816; b=bzPNr+q6oNTht/KeenDsp30YFxz+hPOXkJodzkKkjGN77fGVkM9ZHVLeyLTYY2Z3mj etxaDgKq3ljbsPK03Mi9U3vGxSowojQOzofPi99f/H63v5XcQqq5o44ZiTeyAdAVvzJ7 qeBG5QYAc/q5L+uMQcjcEN1xi45aLUhf5R8kN7GUbfxs8DLXssR7bTpGJN0iwGHZ/B4v szfwpag8zdY+CvcwmU8iCmp+uJvIjdC3RvblEUyOvVXZNDVMUScd6aptmwbyltX4P5bM JFrFxo+ToN9bVEU38fnUHnP5q40AhsEy7NKB2zF+p+nxcRZ9JJFdAUW6yQpRXK+ogrDF jn8w==

Arc-seal: i=2; a=rsa-sha256; t=1566204015; cv=pass; d=google.com; s=arc-20160816; b=PUjv+IpVi4DF1GwjJmLDfuiwcBcJtUk7tT8FPQL9DUEXEylOSM5wo3NZ04WAhz3KbR U8uVW3SdI7tI3OZudzpms3MCT0OlEb/rAPLB2Tlx9sdY7zz8UlEisurxmdtCzojDDarh FPDrlgJQdIIuXB1Z5JJEbmM5jFoGAjAY14I9beaixYT1s1wYTt/Fr/poIx7bUU0gdqRl PCkFRjN2MLc43KXPEZ74DUvRrCWduVlGvxEkSzjphZia0bAMJwasf2nmPdqOydLw3JJB aHW6sf08ovBPHHAURwUhPlMH3xASZaoeGtk5NLN6BreHIT6imrVn4m/McK7dot3ki6IK qBgA==

Arc-seal: i=1; a=rsa-sha256; t=1566204013; cv=none; d=google.com; s=arc-20160816; b=KITIpGy6LvuH90zjkN090/k6G3vuty6qPnXUUzMf364ORFXdI2JEVX6Un33aD4XOso 4V/nM1OHVK+RC9TAxsXXsIZ6hliGo4qkZUjL/5+TRxq7KFYYa3Vl2fCfVa6HtV9LxXSr QAWucmUiwCgMa5NHDEid/wmnH7JopQifFB9UNx5EdnHNYqMLieFIDhMyCKQ7/keiIf+T ZptiObQFvNxBqt/v7OGC/lr3Bkh/5Q9ZmkkE68/H3InQw1eqrk3UmrATr2HYggl9pqnG uTgB2ODlQgDYPhju83e7O8w/dS15CLQ9iaXUiG4s8+kvBqaYn0C79fFZXr0JAYZcp45b f8uA==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=date:from:reply-to:to:in-reply-to:references:subject:mime-version :message-id:x-original-sender:x-original-authentication-results :precedence:mailing-list:list-id:list-post:list-help:list-archive; bh=cvRh/cAKDu7YbNnGPdu12v4k9wuSh0DsEcLOaa+ID8w=; b=N9wfdxujXrCjlzrm+qd1LY/xTyTwlKjRsEBFJDKVZNrzcC6L68pChhfN/n8cueB6kf cSgGWDy2aZdfGEnOVipN6ty9AKRcwTHjqfGvVmnP8aUIO9wy7X2Li7bkL+Zd8R1xp++y nvCymlen7TT19nlqksuuyNkwZJs4CQz0yCe5wv1JBxUL4n5JTtUxLgnZQjGHtbgmtgXV VyCtHpviDxI3A6QAEJvV2Bkk+Io9z+TrjVcr7Ai8DfFyJ5TatxqgCwwcXHpANaes26kY D+hwei8EhzIL83LXKPcMwvXDVvA+Plf2WbzVeH1MXJi1nmR7TuLW0WR/bmipNOw0uCww ukGA==

In-reply-to: <CALQtm8V1+OyRejBQFi=ngeKTSfT7_A6oZjB08vpc_TRN7bUh0w@mail.gmail.com>

List-archive: <https://groups.google.com/a/noaa.gov/group/ferret_users/>

List-help: <https://support.google.com/a/noaa.gov/bin/topic.py?topic=25838>, <mailto:ferret_users+help@noaa.gov>

List-id: <ferret_users.noaa.gov>

List-post: <https://groups.google.com/a/noaa.gov/group/ferret_users/post>, <mailto:ferret_users@noaa.gov>

Mailing-list: list ferret_users@xxxxxxxx; contact ferret_users+owners@xxxxxxxx

References: <CALQtm8V1+OyRejBQFi=ngeKTSfT7_A6oZjB08vpc_TRN7bUh0w@mail.gmail.com>

Reply-to: Dr Kamoru Abiodun Lawal <lawal_ka@xxxxxxxxxxx>

Sender: owner-ferret_users@xxxxxxxx

Dear Linda,

Check this out:

Re: [ferret_users] [FERRET]calculate spatial correlation of 2 variables

I don't think you need the statement "gt=precip[d=1]" because your L from both data sets are L=1:1

Good luck.

------
Kamoru Abiodun LAWAL, PhD

Post-Doctoral Research Fellow

African Climate and Development Initiative (ACDI),

Geological Sciences Building, 6th Floor,

University of Cape Town, Private Mail Bag X3,

Rondebosch, Cape Town, 7701, South Africa.

Website: www.acdi.uct.ac.za

Phones: +27 78 268 5434 (m); +27 21 650 2902 (o)

E-mails: kamoru.lawal@xxxxxxxxxxxxxxxx lawal@xxxxxxxxxxxxxx
.........there's no wrong time for doing the right things..........
<<<<----------------------------------------------------------------------------->>>>

On Monday, 19 August 2019, 08:40:19 GMT+2, Molulaqhooa Maoyi <molulaqhooa.maoyi@xxxxxxxxx> wrote:

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