Re: [R] SVD of a variance matrix

From: Giovanni Petris <GPetris_at_uark.edu>
Date: Tue, 15 Apr 2008 19:06:45 -0500 (CDT)

Hi Ravi,

Thank you for your useful reply. Does the result also hold for variance-covariance matrices that have one or more zero eigenvalues? Do you have a reference to suggest?

Thank you again!

Giovanni

> Date: Tue, 15 Apr 2008 18:14:11 -0400
> From: Ravi Varadhan <rvaradhan_at_jhmi.edu>
> Thread-index: AcifQeEz9B1geo3TQyesYlQGMCSuNgAAWF1QAACa9sA=
>
> Let me correct my reply a bit.
>
> U and V will differ by a factor of (-1) corresponding to negative
> eigenvalues (if any) of a general symmetric A. However, for symmetric
> positive-definite matrices (e.g. variance-covariance matrix), they will be
> identical.
>
> Ravi.
>
> ----------------------------------------------------------------------------
> -------
>
> Ravi Varadhan, Ph.D.
>
> Assistant Professor, The Center on Aging and Health
>
> Division of Geriatric Medicine and Gerontology
>
> Johns Hopkins University
>
> Ph: (410) 502-2619
>
> Fax: (410) 614-9625
>
> Email: rvaradhan_at_jhmi.edu
>
> Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html
>
>
>
> ----------------------------------------------------------------------------
> --------
>
>
> -----Original Message-----
> From: r-help-bounces_at_r-project.org [mailto:r-help-bounces_at_r-project.org] On
> Behalf Of Ravi Varadhan
> Sent: Tuesday, April 15, 2008 6:03 PM
> To: 'Giovanni Petris'; r-help_at_r-project.org
> Subject: Re: [R] SVD of a variance matrix
>
> Yes. SVD of any symmetric (which is, of course, also square) matrix will
> always have U = V. Also, SVD is the same as spectral decomposition, and the
> columns of U and V are the eigenvectors, but the singular values will be the
> absolute value of eigenvalues.
>
> Ravi.
>
> ----------------------------------------------------------------------------
> -------
>
> Ravi Varadhan, Ph.D.
>
> Assistant Professor, The Center on Aging and Health
>
> Division of Geriatric Medicine and Gerontology
>
> Johns Hopkins University
>
> Ph: (410) 502-2619
>
> Fax: (410) 614-9625
>
> Email: rvaradhan_at_jhmi.edu
>
> Webpage: http://www.jhsph.edu/agingandhealth/People/Faculty/Varadhan.html
>
>
>
> ----------------------------------------------------------------------------
> --------
>
> -----Original Message-----
> From: r-help-bounces_at_r-project.org [mailto:r-help-bounces_at_r-project.org] On
> Behalf Of Giovanni Petris
> Sent: Tuesday, April 15, 2008 5:43 PM
> To: r-help_at_r-project.org
> Subject: [R] SVD of a variance matrix
>
>
> Hello!
>
> I suppose this is more a matrix theory question than a question on R,
> but I will give it a try...
>
> I am using La.svd to compute the singular value decomposition (SVD) of
> a variance matrix, i.e., a symmetric nonnegative definite square
> matrix. Let S be my variance matrix, and S = U D V' be its SVD. In my
> numerical experiments I always got U = V. Is this necessarily the
> case? Or I might eventually run into a SVD which has U != V?
>
> Thank you in advance for your insights and pointers.
>
> Giovanni
>
> --
>
> Giovanni Petris <GPetris_at_uark.edu>
> Associate Professor
> Department of Mathematical Sciences
> University of Arkansas - Fayetteville, AR 72701
> Ph: (479) 575-6324, 575-8630 (fax)
> http://definetti.uark.edu/~gpetris/
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Wed 16 Apr 2008 - 00:10:10 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 16 Apr 2008 - 18:30:28 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive