NCL ttest 检验有FillValue

timevigor · 发表于 2016-5-16 21:44:47

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

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

x

请问各位做过ncl里t检验的同学，在计算出的prob值有FillValue是怎么回事？官网上说prob的值应该在0-1之间啊？而且FillValue的值还不少

zhangyw · 发表于 2016-5-17 09:37:19

可以用ismissing函数看看有多少缺省值吧

timevigor · 发表于 2016-5-17 10:08:32

zhangyw 发表于 2016-5-17 09:37
可以用ismissing函数看看有多少缺省值吧

嗯都看了可是prob不应该有缺省值才对吧。probt也是有缺省值，但是算的平均值方差都没有缺省值

zhangyw · 发表于 2016-5-17 13:12:36

timevigor 发表于 2016-5-17 10:08
嗯都看了可是prob不应该有缺省值才对吧。probt也是有缺省值，但是算的平均值方差都没有缺省值

如果直接用函数计算的话他是会跳过缺省值的，你prob是怎么计算的呢

timevigor · 发表于 2016-5-17 14:42:40

zhangyw 发表于 2016-5-17 13:12
如果直接用函数计算的话他是会跳过缺省值的，你prob是怎么计算的呢

probt = ttest(aveX,varX,sX, aveY,varY,sY, iflag, True)
prob=probt(0,:,:)
alpha=100*(1.-prob)
t=abs(probt(1,:,:)
就是按照官网来算的，输出probt检查就发现了缺省值，而平均值和方差没问题。是因为方差为0的关系吗

zhangyw · 发表于 2016-5-17 15:31:52

timevigor 发表于 2016-5-17 14:42
probt = ttest(aveX,varX,sX, aveY,varY,sY, iflag, True)
prob=probt(0,:,:)
alpha=100*(1.-prob)

那就用ismissing过滤掉缺省值好啦，你这如果有缺省值肯定是要报错的呀，那就前面加上一个判断好啦，if(.not.ismissing()) then 就可以啦，我认为。。

我想养猫 · 发表于 2016-8-28 13:20:47

看你用什么数据了，如果是积雪海温这种，有缺省值很正常嘛

潇潇 · 发表于 2016-8-29 10:20:20

楼主，问题解决了吗？有缺省值的话，画图怎么去掉这些缺省值呢？

subtropical · 发表于 2017-2-9 18:11:05

zhangyw 发表于 2016-5-17 15:31
那就用ismissing过滤掉缺省值好啦，你这如果有缺省值肯定是要报错的呀，那就前面加上一个判断好啦，if(.n ...

您好，我在对合成进行检验的时候也遇到问题了，好像是分母（标准差）出现0和缺测，

请问怎么设置，可以说的详细点吗？谢谢了

yeah... · 发表于 2017-4-2 09:57:35

subtropical 发表于 2017-2-9 18:11
您好，我在对合成进行检验的时候也遇到问题了，好像是分母（标准差）出现0和缺测，请问怎么设置，可以说 ...

我在进行的时候遇到提示我说 warning:ttest: encountered 4623 cases where s1 and/or s2 were less than 2. Output set to missing values in these cases
请问这是什么问题呢这和s1 s2《2有什么关系呢？
谢谢谢谢

		自动登录	找回密码
密码			立即注册