【求助帖】东亚夏季风指数回归降水检验异常

鲁滨逊

有大神知道这个怎么回事吗，我在做1950-1971年用东亚夏季风指数回归和相关东亚夏季降水的时候，发现在中国东部的降水始终不能通过95%的检验。降水资料用的是ERA5的再分析月平均降水资料，降水做了去趋势处理。东亚夏季风指数的算法是两个区域夏季850纬向风场区域平均的差值，指数做了去趋势处理和标准化处理，用的风场资料是ncep的月平均再分析风场资料 。                                                                 

这个指数一般要乘以-1

我对之后的年份做回归和相关，中国东部降水就能通过检验。按道理说东亚夏季风和降水的关系是比较显著的呀，有大神能指导一下可能哪里有问题吗？
                                                                             
我想过会不会是自由度的问题，检验的时候用的自由度就是regcoef函数自带的自由度nptxy-2。求大神指点一下。

ZP盛夏光年

程序发一下看看程序有没有问题

来自安卓APP客户端

鲁滨逊

ZP盛夏光年 发表于 2022-3-22 10:54
程序发一下看看程序有没有问题

我添加了附件，感谢您了，附件里用的东亚夏季风指数和降水资料如果需要的话我再添加附件。

鲁滨逊

ZP盛夏光年 发表于 2022-3-22 10:54
程序发一下看看程序有没有问题

1. ;*************************************************
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
   
6. ;*************************************************
   
7. begin
   
8. ;*************************************************
   
9. ; open file and read in data: data are on a gaussian grid
   
10. ;*************************************************
    
11. yrS = 1950
    
12. yrE = 2021
    
13. nt  = yrE - yrS + 1
    
14. f1=addfile("D:\thesis\data\calculate\RAWdata\Removenino34.nc","r")
    
15. IEASM=f1->SDIEASM_IA                                              ;去除了趋势，年代际信号，nino34指数
    
16. f2=addfile("D:\thesis\data\calculate\RAWdata\precIavnoENSO.nc","r")
    
17. precip=f2->precIavnoENSO
    
18. lat=f2->latitude
    
19. lon=f2->longitude
    
20. ny = dimsizes(lat)
    
21. nx = dimsizes(lon)
    
22. ;计算不同时间段内的回归场
    
23.   reg_precip_IEASM_first   = regCoef(IEASM(0:21) , precip(latitude|:, longitude|:, time|0:21)) ;from 1950 to 1971  
    
24.   reg_precip_IEASM_second   = regCoef(IEASM(22:71) , precip(latitude|:, longitude|:, time|22:71)) ;from 1972 to 2021  
    
25.   reg_precip_IEASM_first!0="lat"
    
26.   reg_precip_IEASM_first!1="lon"
    
27.   reg_precip_IEASM_first&lat=lat
    
28.   reg_precip_IEASM_first&lon=lon
    
29.   copy_VarCoords( reg_precip_IEASM_first , reg_precip_IEASM_second)
    
30.   ;print(reg_precip_IEASM_first(40:50,50))
    
31.   ;计算第一时间段内的检验场
    
32. df   = new((/ny*nx/),"float",precip@_FillValue)
    
33. tval = new((/ny*nx/),"float",precip@_FillValue)
    
34. b    = new((/ny*nx/),"float",precip@_FillValue)
    
35. 
    
36.    df   =  reg_precip_IEASM_first@nptxy-2      ;degrees of frerdom
    
37.    tval =  reg_precip_IEASM_first@tval
    
38.    b    =  tval
    
39.    b    =  0.5
    
40.    prob = betainc(df/(df+tval^2),df/2.,b)   
    
41.    prob = (1.-prob)*100
    
42.    precip_first_test     = new((/ny,nx/),"float",precip@_FillValue)
    
43.    do i=0,nx-1
    
44.      do j=0,ny-1
    
45.        precip_first_test(j,i) = prob(nx*j+i)
    
46.      end do
    
47.    end do
    
48.   ;++++++++
    
49. delete(df)
    
50. delete(tval)
    
51. delete(b)
    
52. delete(prob)
    
53. copy_VarCoords( reg_precip_IEASM_first , precip_first_test )
    
54. ;选取通过95%信度检验的场
    
55.   reg_precip_IEASM_first_95   = mask(reg_precip_IEASM_first,(abs(precip_first_test).ge.95),True)              
    
56. ;++++++++
    
57.   copy_VarCoords( reg_precip_IEASM_first , reg_precip_IEASM_first_95 )
    
58.   print(reg_precip_IEASM_first_95(50:70,70))
    
59. ;计算检验第二时间段内的检验场
    
60. df   = new((/ny*nx/),"float",precip@_FillValue)
    
61. tval = new((/ny*nx/),"float",precip@_FillValue)
    
62. b    = new((/ny*nx/),"float",precip@_FillValue)
    
63. 
    
64.    df   =  reg_precip_IEASM_second@nptxy-2      ;degrees of frerdom
    
65.    tval =  reg_precip_IEASM_second@tval
    
66.    b    =  tval
    
67.    b    =  0.5
    
68.    prob = betainc(df/(df+tval^2),df/2.,b)   
    
69.    prob = (1.-prob)*100
    
70.    precip_second_test     = new((/ny,nx/),"float",precip@_FillValue)
    
71.    do i=0,nx-1
    
72.      do j=0,ny-1
    
73.        precip_second_test(j,i) = prob(nx*j+i)
    
74.      end do
    
75.    end do
    
76.   ;++++++++
    
77. delete(df)
    
78. delete(tval)
    
79. delete(b)
    
80. delete(prob)
    
81. copy_VarCoords( reg_precip_IEASM_second , precip_second_test )
    
82. ;选取通过95%信度检验的场
    
83.   reg_precip_IEASM_second_95   = mask(reg_precip_IEASM_second,(abs(precip_second_test).ge.95),True)              
    
84. ;++++++++
    
85.   copy_VarCoords( reg_precip_IEASM_second , reg_precip_IEASM_second_95 )
    
86.   ;print(reg_precip_IEASM_second_95(:,90))
    
87.   ;存入文件
    
88.   outfile        ="./IEASM_precip.nc"
    
89.   system ("rm -f " + outfile)
    
90.   fout=addfile(outfile,"c")
    
91.   filedimdef(fout, "time", -1, True)
    
92.   fout->reg_precip_IEASM_first=reg_precip_IEASM_first
    
93.   fout->reg_precip_IEASM_first_95=reg_precip_IEASM_first_95
    
94.   fout->reg_precip_IEASM_second=reg_precip_IEASM_second
    
95.   fout->reg_precip_IEASM_second_95=reg_precip_IEASM_second_95
    
96.   end

复制代码

ZP盛夏光年

鲁滨逊 发表于 2022-3-23 19:16

你试试IEASM指数标准化一下，然后再回归

鲁滨逊

ZP盛夏光年 发表于 2022-3-24 09:04
你试试IEASM指数标准化一下，然后再回归

我的IEASM就是标准化后的数据，然后做了去除年代际趋势和nino3.4指数的影响。

tender

你好，请问施能老师定义的东亚夏季风指数在哪里下载？可以发一下吗？球球了

鲁滨逊

具体的定义方法可以找一下，然后根据那个方法计算

乌衣年少

tender 发表于 2023-8-21 17:18
你好，请问施能老师定义的东亚夏季风指数在哪里下载？可以发一下吗？球球了

同求，自己算的和施能老师文章里的数据不太一致

乌衣年少

同求，自己算的和施能老师文章里的数据不太一致