Функция dcast, принимающая аргументы из двух переменных значения

Скажем, у меня есть пример dataframe со следующей структурой

cars=c("A","A","A","A", "B","B","B","B", "C","C","C","C","A","A","A","A", "B","B","B","B", "C","C","C","C")
vendor=c("d","e","f","g", "d","e","f","g", "d","e","f","g", "d","e","f","g", "d","e","f","g", "d","e","f","g")
state=c(1,1,1,1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,2,2,2)
PS_mean=c(100, 110, 120, 130, 90, 95, 140, 180, 70, 80, 120, 150, 100, 110, 120, 130, 90, 95, 140, 180, 70, 80, 120, 150)
PS_stdv=c(10, 20, 30, 40, 10, 20, 30, 40, 10, 20, 30, 40, 10, 20, 30, 40, 10, 20, 30, 40, 10, 20, 30, 40)
mycars=data.frame(cars, vendor, state, PS_mean, PS_stdv)

Теперь я хочу применить изменение формы с помощью dcast, как

mycars_cov<-dcast(setDT(mycars[c('cars','state','PS_mean','PS_stdv')]), cars~state, value.var=c("PS_mean", "PS_stdv"), car_PS_var("PS_mean", "PS_stdv"))

Как видите, функция "car_PS_var" определяется пользователем с двумя входами

car_PS_var<- function(x,y){
   x<-as.numeric(x)
   y<-as.numeric(y)
   z=sd(x)*sd(y)/mean(x)
   return(z)
}

Я не знаю, как применить функцию, которая принимает два "value.var" в качестве аргументов и возвращает один. Обычно с помощью dcast вы можете применять функцию только к одной переменной, поэтому car_PS_var("PS_mean", "PS_stdv") не работает

В этой форме R выдаст несколько ошибок, потому что он не может принимать два входа в функции dcast.

Так как я могу правильно сделать это? Если вы предложите любой другой метод R, который выполняет задачу, также хорошо

1 ответ

Решение

Не уверен, что я понял вашу цель, но из моей интерпретации, быстрый и грязный способ состоит в том, чтобы сгруппировать по автомобилям и указать сначала, создать новый столбец, а затем представить новую таблицу данных.

mycars <- as.data.table(mycars)

temp <- mycars[, .(z = car_PS_var(PS_mean, PS_stdv)),
              by = c("cars", "state")]

dcast(temp, cars ~ state)

   cars        1        2
1:    A 1.449275 1.449275
2:    B 4.325825 4.325825
3:    C 4.545340 4.545340
Другие вопросы по тегам