请教EOF的显著性检验的matlab程序

超台十八

yunboyi 发表于 2012-6-16 22:57

                               
登录/注册后可看大图

谁有蒙特卡洛的检验程序?能否共享一下？
不胜感激！

North 的好像不算检验

你用 NCL 吗, 我按照这文章的方法编了一个小程序:
http://journals.ametsoc.org/doi/pdf/10.1175/1520-0493(1982)110%3C0001%3AASTFPC%3E2.0.CO%3B2



NCL SCRIPT:
--------------------------------------------------------------
;;; see Monthly Weather Review Vol 110 No 1, Jan 1982
;;; "A significance test for principal components applied to a ...."

;;; USER MODIFICATION
p = 64                                          ;; observation points
n = 200                                         ;; number of times
rdim = 100                                      ;; number of experiments
cutoff = 95                                     ;; 95 if rdim = 100
;;; END USER MODIFICATION

data = new((/p,n/),"float")
u = new((/p,rdim/),"float")
u_sorted = new((/p,rdim/),"float")

do r=0,rdim-1
    do j=0,p-1                                 
        data(j,:) = random_normal(0.,1.,n)
    end do
    eof = eofcor(data,p)
    sum_eigenvalue = sum(eof@eval)
    u(:,r) = eof@eval/sum_eigenvalue
end do

;;; reordering
do j=0,p-1
    dum1 = u(j,:)
    dum2 = dim_pqsort(dum1,2)
    u_sorted(j,:) = dum1
end do

print(u_sorted(0:4,cutoff-1))     
--------------------------------------------------------------------------------

我不会 matlab, 但这文章用的方法很简单也很容易编程.

问题是算 EOF 只可以是 correlation matrix.

yunboyi

学习了，多谢超台十八 ！

杨韶洋

搜一下吧里都会有，嘿嘿，强大的论坛呀

可乐加盐

无径之林 发表于 2012-5-28 16:21

                               
登录/注册后可看大图

看看我传你的pdf，根据North准则去判断

我看了你这个PDF，里面举例还真不错。请问 还有其他统计方法的举例分析没有呢 比较感兴趣。如果有能否分享下。嘻嘻

流年忘了彼此

求matlab 程序

ayzqs

North准则检验,我是来学习这个的

哞哞三年半

无径之林 发表于 2012-5-28 16:21

                               
登录/注册后可看大图

看看我传你的pdf，根据North准则去判断

useful，thank you!

kongfeng0824

yunboyi 发表于 2012-6-16 22:57

                               
登录/注册后可看大图

谁有蒙特卡洛的检验程序?能否共享一下？
不胜感激！

可以封你想一下蒙特卡洛程序吗？

kongfeng0824

yunboyi 发表于 2012-6-16 22:57

                               
登录/注册后可看大图

谁有蒙特卡洛的检验程序?能否共享一下？
不胜感激！

不知道你是否有蒙特卡洛检验程序？

lantian

又学习了一些东西，谢谢