Re: [ferret_users] variance Calculation

To: USHNANSHU DUTTA <ushnanshudutta@xxxxxxxxx>

Subject: Re: [ferret_users] variance Calculation

From: Ryo Furue <furue@xxxxxxxxxx>

Date: Thu, 4 Oct 2018 12:09:12 +0900

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@hawaii-edu.20150623.gappssmtp.com header.s=20150623 header.b="xBx2eu/E"; arc=pass (i=2 dkim=pass dkdomain=hawaii-edu.20150623.gappssmtp.com); spf=softfail (google.com: domain of transitioning furue@xxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=furue@xxxxxxxxxx

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@hawaii-edu.20150623.gappssmtp.com header.s=20150623 header.b="xBx2eu/E"; arc=pass (i=2 dkim=pass dkdomain=hawaii-edu.20150623.gappssmtp.com); spf=softfail (google.com: domain of transitioning furue@xxxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=furue@xxxxxxxxxx

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@hawaii-edu.20150623.gappssmtp.com header.s=20150623 header.b="xBx2eu/E"; arc=pass (i=1 spf=pass spfdomain=hawaii.edu dkim=pass dkdomain=hawaii-edu.20150623.gappssmtp.com); spf=softfail (google.com: domain of transitioning furue@xxxxxxxxxx does not designate 65.127.216.227 as permitted sender) smtp.mailfrom=furue@xxxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@hawaii-edu.20150623.gappssmtp.com header.s=20150623 header.b="xBx2eu/E"; spf=pass (google.com: domain of furue@xxxxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=furue@xxxxxxxxxx

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:cc:to:subject:date:from:in-reply-to:references :mime-version:sender:dkim-signature; bh=vlKk/2AVLM5up7n5p70FBqx2iyTlhux5a9PH71L9Tmg=; b=GUA32ADlIarqVJpWWxPs9T0Wpj1LCQUyJOsmBor+FfYfXaoc6/iAn+KcXGkSbgVrwv nctXq7iHRyJz9NZMZH99KzpLRFx64Rfxt69FrocurxGrDxwQnow0IcS6W3wIi3T1jG2S YCfw2084BZCIfxnVcRDak2vyWZnK/xL1npzCtqtpD/9qAQ/tVfgCLMS+GnDr51qDldoc E7sXehSfnRN3CSrmhEkBJC6aMLthIwbRS/bYWMjvnsQdLl555aBaIrFqEZ6jdscnFAth Foxw6rC0vBUmMulIrYoXaKoPrILyh90V74uexw+WnpwHSJaN2RwPwlQ1bFNOhA1cjEFn 8jlQ==

Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=message-id:cc:to:subject:date:from:in-reply-to:references :mime-version:dkim-signature; bh=vlKk/2AVLM5up7n5p70FBqx2iyTlhux5a9PH71L9Tmg=; b=dgPGaUTQDIbT4YKlXuEA9LyWnHvPpI0YKTrpKQMdKLhW+QlUDjBPapie0cSNlTUS3i u8xO+r7jlQ0KRUEEuhkGCJMkbIOQE7v9ndhHPvTV+kwRv/O9pspNouXeNlVeINBD9PF9 GOqVoKtnIgLH0xyXbObBVq1TaBbGEfwmni6lm2c1ZdOn/TyumgKdQRXA0A7v4lVWXIS8 MbEAIVg/pFdK+BuLvI/Ad6IpxpGPoZq5CxpOfFGxprr77SqWpOxsUu2HiKpw24hoT86t hW31B7W16qWDqZ1+i4e5Mxk8pTef/jilIa4SjWnTuxyRPpG5Xa1XJH//bhlp2XJOGf7Z +ylQ==

Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=cc:to:subject:message-id:date:from:in-reply-to:references :mime-version:dkim-signature; bh=vlKk/2AVLM5up7n5p70FBqx2iyTlhux5a9PH71L9Tmg=; b=D3TEAAQcGqOJRRDPVAp+SpaUhbizVqstsiiq9P0c2pLroLao59aBb1mg8g1R/exnz5 hdzPMqom4RLGxw5tY7xwWRQ0XTIysJgsrG8HnREhMkoufdbgLc9CkFxjYd+3oDbWf2ic 5qQepkAXXd5bK46971oV5ii3y4rh86c1k10gDuYRCTe1f+RrKH5eIclonAFyrfS4C45R SDzRkuOST9bwi2YbFNDCUReBYCbJG+8mwdjAjIguIcxdE+6KspHhf8V4c4UTwB+fzm1i HBfLTV/HNsIMEn/v/H1ct9mWfG6f0ApKJmGNggxuAHMoLmVKM/qFyGbI7e5VZPVwjnd7 TwHA==

Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=cc:to:subject:message-id:date:from:in-reply-to:references :mime-version:dkim-signature; bh=vlKk/2AVLM5up7n5p70FBqx2iyTlhux5a9PH71L9Tmg=; b=lzF5oU8gcO5O/OIMh/AJIueWrG+fK5ytTmf4HcuMi9Vxr/9/e8Y+i8hjTQsVpNrp0C FhvUuW1fM1PBlDgF1Crdq1VAudHlRUrN4L18xLoyfqjM1G/juhxT/DMLhqirped+xUM1 Zp0atFignJeQv2cnN9oFxAqKU+JDpgg5hcEafYN1Ftjq+W8WNztpr86NtvwKP7HWbgI+ OB/LjDkZssZb4zRhbAwJghNRg/jNVBIqsXys27M/4bPu39ytCjjYH+zR/9v2q3Idc/sc BKYbc/BdqqNhbsMKhIfMoTp+qkjnG4DqAHyuvxVsf27k24xMH1nZoXlzrJr2XeP9rQJS 5Hdw==

Arc-seal: i=4; a=rsa-sha256; t=1538622570; cv=pass; d=google.com; s=arc-20160816; b=jYlo/GNxLUhMYBVhhUQLTLYJ8PhtW+jGdPnVK7CtePCvLEMc3Wcon2g8bDokcrODwj ninMgZti+GDJEU7QkNMi7Wbf6YqqFep37jBUR9kDQf0kIkppeoSkBcxxg8a9hdUxVCq+ Jgk1yxoRTkrNsbGCNm3e93M9QvIgOXiMXlzAx0SctNOjgV2tppYiEQHB3ZQcZDb1sYT2 tu+cCLfVW7nL+dSAZ9OuwvpRD8Qp5lz33+WLPzm/aIiPwE3FveYUM1uTHJ+7PhBxXCsT eis/RbPO3acL1Uu/kfHz7hTAPtICRlPXkgSQeigXZo93G2d3j1zlnOk7zGdXjl8wSZ5V q4HQ==

Arc-seal: i=3; a=rsa-sha256; t=1538622570; cv=pass; d=google.com; s=arc-20160816; b=yinMd131BKUYkyieqCwVVrz6+Eiq4238I7rKHp2JuumqftapKXZaiiAzS8DO7luccc vawOtg+nnv2oQbXXmblAKlQIO9Yqoleu6UQ1IWKXOZH1glc8gZdp8fCmnOfQWWKR5NXN wbLse20OSDrXPZAM8ihEJcNIRbEvwhEwvCdwSBn1oCKSzscVLCciRyGYB/Ja0JybSRkc PbP9BaOKBlZ4uesV4UC7k3q03Frn9VSEwM7rFglIqX1UcREBFsOuE3Go1rdYXHxSP9Pn jKMjOgqbTEvdf52NwqgAQxavVmi1wFQJtktkf44iz0JleRt5+2I4Tv/zGpPc1cEIUlAL Xhmw==

Arc-seal: i=2; a=rsa-sha256; t=1538622568; cv=pass; d=google.com; s=arc-20160816; b=im+8p9HYGfYzOqyNYPreNbmGM7xmo8bHO6exJEzFfX14ujw+HfcxxDvWIO2/OkZ773 5zq9Vewp2Ks/pmPf8vbTySRmIt6IHwLGLScKnzoHzaIrMpFqI8dDEJnMHimJrEAqoT/F YJOQGtegsv9CMt667WdZS0FzRTZcjO0rcPOvAigkw5Qv2ZGHYlz4N9/TvlJBNdD7HMKC 1ZzZUcM43rEUy3SEPoglCLN3Ui6frtvZ2v9SyNMkKNlKc/BEMfhVyi6Y/zJzfuwXMffK pywaKoDwmJ6++icbmO0rVmPkPuCD1ccX+ttJkg50n7oOxvFzUqEvvejTznOG7wtVWUnc hMEg==

Arc-seal: i=1; a=rsa-sha256; t=1538622564; cv=none; d=google.com; s=arc-20160816; b=XVzO23J+E0etAIoxhAHJ7ZuTUgggChnOz7qU7bWi/oqT5d5w+Bbtk92wPZ/Bi83IMZ Wo5rx1Te3XmozyuNYpTTHR4e+8b0xnExRzM96VUJxM8vGsVQp/n2s9rVXin2KxZTLDBV X4ck5ul0L/C7kiV27VEI9XqQzQav8lh6y3m5jn6XQrHN8URCc5JmZ8R6nYzNkDXLFvEP 36xhskS+xwprS1SU/9goiSU9IHvkq9u5/aGyjLWgzKS20kJuvh2vZfr+Lp6XvsI7lsrh E7ruZgZ+EH/ge6WWS+qgWtAvaisr70IgDCe8iTKuVKvj1lLcpQ7msbu0WheWxSvFeq1n 3w8Q==

Cc: ferret users <ferret_users@xxxxxxxx>

Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa.gov; s=google; h=sender:mime-version:references:in-reply-to:from:date:subject:to:cc :message-id:x-original-sender:x-original-authentication-results :precedence:mailing-list:list-id:list-post:list-help:list-archive; bh=vlKk/2AVLM5up7n5p70FBqx2iyTlhux5a9PH71L9Tmg=; b=pBTQymzE5K6B1fewHeNrii75bMJZBUAVWPcg5rdZtKntayCAZY52IMRHM0lMdsED1s dBYo/khxUhuru3uzukhUaGDbgbAx51GRM1FnL3F/RQL+Qg0EaGeIREEp33XGZhd36TO8 waF0Bk2dm+upR4dUY3sv0/ty/i/rTjYfd14QKETwgzjE7UWK4OlTNizXbug96awVhtsL 3Lig+MBE+RcOMnfywbX+S28ViR6Vb4iTGMbtRz439/8uJPu8P46VUy7JOWTHpFFCk38X gMatHFr4NVmWszBuiMbGWs7dCvg/rcynP++klxdDeRfbaV0ZoiyjFie6sQrhBHhdFhRl O4ww==

In-reply-to: <CAC67f-Z1XOEUZk-3uARZ_ptjF9eapYGqh+4GBNB2MXaCzjjbfQ@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: <CAC67f-Z1XOEUZk-3uARZ_ptjF9eapYGqh+4GBNB2MXaCzjjbfQ@mail.gmail.com>

Sender: owner-ferret_users@xxxxxxxx

On Wed, Oct 3, 2018 at 9:02 PM USHNANSHU DUTTA <ushnanshudutta@xxxxxxxxx> wrote:

Is it possible to compute variance for each year in a single scrpt in ferret?

Here Andrew shows a method to calculate a "running correlation":

https://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_2018/msg00306.html

The method includes (or implies) how to calculate the running variance. (See my response to Andrew.)

The value you want is the running variance at the center of each year. (I haven't thought about how to extend it to cases where the length of the year isn't constant.)

Hope this helps,

Ryo