<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Feb 10, 2014 at 9:08 PM, alex <span dir="ltr"><<a href="mailto:argriffi@ncsu.edu" target="_blank">argriffi@ncsu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">On Mon, Feb 10, 2014 at 2:03 PM, eat <<a href="mailto:e.antero.tammi@gmail.com">e.antero.tammi@gmail.com</a>> wrote:<br>

> Rhetorical or not, but FWIW I'll prefer to take singular value decomposition<br>
> (u, s, vt= svd(x)) and then based on the singular values s I'll estimate a<br>
> "numerically feasible rank" r. Thus the diagonal of such hat matrix would be<br>
> (u[:, :r]** 2).sum(1).<br>
<br>
</div>It's a small detail but you probably want svd(x, full_matrices=False)<br>
to avoid anything NxN.<br></blockquote><div>Indeed.</div><div><br></div><div>Thanks,</div><div>-eat </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="HOEnZb"><div class="h5">_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</div></div></blockquote></div><br></div></div>