similar to: NA in eigen()

Displaying 20 results from an estimated 4000 matches similar to: "NA in eigen()"

2009 Feb 11
2
error in my previous message
i'm sorry. i had an error in my previous code because i left out a letter in the rownames. while fixing that, i also found a solution. so i'm sorry for the confusion. below is my fix. temp2 <- matrix(rnorm(10),nc=1,nrow=10) rownames(temp2) <-
2009 Feb 11
2
sorting a matrix by the column
this is a bad question but I can't figure it out and i've tried. if i sort the 2 column matrix , temp1, by the first column, then things work as expected. But, if I sort the 1 column matrix, temp2, then it gets turned coerced to a vector. I realize that I need to use drop=FALSE but i've put it in a few different places with no success. Thanks. temp1 <-
2006 Jan 04
3
matrix math
I am using R 2.1.1 in an windows XP environment. I have 2 dataframes, temp1 and temp2. Each dataframe has 20 variables (“cocolumns") and 525 observations (“rows”). All variables are numeric. I want to create a new dataframe that also has 20 columns and 525 rows. The values in this dataframe should be the sum of the 2 other dataframe. (i.e. temp1$column
2005 May 01
2
eigen() may fail for some symmetric matrices, affects mvrnorm()
Hi all, Recently our statistics students noticed that their Gibbs samplers were crashing due to some NaNs in some parameters. The NaNs came from mvrnorm (Ripley & Venables' MASS package multivariate normal sampling function) and with some more investigation it turned out that they were generated by function eigen, the eigenvalue computing function. The problem did not seem to happen
2007 Mar 18
1
eigen returns NAs from a real matrix
Hi, All: Attached please find a symmetric, indefinite matrix for which 'eigen(...)$vectors' included NAs: > load("eigenBug.Rdata") > sum(is.na(eigen(eigenBug)$vectors)) [1] 5670 > sessioninfo() Error: could not find function "sessioninfo" > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United
2012 Oct 31
2
Aggregate Table Data into Cell Frequencies
R-help - I have this set of aggregated tables (sample data below via dput()). And I would like to have delayValue as the column variables with the "temp" (temp1, temp2, temp3) values as the row variables. However I would like to have the temp variables *aggregated into single rows* so that I have the frequency ("Freq" | counts) of each time each "delayValue" occurs
2006 Mar 28
2
R crashes during 'eigen'
Hi all, Hi, When I want to compute the eigenvalues & eigenvectors of a specific matrix, R crashes (i.e. it stops responding to any input). I've tried it with different versions of R (2.1.1, 2.2.0, 2.2.1) - all with crashing as result. What I did before the crash was: M <- as.matrix(read.table("thematrix",header=T)) eigen(M) If, instead of eigen(M), I use eigen(M,
2010 Nov 19
3
Sweave Dynamic Graph Question
i have a time Series of IBM closing px from 1/1/2000 to today I want to graph the time serie by dividing the graph by year and month all the monthly graphs with the same year will go to one page. so from 1/1/2000 to 11/19/2010. i will have 11 pages, and each page will have 12 graphs (jan to dec) except for 2010. I am able to do it in R, but when i use sweave, I can only print the last page.
2009 May 28
2
Replace is leaking?
Okay, someone explain this behaviour to me: Browse[1]> replace(rep(0, 4000), temp1[12] , temp2[12])[3925] [1] 0.4462404 Browse[1]> temp1[12] [1] 3926 Browse[1]> temp2[12] [1] 0.4462404 Browse[1]> replace(rep(0, 4000), 3926 , temp2[12])[3925] [1] 0 For some reason, R seems to shift indices along when doing this replacement. Has anyone encountered this bug before? It seems to crop up
2011 Mar 09
2
Anomaly with unique and match
I stumbled onto this working on an update to coxph. The last 6 lines below are the question, the rest create a test data set. tmt585% R R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-unknown-linux-gnu (64-bit) # Lines of code from survival/tests/singtest.R > library(survival) Loading required package: splines
2008 Aug 26
1
parse and eval character vector
Dear R-help, I have a character vector, some elements will be numeric, some not, and some even empty. E.g.: temp1 <- c("abcd"," 2 ","") I'm only interested in the numeric elements, the rest I can just throw away. It is easy enough to loop through the vector: temp <- try(eval(parse(text=temp1[1])), silent=TRUE); class(temp) # try-error temp <-
2003 Apr 18
1
Problem with eigen() and LAPACK
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
2004 Jul 27
4
Problems with Lapack's eigen() on 64-bit
I'm only now realizing that we have severe problems with R on our AMD 'Opteron' and 'Athlon64' clients running Redhat Enterprise with all 64-bit libraries (AFAICS). The Lapack problem happens for R-patched and R-devel both on the Opteron and the Athlon64. Here are platform details: o "gcc -v" and "g77 -v" both end with the line gcc version 3.2.3
2005 Nov 23
2
vector of permutated products
Given an x-vector with, say, 3 elements, I would like to compute the following vector of permutated products (1-x1)*(1-x2)*(1-x3) (1-x1)*(1-x2)*x3 (1-x1)*x2*(1-x3) x1*(1-x2)*(1-x3) (1-x1)*x2*x3 x1*(1-x2)*x3 x1*x2*(1-x3) x1*x2*x3 Now, I already have the correctly sorted matrix of permutations! So, the input looks something like: #input x<-c(0.3,0.1,0.2) Nx<-length(x) Ncomb<-2^Nx
2011 Mar 29
2
List extraction
I have created a list of tables with the same columns but different number of row. Example (actual list has ~200 elements): > temp1<- data.frame(ID=c("Herb","Shrub"),stat=c(4,5),pvalue=c(.03,.04)) > temp2<- data.frame(ID=c("Herb","Shrub", > "Tree"),stat=c(12,15,13),pvalue=c(.2,0.4,.3)) > L<-list(a=temp1,b=temp2) > L $a
2005 Dec 29
1
Problems with calloc function.
Hi all, I have a C code in Linux, it has 7 pointers and compile e run OK, but when I run in R happens problems with calloc function, it returns NULL. ############################################### > int *temp1,*temp2,*temp3,*temp4; temp1 = (int *)calloc(col,sizeof(int)); if(temp1 == NULL){ printf("\n\n No Memory1!"); exit(1); } temp2 = (int *)calloc(col,sizeof(int));
2006 Oct 16
1
Bugs with partial name matching during partial replacement (PR#9299)
This is a rather interesting, but I don't think it is a bug - it is just things that "you are not supposed to do"... you are assuming a certain evaluation order of the 4 "$" operators in " D$ABC[D$M] = D$V[D$M] " as in: temp1 <- D$M # 2nd and 4th temp2 <- D$V[temp1] # 3rd D$ABC[temp1] = temp2 # 1st What R did was this:
2011 May 14
1
Identify Objects that end with .f (and all caps)
Dear R Helpers, I am trying to find a way to identify all the objects in my environment that are all caps and then end with .f. I can do the all caps part pretty easily, but I have tried a number of variations on the \ and can't get a recognition of that operator. As a simple example A.f<-"foo1" AA.f<-"foo2" aa.f<-"foo3" A.a<-"foo4" ls()
2008 Mar 11
2
persp question
someone sent in a question earlier about doing something in 3D so i took a stab at it purely for educational purposes ( i'm not even sure that I understood the question actually ). Unfortunately, persp gives me an error that I don't understand because it says "object y not found". I'm sending y in as a parameter to persp similar to what ?persp shows in one of oits examples
1997 Aug 25
2
R-alpha: eigen and batch
Batch: Putting q(save=F) at the end of my file does not work in my context because I can no longer source the file without quitting. I have that quit statement in my .First so that I always quit that way interactively. The problem is that it is ignored in batch. eigen: The crash occurs on my 586 running Red Hat Linux 2.0.27 but not on my son's 486 running SLackware Linix 2.0.29. We both