NCL用ttest画显著性区域

mangooooo

前面读数据的部分省略了,,sst是我要用的部分年份的海温,ano是其中的几年的海温,sst_djf_avg和ano_ave是求了冬季平均.
;t-test
siglvl = 0.05
var_sst = dim_variance_n_Wrap(sst,0)
var_ano = dim_variance_n_Wrap(ano,0)
sX = dimsizes(sst)
sY = dimsizes(ano)
iflag = False
prob = ttest(sst_DJF_avg, var_sst,sX, ano_ave, var_ano, sY, iflag, False)
copy_VarMeta(ano_ave, prob)

data = ano_ave - sst_DJF_avg
copy_VarMeta(ano_ave, data)


;plot
wtype = "png"
wtype@wkWidth = 2500
wtype@wkHeight = 2500
wks = gsn_open_wks(wtype, "sst_ano")
gsn_define_colormap(wks, "rainbow")

res = True

res@cnFillOn = True
res@tiMainString = "sst ano"


plot = gsn_csm_contour_map(wks,prob,res)


end


报错是s1应该是一个标量或者与ave1,var1具有相同的维数.

jinwithsky

我遇到跟你相同的问题，还没解决。我贴一下我的程序
dimxy = dimsizes(sp2)
    mtim = dimxy(0)
    mlev = dimxy(1)
    mlat = dimxy(2)

    xtmp = sp2(latitude|:, longitude|:, time|:)
    ytmp = sp1(latitude|:, longitude|:, time|:)

    xave = dim_avg(xtmp)
    yave = dim_avg(ytmp)
    xvar = dim_variance(xtmp)
    yvar = dim_variance(ytmp)

    sigr = 0.05
    xeqv = equiv_sample_size(xtmp, sigr, 0)
    yeqv = equiv_sample_size(ytmp, sigr, 0)

    xn = wgt_areaave(xeqv, 121, 720, 0)
    yn = wgt_areaave(yeqv, 121, 720, 0)

    iflag = False
    pro  =ttest(xave, xvar, xn, yave, yvar, yn, iflag, False)

    sp5 = xave - yave

TB~suo

为什么不用student_t？  S1和S2不是要求样本数量吗？你的dimsizis(sst)又是几维数？

逝风丝痕

可以将sX、sY扩展成和对应变量场维数一致的空间场，只不过数值代表每个格点的样本数

常山赵子龙

逝风丝痕 发表于 2019-5-9 22:28
可以将sX、sY扩展成和对应变量场维数一致的空间场，只不过数值代表每个格点的样本数

谢谢你，原来是这样啊。

Kminl

你好，请问一下sX sY要如何选取呢