Re: [ferret_users] cumulative distribution curve

To: Ryo Furue <furue@xxxxxxxxxx>

Subject: Re: [ferret_users] cumulative distribution curve

From: gopika suresh <gopikasurvay@xxxxxxxxx>

Date: Mon, 2 Dec 2019 11:52:09 +0530

Arc-authentication-results: i=4; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=amArzn4X; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=softfail (google.com: domain of transitioning gopikasurvay@xxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=gopikasurvay@xxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=gmail.com

Arc-authentication-results: i=3; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=amArzn4X; arc=pass (i=2 dkim=pass dkdomain=gmail.com); spf=softfail (google.com: domain of transitioning gopikasurvay@xxxxxxxxx does not designate 52.222.5.176 as permitted sender) smtp.mailfrom=gopikasurvay@xxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=gmail.com

Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=amArzn4X; arc=pass (i=1 spf=pass spfdomain=gmail.com dkim=pass dkdomain=gmail.com dmarc=pass fromdomain=gmail.com); spf=softfail (google.com: domain of transitioning gopikasurvay@xxxxxxxxx does not designate 65.153.203.175 as permitted sender) smtp.mailfrom=gopikasurvay@xxxxxxxxx

Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@gmail.com header.s=20161025 header.b=amArzn4X; spf=pass (google.com: domain of gopikasurvay@xxxxxxxxx designates 209.85.220.41 as permitted sender) smtp.mailfrom=gopikasurvay@xxxxxxxxx; dmarc=pass (p=NONE sp=QUARANTINE dis=NONE) header.from=gmail.com

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=EwoelNv3Nwv/yobQc2D6ojSiQEEdVDIiJtDemjyON2w=; b=JjqHsUGz5ILnS0I1VQ5BmOz6IT7s2tv/fHZxsL8jTI37zg8rT42En9vBYYsQNreyjA MlplIas/TaAcSZfwvaw53+TOVQ9amXu5bK491LO7by3cPSzJxOqWxkr0x/P7wCn0Mbgj THkFaoIEkZSx2erqulyHptgRrzD4DMUaC6MWq8hJv3NSlx5eL9JdNPzKvdFAqh91ekry PGcG+qNsPpKrNJvx73gve+nrvCHjCQ9FWzInmTyh2YoC4hf7QONzUe/9lqzB+wlego4g +E2+Xr2uWmY/l/UGBoowIMHOWf8fem/ZyOPobnbIhiGO0A2+IQHugmh4Z1FBqXzQJF3m xlOw==

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=EwoelNv3Nwv/yobQc2D6ojSiQEEdVDIiJtDemjyON2w=; b=k77ZXsqXCE0aTOJCBtNumiG/11uBq1f22Ra38IEQXULR6PLDFXNZFEAo/Dxl3FRaEk pZiSUEGXi/BnkqOyyBEfaEFibEpEomLuFysfXKxDyULUPQmjQherEAqfeJqBdqOjQJzS omat3ua0PUUKnrPmGg+Wj4CcIeyvGtvWwQ2S7+7F2pXDREQbRHC/lgUeUZqdMy7S2Cjs 6exIKvqUsfEO12P0wXSdeJhQjAv3mYv84sqsWiuIiFaDU2BWcjDgnRaBIgJvZQFlDYJg PRXZCcU+rthH2FTcZ3XMD7nqID/9YyjX8nKcT8n8cO0HoHm6RK8eBLU+ygd9YxY65JSz MX3Q==

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=EwoelNv3Nwv/yobQc2D6ojSiQEEdVDIiJtDemjyON2w=; b=VcfnOqsLp2a5hdXIUJQIz7hp36z8iRXL5zFubAzW5vdx5KAdq+r4/OXPOo3llg4yZi T8mQg1pjLc8PkJbbyB9yDhyPVN/x8w23wp1K/eC8a3dw0hfQbYrmaPdwoL8Zl5S77MKK 8kkRjXXhBE2dUsPYsG4WVtTsL+wPdaqidqkzdSRA6CgJuvNVdHtBBTVe1lOrMMwM8H76 j+C+hTCCUIls3LJS6ZQ7brPvi1muzMcGNKKhQNHsfn2IwHa4R+T2iMedo0H2WmimdAIw TqqaEepGlylgU35AkXf7dMr1NmXNZzQPZZ9w7OQxovzBkjBroZs5yP519V6JOsmeAuju 4Pgw==

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=EwoelNv3Nwv/yobQc2D6ojSiQEEdVDIiJtDemjyON2w=; b=lqPSu0n8RUMLKIR+Mz9BVbXsqnNiQfc9r53KFTOuyi12HxlV/MgeGFQ7xNVEsIbrRK kSuHWmhACHGNLxGofYiEQ9Q2vWUXXRJq7gNqi9+xnlYmMwB0QjdmM40G1+rx61DiLtTL vN67/q3ORk7xrl0XKlDZBxhdXk4icG7WCOmaQXtES1AIixqap8XzChcbp1SbJBS+Fz4u or+BM4OPNRccBOHK4zA90Dik0W9p8DkGKvbcdnArQjcD0aAA9Ehz0X0oo52cOAbcUdky bHapxbU+NOCZjT8resCzgZAiAe3RCMlJze1Q3u1MZfjNyQ03Q9KAKnTtctsg517q1y0i jOXw==

Arc-seal: i=4; a=rsa-sha256; t=1575267779; cv=pass; d=google.com; s=arc-20160816; b=VaZ1i2LREu29QiKpDj5+a5wNqM6V6hMo8/TJHJbio10sz3IDU4xWqY30QMQM0rM4KC RDT8BwHUHqzhHqpWDg/KS36WyScUACYchm802fyAgzwOtsDayoBldHzjCCpojzFxh9ca lml+CvstBK8x7g07SUvYs1SmzJhRcdl1UzMXnn2dMU1IhsHsXSG+EMitvxheQ7dn3VJ3 +4uk6Cpc/NDp3EwkyV2bwvX2pXvyoeyCauheoFvdpovjyKZq9GiFiE/0m87PudWys1Xv /3XkMXvcNU9J9643Lfcc3rm97GzT+a5VIRWm1/++2TeCKwa1J8TXVTFfHFAWTN2op+aL G/dw==

Arc-seal: i=3; a=rsa-sha256; t=1575267778; cv=pass; d=google.com; s=arc-20160816; b=Y39+k7FYVU1fnrPUoM2Il/e9X9ZnYcwNwPIcAYw+cR96iirfaAhVxfp/4ozMV3kLhg emL1iNC/gtjQWv/swn8rvL9mNIJwnLAPKqo1fA2+QNt4zYmLoJYcWbtBrnL39sPiwJyQ gYQ5GvprqEyXHTzW/Mfx9cwUWokIvn7ZWNb6OTamh2qLTxYCrbwHbTjX24jkR2tv2TR4 WvMuCgTTUnrqXhrpDF76X+J2tYSWJjmdVJSCtlFgq8E9UBNSMSfojDmRE4TyAewwbIBQ ku8AxiXGwdVchfsR6/I+FNJA0/VChVkOmyJjwr7mouiE05PckMxkMOGfT4sA4cD/r3Z+ 3gaw==

Arc-seal: i=2; a=rsa-sha256; t=1575267770; cv=pass; d=google.com; s=arc-20160816; b=WNz1oQ4K6lEKnrpG7m/ElhPAKrR/a+DbTEVLPnGmLy5hNCLLnMWHD5tIXLgKjJo04M HULc3MwwYxrfKrOcp9RWEOgyW6JBklFb+XfJcv5/uIgrb8004HeoMbQyvzaBcjtF3nIz H5aQHqjDwY4E5X8EUHsfnvy+4joW7OS508dp9pr/ixH8ZvVk4hhQEFY8ou3+wIaTxNCK rZmkCw2UA/Gmg25LqtqMNmgzxZu3w9EAKBZDhCe+LpHLac4RGkD9NjyoaFYGQTEHO+mS VUYgRxfTCF2T6TEDyPuForJRmUBcMNSNW9L/HkwOagIpVcBaSt9S93N3iVJyhZEuAkLJ iNEQ==

Arc-seal: i=1; a=rsa-sha256; t=1575267766; cv=none; d=google.com; s=arc-20160816; b=zxOpWA+ebAnnE/EuCD84XL9MCKRM+moVbjV6gaV0Xzo96Z7vAByL5SDqIViacvbH6W THSnSNxZuPio+JXefPK/eY+siGS6fGOPbQUnzQTBAhuTAJ25kNrcYK0hCmfAO2671dnf QKEnuhf0jPj/TSjQ8GOYAv+37orvlH7bAhdPw7myZUV4M04fmXcvc2CuQ711VqrHbACu fNoknzo7aFK3PrPi6E1ZdkQ7heaqD5MJLe5Yo7axIv03GSC2AjAeVTMgdcwd+elSu5fs bls/3GFjiMX80vXoH2nr0YQWh7BYra8R0GRUoqulpF03y6QHv2Uv4kI0xtZLrxBY6dSU 8k4Q==

Cc: "ferret_users@xxxxxxxx" <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=EwoelNv3Nwv/yobQc2D6ojSiQEEdVDIiJtDemjyON2w=; b=BSybxcYc6JY3HMTQydtuj7oNvbR0OSA5mliTDFc4UL800FA8MaziJeLgbAJf1KGVDF vWAnhUNNjIVAyMHPALEC6dANUGQgYJF4+qBYNAKgx9x6NVby3DymeD6yXoBO6MhmOxL2 5neqpk4m4g5x0j8S05PTpZ4kow0hGi8egi1PcihP8F+3wFst7JJ/B72ISF3/lKCwDwFi cPFPsL/6KhpHbYakTNbs/qRf61Dl6uICAYcJCfF2tcJ2pBi32++mrXBcf+y4NBz6f64F 3g40/euNXIaQJXnSSvi46Afw7RDa/XRkP7y8oyw7F+dCx4CbvLYeeicI4vSgvborEgGI Aw7w==

In-reply-to: <CAEwQpaK3GzHUWXPxOeEnJzg724+VR_1VG5zKG34yyG7s6LjY9g@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: <CAC4D6+9bAV2QsYtJz0zamNiaBhDKj6GG5D4g7Wru=3vQPtAj-g@mail.gmail.com> <CAC4D6+-J0tCq4fUVH=r_5LrzFQNKiQ3Ex7jFAqgLs=gR59o9ag@mail.gmail.com> <CAC4D6+8ak3ohdKPguxzceGG1MupvYmMBo1McXD7gsHbKPgOSUQ@mail.gmail.com> <CAEwQpaK3GzHUWXPxOeEnJzg724+VR_1VG5zKG34yyG7s6LjY9g@mail.gmail.com>

Sender: owner-ferret_users@xxxxxxxx

Dear Ryan

thank you for the reply, but I think I will explain the problem little more specifically

here i am having 36 points in my x-axis (attached as hval)

and I have tried to find out the frequency per year (1900-2100 which means 201yrs) by running the below command. But I am getting values for the years which don't have any occurrence just the same as the previous year frequency (attachment frequency.txt). I want those years with zero frequency as the number zero only not the value from previous year.

code:

DEFINE AXIS/X occurences = hval[x=1:33]

LET pdf = I[gx=occurences]

DEFINE AXIS/X=1900:2100:1 my_intervals

let frequency=pdf[gx=my_intervals@sum]

thanks

gopika

On Mon, Dec 2, 2019 at 11:35 AM Ryo Furue <furue@xxxxxxxxxx> wrote:

Dear Gopika,
On Sun, Dec 1, 2019 at 11:32 AM gopika suresh <gopikasurvay@xxxxxxxxx> wrote:
If some value is there in any year then it will be increasing otherwise straight line.
Basically, the @RSUM does that, if I understand your requirements correctly:
yes? let events = {1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1}
yes? plot/symbols/lines events, events[i=@RSUM]

The problem is that Ferret defines RSUM on the "wrong" gridpoints. You want to shift it to the left by half gridspacing. To do so would require a bit of complicated manipulation . . . and you also want the initial value of the cumulative sum to be zero. . . . (Personally, I would write the values to a text file, read it by a program and compute the cumulative sum, and plot the result with gnuplot.)

Cheers,
Ryo

Gopz

VARIABLE : PDF[GX=MY_INTERVALS@SUM] SUBSET : 201 points (X) 1900 / 1: .... 1901 / 2: .... 1902 / 3: .... 1903 / 4: .... 1904 / 5: .... 1905 / 6: .... 1906 / 7: .... 1907 / 8: .... 1908 / 9: .... 1909 / 10: .... 1910 / 11: .... 1911 / 12: .... 1912 / 13: .... 1913 / 14: .... 1914 / 15: .... 1915 / 16: .... 1916 / 17: .... 1917 / 18: .... 1918 / 19: .... 1919 / 20: .... 1920 / 21: .... 1921 / 22: .... 1922 / 23: .... 1923 / 24: .... 1924 / 25: .... 1925 / 26: .... 1926 / 27: .... 1927 / 28: .... 1928 / 29: .... 1929 / 30: .... 1930 / 31: .... 1931 / 32: .... 1932 / 33: .... 1933 / 34: 0.20 1934 / 35: 0.20 1935 / 36: 0.20 1936 / 37: 0.20 1937 / 38: 0.20 1938 / 39: 0.16 1939 / 40: 0.16 1940 / 41: 0.16 1941 / 42: 0.16 1942 / 43: 0.16 1943 / 44: 0.16 1944 / 45: 0.16 1945 / 46: 0.16 1946 / 47: 0.16 1947 / 48: 0.16 1948 / 49: 0.16 1949 / 50: 0.16 1950 / 51: 0.15 1951 / 52: 0.13 1952 / 53: 0.13 1953 / 54: 0.13 1954 / 55: 0.13 1955 / 56: 0.13 1956 / 57: 0.13 1957 / 58: 0.13 1958 / 59: 0.13 1959 / 60: 0.13 1960 / 61: 0.13 1961 / 62: 0.13 1962 / 63: 0.13 1963 / 64: 0.13 1964 / 65: 0.13 1965 / 66: 0.13 1966 / 67: 0.13 1967 / 68: 0.13 1968 / 69: 0.13 1969 / 70: 0.13 1970 / 71: 0.13 1971 / 72: 0.13 1972 / 73: 0.13 1973 / 74: 0.27 1974 / 75: 0.27 1975 / 76: 0.27 1976 / 77: 0.27 1977 / 78: 0.27 1978 / 79: 0.27 1979 / 80: 0.27 1980 / 81: 0.27 1981 / 82: 0.27 1982 / 83: 0.27 1983 / 84: 0.27 1984 / 85: 0.27 1985 / 86: 0.27 1986 / 87: 0.27 1987 / 88: 0.27 1988 / 89: 2.00 1989 / 90: 2.00 1990 / 91: 2.20 1991 / 92: 2.40 1992 / 93: 2.40 1993 / 94: 0.70 1994 / 95: 0.70 1995 / 96: 0.70 1996 / 97: 0.70 1997 / 98: 0.70 1998 / 99: 0.70 1999 / 100: 0.70 2000 / 101: 0.70 2001 / 102: 0.70 2002 / 103: 0.70 2003 / 104: 1.07 2004 / 105: 1.07 2005 / 106: 1.07 2006 / 107: 1.07 2007 / 108: 1.07 2008 / 109: 1.07 2009 / 110: 1.07 2010 / 111: 2.33 2011 / 112: 3.60 2012 / 113: 3.60 2013 / 114: 4.00 2014 / 115: 4.00 2015 / 116: 4.20 2016 / 117: 4.40 2017 / 118: 4.40 2018 / 119: 4.80 2019 / 120: 4.80 2020 / 121: 16.80 2021 / 122: 2.80 2022 / 123: 2.80 2023 / 124: 2.80 2024 / 125: 2.80 2025 / 126: 2.90 2026 / 127: 3.00 2027 / 128: 3.00 2028 / 129: 3.00 2029 / 130: 3.00 2030 / 131: 38.10 2031 / 132: 7.20 2032 / 133: 7.20 2033 / 134: 7.60 2034 / 135: 7.60 2035 / 136: 7.80 2036 / 137: 8.00 2037 / 138: 8.00 2038 / 139: 4.20 2039 / 140: 4.20 2040 / 141: 4.20 2041 / 142: 4.20 2042 / 143: 4.20 2043 / 144: 4.40 2044 / 145: 4.40 2045 / 146: 4.40 2046 / 147: 4.40 2047 / 148: 4.40 2048 / 149: 9.20 2049 / 150: 9.20 2050 / 151: 9.40 2051 / 152: 9.60 2052 / 153: 9.60 2053 / 154: 10.00 2054 / 155: 10.00 2055 / 156: 60.80 2056 / 157: 5.60 2057 / 158: 5.60 2058 / 159: 5.60 2059 / 160: 5.60 2060 / 161: 5.70 2061 / 162: 5.80 2062 / 163: 5.80 2063 / 164: 5.80 2064 / 165: 5.80 2065 / 166: 33.79 2066 / 167: 1.77 2067 / 168: 1.77 2068 / 169: 1.77 2069 / 170: 1.77 2070 / 171: 1.77 2071 / 172: 1.77 2072 / 173: 1.77 2073 / 174: 1.77 2074 / 175: 1.77 2075 / 176: 1.77 2076 / 177: 1.77 2077 / 178: 1.77 2078 / 179: 1.77 2079 / 180: 1.77 2080 / 181: 1.77 2081 / 182: 1.77 2082 / 183: 1.77 2083 / 184: 1.83 2084 / 185: 1.83 2085 / 186: 1.83 2086 / 187: 1.83 2087 / 188: 1.83 2088 / 189: 1.83 2089 / 190: 1.83 2090 / 191: 1.83 2091 / 192: 1.83 2092 / 193: 1.83 2093 / 194: 1.83 2094 / 195: 1.83 2095 / 196: 1.83 2096 / 197: 1.83 2097 / 198: 1.83 2098 / 199: 1.83 2099 / 200: 1.83 2100 / 201: 33.91

Attachment: hval.png
Description: PNG image