爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 16878|回复: 5

[作图] 回归分析画垂直剖面图中显著性检验的问题

[复制链接]
发表于 2021-7-31 16:23:36 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
作海冰与高度场的回归问题,需要画120°E上的垂直剖面图,填色图是回归出来的高度场异常,显著性检验的打点部分为什么只在30°N以南呢?(显著性检验的部分程序在画经纬度图的时候都是对的,放到垂直剖面图里就出现问题了),部分程序如下,请各位大佬帮忙答疑解惑~蟹蟹~
计算部分:     
     rc0_1   = regCoef_n(iceindex0_1djf,y_new1,0,0)
     rc0_1 = rc0_1*10
     copy_VarCoords(hgt0_1,rc0_1)
     tval0_1  = onedtond(rc0_1@tval , dimsizes(rc0_1))
     df0_1    = onedtond(rc0_1@nptxy, dimsizes(rc0_1)) - 2
     b     = tval0_1
     b     = 0.5
     prob_rc0_1 = betainc(df0_1/(df0_1+tval0_1^2),df0_1/2.0,b)
     prob_rc0_1 = prob_rc0_1*0.9/0.95
     copy_VarCoords(hgt0_1,prob_rc0_1)
     ;printVarSummary(prob_rc0_1)
     prob_rc0_1!0 ="level"
     prob_rc0_1!1 ="lat"
     prob_rc0_1&level = y_new1&level
     prob_rc0_1&lat   = y_new1&lat

画图部分:
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.1/)                       ;; set to significance level
     res2@cnFillPatterns     = (/17,-1/)
     res2@cnFillColors       = (/1,0/)
     res2@cnFillScaleF       = 1        ; add extra density
     res2@cnFillDotSizeF     = 0.003
     res2@gsnAddCyclic       = True
     res2@gsnLeftString      = ""
     res2@gsnLeftString      = ""
     res2@gsnRightString     = ""


    res@gsnLeftString = "(a) 1979-1997 Dec"
     base(0) = gsn_csm_pres_hgt(wks,rc0_12,res)
     plot(0) = gsn_csm_pres_hgt(wks,prob_rc0_12,res2)
     overlay(base(0),plot(0))


SIC regcoef hgt[1979-1997,120E].png
密码修改失败请联系微信:mofangbao
发表于 2021-9-26 22:06:38 | 显示全部楼层
请问楼主解决了嘛,我对剖面图做检验也遇到这样的问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

李钰岚 该用户已被删除
发表于 2021-10-13 22:15:38 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-3-17 13:52:07 | 显示全部楼层
有人解决这个问题了吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-3-17 15:39:19 | 显示全部楼层
regcoef不是求回归系数的吗,请问后面应该怎么回归到距平上去呀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-3-17 15:41:40 | 显示全部楼层
  请问prob_rc0_1 = prob_rc0_1*0.9/0.95这个是什么意思呀,为什么是0.9和0.95呀
还有这个rc0_1 = rc0_1*10为什么要乘以10呀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表