Arhn - архитектура программирования

dplyr: Как рассчитать кратность изменения внутри группы на основе значений в другом столбце

Мои текущие данные примерно имеют следующий шаблон:

Tree   Fertilized   Region   Fruits

apple  lightly      sunny    100
apple  lightly      dark     50
apple  heavily      sunny    300
apple  heavily      dark     200
pear   lightly      sunny    150
pear   lightly      dark     200
pear   heavily      sunny    300
pear   heavily      dark     150

Здесь я хочу рассчитать (как часть более крупной функции) кратное изменение размещения дерева на солнечном месте по сравнению с темным в каждой комбинации количества удобрений и типа дерева (например, 2-кратное изменение для слегка удобренных яблони):

df%<>%
  group_by(Tree,Fertilized) %>% 
  summarise(!!paste0("fold_change_", quote(Fruits)) := .[Region == "sunny","Fruits"]/.[type == "dark","Fruits"])

Однако я получаю сообщение об ошибке, говорящее о том, что столбца Fruits не существует. У кого-нибудь есть предложение о том, как заставить это работать? Я предполагаю, что решение заключается в небольшой настройке синтаксиса, но я не могу найти его ни сам, ни в Интернете.

Фактический набор данных имеет гораздо больше типов деревьев и параметров, таких как Fruits, поэтому я выбрал структуру конвейера и динамическую маркировку столбцов (!!paste0(), :=), которые могут быть актуальными или нерелевантными для решения этой проблемы.

Заранее спасибо всем, кто попытается помочь!

Привет, Роб

r data-wrangling dplyr
31.08.2020

  • так как проблема может исходить от вашего quosure/quote/!! использования, можете ли вы предоставить минимальный воспроизводимый пример, показывающий, как ваш код встроен в функцию? 31.08.2020

Ответы:


1

Я бы использовал операцию группировки:

library(data.table)
library(dplyr)


f <- tempfile()
writeLines("
Tree,  Fertilized,  Region,  Fruits,
apple, lightly, sunny, 100,
apple, lightly, dark, 50,
apple, heavily, sunny, 300,
apple, heavily, dark, 200,
pear, lightly, sunny, 150,
pear, lightly, dark, 200,
pear, heavily, sunny, 300,
pear, heavily, dark, 150
", f)
dat <- read.csv(f)

Таблица данных

dat <- data.table(dat)

dat[order(Region), .(fold_change = Fruits[2] / Fruits[1]), by=.(Tree, Fertilized)]
#>     Tree Fertilized fold_change
#> 1: apple    lightly        2.00
#> 2: apple    heavily        1.50
#> 3:  pear    lightly        0.75
#> 4:  pear    heavily        2.00

tidyverse

dat %>% 
  arrange(Region) %>%
  group_by(Tree, Fertilized)  %>%
        summarize(fold_change = Fruits[2] / Fruits[1])
#> `summarise()` regrouping output by 'Tree' (override with `.groups` argument)
#> # A tibble: 4 x 3
#> # Groups:   Tree [2]
#>   Tree  Fertilized fold_change
#>   <chr> <chr>            <dbl>
#> 1 apple " heavily"        1.5 
#> 2 apple " lightly"        2   
#> 3 pear  " heavily"        2   
#> 4 pear  " lightly"        0.75
31.08.2020
  • Большое спасибо! Такой подход делает свое дело. Возможно, мой первоначальный план доступа к значениям по значению, а не по индексу, на самом деле не работает. 01.09.2020
    • Новые материалы
    Коллекции публикаций по глубокому обучению
    Последние пару месяцев я создавал коллекции последних академических публикаций по различным подполям глубокого обучения в моем блоге https://amundtveit.com - эта публикация дает обзор 25..
    Представляем: Pepita
    Фреймворк JavaScript с открытым исходным кодом Я знаю, что недостатка в фреймворках JavaScript нет. Но я просто не мог остановиться. Я хотел написать что-то сам, со своими собственными..
    Советы по коду Laravel #2
    1-) Найти // You can specify the columns you need // in when you use the find method on a model User::find(‘id’, [‘email’,’name’]); // You can increment or decrement // a field in..
    Работа с временными рядами спутниковых изображений, часть 3 (аналитика данных)
    Анализ временных рядов спутниковых изображений для данных наблюдений за большой Землей (arXiv) Автор: Рольф Симоэс , Жильберто Камара , Жильберто Кейрос , Фелипе Соуза , Педро Р. Андраде ,..
    3 способа решить квадратное уравнение (3-й мой любимый) -
    1. Методом факторизации — 2. Используя квадратичную формулу — 3. Заполнив квадрат — Давайте поймем это, решив это простое уравнение: Мы пытаемся сделать LHS,..
    Создание VR-миров с A-Frame
    Виртуальная реальность (и дополненная реальность) стали главными модными терминами в образовательных технологиях. С недорогими VR-гарнитурами, такими как Google Cardboard , и использованием..
    Демистификация рекурсии
    КОДЕКС Демистификация рекурсии Упрощенная концепция ошеломляющей О чем весь этот шум? Рекурсия, кажется, единственная тема, от которой у каждого начинающего студента-информатика..