У меня вопрос про subset. Как бы создать список из нескольких сабсетов из объекта phyloseq? Функция для сабсетов для таких объектов основана на обычном subset() для датафреймо...
# Тут несколько проблем: # (1) as.data.frame как таковой не нужен, без него тоже будет работать. # (2) return() как таковой не нужен, Хедли рекомендует его вставлять только ко...
У меня наверняка дурацкий вопрос. Как можно сопоставить две матрицы по именам строк и столбцов, чтобы они стали одного размера, а недостающие ячейки заполнились NA? # было A <...
А точно нет никакого такого хитрого джойна? А то у меня монстр получился... rework <- function(a = A, b = B){ add_row_filler <- function(a, b){ not_in_a <- setdiff(rowna...
Вдогонку. Пост хорош, как выжимка, но для начинающих рекомендую Richard McElreath в оригинале. Кроме книжки Statistical Rethinking и сопровождающего ее пакета rethinking, у не...
Скажите, господа, а существует ли метод, позволяющий визуализировать одновременно несколько таблиц частот? Такой как Correspondence Analysis, но с учетом блочной структуры сбо...
Не воспроизвести ошибку. Аналогичный код на искусственных данных работает. Что-то с данными? Или с версией ggplot2? library(dplyr) library(ggplot2) library(scales) set.seed(42...
Товарищи, кто-нибудь понимает в копирайте? - Пакет R распространяется по GNU GPL, значит ли это, что данные в нем тоже по этой лицензии распространяются, если не сказано иное?...
А важно именно доли моделировать? Можно, например, моделировать число успехов для каждого субъекта, а общее число испытаний взять в качестве оффсета. exp(коэффициенты) этой мо...
Интересно, похожи ли наши разногласия по поводу первой формулировки на то, что Эндрю Гельман описал в этом посте? Гельман считает, что упоминание в явном виде значений границ ...
Не хватает вариантов. Где же -> для олдскульных пользователей консоли? И потом <- и = не исключают, а дополняют друг друга, только юзкейсы разные)
А с примером-то работает? У меня там закрывающая скобка была пропущена, кстати. your_dfr <- data.frame(year = 1:3, m1 = c(2, 3, 3), m2 = c(3, 6, 8), m3 = c(4, NA, 12)) your_d...
А это не моветон, все пихать в один пайп - от открытия данных до картинки?
А почему бы не сохранить результаты в переменную и ее проверить?
А почему, кстати, не биномиальная GLM? В смысле glm (cbind(#у, #н) ~ предикторы, family = binomial) Она все кейсы взвесит по числу наблюдений.
М.б. просто усеченное среднее?
Подскажите как по-русски сказать studentized range distribution?
А почему такая разница между машинами?
Нет такого для Stepik?
Что значит сравнить и какая у них структура? Нужна более четкая формулировка задачи.