[Thread Prev][Thread Next][Index]
[ferret_users] lag correlations: time-lag products for all lags
- To: ferret_users <ferret_users@xxxxxxxx>
- Subject: [ferret_users] lag correlations: time-lag products for all lags
- From: William Kessler <william.s.kessler@xxxxxxxx>
- Date: Sat, 22 Jul 2017 10:24:14 -0700
- Arc-authentication-results: i=2; mx.google.com; dkim=pass header.i=@noaa-gov.20150623.gappssmtp.com header.b=StsVB0eu; spf=pass (google.com: domain of william.s.kessler@xxxxxxxx designates 2607:f8b0:400e:c05::22f as permitted sender) smtp.mailfrom=william.s.kessler@xxxxxxxx
- Arc-authentication-results: i=1; mx.google.com; dkim=pass header.i=@noaa-gov.20150623.gappssmtp.com header.b=StsVB0eu; spf=pass (google.com: domain of william.s.kessler@xxxxxxxx designates 2607:f8b0:400e:c05::22f as permitted sender) smtp.mailfrom=william.s.kessler@xxxxxxxx
- 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 :mime-version:to:cc:message-id:date:subject :content-transfer-encoding:from:arc-authentication-results :arc-message-signature:dkim-signature:arc-authentication-results; bh=1p7o1I+OMZ+F+HrBOR7SwaxA9FLKsz3HNlSuHD429Y8=; b=CUGN4vsZg5YFBtVXbXxFzoTmouKgoMKvT+khhM9mt8yG+PaduFeiLNZENd8Yvv0hlv +fKXZwXiOopMVBJQEtFuKQwjScTnwE9b/NCPPqd8ZM5P21sSzIPVZYZLk6OtS7c67riq /Za+dnTEnX1Df+g18xC01zKwThn86O/p914rmp70gxBT17vRmH/k1kvYOtMukK0csbrO vGAPiLYrixINFkzLL2SsJHVf0JwCICAHGiNvItk1i/3ImruIKli8Y0vG1J3U/b/CJ3NZ jtJqfcMn0f8OBzdjB5O0V4loBfwcZNRuvhaGnBPjq8U2qIGHx8Y1mJvffdpJ3aFE/jZS Ucvw==
- Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=mime-version:to:cc:message-id:date:subject :content-transfer-encoding:from:dkim-signature :arc-authentication-results; bh=1p7o1I+OMZ+F+HrBOR7SwaxA9FLKsz3HNlSuHD429Y8=; b=S+7e5+YtlAnH427VVwT78bvVl+DwE3g+Uf9OE52ywCTm7AMCDMHfZOz5fwZZq0EUii 00ZVIwAq2aWwYvuxz3N+xzw5PPQ0B4NagGxTnQVK8Az+UFw/JjVXrsxi6K8yII3h+gis BQLK7Xl4YVcFrr0eE8x+TohslWZ2rz/5SVDZuVxyo/Ev1FoMcDNGLNRDyycmuKbA4q38 55q1YWKpxG4del2qo/2nOi8IPRLLm9sqcRcBBq7kGuN1qh32tKbU0jX3qnViSbqhhWeg mwohZutRGQsRQP/ogfxYsaze9zQStaSIbuAci27UPeItlF3iTwy3/767BOKQHW6udj0n Psjw==
- Arc-seal: i=2; a=rsa-sha256; t=1500744258; cv=pass; d=google.com; s=arc-20160816; b=Vl71eUAfwcRMgBrfrdmmP9aYp74BZDa44aTwPIg81GFS6Sb8IbzDddIYedHEJhRJGo cEW9Xh/7HD3X+ROX+SRBwbZEHiz6/XVG0RH+E5qZ92zXCDtYrMS8qQpzeQ/2vH9ODs+g jLYgOkg/Dxe42OuOY+dT/zWwO0PBrjqsqX/gGtgyIzqt6NceyxFSlWgLmkMUAQ5nIg6h rm7qqqK3flwMdTZazBhMDljvGvgvxBnsuA7JxvBPrxuOmqc/JCVDw1oNAcvjEYw/3e7z bnCpe+cuHUA0PaqKdC/Rw772ASkHOjrryDIyCp9L5PkIvO5DCLvpRaqQrQtsnmLFLSb9 /xFA==
- Arc-seal: i=1; a=rsa-sha256; t=1500744257; cv=none; d=google.com; s=arc-20160816; b=ocRc2whIFQ30w6GUi16EahuqgY067H+wcOLKK+m9pwjPp6w4MVyb9fKIA4rSBWH5TS XrbdruNIfIV2vNevbJvNlxV14b0G0Lbnri1nt92G2J42ZbcscpNPi/2FPCP6C9Ub5j+p Q6TPrIif3AFt6eD1qmFzp3yNdZ2rweRzoHv7Wo2Mbj/hOakrjTfTB9CY1YfaAJxx6xgY m5dbc4DkqOEO8yjd2aB7D9Oy4BEnroT0NcrsJsHXsq2GzrRBdTDZjl7Maq4d2mdT22e6 BXii4YlDzhoU8gYmpIScxKT01E1kl3n3+sCz0d8E5l2ZXgFAq+SJKxYAbilK4AzoZYjr UCOQ==
- Cc: "William S. Kessler" <William.S.Kessler@xxxxxxxx>
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=noaa-gov.20150623.gappssmtp.com; s=20150623; h=from:content-transfer-encoding:subject:date:message-id:cc:to :mime-version:x-original-sender:x-original-authentication-results :precedence:mailing-list:list-id:list-post:list-help:list-archive; bh=1p7o1I+OMZ+F+HrBOR7SwaxA9FLKsz3HNlSuHD429Y8=; b=YpEjxRl2jlEf3HmW/xTlCgoqBbu3L+1xd9Sk+1jUm7N6Ft16bEFAZn/lg8WUTs2IFg RMg81kYAvtWTwRxB8Pa8EYttW7GI4+fLdvH7ry37KFmmIfE9LFXQnr5TLzOfClEUnAay OjuCn5+CCF5hrkVpZT+I7Hzoatlha0RCeyOCjhg1b4F3C4fmDUEbkLNCekiTLmCNhKzz fj7+sC9wZpv/eE9t/5NtyBnXn0zMYgj2gq3K3P+Ona/3N/pG/gbUGxCfLCt7kxy3AdKX SpyObR72/nLzHUO+4Mx+sYNaoi34lxwrJKrfoA5SdDJvzRWpq5f/VPswZVE4cd8rSgEb cpIw==
- 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
Hi Ferreters -
I want to construct a lag-correlation matrix, ideally without a REPEAT loop (there are thousands of points in my time series).
Thus I need to find the time-lagged products, including all lags within a range.
Suppose I have a time series of values DATA[L=1:LL], and the lag range is from L-R to L+R.
If I put the lag range on a M axis, _E=1:2*R+1, each element would be the sum of lagged products:
LAG[M=1] is the sum of: DATA*DATA (over all L)
LAG[M=2] is the sum of: DATA*DATA[L=@shf:1]
LAG[M=3] is the sum of: DATA*DATA[L=@shf:2]
...
LAG[M=2*R+1] is the sum of: DATA*DATA[L=@shf:2*R]
For the moment, let's not worry about the ones that run off the end of the L-axis.
Can I do this without a massive REPEAT loop? For instance, is there a way to promote the time series DATA to 2-d where the 2nd dimension is the time-shifted terms as above?
Billy K
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
Ferret
Privacy Policy | Disclaimer | Accessibility Statement