<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>