求两个不同资料中的变量相关出问题，求助

subtropical

我用的是两个不同的资料，风场和海温场，我要求的是一个区域的风场强度与海温场的相关情况，我稍作了处理 ，也就是一个是一维的和二维的，然后参照官网的例子和函数自己写了一个脚本，截取一部分，但是出错了，求助，怎么解决。。

曦尘

简单点说，直接填个值，是总年份减去2.

曦尘

Nr代表的是某个格点的自由度，一般情况下，我们认为时间序列的自由度是其长度减去2，比如你的数据是30年，填28就可以了，这里把每个格点的自由度都认为是一样的，所以就是个标量scalar.

曦尘

用错函数了，线性相关系数检验用rtest，https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtml

曦尘

嗯，很不错哟~

subtropical

已解决，总结一下：1.提取在源文件中直接写lat,lon,2参考官网例子时，复制坐标信息写错了，应该是copy_VarCoors.

subtropical

曦尘 发表于 2016-3-17 22:33
嗯，很不错哟~

{:loveliness:}，谢谢回复，多有打扰

subtropical

曦尘 发表于 2016-3-17 22:33
嗯，很不错哟~

请问求完一个一维的变量（time）和场变量(time,lat,lon)相关后，用ttest做信度检验，出错，一维的和二维的应该怎么做信度检验呢，大神不好意思，又来打扰你了

曦尘

或者你直接根据遵从自由度为n-2时，α=0.05显著水平下的ta值（查表，书上或网上都可以），根据公式R=sqrt(ta**2/(n-2+ta**2)) 求出r的临界值，超过这个临界值得就是通过显著性检验的。

subtropical

曦尘 发表于 2016-3-18 16:50
用错函数了，线性相关系数检验用rtest，https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtm ...

太感谢太感谢太感谢您了，菜鸟一只，能得到大神帮助，真是感激不尽，我一定好好学习下。{:eb511:}{:eb511:}

曦尘

小事，不客气哈…

subtropical

曦尘 发表于 2016-3-18 16:50
用错函数了，线性相关系数检验用rtest，https://www.ncl.ucar.edu/Document/Functions/Built-in/rtest.shtm ...

我用了prob  = rtest(ccr, Nr, 0)
但是因为Nr必须与ccr维数一样，所以我用1： prob = rtest(ccr,dimsizes(y(0,:,:)),0)；
2：prob = rtest(ccr,dimsizes(ccr),0)；都出错
fatal:rtest: If Nr is not a scalar, then it must have the same number of dimensions as r
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 190 in file vorSSTA29elninodjfttest.ncl
而且我在算相关的时候已经把信息赋给ccr了，如
ccr = escorc_n(vorarea(time|:),y(time|:,lat|:,lon|:),0,0)
copy_VarAtts(y(0,:,:),ccr)
  copy_VarCoords(y(0,:,:),ccr)
输出变量ccr也是对的。为什么呢？大神，求助
Variable: ccr
Type: float
Total Size: 58320 bytes
            14580 values
Number of Dimensions: 2
Dimensions and sizes:   [lat | 81] x [lon | 180]
Coordinates:
            lat: [-80..80]
            lon: [ 0..358]
Number Of Attributes: 13
  missing_value :       -9.96921e+36
  long_name :   DJF: Monthly Means of Sea Surface Temperature
  units :       degC
  var_desc :    Sea Surface Temperature
  level_desc :  Surface
  statistic :   Mean
  actual_range :        ( -1.8, 33.95 )
  valid_range : ( -5, 40 )
  dataset :     NOAA Extended Reconstructed SST V4
  parent_stat : Individual Values
  NMO : 0
  time :        21914
  _FillValue :  -9.96921e+36