<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="m_-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.<wbr>com/decuslib20-02/01/decus/20-<wbr>0026/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">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><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>