apply()并计算所有dataframe行第一行的比例
apply()并计算所有dataframe行第一行的比例
提问于 2011-04-06 10:18:42
我有一个数据框，如下所示，列出了按车辆类型划分的受伤人数：
trqldnum <- data.frame(motorveh=c(796,912,908,880,941,966,989,984),
motorcyc=c(257,295,326,313,403,389,474,496),
bicyc=c(109,127,125,137,172,146,173,178))
trqldnum
# motorveh motorcyc bicyc
#1 796 257 109
#2 912 295 127
#3 908 326 125
#4 880 313 137
#5 941 403 172
#6 966 389 146
#7 989 474 173
#8 984 496 178目前，我正在计算每个车辆类型的第一行的比例，使用：
trqldprop <- t(apply(trqldnum,1,function(x) {
x/c(trqldnum[1,1],trqldnum[1,2],trqldnum[1,3])
}))
trqldprop
# motorveh motorcyc bicyc
#1 1.000000 1.000000 1.000000
#2 1.145729 1.147860 1.165138
#3 1.140704 1.268482 1.146789
#4 1.105528 1.217899 1.256881
#5 1.182161 1.568093 1.577982
#6 1.213568 1.513619 1.339450
#7 1.242462 1.844358 1.587156
#8 1.236181 1.929961 1.633028这看起来有点丑陋，如果数据改变了形状，我需要手动更改函数的分母。如果我尝试在apply()语句中仅使用以下内容，我最终会得到一个列表列表中的输出。
function(x) x/c(trqldnum[1,])我更喜欢以上面的数据帧结果结束，但我只是在试图弄清楚它。
回答 4
Stack Overflow用户
回答已采纳
发布于 2011-04-06 23:54:04
怎么样
sweep(trqldnum,2,unlist(trqldnum[1,]),"/")unlist需要将数据帧的第一行转换为可以扫描的向量……
Stack Overflow用户
发布于 2011-04-06 10:41:34
将数据帧转换为矩阵并使用矩阵运算：
m <- as.matrix(trqldnum)
trqldprop <- as.data.frame( t(t(m)/m[1,]) )
> trqldprop
motorveh motorcyc bicyc
1 1.000000 1.000000 1.000000
2 1.145729 1.147860 1.165138
3 1.140704 1.268482 1.146789
4 1.105528 1.217899 1.256881
5 1.182161 1.568093 1.577982
6 1.213568 1.513619 1.339450
7 1.242462 1.844358 1.587156
8 1.236181 1.929961 1.633028请注意，我们需要转置矩阵(请参阅t(m))，因为当您将矩阵除以向量时，操作是按列完成的。
Stack Overflow用户
发布于 2011-04-06 11:35:14
我喜欢这些任务的plyr，因为它们允许您指定输出的格式。您可以将其转换为一个函数，该函数可以轻松地扩展到更多的列和不同的划分基准面。
FUN <- function(dat, baseRow = 1){
require(plyr)
divisors <- dat[baseRow ,]
adply(dat, 1, function(x) x / divisors)
}
FUN(trqldnum, 1)
motorveh motorcyc bicyc
1 1.000000 1.000000 1.000000
2 1.145729 1.147860 1.165138
3 1.140704 1.268482 1.146789
4 1.105528 1.217899 1.256881
5 1.182161 1.568093 1.577982
6 1.213568 1.513619 1.339450
7 1.242462 1.844358 1.587156
8 1.236181 1.929961 1.633028页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/5560772
复制相关文章
相似问题
腾讯云开发者
