- 积分
- 781
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-7-5
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
RT,NCL在计算相关系数的显著性水平时有一个方便的函数rtest(),直接一条命令就可以解决。但是返回的值都是[0,1]范围内的显著性水平a,如果想区别正负相关区,就需要对结果进行重新处理下。一开始用prob=rtest(ccr,N,0) ;ccr是相关系数
prob=where(ccr.lt.0,-prob,prob)
这样处理后画出的图多出了一条类似0线的东西,仔细看会发现在线是由正负颜色的线叠加起来的。
为了避免这种情况,可以采用下面的小方法
prob=rtest(ccr,N,0) ;ccr是相关系数
prob=1.0-prob
prob=where(ccr.lt.0,-prob,prob) ;这样的话画图的时候就可以以-0.99 , -0.95 , - 0.9 , 0.9 , 0.95 , 0.99来作为levels了,
而不是-0.1 , -0.05 , -0.01 , 0 , 0.01 , 0.05 , 0.1
图形效果如下
多出了一条线
线去掉了
|
评分
-
查看全部评分
|