[R] how to use a function in aggregate which accepts matrix and outputs matrix?
jim holtman
jholtman at gmail.com
Thu Apr 17 21:16:41 CEST 2008
Does something like this work for you? It is using 'lapply' with the
indices of the rows:
> x <- matrix(c( 'c1' , 'r6', '150', 'c1' , 'r4' ,'70' ,'c1' , 'r2' ,'20',
+ 'c1' , 'r5' ,'90', 'c2' ,'r2' ,'20', 'c3' , 'r1' ,'10'),
byrow=TRUE, ncol=3)
> # use lapply
> result <- lapply(split(seq(nrow(x)), x[,1]), function(.rows){
+ c(x[.rows[1],1], length(.rows), mean(as.numeric(x[.rows,3]))) #
return something
+ })
> result
$c1
[1] "c1" "4" "82.5"
$c2
[1] "c2" "1" "20"
$c3
[1] "c3" "1" "10"
> do.call(rbind, result)
[,1] [,2] [,3]
c1 "c1" "4" "82.5"
c2 "c2" "1" "20"
c3 "c3" "1" "10"
2008/4/17 zhihuali <lzhtom at hotmail.com>:
>
> Dear netters, suppose I have a matrix X [1,] 'c1' 'r6' '150'[2,] 'c1' 'r4' '70'[3,] 'c1' 'r2' '20'[4,] 'c1' 'r5' '90'[5,] 'c2' 'r2' '20'[6,] 'c3' 'r1' '10'I want to apply some funciton to groups of rows by the first column.If the function is just to calculate the average X[,3], it will be easy: aggregate(as.numeric(X[,3]),by=list(X[,1]),mean)But the function I want to use is more complicated. It will take as input a matrix(X[rows in the same group,c(2,3)], do some computation, and output another matrixwith the dimension (y,3), where y depends on the input. And I'd like the resultto be a rbind of each of the subset outputs. aggregate can not do that because the function is supposed to take vectors and output scalars.How can I apply the complicated function to groups of a matrix? Thanks! Zhihua Li
> _________________________________________________________________
> Windows Live Photo gallery 数码相机的超级伴侣，轻松管理和编辑照片，还能制作全景美图！
> http://get.live.cn/product/photo.html
> [[alternative HTML version deleted]]
>
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem you are trying to solve?
More information about the R-help
mailing list