АК
А под бубунту надо подождать 4.0, да? Когда будет?
Size: a a a
2020 April 25
A
Это видел. Я хочу не тестить, а так чтоб вжух и работало. До 20.04 обновился благополучно на одной пеке
A
получилось!
e
library(data.table)
mtcars <- as.data.table(mtcars)
mtcars[, .(p_value = shapiro.test(mpg)$p.value), by = .(am, vs)]
f <- function(dt, x, gr) {
dt[, .(p_value = shapiro.test(get(x))$p.value), by = gr]
}
f(mtcars, "mpg", c("am", "vs"))
вот аналог с
data.table. Почему-то выходит проще, чем в "простом" tidyverseСпасибо)
e
foo<-function(a,...){a %>% group_by_(...) %>% summarise(p_value = shapiro.test(mpg)$p.value)}
mtcars %>% foo(.,'am','vs')
И все работает
mtcars %>% foo(.,'am','vs')
И все работает
Вот в такой постановке да, когда вы уверены, что используете mtcars. А если у вас только рандомный датафрейм, и известно, что нужно найти нормальном ь первого количественнлго столбца по сочетанию бинарных факторов второго и третьего столбца?
АК
получилось!
Это временно. Должна быть нормальная репа со временем.
IS
este
Вот в такой постановке да, когда вы уверены, что используете mtcars. А если у вас только рандомный датафрейм, и известно, что нужно найти нормальном ь первого количественнлго столбца по сочетанию бинарных факторов второго и третьего столбца?
есть вот такое половинчатое решение, раздельная группировка и суммарайз с шапиро
my_summarise2 <- function(df, expr) {
expr <- enquo(expr)
summarise(df,
p_value = shapiro.test(!! expr)$p.value
)}
}
mtcars %>% group_by(am,vs) %>% my_summarise2(mpg)
Я постараюсь объединить в одну
my_summarise2 <- function(df, expr) {
expr <- enquo(expr)
summarise(df,
p_value = shapiro.test(!! expr)$p.value
)}
}
mtcars %>% group_by(am,vs) %>% my_summarise2(mpg)
Я постараюсь объединить в одну
АК
есть вот такое половинчатое решение, раздельная группировка и суммарайз с шапиро
my_summarise2 <- function(df, expr) {
expr <- enquo(expr)
summarise(df,
p_value = shapiro.test(!! expr)$p.value
)}
}
mtcars %>% group_by(am,vs) %>% my_summarise2(mpg)
Я постараюсь объединить в одну
my_summarise2 <- function(df, expr) {
expr <- enquo(expr)
summarise(df,
p_value = shapiro.test(!! expr)$p.value
)}
}
mtcars %>% group_by(am,vs) %>% my_summarise2(mpg)
Я постараюсь объединить в одну
В
dplyr когда-то ещё был do.e
есть вот такое половинчатое решение, раздельная группировка и суммарайз с шапиро
my_summarise2 <- function(df, expr) {
expr <- enquo(expr)
summarise(df,
p_value = shapiro.test(!! expr)$p.value
)}
}
mtcars %>% group_by(am,vs) %>% my_summarise2(mpg)
Я постараюсь объединить в одну
my_summarise2 <- function(df, expr) {
expr <- enquo(expr)
summarise(df,
p_value = shapiro.test(!! expr)$p.value
)}
}
mtcars %>% group_by(am,vs) %>% my_summarise2(mpg)
Я постараюсь объединить в одну
Спасибо!
IS
В
dplyr когда-то ещё был do.было такое) Ну и в практике - если бы я для себя писал - наверное оставил бы группирову нетронутой, указывать в кавычках поля для группировки в одной функции - как то не dplyr way)но чисто из принципа еще поковыряю
e
было такое) Ну и в практике - если бы я для себя писал - наверное оставил бы группирову нетронутой, указывать в кавычках поля для группировки в одной функции - как то не dplyr way)но чисто из принципа еще поковыряю
А хотите еще деталь?)
e
В функцию подается ТОЛЬКО датасет. то, что группировка по 2,3 столбцам, а шапиро - p значение - по первым - можно захардкодить
e
один аргумент
IS
А ну хардкодить то много ума не надо
foo1<-function(a){
f_shap <- function(x){st = shapiro.test(x); return(st$p.value)}
a %>% group_by_at(2:3) %>% summarise_at(vars(1),f_shap ) %>% rename_at(vars(3),
~str_replace(.,pattern = '.*',replacement = 'p_val')) }
mtcars %>% select(mpg,vs,am,everything()) %>% foo1()
foo1<-function(a){
f_shap <- function(x){st = shapiro.test(x); return(st$p.value)}
a %>% group_by_at(2:3) %>% summarise_at(vars(1),f_shap ) %>% rename_at(vars(3),
~str_replace(.,pattern = '.*',replacement = 'p_val')) }
mtcars %>% select(mpg,vs,am,everything()) %>% foo1()
IS
А извняюсь - группировка по 1 и 2 - шапиро по третьему - ну поменять не сложно
e
Круто))) спасибо)
e
А как бы сюда впихнуть еще переименование столбца, который тестируется на р_val?
IS
Добавил ренейминг И поставил группироввку корректную в сообщении выше
e
Добавил ренейминг И поставил группироввку корректную в сообщении выше
Спасибо большое!