similar to: EIGEN VECTOR PROBLEM

Please excuse what is obviously a trivial matter... I have a large 3-d array. I wish to set the third dimension (z) to NA if there are any NA values in the first two dimnesions (xy). So, given array foo: foo <- array(data = NA, dim = c(5,5,3)) foo[,,1] <- matrix(rnorm(25), 5, 5) foo[,,2] <- matrix(rnorm(25), 5, 5) foo[,,3] <- matrix(rnorm(25), 5, 5) I'll set two elements

hello everyone. Look at the following R idiom: a <- array(1:30,c(3,5,2)) M <- (matrix(1:15,c(3,5)) %% 4) < 2 a[M,] <- 0 Now, I think that "a[M,]" has an unambiguous meaning (to a human). However, the last line doesn't work as desired, but I expected it to...and it recently took me an indecent amount of time to debug an analogous case. Just to be explicit, I would

On Wed, Mar 18, 2020 at 01:46:14PM -0400, habib dan aouta wrote: > Hello Richard, > > Hope you are doing well. My name is Habib and I am current student > at the University of North Carolina at Charlotte(U.S). I am > currently following your Libnbd client fuzzing tutorial from the > Wordpress articles

[I thought I'd submitted this bug report some time ago, but it's never showed up on the bug tracking system, so I'm submitting again.] qr.solve() gives incorrect results when dealing with complex matrices or with qr objects that have been computed with LAPACK=TRUE, whenever the b argument has more than one column. This bug flows from qr.coef(), which has a similar problem. I believe

Hi everyone, It's a very minor point, but could we ensure that eigen and La.eigen return a *matrix* for the "vectors" component of the list by including a "drop = FALSE", as specified in the help file, ie put list(values = z$values[ord], vectors = if (!only.values) z$vectors[, ord, drop = FALSE]) Thanks, Jonathan. --please do not edit the information

Dear all, It is not a R related problem rather than statistical/mathematical. However I am posting this query hoping that anyone can help me on this matter. My problem is to get the Geometrical Interpretation of Eigen value and Eigen vector of any square matrix. Can anyone give me a light on it? Thanks and regards, Arun [[alternative HTML version deleted]]

Dear R-users, I'm estimating a two- dimensional state-space model using the FKF package. The resulting log likelihood function is maximized using auglag from the Alabama package. The procedure works well for a subset of my data, but if I try to use the entire data set I get the following error message. Error in eigen(a$hessian) : infinite or missing values in 'x' What's even

Hi all, when testing the new improvements in the new 1.7.0-version I stumbled over the following: >eigen(matrix(c(0,.3,2,.9),2,2)) Error in eigen(matrix(c(0,.3,2,.9),2,2)) : LAPACK routine DGEEV gave error code -13 >eigen(matrix(c(0,.3,2,.9),2,2),EISPACK=TRUE) $values [1] 1.3458236 -0.4458236 $vectors [,1] [,2] [1,] -1.1436890 -0.9760443 [2,] -0.7696018

Hi. I was comparing the components from princomp's loadings and the eigen given the same input. I found that the sign of componenets (+/-) are opposite between the two components (from princmop and eigen) but the magnitudes are identical. Why? Thanks! [[alternative HTML version deleted]]

Dear all R users, Can anyone tell me to calculate Eigen value of any real symmetric matrix which algorithm R uses? Is it Jacobi method ? If not is it possible to get explicit algorithm for calculating it? Thanks and regards, Arun [[alternative HTML version deleted]]

Hi, Can anyone please point to how to decompose BCCB (Block-Circulant-Circulant-Block) matrices? I am interested in the derivations: I do know that this can be numerically done using 2-dimensional FFTs. Many thanks and best wishes!

Hi all. I got a question about eigenvector. I've tried input a symmetric matrix to both R (using eigen function) and minitab, but the result is really different. Can anyone tell me what's wrong with that? Thanks. [[alternative HTML version deleted]]

Dear all, I am running an Ecological Niche Factor Analysis in R but I am stuck with a problem. As soon as I ask the software to compute the ENFA, I get the following error message: Error in eigen(Se) : infinite or missing values in 'x' Does anyone know what could be wrong. I do not have missing values in my dataset. Please follow the codes I used for my analysis so far: library

Dear, My problems are; 1. Simulate a P-variate multivariate data set (N by P) 2. Draw samples of size n from N (with or without replacement) 3. Obtain eigen values of the variance-covariance matrix of each sample drawn in (2) above. 4. print out the matrix of eigen values as well as the identity of sample generating them. Thanks ---------------------------------

Hi everyone, I am using eigen() to extract the 2 major eigenpairs from a large real square symmetric matrix. The procedure is already rather efficient, but becomes somehow slow for real time needs with moderately large matrices (few thousand lines). The R implementation statically extracts all eigenvalues (and optionally associated eigenvectors). I heard about optimizations of the eigen

Sorry if this bug is fixed in the meantime, but what is this?? > m [,1] [,2] [,3] [1,] 1 2 1 [2,] 0 2 0 [3,] 1 2 1 > eigen(m) $values [1] 2 2 0 $vectors [,1] [,2] [,3] [1,] -4.194304e+06 0.7071068 -0.7071068 [2,] -4.656613e-10 0.0000000 0.0000000 [3,] -4.194304e+06 0.7071068 0.7071068 > eigen(m, symmetric=T) $values [1]

eigen() seems to work for symmetric matrices only. This is out of sync with the help file. > trpr.37 0 1 2 3 4 0 1.00000000 0.0000000 0.0000000 0.0000000 0.0000000 1 0.44444444 0.5555556 0.0000000 0.0000000 0.0000000 2 0.02439024 0.2439024 0.7317073 0.0000000 0.0000000 3 0.00000000 0.0000000 0.2307692 0.7692308 0.0000000 4 0.00000000 0.0000000

I don't have time to track this down right now -- if no-one does it in the next few days before I finish the current piece of work I will try to find it myself. R gets the wrong answer for the eigensystem of a (particular) 4x4 matrix (SunOS 4.1.4). Running these commands: x _ matrix( c(-10, 0 ,0 ,0, 0 ,-0.3 ,0 ,0, 1 , 0.03 ,-0.594,0.0020,

eigen(matrix(rep(1,9),ncol=3),only.values=T) works repeatedly but eventually (after 3 or 4 goes) gave a core dump of 11mb (I think only the first answer was right). if I type gc() immediately after the first call, I get a core dump right away (this might possibly be related to my huge core dumps with nlm) eigen(matrix(rep(1,9),ncol=3)) core dumps immediately after the first call Jim

Eigen seems to be giving different results in 0.64 : eigen(matrix( c(0.00000000, 0.0000000, 0.00000000, 1.000000000, 0.01262442, -0.1141633, 0.04989433, 0.006112607, 0.04781087, 0.4177869, -0.15569656, 0.074242237, 0.20167331, 0.2531109, -0.04973573, 0.034046467) ,4,4))$values R: [1] 0.35519256 -0.29294606 -0.26637369 -0.03168619 S: [1] 0.4689288205 -0.4246063863 -0.2810596564