Re: [ferret_users] Linear regression on a vs plot

To: ferret_users@xxxxxxxx

Subject: Re: [ferret_users] Linear regression on a vs plot

From: "Ansley C. Manke" <ansley.b.manke@xxxxxxxx>

Date: Mon, 29 Jul 2019 10:28:48 -0700

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@noaa.gov header.s=google header.b=aCg7MltZ; spf=pass (google.com: domain of ansley.b.manke@xxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=ansley.b.manke@xxxxxxxx; dmarc=pass (p=REJECT sp=NONE dis=NONE) header.from=noaa.gov

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@noaa.gov header.s=google header.b=aCg7MltZ; spf=pass (google.com: domain of ansley.b.manke@xxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=ansley.b.manke@xxxxxxxx; dmarc=pass (p=REJECT sp=NONE dis=NONE) header.from=noaa.gov

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=list-archive:list-help:list-post:list-id:mailing-list:precedence :content-language:in-reply-to:mime-version:user-agent:date :message-id:from:references:to:subject:dkim-signature; bh=KncZw/KVZLanAtIy+NXUMsOgszv5TM5e2rVFx/7v7oA=; b=Kot9R8asGCFNEUq6X7K9ixaL+qP0of8N2EBxfkBgqvDLFKX36RmHZNdR711GT3A8HC uZ7qPAgMCBqpGDGXi5+vFHFuN0kAdy54TPri0xyOQCqmsjNR+my/7X/4Ff9NeVjYFDiT 1ewo/lNuH+aplpF0fj3RAWxySKej4IsqC86xsUoASoS7J+IYPgJWD2HOQl4ygOe41k2i YtiP6+EpVP4SRgN0scuNds3HaTiFknYHoWQIhbDxhh6NvaURZwi/RNPKZ4oDuk4iRlNN irpepgsiAD0mguhIznnBpdGMXno6GddqJonyWJ5W20HPGPqr6orMXzgFGbETD2c2t7AA Zydw==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=content-language:in-reply-to:mime-version:user-agent:date :message-id:from:references:to:subject:dkim-signature; bh=KncZw/KVZLanAtIy+NXUMsOgszv5TM5e2rVFx/7v7oA=; b=EwCwKGVzCQ97RMF/pw52ELFqYLlr8HKg0/Jyat1snbbgYFavr1L/viS0n9rkN4ioFW kLFu8ELvDXdUH/qfiFqtbAapAGMBYV6zY6MPsG8qRhW5xufIMp7If0O+sCQdTkMDtgqE BWYtJKRr4MKI9iVsXSaivqgPlu50KZzWUIr98plYuZKQkXAsAVbjaqFEznSEbovFLLSn nzwJZULtifoC8UMDKSWd/rCHXBNEpxtJCWieuaepjrCH2PD6U7H41S3aZi8TbDKa235n pzA/KKNZtvRpLsBBY69yVwmkNMCkJhKE+jUJm9z519I/t1wQGvCEuRcjU4RLpUfn2lza qL/A==

Arc-seal: i=2; a=rsa-sha256; t=1564421331; cv=pass; d=google.com; s=arc-20160816; b=bjbBubFd8tZvRfzPeaEH42wjSpDyK87jWYbkM/h4KjR+U/sM2o7ws2A9XTRnE9+RBi /CNHRW/3OQ0gK1UECevsikAsJrzpZn6tmGh5RLuB7O43f0krhYTqYDYY2/9G9MBwArGA yGAA0Q7A6EOrWLvZYI+Olm8EdmPeyl6eBaurDJA1aHwXsFqvw6JfJ0jt0RLkFI9AdbnF CfujeHYG75jZOQwsWZwtv1d4DheFWX6RV5zn38yyegenvMwOoKW1dr0pc6H9TLOvWuun AjMsJ5T4n/9e+4ZC48jJWe/5Yr5xOTQD3OscTVWtMrty+mSu31T5hHKAqSJf95nydYUG tN6w==

Arc-seal: i=1; a=rsa-sha256; t=1564421330; cv=none; d=google.com; s=arc-20160816; b=hZcuywP68o4nv8xW80W1OEcribkIN89qZ15NpnXQNp0IlXYS0G2FDn1GRWYFMc5T94 dCqdZmbGh0oU9smkSd+eeSacAK5E2bC5DYlivxVPt9D3tvn16IOaeSsY8B5gJHoZ7bQB UeEVB//AS82QvUu3C1LDvkYLQmiyAI1JG7emTr3huwIln2K0DAAOz5wzzL4z+UhMvzPx /QvxFg8fevjYbvw7bwHZnBDe2IVxJmipwD+s3VsLJCwzjhJ1duKRwqkKAhJ0LakWIJtK v3V8L4C+LdFi/MChaRZPjnAzyvbvY/dAf+DOSOzr9V3iv7GIwTwRi+/a9rnEqguRnxJ5 RUHQ==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=subject:to:references:from:message-id:date:user-agent:mime-version :in-reply-to:content-language:x-original-sender :x-original-authentication-results:precedence:mailing-list:list-id :list-post:list-help:list-archive; bh=KncZw/KVZLanAtIy+NXUMsOgszv5TM5e2rVFx/7v7oA=; b=OFVDwyIVrS2GVwNA8LP6UwKeWMbF3dPwn+TvwRVe2qo0YWpMd7rmcYATPDIN92QPwL cBoKn7ggOIGypT44sZUAvFT952EMBh8nx7uxjupx1r3OKXNt8Z5E3gbWgZT6Us3rkjfh Haiwid3/M0uoFIEln+rxa1kOVpx3gm08/MuHG06gUJGyvT0QH5qy4sDKjvhBDUJfq/// 1VMYnx9i/Uoy839+MZEyC7RqhXuqOU2gMYpm1L7j6gNNcJhXL+vVFXwJcKKvgbuL4UT+ yQsGFocYQZjA7wbFS1fyFvzGtQSwYGaSVZtMb1kMPHrn9louwl80TV1rXYdUsL5Svebc MA4g==

In-reply-to: <AM0PR06MB5908067A7D48BD43A9F6FB7FB7DD0@AM0PR06MB5908.eurprd06.prod.outlook.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: <AM0PR06MB5908067A7D48BD43A9F6FB7FB7DD0@AM0PR06MB5908.eurprd06.prod.outlook.com>

Sender: owner-ferret_users@xxxxxxxx

User-agent: Mozilla/5.0 (Windows NT 6.1; WOW64; rv:60.0) Gecko/20100101 Thunderbird/60.8.0

Hi Emma,

You need to define P and Q using expressions that include the averaging so that they are 1-D Time variables. Then the regresst script will define qhat correctly.

Yes? Plot/vs mld[i=@ave,j=@ave],sst[I=@ave,j=@ave]

Yes? Let p = sst[I=@ave,j=@ave]

Yes? Let q = mld[i=@ave,j=@ave]

Yes? Go regresst

Yes? plot/vs/line/over p, qhat

On 7/29/2019 2:18 AM, Emma Neudegg wrote:

Hi,

I am relatively new to ferret and am struggling to plot a regression line on a vs plot. I am using annual mixed layer depth and sea surface temperature data from 1970 – 2099.

This is what im currently doing, following a previous post https://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_98/msg00121.html –

Yes? Plot/vs mld[i=@ave,j=@ave],sst[I=@ave,j=@ave]

Yes? Let p = sst

Yes? Let q = mld

Yes? Go regresst

Yes? plot/vs/line/over p[I=@ave,j=@ave],qhat[I=@ave,j=@ave]

However plotting ‘p,qhat’ doesnt give me a single line and also doesn’t look correct - shown in the picture I have attached. I would also like to obtain the rsquared value for this plot and am unsure how to.

Many thanks,

Emma

--
Ansley Manke
NOAA/PMEL Science Data Integration Group
7600 Sand Point Way NE
206-526-6246

Attachment: regresst.png
Description: PNG image