我有以下数据框
dat <- data.frame(x=c(1,2,3,3,2,1), y=c(3,4,4,5,2,5))
现在我想得到第三列,将y行值除以聚合y值(基于列x中的唯一值)。那么,我得到第1行如下:1,3,0.375; 0.375计算为3 /(5 + 3)。
我对R比较陌生,我希望你能帮助我。谢谢!
答案 0 :(得分:13)
有多种方法可以解决这个问题,这里有一个
with(dat, ave(y, x, FUN = function(x) x/sum(x)))
## [1] 0.3750000 0.6666667 0.4444444 0.5555556 0.3333333 0.6250000
这是另一种可能性
library(data.table)
setDT(dat)[, z := y/sum(y), by = x]
dat
# x y z
# 1: 1 3 0.3750000
# 2: 2 4 0.6666667
# 3: 3 4 0.4444444
# 4: 3 5 0.5555556
# 5: 2 2 0.3333333
# 6: 1 5 0.6250000
这是第三个
library(dplyr)
dat %>%
group_by(x) %>%
mutate(z = y/sum(y))
# Source: local data frame [6 x 3]
# Groups: x
#
# x y z
# 1 1 3 0.3750000
# 2 2 4 0.6666667
# 3 3 4 0.4444444
# 4 3 5 0.5555556
# 5 2 2 0.3333333
# 6 1 5 0.6250000
答案 1 :(得分:1)
当然,人们在SQL中思考的方式,在这种情况下非常冗长,但很好地推广了各种其他类似的问题:
library(sqldf)
dat <- sqldf("
with sums as (
select
x
,sum(y) as sy
from dat
group by x
)
select
d.x
,d.y
,d.y/s.sy as z
from dat d
inner join sums s
on d.x = s.x
")
答案 2 :(得分:1)
以下是一些基础R解决方案:
1)prop.table 使用prop.table
基础ave
功能,如下所示:
transform(dat, z = ave(y, x, FUN = prop.table))
,并提供:
x y z
1 1 3 0.3750000
2 2 4 0.6666667
3 3 4 0.4444444
4 3 5 0.5555556
5 2 2 0.3333333
6 1 5 0.6250000
2)总和这也有效:
transform(dat, z = y / ave(y, x, FUN = sum))