关于滑动相关函数（run_cor）的返回值的坐标属性

LiuW

run_cor[size=13.3333px](x,y,time,5)，参考官网的函数说明，我将x,y,time变量的Coordinates属性的范围均设置为1438-1473，但是返回值的Coordinates属性范围是在2002-2033，返回值的size倒是没问题，但是时间坐标怎么差这么远？求解

图中cor10是返回值，其他两个x和time变量，y也是一样的

LiuW

经过学习，发现RUN_COR函数中time变量的单位影响返回值的坐标属性，但是还是存在问题，按照官网给的介绍和计算过程的详细代码，计算11年滑动相关时，第i年的返回值应该是两气候要素的i至i+11量的相关系数，也就是最后画XY曲线图，最后10年应该没有，但是我的返回值看上去是两气候要素i-5至i+5年的相关系数，画XY图前后各少五年，到底哪个变量的属性设置引起这样的变化？向有了解的大佬求助。

小悦悦

请问您的time的单位属性设置的什么？我的cor出来的时间序列都是2002-2031，太奇怪了。

LiuW

小悦悦 发表于 2020-5-5 12:17
请问您的time的单位属性设置的什么？我的cor出来的时间序列都是2002-2031，太奇怪了。

time的单位属性必须要设置为ut_calendar（cd_calendar）能够识别的单位，我的是年际的数据，就设置了“years since .........”，好像设置的识别不了就会从2000开始，你的是五年滑动吗？

小悦悦

我做的11年滑动。我看了一下，这个函数做的结果和我自己写的程序做滑动过后然后用escorc做相关的结果一样，但是我还是不知道这样做滑动相关的公式是否和林学椿老师提出来的做法一样，而且我看做滑动相关也有两种算法吧，纠结。

LiuW

小悦悦 发表于 2020-5-8 15:53
我做的11年滑动。我看了一下，这个函数做的结果和我自己写的程序做滑动过后然后用escorc做相关的结果一样， ...

按照NCL官网的函数提供者的代码，这就是以i至i+n-1（n为滑动点数）来计算i时刻的相关系数，我数理基础比较弱，我理解林老师1978年文献中提到的滑动系数计算方法和这个是一样，下面附上面提到的代码，仅供参考。
Here I am contributing a function to calculate running/moving
correlation..............
; Function to calculate running/moving correlation..
undef ("run_cor")
function run_cor(x:numeric,y:numeric,time:numeric,wSize:integer)
;**********************************************************************
; x and y: numeric array of 1D [eg: x(time) and y(time)]
; x and y should have same size
; time: numeric array of 1D that represent time coordinate of x or y
; time coordinate must have recognized by ut_calendar
; wSize is the scalar integer that specifies window size.[e.g 11 yr, 15
year, 21 year, etc]
;Sample usage: x(time) and y(time).
; cor_11=run_cor(x,y,time,11)
; cor_15=run_cor(x,y,time,15)
; cor_21=run_cor(x,y,time,21)
;************************************************************************
local aa,bb,i,cor,k,N
begin
N=dimsizes(x)
cor=new((/(N-(wSize-1))/),float)
date=new((/(N-(wSize-1))/),integer)
utc_date=ut_calendar(time,0)
year = floattointeger(utc_date(:,0))
k=(wSize-1)/2
do j=0,N-(wSize)
date(j)=year(k)+j
end do
do i=0,(N-wSize)
aa=x(i:i+(wSize-1))
bb=y(i:i+(wSize-1))
cor(i)=escorc(aa,bb)
end do
cor!0="time"
cor&time=date
return(cor)
end

小悦悦

LiuW 发表于 2020-5-8 21:45
按照NCL官网的函数提供者的代码，这就是以i至i+n-1（n为滑动点数）来计算i时刻的相关系数，我数理基础比 ...

非常感谢！

qxlong

请问做了滑动相关后，您有做显著性检验吗？是用公式计算有效自由度来做检验的吗？

LiuW

qxlong 发表于 2021-3-16 10:14
请问做了滑动相关后，您有做显著性检验吗？是用公式计算有效自由度来做检验的吗？

我没做，如果要做是需要用公式计算出有效自由度。

dh5002

楼主可以贴上完整代码吗