[ferret_users] Linear regression on a vs plot

To: "ferret_users@xxxxxxxx" <ferret_users@xxxxxxxx>

Subject: [ferret_users] Linear regression on a vs plot

From: Emma Neudegg <fq18727@xxxxxxxxxxxxx>

Date: Mon, 29 Jul 2019 09:18:27 +0000

Accept-language: en-GB, en-US

Arc-authentication-results: i=5; mx.google.com; dkim=pass header.i=@bristol.ac.uk header.s=selector2 header.b=eDIqk0QC; arc=pass (i=3 dkim=pass dkdomain=bristol.ac.uk); spf=softfail (google.com: domain of transitioning fq18727@xxxxxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=fq18727@xxxxxxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=bristol.ac.uk

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@bristol.ac.uk header.s=selector2 header.b=eDIqk0QC; arc=pass (i=3 dkim=pass dkdomain=bristol.ac.uk); spf=softfail (google.com: domain of transitioning fq18727@xxxxxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=fq18727@xxxxxxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=bristol.ac.uk

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@bristol.ac.uk header.s=selector2 header.b=eDIqk0QC; arc=pass (i=2 spf=pass spfdomain=bristol.ac.uk dkim=pass dkdomain=bristol.ac.uk dmarc=pass fromdomain=bristol.ac.uk); spf=softfail (google.com: domain of transitioning fq18727@xxxxxxxxxxxxx does not designate 65.153.203.168 as permitted sender) smtp.mailfrom=fq18727@xxxxxxxxxxxxx

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@bristol.ac.uk header.s=selector2 header.b=eDIqk0QC; arc=pass (i=1 spf=pass spfdomain=bristol.ac.uk dkim=pass dkdomain=bristol.ac.uk dmarc=pass fromdomain=bristol.ac.uk); spf=pass (google.com: domain of fq18727@xxxxxxxxxxxxx designates 40.107.14.90 as permitted sender) smtp.mailfrom=fq18727@xxxxxxxxxxxxx; dmarc=pass (p=NONE sp=NONE dis=NONE) header.from=bristol.ac.uk

Arc-authentication-results: i=1; mx.microsoft.com 1;spf=pass smtp.mailfrom=bristol.ac.uk;dmarc=pass action=none header.from=bristol.ac.uk;dkim=pass header.d=bristol.ac.uk;arc=none

Arc-message-signature: i=5; 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:content-language:accept-language:date :thread-index:thread-topic:subject:to:from:sender:dkim-signature; bh=v0BZ/ZjKEGZQJ696hRQwFfSve+IoszrxYwzzCvruyS8=; b=wwakaYCoH5zRwRFhSJhmV8Ium7jM0eHsW0/eyjOsGCaML81v/zhFlV/oenFgbWxRrE oUlsGkPjwJX+eZyD70/32vV3MdW6RMdU0FBGMXqjOZhTck9iickXZHpgnR7RtdUhtfmx g0AzFkpHd5K2Vf9funBfJDAlWnJkhVt5hQhPHbRx9NYx6Puc38UQbNfKgoe7ZrLzxkTn RuaKFuTF2NYVc3HgGJhETA64kiBoFtyBanIEykwYOxcJEyplw98S65DqN5q3MXI4ppjL 3GUmLfsBeB+xr8Vu2IA8Ln5CSQdim1jbZyMnL3w7F5KCs/jr1JAknHwY4IgA7reA2AwU 3zZA==

Arc-message-signature: i=4; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:mime-version:content-language:accept-language:date :thread-index:thread-topic:subject:to:from:dkim-signature; bh=v0BZ/ZjKEGZQJ696hRQwFfSve+IoszrxYwzzCvruyS8=; b=fpRfaHXg6FK5D5LOW7yRdrS1Za6kz0o/W2sOxLWl2G3c4V39bFBhvxg5mvTv0+Ehcv hlK8QCwRegilFE2o7xsziHCH+pO4tB5MCKT/lqTQAva2cTcAd9CTBbScfkaPpNG/KB+M 7UzQYzmb7lUBojMxZ2HsCamqRGEj9/X2wpS7JAaZ6AqGAhM9euyq/X5ytJUz8cFoaPLw KYZ/tyahgNFGePwYOh1Y+ZFW9kllY5Ot9PsM4TDBP62eU+cy95Q2ET8XjE+ap2uMnYIT boO5b28uNouEE5Vs9kN29sWxXnYJ5n6OyR4/ThhYcIRZSMFeEghV0TiLGK2hvul3gXJ5 B8cQ==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:message-id:date :thread-index:thread-topic:subject:to:from:dkim-signature; bh=v0BZ/ZjKEGZQJ696hRQwFfSve+IoszrxYwzzCvruyS8=; b=hjsIXrHSgl1EDhtHGy/EJfEGLRR8WfoO2IUG1bdelfysOwqz9Yhrd5938Gi67EvCIv IWgaLGAOlY6u0i1XzK68xrYR1BBbXWiEFClV0F6dUWuxID/VkPlUBKAf8YZ2rxZObZtN XKHdfqKCE0+Cm6zJbBGoXFQrwKF76wmWrzGzD4UDvzvmbzZz6++TlCEAykPyLwRX57Bc XugLx26ibhkcSIuG2QCRbeXg0VNsbk7/Dc9EB82tdLmJTDq9JdTunoWKOAQA7YRuGm38 2RYzimY+EoUJABEyOEnTciCIYJPhOC2jn1bvN0M+1IoACAfmykWkzVOwGY+l70RmXZsH knmg==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:content-language:accept-language:message-id:date :thread-index:thread-topic:subject:to:from:dkim-signature; bh=v0BZ/ZjKEGZQJ696hRQwFfSve+IoszrxYwzzCvruyS8=; b=dusmIas6hY4mkdGpx3xN05ofLgI2NqaBattlF0jGm8dzuHjEiHm4JlBBnJxMEKdPdy qPzmzKg0q80cvm+BhwDUt2DKQLxG0oiBIUedq/2P3UtYZaFrs0iUN9cos5THnVtfdbqL DtwsE1+VrduqX5H1HM9gHFm8s2bLIBdup0vPreBloAd20B0VGGblxfLWwDFOw9yzzLHs xebr1X0Mw9o9TyfQN/N0su/QXfm/DiCmipQ7nn3WgjglJ2PpuMsN7VeG06nE0LecE0ZS HUFjDocNYd0fO3uOrO1uKxyb5JvU4cCfnXG4y5fZYnkz/Sz+ArJtD2/EVB8YK/s/UOYD dGlQ==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector9901; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-SenderADCheck; bh=v0BZ/ZjKEGZQJ696hRQwFfSve+IoszrxYwzzCvruyS8=; b=mI4CwHfKtaEwRepf5805RsJDEcGpos/a4SEdj8MDP0824uoCtEPlLL4Fp9aCOjN+Diyterx7K+XddA34O8Qd7HJ6TbdBn9hwi5ddEG4jRGub1MigX9EGX+8asa7gaHflDfnXfXlFgySafNkXevFHDjW7HiYYnRPI2nkmAPkOs5kjD4JV4QX0EdLwGYx+QVhBxxTwVDMaE9/o8DnyYLd33suD2Wx+z4FXWJi1eK7ZVIJ7F92GazcfJOL0jnnIztDwDfePUwZPZpiM0hB8B/AAWTpFP+vUn8mzw7ueBVTcTMVAzaOlXWIhuZ5zU38m36iJMQegClnOakuXoe97VFimyw==

Arc-seal: i=5; a=rsa-sha256; t=1564391915; cv=pass; d=google.com; s=arc-20160816; b=BOMDq0UfLohizVo5w1jiMxe1HaYmeInIr1VkyL+AsTwgJ5d12NeCU87ddy+PLTKaNE x3f1KxAf89t42peEb+ePzyc9pm3VPunIBV4KyytiY40n1ivMK7+DVGVHZbSPusZB/+tz ZNaatBupa54MQxUVVseiGyKm01Bxy9XqxNS5gDTD20MonZqTLyEcei8EIyG1bXvln8CK KQvVeaTUiTm21JU7iNhMhT7I+GDxn9Nk3TvaRtbKz8Ff/7dMdjelgIY1kAgEASbkMNV3 T92cSZeN+3dx9buG88pVC8FhfcEuH7/pvj3jkBGouVnhupUNBLNSxdPI4eAbUbJx2rax XQ6A==

Arc-seal: i=4; a=rsa-sha256; t=1564391914; cv=pass; d=google.com; s=arc-20160816; b=QXJpcSpnwcMFP6w6KyCNLAmOMn+1QgupnAYpwMriFZix3beIe16MnK574M4aKsCS2e f1JdolrLyvAwAcqhKp4bVCLJSmLzuvVi0U5EIC3lTSuhSXrQxrYjwg1PwrcYkC0rsG4l prfRK2HdcGo4ygZ9vLzQipDWuM92kOcmBowR1Uulr2+3lxNHiJTwSESQ6e23SDBm06jo DAcjVU3hjitk5e/oL1E8oirg0ADMl6EkM39O0RG4JJepUjUJE/EGgl/g2jbFKnezPvJG VEd1AJuog+I0IVnoQraetiXt8K8xMRhlcAXAA8Jpvz9CZcTjUfYKqV2zDwA84o3m3V0J FALA==

Arc-seal: i=3; a=rsa-sha256; t=1564391911; cv=pass; d=google.com; s=arc-20160816; b=vgW7ya1bMoAWXozUm6/iDubDJ6uo5ZcK1jloESzfTCuptYZRbACPOvTA/+OaQN6mYE cVfh2AOAhUMYTFLHBUtYLb+V7zLB7C2fNW3Mm0a4msGkb4qPRIT2KPsuldXsYda8VgQG 7vHPiRqcutUYtXc9oOgtzBUKlRBR5K48TglzolHBsp3H2rT+8VaIP/YX2ElkvFZsbPcN Qp1XhdbSY4mVwzNjSjoQFFg789DGldelB2gQuWzU2ILlaF1rCsGCfmMzdG09FR1NQfal 6ghvX+W4QtzMFGQQrnH8T+UlFaski1ERQOu/OFHty2tK2r8PPfWsjkaIAt3qLnrVqjNr wvSw==

Arc-seal: i=2; a=rsa-sha256; t=1564391908; cv=pass; d=google.com; s=arc-20160816; b=HY09FXlPPYB+Ptwg4V+ShALfmM1P8xB9EeRRxzqnYmnYMJv00QakC0blyCb1LhJLla B34EFltbZ+QJ33xKcfvKAOZDlhyK0fzuI5GiHCsRbQvKs2Lr7+0rdDX9Jz/AhSv1qZtM miI8ykExDS4zlVOH37o0sjeow26MzaFuEupDb7Ot/mgBnbZU9t9lGrJddaTYuUZI0kqj 1iEj2SABTUMeIcQxsXKXsHxfiQei4HXh5u4ymvLmzgkJCEcnBXmhmi+NSuwvkGGotuc/ 1JOzOJvCuoJqFCraS8f4Yh5IZ3k3GwGiNFWtIi8Gmf82e9ghMLpixWfnBDJ+hCAbEkK8 mPRQ==

Arc-seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=OqKQOUqDsTPND9JsXOyqC3WEkn9PuUb1asMN36WbcX2yiCVrAftlwu4OPD29wD5aDz8RH6WwQB8noIY7RFbWahofgXN5PCxks3rqU2+El9lKY+LASirqzKu/xykLqpKh8bPp/cwAWxllv+SCg84lXXVzRk8say1gruR6YCQADI4nQXjqtFauZm3wjTUQAf15yUpgUTzig5ZKmNkV8uzDHns2TeQwD6BSFpKDuYV4pda/f0/TlP/9UWHeWwdQFX6pokXgixMj5zPijfLv2HeJ0T0bWSikg+gl5Fo4MjUkQgb8sg5wi8IBkXjly0AeXHbd+fj6DHV8i1A3qezaJ5lwLQ==

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:from:to:subject:thread-topic:thread-index:date :accept-language:content-language:mime-version:message-id :x-original-sender:x-original-authentication-results:precedence :mailing-list:list-id:list-post:list-help:list-archive; bh=v0BZ/ZjKEGZQJ696hRQwFfSve+IoszrxYwzzCvruyS8=; b=nI93hpIHv1fQ3ttRoy638tUZIChF64GRiA1dG5+pi8/RP+mIpLnOVgnnC99532k3B3 JleHsPnzPf25v2osUwc4B0MgNGTF3yP5+i3KVNoPSBI6NjQ5O1wpCQxQjhNo+DnqhSd9 hvqgLg397VzTPpllKiL1x9ZG0EW9ApMHvPcBp15j2nMxdSL6c89zcvGdxDs4rFSrdgT3 3N7LAQvj90QQzBzkA9YWQ/8DW9Mzvuzld09VT0JYVFq8oyQSjXBlNHX+GeJ7cNUE+UrD M4lnauuNOnR+Q9PK5SnHpVStngBB5Ib2wsL+QqdK83s7JTLUhnqZf4EfiQehTxgyWdO+ aoQw==

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

Sender: owner-ferret_users@xxxxxxxx

Thread-index: AQHVRewNFLThSV0bTEuIOc8OejbC8g==

Thread-topic: Linear regression on a vs plot

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