关于pc1回归850百帕高度场后的显著性检验问题

北岛夏屿

①我先是利用回归后返回的tval来进行检验，主要步骤如下：;t

tval = onedtond(rchgt@tval , dimsizes(rchgt))
df   = onedtond(rchgt@nptxy, dimsizes(rchgt)) - 2
b = tval
b = 0.5
prob_rchgt = betainc(df/(df+tval^2),df/2.0,b)


;打点
res2@cnLevelSelectionMode ="ExplicitLevels"
res2@cnLevels = (/0.1/)                    
res2@cnFillPatterns = (/17,-1/)    ;点的样式
res2@cnFillColors = (/1,-1/)       ;点的颜色


得到如第一张图的结果，我感觉不太好，因为
第一张图中，30W左右的赤道区域不知道为什么会有密集的打点区


②所以自己写公式又算了一下，公式为：t =rc*sqrt(df)/sqrt(1-rc^2)，主要步骤如下：
df   = onedtond(rchgt@nptxy, dimsizes(rchgt)) - 2
tval1=new((/97,288/),"float")
do i=0,96
    do j=0,287
        tval1(i,j)=rchgt(i,j)*sqrt(df(i,j))/sqrt(1-rchgt(i,j)*rchgt(i,j))
    end do
end do
b = tval1
b = 0.5
prob_rchgt = betainc(df/(df+tval1^2),df/2.0,b)


;打点
res2@cnLevelSelectionMode ="ExplicitLevels"
res2@cnLevels = (/0.1/)                    
res2@cnFillPatterns = (/17,-1/)    ;点的样式
res2@cnFillColors = (/1,-1/)       ;点的颜色

结果如第二张图所示。
第二张图中，像是画反了一样，所有高值区都没打上点，反而是其他区域打上了点，完全反了过来。想请教一下大家为什么会有这种情况。

李钰岚

请问楼主这个问题解决了吗？

lz20200901@

mark一下，后期要用