EOF模态显著性检验方面的疑惑

麦田_smile

w1215lin 发表于 2013-9-18 22:03

                               
登录/注册后可看大图

楼主，我看了下，因为你这个特征值是按从大到小的顺序排列的，所以对我图片里提到的公式应该稍微修改一下， ...

不好意思 我没细看你的图片 我的公式的确不是你图中给出的不等式 我用的就是你刚刚说的那个不等式，不过我上面的确错了，应该是前18个模态显著。但是我的问题是这么多模态显著，我在研究时究竟应该选择前几个呢？

w1215lin

麦田_smile 发表于 2013-9-19 09:45

                               
登录/注册后可看大图

不好意思 我没细看你的图片 我的公式的确不是你图中给出的不等式 我用的就是你刚刚说的那个不等式，不过我 ...

你再用蒙特卡洛检验试试

麦田_smile

w1215lin 发表于 2013-9-19 13:34

                               
登录/注册后可看大图

你再用蒙特卡洛检验试试

请问你有相关程序吗？

w1215lin

没有，我都是用matlab自己写的，还不知道对不对，看你们都用fortran的。论坛上有个帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=13501，里面有你要的程序。

麦田_smile

w1215lin 发表于 2013-9-19 14:49

                               
登录/注册后可看大图

没有，我都是用matlab自己写的，还不知道对不对，看你们都用fortran的。论坛上有个帖子http://bbs.06climat ...

我用的就是这个 0.05的信度出来还是很多 我想提高下信度 但还不知道怎么改参数 在琢磨呢

lqouc

麦田_smile 发表于 2013-9-19 15:59

                               
登录/注册后可看大图

我用的就是这个 0.05的信度出来还是很多 我想提高下信度 但还不知道怎么改参数 在琢磨呢

看你俩聊天都给我笑坏了，之前不是告诉你方法了，自己尝试下嘛。
不过你的问题根本就不在于这个信度怎么设置，你就是设置成0.0001也可能还是有10个模态通过检验。
这两种检验只是为了证明你想分析的那个模态是独立的，不是随机误差造成的。
那么一般情况下EOF分析很少有选用超过前三个模态的，原因在于后面的方差贡献太小，我们这个方法就是要抓住成分，你弄得太多就失去意义了。而且EOF只是个开始，最关键是你对主模态的物理解释，设想你要分析20个模态，那不得写本书出来啊。
我觉得你的问题就在于对自己研究的对象没有一个科学的假设和判断，在做EOF之前你就应该知道大概会是什么结果，因为这一领域肯定有人进行过相关研究。根据前人的结果，选取你要分析的主成分。
比如你要研究将岁的年际变化，那就先把数据进行一个年平均和去趋势，得到的结果进行EOF，那么前几个模态可能有年代际和年际的，你要根据pc的波动特征判断。一般是选取前四个模态，看哪个模态是有年际变化周期，然后对其进行深入分析。
说这么多不知道你能不能明白，我觉得吧，不要太在意程序怎么编，数据怎么算，多看看文献，找到正确的思路最重要。

rceclx

lqouc 发表于 2013-9-19 16:35

                               
登录/注册后可看大图

看你俩聊天都给我笑坏了，之前不是告诉你方法了，自己尝试下嘛。
不过你的问题根本就不在于这个信度怎么 ...

版主讲得不错啊！

大太阳

有道理呀，学习了

meredithxf

楼主，我看到有文献写那个n是自由度，你确定是那个测量的次数么？

tsl08

lqouc 发表于 2013-9-19 16:35
看你俩聊天都给我笑坏了，之前不是告诉你方法了，自己尝试下嘛。
不过你的问题根本就不在于这个信度怎么 ...

“我觉得你的问题就在于对自己研究的对象没有一个科学的假设和判断，在做EOF之前你就应该知道大概会是什么结果。
我觉得吧，不要太在意程序怎么编，数据怎么算，多看看文献，找到正确的思路最重要。”
这两句话说的太好了~