Skip to content
Prev 65985 / 398502 Next

Summing up matrices in a list

try this:

matSums <- function(lis, na.rm=FALSE){
    if(!is.list(lis) || !all(sapply(lis, is.matrix))) stop("'lis' must 
be a list containing 2-dimensional arrays")
    dims <- sapply(lis, dim)
    n <- dims[1, 1]
    p <- dims[2, 1]
    if(!all(n == dims[1, ]) || !all(p == dims[2, ])) stop("the 
matrices must have the same dimensions")
    mat <- matrix(unlist(lis), n * p, length(lis))
    matrix(rowSums(mat, na.rm=na.rm), n, p)
}
#########
mylist <- list(matrix(1:6, 2), matrix(7:12, 2))
matSums(mylist)


I hope it helps.

Best,
Dimitris

----
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven

Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
Web: http://www.med.kuleuven.ac.be/biostat/
     http://www.student.kuleuven.ac.be/~m0390867/dimitris.htm


----- Original Message ----- 
From: "Vicky Landsman" <msvika at mscc.huji.ac.il>
To: <r-help at stat.math.ethz.ch>
Sent: Wednesday, March 16, 2005 5:21 PM
Subject: [R] Summing up matrices in a list