低通滤波之后的回归分析的显著性检验，出现了填色区域都过检的情况！如何修改呢

JGD飘

;----年代际滤波------
ihp  = 0            ; low pass 得到低频信号，
nwt  = 49
fca  = 1./120
fcb  = -999
sigma= 1.0
wgt  =  filwgts_lanczos(nwt, ihp, fca, fcb, sigma )
pc  := wgt_runave_Wrap(pc, wgt, 0)
slp := wgt_runave_n_Wrap(slp,wgt,0,0)
;---------------
yy     = regCoef(pc, slp({lat|:},{lon|:},{time|:}))
copy_VarCoords(slp(0,:,:),yy)
;---------
p      = student_t(yy@tval,yy@nptxy-2)
p_yy   = onedtond(p,(/dims(1),dims(2)/))
copy_VarCoords(slp(0,:,:),p_yy)
printMinMax(p,False)
print(avg(p_yy(:,23)))  
;-------------------------
最后一张是打点的代码，过了0.01都打上点
但是出图的结果有问题，有朋友帮忙看看吗！

学生党少华

滤波之后有效自由度都变了，不能按照以前的来检验了吧

JGD飘

学生党少华 发表于 2021-11-9 12:46
滤波之后有效自由度都变了，不能按照以前的来检验了吧

p      = student_t(yy@tval,yy@nptxy-2-49)
在自由度这里减去损失的值，得到的图也是差不多的

lqouc

JGD飘 发表于 2021-11-9 14:27
p      = student_t(yy@tval,yy@nptxy-2-49)
在自由度这里减去损失的值，得到的图也是差不多的

自由度减多少你是怎么算的？如果是通过计算有效样本量得到的结果，那就只能是这样了，当然也可以提高显著性水平来减少打点范围。

JGD飘

lqouc 发表于 2021-11-9 15:13
自由度减多少你是怎么算的？如果是通过计算有效样本量得到的结果，那就只能是这样了，当然也可以提高显著 ...

如何计算有效自由度呢，我用了家园上的公式，计算出来是负值

lqouc

JGD飘 发表于 2021-11-10 17:36
如何计算有效自由度呢，我用了家园上的公式，计算出来是负值

Metz W 1991 Optimal relationship of large-scale flow patterns and the barotropic feedback due to high-frequency eddies J. Atmos. Sci. 48 1141-1159这篇论文的附录A有介绍，按照公式编程计算一下就行

学气象的小弱鸡

你好大佬，请问下您那个权重系数49是依据什么来的啊？

zyy542071902

想请教一下楼主，权重系数应该怎么确定啊？

gjlxr

真厉害，可惜我搞气象多年，还是个大老白

lz20200901@

mark一下，后期学习了