Вычислите смежные пропорции
У меня есть кадр данных:
x <- data.frame(id = letters[1:3], val0 = 1:3, val1 = 4:6, val2 = 7:9)
# id val0 val1 val2
# 1 a 1 4 7
# 2 b 2 5 8
# 3 c 3 6 9
В каждой строке я хочу рассчитать соответствующие пропорции (отношение) для каждого значения. Например. для значения в столбце "val0", я хочу рассчитать весовые значения val0/(val0 + val1 + val2).
Требуемый вывод:
id val0 val1 val2
1 a 0.083 0.33 0.583
2 b 0.133 0.33 0.533
3 c 0.167 0.33 0.5
Может ли кто-нибудь сказать мне, что это лучший способ сделать это? Здесь это всего три столбца, но может быть много столбцов.
Ответы
Ответ 1
И еще одна альтернатива (хотя это в основном красивая версия sweep)... prop.table:
> cbind(x[1], prop.table(as.matrix(x[-1]), margin = 1))
id val0 val1 val2
1 a 0.08333333 0.3333333 0.5833333
2 b 0.13333333 0.3333333 0.5333333
3 c 0.16666667 0.3333333 0.5000000
В разделе "описание" файла справки в ?prop.table:
Это действительно sweep(x, margin, margin.table(x, margin), "/") для новичков, за исключением того, что если margin имеет длину ноль, тогда получается x/sum (x).
Итак, вы можете видеть, что внизу, это действительно очень похоже на решение @Jilber.
И... это приятно для разработчиков R, чтобы быть внимательными к нам новичками, не так ли?:)
Ответ 2
следует выполнить трюк
cbind(id = x[, 1], x[, -1]/rowSums(x[, -1]))
## id val0 val1 val2
## 1 a 0.08333333 0.3333333 0.5833333
## 2 b 0.13333333 0.3333333 0.5333333
## 3 c 0.16666667 0.3333333 0.5000000
Ответ 3
Другая альтернатива, использующая sweep
sweep(x[,-1], 1, rowSums(x[,-1]), FUN="/")
val0 val1 val2
1 0.08333333 0.3333333 0.5833333
2 0.13333333 0.3333333 0.5333333
3 0.16666667 0.3333333 0.5000000
Ответ 4
Функция ns_to_percents из пакета janitor выполняет следующее:
library(janitor)
ns_to_percents(x)
id val0 val1 val2
1 a 0.08333333 0.3333333 0.5833333
2 b 0.13333333 0.3333333 0.5333333
3 c 0.16666667 0.3333333 0.5000000
Это эквивалентно ns_to_percents(x, denom = "row"), хотя "row" является аргументом по умолчанию, поэтому в этом примере не требуется.
Если вы показываете результат, вы можете предпочесть janitor::adorn_crosstab.
Отказ от ответственности: я создал пакет janitor, но считаю целесообразным опубликовать его; функция была создана для выполнения именно этой задачи, делая код более понятным для чтения, и пакет можно установить из CRAN.
