关于导师的话，求相关，我不明白是什么意思

xiaohundanpipi

    导师出差前说用一个时间序列和一个空间场做相关系数的分布图，这个我明白，就是用grads中的tcorr函数实现。    这个函数要求时间序列是时间的一维函数，空间场是一个时次的平面，也就是二维的。

    但是他说要把这个序列和前一年的冬季空间场做个相关系数分布图，再跟当年的春季和夏季分别再做相关系数分布图，看看到底谁对这个时间序列的影响大。我不明白的是，这个时间序列本来就是很多年的数据平均出一年一个点，有这么多年，要跟前一年的冬季作相关，哪一年的前一年？

    导师出差就要回来了，就差一步了，图还没出，好恐慌，请各位帮忙。

xiaohundanpipi

难道空间场要取一个平均场吗？

xiaohundanpipi

难道空间场要取一个平均场吗？

xiaohundanpipi

我明白了，这个空间场是三维的，也是随时间变化的。

言深深

xiaohundanpipi 发表于 2012-9-23 16:07

                               
登录/注册后可看大图

我明白了，这个空间场是三维的，也是随时间变化的。

哈哈，看出来你很捉急啊····
问题主要有两个：
1，同期相关，就是任意两条序列计算相关系数，你的第一个问算的就是这个，应该将你的序列和场上面的每一个个点算一个相关系数，然后画出来每一个个点上的系数值，看分布情况；
2，实际上和前冬算相关是算的滞后相关，两条序列一样长，不过一个是当前的，一个是几个月以前的，相关性好说明变化趋势一致，表面影响较大

xiaohundanpipi

言深深 发表于 2012-9-23 16:56

                               
登录/注册后可看大图

哈哈，看出来你很捉急啊····
问题主要有两个：
1，同期相关，就是任意两条序列计算相关系数，你的第 ...

嗯嗯嗯嗯，我也终于想明白是这个意思了，主要是对grads那个函数不了解，而且基本功不扎实，还得继续努力呀。
哈哈，谢谢你啦~

言深深

xiaohundanpipi 发表于 2012-9-23 17:42

                               
登录/注册后可看大图

嗯嗯嗯嗯，我也终于想明白是这个意思了，主要是对grads那个函数不了解，而且基本功不扎实，还得继续努力呀 ...

呵呵，不客气

jinglj2001

能给出GS语句吗

xiaohundanpipi

jinglj2001 发表于 2012-9-23 20:15

                               
登录/注册后可看大图

能给出GS语句吗

sorry啊，我还写不出来。
我有的数据是一个时间序列和一个平面。
我发现grads同时open两个ctl时要求维数都是一至的才可以进行运算，现在这两个数据意味着我不能同时读进grads里，我不知道怎么办了
这是grads官网上对这个函数的描述，希望对你有帮助。
This example calculates the temporal correlation between sea level pressure and the defined variable elnino, an areal average of surface temperature in the equatorial Pacific.
set x 1
set y 1
set z 1
set t 1 100
define elnino = aave(ts,lon=-160,lon=-80,lat=-10,lat=10)
set lon -180 180
set lat -90 90
set z 1
set t 1
d tcorr(elnino, slp, t=1, t=100)

xiaohundanpipi

jinglj2001 发表于 2012-9-23 20:15

                               
登录/注册后可看大图

能给出GS语句吗

我把时间序列做成平面三维的，带进刚才的函数，就能画出来了。
gs文件是：
'open tt.ctl'
'open hgt.ctl'
'set z 1'
'set t 1 31'
'define cc=aave(tt.1,x=1,x=35,y=1,y=35)'
'set t 1'
'd tcorr(cc,hgt.2,t=1,t=31)'


大致就是这样了，做时间序列与平面的相关系数分布图，还是很简单的，要牢记grads打开多个文件并处理的时候，数据维数必须一致。