回归分析的显著性检验问题

李小毛123

各位大神，下面是我做回归检验的部分程序：
rc = regCoef(x, SST(level|:,lat|:,lon|:,shijian|:) )   
rcc     = rc(:,:,0)                             
rcc!0   = "level"    ; name dimensions
rcc!1   = "lat"
rcc&level = SST&level
rcc&lat = SST&lat  
tval = onedtond(rcc@tval , dimsizes(rcc))  
printVarSummary(rcc)
tval@FillValue=  -9.96921e+36
tval!0   = "level"    ; name dimensions
tval!1   = "lat"
tval&level = SST&level
tval&lat = SST&lat   
tval = onedtond(rcc@tval , dimsizes(rcc))   ;t-statistic
df   = onedtond(rcc@nptxy, dimsizes(rcc)) - 2  ;自由度
b = tval    ; b must be same size as tval (and df)
b = 0.5
prob =  betainc(df/(df+tval^2),df/2.0,b)
prob!0   = "level"    ; name dimensions
prob!1   = "lat"
prob&level = SST&level
prob&lat = SST&lat   ; assign coordinate values to named dimensions
res2                    = True
res2@gsnDraw            = False
res2@gsnFrame           = False
res2@cnFillOn           = True
res2@cnLinesOn          = False
res2@cnLineLabelsOn     = False
res2@cnInfoLabelOn      = False
res2@lbLabelBarOn       = False
res2@cnMonoFillPattern  = False
res2@cnLevelSelectionMode = "ExplicitLevels"
res2@cnLevels           = (/0.05/)                       ;; set to significance level
res2@cnFillPatterns     = (/17,-1/)
res2@cnFillColors       = (/1,0/)
res2@cnFillScaleF       = 1        ; add extra density
res2@gsnAddCyclic         = True
res2@gsnLeftString      = ""
res2@gsnLeftString      = ""
res2@gsnRightString     = ""   
plots_11= gsn_csm_contour(wks,prob,res2)
overlay(plot,plots_11)
draw(plot)
frame(wks)
思路比较常规，就是官网上给出的，但是不知道为什么画出来的图通过显著性区域的地方，打点位置明显有问题，只在700hpa 那块出现，是哪一步出问题了呢？？

李小毛123

龙猫猫 发表于 2017-11-2 13:08
不错的东西  持续关注

有木有人回复