<div dir="ltr"><div>'eof2data' will construct fields based on the rotated components. The code is just 'do loops' <br><br></div><div>BUT...<br></div><div><br>--<br></div><div>Classic EOFs are orthogonal (independent). The contributions of each component can be added. The associated principle components are not correlated.<br></div><div><br></div><div>However ..<br><br></div><div>If you rotate a subset via varimax, the subset will be orthogonal to one another within the subset BUT **not** to the entire field.<br></div><div>You are 'mixing' the subset eofs. The principle components associated with the patterns **are** correlated<br><br>---<br><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 20, 2015 at 4:05 AM, Ipsita Putatunda <span dir="ltr"><<a href="mailto:ipsita.putatunda@gmail.com" target="_blank">ipsita.putatunda@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi all,<br></div> Am doing varimax roatation of first three EOFs with my data using "eofunc_varimax_Wrap" function, and then trying to reconstruct the plot using "eof2data" function. Though my first 3 REOFs is contributing more than 98% of variance but still am not getting similarity in my reconstructed plot with the original one. I doubt if am making any mistake. Part of my script has been provided herewith. Any help will be appreciated.<br><br><br></div>.........[snip]...<br><div><div> xo=LHo(lev|:,time|:)<br> optEOF =True<br> eofo=eofunc_Wrap(xo,neof,optEOF)<br> eof_tso=eofunc_ts(xo,eofo,False)<br> printVarSummary(eofo)<br> eof_roto=eofunc_varimax_Wrap(eofo,1)<br> eofunc_varimax_reorder(eof_roto)<br> printVarSummary(eof_roto)<br> do n=0,neof-1<br> eof_tso(n,:)=eof_tso(n,:)+eof_tso@ts_mean(n)<br> end do<br> printVarSummary(eof_tso)<br>;=====================================================<br>; Reconstructed Mean <br>;===================================================<br> oo=eof2data(eof_roto(0:2,:),eof_tso(0:2,:))<br> ; oo=eof2data(eof_roto,eof_tso)<br><br></div><div> [snip]<br><br><br></div><div>Thanks,<br></div><div>Ipsita<br></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>