python 做REOF

liubr11

感觉论坛里REOF的代码不少，但是教程和心得严重不足呀，新手根本不会用。我一直都坚持用python处理数据，希望python也有REOF的函数和代码。我仔细看了文章教程步骤和疑问如下：1、对于降水资料阵X,求资料阵X的相关阵的特征值及其特征向量，构成矩阵V和R，则因子荷载阵A=VR**0.5（问题：相关阵应该是C=XX'吧）
2、主分量矩阵为FR**0.5'V'X
3、确定显著的主分量个数（怎么确定）
4、对因子荷载阵􀀁 进行极大方差转动,最后求出F*（这个最不懂）􀀁另外，下面这个函数是不是极大方差转动呀，具体怎么用，

def varimax(Phi, gamma = 1.0, q = 20, tol = 1e-6):    from scipy import eye, asarray, dot, sum, svd    p,k = Phi.shape    R = eye(k)    d=0    for i in xrange(q):        d_old = d        Lambda = dot(Phi, R)        u,s,vh = svd(dot(Phi.T,asarray(Lambda)**3 - (gamma/p) * dot(Lambda, diag(diag(dot(Lambda.T,Lambda))))))        R = dot(u,vh)        d = sum(s)        if d_old!=0 and d/d_old < 1 + tol: break    return dot(Phi, R)

总之还请大家多多指点，我的同学学长都不做这个，问老师又叫我自己找，还望大家多多帮助：）

liubr11

抱歉粘贴出错了，代码在http://en.wikipedia.org/wiki/Talk%3aVarimax_rotation里面

liubr11

我的理解是varimax函数里的phi就是A，返回的就是A*,欢迎指教！

po_po1

没用过eof分析，所以对那个算法不了解

平流层的萝卜

http://fisica.ehu.es/pyclimate/DownloadPage   楼主可以了解一下这个扩展包，有EOF的程序

fermidq401

这种纯数字的计算用R做很轻松啊，之后可以在python中调用R。

找不回密码了QAQ

没看懂！我也想知道怎么做啊

dragonloon

学习了{:eb502:}{:eb502:}

xidaly

多谢分享。

kust

{:eb502:}