<div dir="ltr">Thanks for clarifying! But it is surprising that I couldn't find anything so far on the  Barnett and Preisendorfer CCA. I hope it is not because of some technical intricacies that I'm currently not aware of <span style="font-family:HelveticaNeue,"Helvetica Neue",Helvetica,Arial,sans-serif">⊙﹏⊙ Nevertheless, extra gain if we can get cancor to work after all! Again, thank you for looking into this issue.</span><div><font face="HelveticaNeue, Helvetica Neue, Helvetica, Arial, sans-serif"><br></font></div><div><font face="HelveticaNeue, Helvetica Neue, Helvetica, Arial, sans-serif">Cheers,</font></div><div><font face="HelveticaNeue, Helvetica Neue, Helvetica, Arial, sans-serif">Xiaoning<br></font><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 13, 2017 at 10:28 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div><div><div><span style="color:rgb(255,0,0)"><b>I have no comment on the Barnett & Priesendorfer canonical correlation.</b><br><b>Somebody else will have to comment.</b></span><br><br>===<br>re: "... the built-in “cancor” (which fails often) ..."<br><br>This is true. There is an issue with NCL's 'cancor function. <br>In fact, there is a 'trouble ticket' on this from **2010**: <br>      JIRA 589: Possible bug in CANCOR<br><br></div><div>My recollection is that It works correctly for "smaller" array sizes but fails when some critical size is reached.<br></div><div>Obviously,  investigating this issue has been put on the back-burner. <br></div><div>-----<br><br>FYI: NCL's canonical correlation function uses a 'double precision' version of the f77 'SUBROUTINE CANOR' from the IBM Scientific Software Package <span class="gmail-m_-4128378100544233116m_-5721434046672628824error">[SPP]</span>.scientific library. <br><br>  <a href="http://pdp-10.trailing-edge.com/decuslib20-02/01/decus/20-0026/canor.ssp.html" target="_blank">http://pdp-10.trailing-edge.co<wbr>m/decuslib20-02/01/decus/20-00<wbr>26/canor.ssp.html</a><br><br>----<br></div>Actually, NCL has two hidden 'interfaces' between the NCL code and the f77 subroutine. :<br><br></div>  (a) a f77subroutine (written by me) <br></div>  (b) a C-language interface between the NCL code and the fortran interface code (a)<br>----<br></div><div>**Maybe**, I will take a look<br><br></div><div>Cheers<br></div><div><br><br><br><br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="gmail-h5">On Sun, Nov 12, 2017 at 12:53 PM, Xiaoning Wu <span dir="ltr"><<a href="mailto:xiaoning.wu.1@stonybrook.edu" target="_blank">xiaoning.wu.1@stonybrook.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-h5"><div dir="ltr">Hello,<div><br></div><div>On a statistical application of NCL -  has anyone tried coding up a script to perform the Barnett and Preisendorfer approach of Canonical Correlation Analysis (BPCCA) between two datasets?</div><div><br></div><div>Compared to the built-in “cancor” (which fails often), the trick with BPCCA is to  decompose the data with PCA/EOF first, and the subsequent linear algebra would involve the dominant modes only, thus bypassing direct inversion of the original (large) matrices which is presumably where the classical approach usually gets stuck. The method looks pretty manageable, with existing matlab scripts to refer to, but I'm trying to avoid reinventing the wheel here assuming it must have been done before in NCL. Any thoughts would be very much appreciated!</div><div><br></div><div>Thank you for your time. Have a great week!</div><div><br></div><div>Best regards,</div><div>Xiaoning  </div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div></div></div>