similar to: Matrix library, CHOLMOD error: problem too large

Hi: I was trying to use apply on a sparse matrix from package Matrix, and I get the error: Error in asMethod(object) : Cholmod error 'problem too large' at file ../Core/cholmod_dense.c, line 106 Is there a way to apply a function to all the rows without bumping into this problem? Here is a simplified example: > dim(sm) [1] 72913 43052 > class(sm) [1] "dgCMatrix"

Hi all! I would like to import a matlab or gauss code to R. Could you help me? Bye, Sebasti?n. 2007/6/23, r-help-request en stat.math.ethz.ch <r-help-request en stat.math.ethz.ch>: > Send R-help mailing list submissions to > r-help en stat.math.ethz.ch > > To subscribe or unsubscribe via the World Wide Web, visit >

Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod_rowfac.c, line 432 My data are counts of microbe colony forming units (CFUs) collected from

Hi, I am getting a warning message when I am fitting a generalized linear mixed model (m1.2 below). CHOLMOD warning: matrix not positive definite Error in objective(.par, ...) : Cholmod error `matrix not positive definite' at file:../Supernodal/t_cholmod_super_numeric.c, line 614 Any idea? Thanks for your help, Reza > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32

Dear R-users, I'm new to R, so my apologies if this question doesn't make sense. I've tried the following model in lmer, and it works perfectly: model<-lmer(aphids~densroot+zone+(1|zone/site), family=quasipoisson) But if I try the exact same model with a different variable, totmas, the model looks as follows: model<-lmer(aphids~totmas+zone+(1|zone/site), family=quasipoisson)

Hi, I need to multiply all columns in a matrix so something like apply(x,2,sum), but using multiplication should do. I have tried apply(x,2,"*") I know this must be trivial, but I get: Error in FUN(newX[, i], ...) : invalid unary operator The help for apply states that unary operators must be quoted. I tried single quotes too, with the same results. Thanks, -Jose -- Jose

Hi All, I recently found gmane http://gmane.org/ It's a system to covert mail to news and back, with the nice property of keeping a searchable archive... Very convenient if you are subscribed to many lists and don't want to have your mail box cluttered. I use it to read several mailing lists already, but R is not available there. I wonder if the admins know about gmane and if they

Hi, Is there any simple solution to get ranks in descending order? Example, a <- c(10, 98, 98, 98, 99, 100) r <- rank(a, ties.method="average") produces 1 3 3 3 5 6 I would want this instead: 6 5 3 3 3 1 Note that reversing r doesn't work but in small examples. Thanks, -Jose -- jquesada at andrew.cmu.edu Research associate http://lsa.colorado.edu/~quesadaj Dept. of

Hi, Has anyone run any R vs Python (numpy) tests? I'd love to see what the differences performance-wise are, specially handling large sparse matrices. Since both rely on external C code, there might not be much of a difference. If you know and use both languages, what are the main differences and what made you stick to one over another? I also noticed that there are strong libraries

Eric ¿Que dataset utilizo? Por curiosidad probé su código, pero me salen errores. Copio y pego todo como me sale en la consola. > # https://raw.githubusercontent.com/datasets/covid-19/master/data/time-series-19-covid-combined.csv > datos <- read.csv("C:/Users/HP/Downloads/time-series-19-covid-combined.csv") > > library(data.table) > library(ggplot2) >

Muchas gracias a todos por su ayuda. Hace un rato antes de abrir el correo, he encontrado dos soluciones a mi problema de como calcular el porcentaje de variación. La primera es usando el paquete dplyr: https://stackoverflow.com/questions/48196552/calculate-percentage-change-in-r-using-dplyr/48196871 La segunda usando la función PercChange del paquete DataCombine

Hi R usrs, I want to build R-2.10.0 on HP-UX, but I got following error message: ld: Unsatisfied symbol "zgemm" in file CHOLMOD.a[cholmod_l_super_numeric.o] ld: Unsatisfied symbol "zgemv" in file CHOLMOD.a[cholmod_l_super_solve.o] ld: Unsatisfied symbol "zherk" in file CHOLMOD.a[cholmod_l_super_numeric.o] ld: Unsatisfied symbol "ztrsm" in file

Hi, I got the following error. Sorry but this time I couldn't reproduce it with a simple chunk of code: .TM.repl.i.2col(): drop 'matrix' case ... Error in .nextMethod(x = x, i = i, j = j) : 'i' has no integer column number should never happen; please report In addition: Warning messages: 1: Ambiguous method selection for "%*%", target

(Extremely sorry, disregard previous email as I hit send before pasting the latest version of the example; this one is smaller too) Dear R users, I want to apply a function that takes two vectors as input to all pairs (combinations (nrow(X), 2))of matrix rows in a matrix. I know that ideally, one should avoid loops in R, but after reading the docs for do.call, apply, etc, I still don't know

(Extremely sorry, disregard previous email as I hit send before pasting the latest version of the example; this one is smaller too) Dear R users, I want to apply a function that takes two vectors as input to all pairs (combinations (nrow(X), 2))of matrix rows in a matrix. I know that ideally, one should avoid loops in R, but after reading the docs for do.call, apply, etc, I still don't know

Gracias Ruben, muy interesante. Algo así estuve buscando, por ciudades y a nivel internacional, pero no lo encontré. Me bastaría con el número de fallecimientos por día, de cuantas más ciudades del mundo mejor. Si alguno supiera de algo así, le agradecería que me lo comunicase. Un saludo, Manuel El mar., 24 mar. 2020 a las 21:04, Rubén Fernández Casal (< rubenfcasal en gmail.com>) escribió:

Hi, I'm finding that readMM() cannot read a file written with writeMM(). Example: library(Matrix) a = Matrix(c(1,0,3,0,0,5), 10, 10) a = as(a, "CsparseMatrix") writeMM(a, "kk.mm") b = readMM("kk.mm") Error in validObject(.Object) : invalid class "dgTMatrix" object: all row indices must be between 0 and nrow-1 Thoughts? Thanks, -Jose -- Jose

Hi, I'm attending summer School at UCLA (IPAM) on "probabilistics models of cognition". I have been an R-user since v. 1.4.1, but was trained in the frequentist tradition (as most psychologists!). I found that all faculty here use matlab and Murphy's bayes net toolbox. I have not had the need to use matlab before, and would love to stick to R for graphics models and

Hello, The following piece of code works fine in R.2.10 (ubuntu): switch(distr,         normal    = {if (is.infinite(param["desv"]))                      n <- c(n,"La desv. estándar no puede ser Inf.")                    if (param["desv"]<0)                      n <- c(n,"La desv. estándar no puede ser <0.")                    },        

Dear R users, I need to normalize a bunch of row vectors. At a certain point I need to divide a matrix by a vector of norms. I find that the behavior of Matrix objects differs from normal matrix objects. Example the following code examples differ only in xnormed changing from normal to Matrix object: x = matrix(1:12,3,4) x = as(x, "CsparseMatrix") xnorms = sqrt(colSums(x^2))