[R] GSVD and qr questions

[Richard M. Heiberger]

off line
I looked at the matlab page and the wikipedia page.
Wikipedia says there are two different things with the gsvd name.

So I guess we need to see a tiny example of your starting matrix
and the decomposition that the matlab function you call returns to you.
I guess you will have to send the "reproducible example" to the list.

Rich

On Tue, May 5, 2015 at 12:14 AM, Roy Mendelssohn - NOAA Federal
<roy.mendelssohn at noaa.gov> wrote:
> Thanks.  But to be honest I am outside of what I really know. so I don’t know if what you describe is what I want  I know what the Matlab function I have calls and what it does, and I am trying to mimic its functionality.  So I am trying to see if there is equivalent to the Matlab qr(A,0) call, that way I know I am getting the same result.
>
> -Roy M.
>
>
>> On May 4, 2015, at 9:10 PM, Richard M. Heiberger <rmh at temple.edu> wrote:
>>
>> Roy,
>>
>> Look at the qr.R and qr.Q functions.  With the complete=TRUE argument
>> they bring their matrices up to full size,
>> by constructing an arbitrary orthogonal completion.
>> Does this provide what you are looking for?
>>
>> Rich
>>
>> On Mon, May 4, 2015 at 11:56 PM, Roy Mendelssohn - NOAA Federal
>> <roy.mendelssohn at noaa.gov> wrote:
>>> Hi All:
>>>
>>> I posted some of this earlier to the r-sig-mac because at that time I got a Mac specific error, but now I have a more general question.  For a canonical variate analysis (CVA) used in some system identification routines, I need to calculate a generalized svd (GSVD).  On the web there are several scripts that are suppose to work with the LAPACK implementation of GSVD, but when I try them with my matrices I get an index error.  I can save the two arrays using R.matlab and read the matrices into Matlab, and the GSVD function in Matlab works fine for the same arrays.  I get the same error trying to use the GSVD routine in the PEIP package in R.
>>>
>>> I have found a Matlab routine that calculates the GSVD using only QR, CHOL and SVD calls, but some of this is outside of what I know and I am trying to translate the routine to R.  I have two questions:
>>>
>>> 1.  Does anyone know of a working routine for GSVD in R  (and I looked in ‘sos’ that is how I found the one in PEIP)?
>>>
>>> 2. In the Matlab function I am trying to translate, they call for the QR decomposition  of the matrix A as qr(A,0).  The second argument (0) according to Matlab does the economy QR decomposition:
>>>
>>>> qr(A,0) produces the economy-size decomposition. If m > n, only the first n columns of Q and the first n rows of R are computed. If m<=n, this is the same as [Q,R] = qr(A).
>>>
>>> Reading over the help for the”qr”  function in R it is unclear that there is a similar way to do an “economy” decomposition with one of the optional arguments.  I am probably being brain-dead, but I don’t see how to do that in the R function.
>>>
>>> Thanks for any help.
>>>
>>> -Roy M.
>>>
>>> **********************
>>> "The contents of this message do not reflect any position of the U.S. Government or NOAA."
>>> **********************
>>> Roy Mendelssohn
>>> Supervisory Operations Research Analyst
>>> NOAA/NMFS
>>> Environmental Research Division
>>> Southwest Fisheries Science Center
>>> ***Note new address and phone***
>>> 110 Shaffer Road
>>> Santa Cruz, CA 95060
>>> Phone: (831)-420-3666
>>> Fax: (831) 420-3980
>>> e-mail: Roy.Mendelssohn at noaa.gov www: http://www.pfeg.noaa.gov/
>>>
>>> "Old age and treachery will overcome youth and skill."
>>> "From those who have been given much, much will be expected"
>>> "the arc of the moral universe is long, but it bends toward justice" -MLK Jr.
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> 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.
>
> **********************
> "The contents of this message do not reflect any position of the U.S. Government or NOAA."
> **********************
> Roy Mendelssohn
> Supervisory Operations Research Analyst
> NOAA/NMFS
> Environmental Research Division
> Southwest Fisheries Science Center
> ***Note new address and phone***
> 110 Shaffer Road
> Santa Cruz, CA 95060
> Phone: (831)-420-3666
> Fax: (831) 420-3980
> e-mail: Roy.Mendelssohn at noaa.gov www: http://www.pfeg.noaa.gov/
>
> "Old age and treachery will overcome youth and skill."
> "From those who have been given much, much will be expected"
> "the arc of the moral universe is long, but it bends toward justice" -MLK Jr.
>