Displaying 20 results from an estimated 3000 matches similar to: "Rank of matrix"
2018 Jan 22
3
Inconsistent rank in qr()
Le 22/01/2018 ? 17:40, Keith O'Hara a ?crit?:
> This behavior is noted in the qr documentation, no?
>
> rank - the rank of x as computed by the decomposition(*): always full rank in the LAPACK case.
For a me a "full rank matrix" is a matrix the rank of which is indeed min(nrow(A), ncol(A))
but here the meaning of "always is full rank" is somewhat confusing. Does it
2018 Jan 23
1
Inconsistent rank in qr()
Le 23/01/2018 ? 08:47, Martin Maechler a ?crit?:
>>>>>> Serguei Sokol <sokol at insa-toulouse.fr>
>>>>>> on Mon, 22 Jan 2018 17:57:47 +0100 writes:
> > Le 22/01/2018 ? 17:40, Keith O'Hara a ?crit?:
> >> This behavior is noted in the qr documentation, no?
> >>
> >> rank - the rank of x as
2011 Mar 31
1
rank of Matrix
Dear list,
Can anyone tell me how to obtain the rank of a sparse Matrix, for
example from package Matrix (class dgCMatrix)? Here is an example of
QR decomposition of a sparse matrix (from the sparseQR class help).
library(Matrix)
data(KNex)
mm <- KNex$mm
str(mmQR <- qr(mm))
Similarly, using the functions/classes from the relatively new
MatrixModels package:
library(MatrixModels)
2017 Jun 22
1
Unexpected behaviour of base::qr()$rank
2017-06-22 20:31 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>:
>
>
> On 22.06.2017 20:09, I?aki ?car wrote:
>>
>> 2017-06-22 19:49 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>:
>>>
>>> On 22.06.2017 17:11, Bernd Funovits wrote:
>>>>
>>>>
>>>> Hello,
>>>>
>>>> I
2013 Oct 17
1
pamer.fnc y la nueva versión de R
Hola buenas noches,
tengo un problema bastante gordo. ¿A alguno le ha dejado de funcionar las
funciones pamer.fnc y mcp.fnc con la nueva versión de R? La semana pasada
formatee el ordenador y ahora scripts antiguos no funcionan. La cuestión es
que me precupa que no funcione el ejemplo de tutorial del autor.
Os dejo un script que debería de funcionar y no lo hace
2017 Nov 20
2
package check fail on Windows-release only?
I mistakenly left a write in "/tmp" in the rockchalk package (version
1.8.109) that I uploaded last Friday. Kurt H wrote and asked me to fix
today.
While uploading a new one, I became aware of a problem I had not seen.
The version I uploaded last Friday, 1.8.109, has OK status on all
platforms except r-release-windows-ix86+x86_64. I get OK on
oldrel-windows and also on devel-windows.
2018 Jan 23
0
Inconsistent rank in qr()
>>>>> Serguei Sokol <sokol at insa-toulouse.fr>
>>>>> on Mon, 22 Jan 2018 17:57:47 +0100 writes:
> Le 22/01/2018 ? 17:40, Keith O'Hara a ?crit?:
>> This behavior is noted in the qr documentation, no?
>>
>> rank - the rank of x as computed by the decomposition(*): always full rank in the LAPACK case.
> For a
2018 Jan 22
2
Inconsistent rank in qr()
Hi,
I have noticed different rank values calculated by qr() depending on
LAPACK parameter. When it is FALSE (default) a true rank is estimated and returned.
Unfortunately, when LAPACK is set to TRUE, the min(nrow(A), ncol(A)) is returned
which is only occasionally a true rank.
Would not it be more consistent to replace the rank in the latter case by something
based on the following pseudo code ?
2017 Jun 22
2
Unexpected behaviour of base::qr()$rank
2017-06-22 19:49 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>:
> On 22.06.2017 17:11, Bernd Funovits wrote:
>>
>> Hello,
>>
>> I experienced some unexpected behaviour while determining the rank of matrices (sometimes 1x1 matrices):
>> base::qr(matrix(1e-20))$rank returns 1 (incorrect)
>> base::qr(diag(c(1, 1e-20)))$rank returns 2
2017 Jun 22
2
Unexpected behaviour of base::qr()$rank
Hello,
I experienced some unexpected behaviour while determining the rank of
matrices (sometimes 1x1 matrices):
base::qr(matrix(1e-20))$rank returns 1 (incorrect)
base::qr(diag(c(1, 1e-20)))$rank returns 2 (incorrect)
Best regards,
Bernd
> R.version
_
platform x86_64-w64-mingw32
arch x86_64
os mingw32
system x86_64, mingw32
status
major
2011 Aug 16
2
generalized inverse using matinv (Design)
i am trying to use matinv from the Design package
to compute the generalized inverse of the normal equations
of a 3x3 design via the sweep operator.
That is, for the linear model
y = ? + x1 + x2 + x1*x2
where x1, x2 are 3-level factors and dummy coding is being used
the matrix to be inverted is
X'X =
9 3 3 3 3 3 3 1 1 1 1 1 1 1 1 1
3 3 0 0 1 1 1 1 0 0 1 0 0 1 0 0
3 0 3 0 1 1 1 0 1 0 0 1
2010 Jul 06
1
Error message using mi() in mi package
Hello!
I get the following message when I run the mi() function from the mi
package.
Error while imputing variable: c3 , model: mi.polr
Error in eval(expr, envir, enclos) : could not find function
"c14ordered"
Here's the situation:
I am running R v. 2.9.2 on Mac OSX v. 10.5.8. I am trying to impute
missing data in a data set that I've trimmed down to 302 variables.
1998 Feb 24
3
R-beta: qr(A)$rank
I would like to find the rank of a matrix using R. Looking through a list
of commands it seems that the command qr should do the job. Checking a
simple example gave me the wrong answer. Here is the example:
> A <- c(1,2,3,4,5,6)
> A <- matrix(A, nrow=2, ncol=3)
> A
[,1] [,2] [,3]
[1,] 1 3 5
[2,] 2 4 6
> qr(A)$rank
[1] 3
The correct answer is two.
I
1998 Feb 24
3
R-beta: qr(A)$rank
I would like to find the rank of a matrix using R. Looking through a list
of commands it seems that the command qr should do the job. Checking a
simple example gave me the wrong answer. Here is the example:
> A <- c(1,2,3,4,5,6)
> A <- matrix(A, nrow=2, ncol=3)
> A
[,1] [,2] [,3]
[1,] 1 3 5
[2,] 2 4 6
> qr(A)$rank
[1] 3
The correct answer is two.
I
2003 Nov 22
3
summary.manova and rank deficiency
Hi all,
I have received the following error from summary.manova:
Error in summary.manova(manova.test, test = "Pillai") :
residuals have rank 36 < 64
The data is simulated data for 64 variables. The design is a 2*2 factorial with 10 replicates per treatment. Looking at the code for summary.manova, the error involves a problem with qr(). Does anyone have a suggestion as to how to
2009 Mar 13
1
lsfit w/ rank-deficient x
Dear R-devel,
It seems that lsfit incorrectly reports coefficients when the input matrix 'x' is rank-deficient, see the example below:
## here values of 'b' and 'c' are incorrectly swapped
> x <- cbind(a=rnorm(100), b=0, c=rnorm(100)); y <- rnorm(100); lsfit(x, y)$coef
Intercept a b c
-0.0227787 0.1042860 -0.1729261 0.0000000
Warning
2004 Jul 01
1
QR decomposition and rank of a matrix
In summary.manova the qr decomposition of a NxN
matrix
is calculated and for some cases is giving me
a rank < N.
However, following suggestions of professor Ripley to
calculate the rank of a Matrix
On 7 Jun 2002, Brian Ripley wrote:
> For a more reliable answer, look at the SVD
> (function svd) and look at the
> singular values. For example (from lda.default)
X.s <-
2007 Apr 06
4
Computing the rank of a matrix.
Hi! Maybe this is a silly question, but I need the column rank
(http://en.wikipedia.org/wiki/Rank_matrix) of a matrix and R function
'rank()' only gives me the ordering of the elements of my matrix.
How can I compute the column rank of a matrix? Is there not an R
equivalent to Matlab's 'rank()'?
I've been browsing for a time now and I can't find anything, so any
2005 May 04
4
rank of a matrix
how do I check the rank of a matrix ?
say
A= 1 0 0
0 1 0
then rank(A)=2
what is this function?
thanks
I did try help.search("rank"), but all the returned help information
seem irrelevant to what I want.
I would like to know how people search for help information like this.
rank(base) Sample Ranks
SignRank(stats) Distribution of the
2010 Jul 14
5
Matrix Size
hi -
i just started using R as i am trying to figure out how perform a linear
regression on a huge matrix.
i am sure this topic has passed through the email list before but could
not find anything in the archives.
i have a matrix that is 2,000,000 x 170,000 the values right now are
arbitray.
i try to allocate this on a x86_64 machine with 16G of ram and i get the
following:
> x <-