Контекст состоит в том, что у меня есть много временных рядов для многих id
и многих band
, и я включил выборку из девяти id
и двух band
. Здесь мы видим, что я могу легко построить временной ряд для всех id
s:
library(tidyverse)
df <- structure(list(id = c(1001L, 1001L, 1001L, 1001L, 1001L, 1001L, 1001L, 1001L, 1001L, 1001L, 1002L, 1002L, 1002L, 1002L, 1002L, 1002L, 1002L, 1002L, 1002L, 1002L, 1004L, 1004L, 1004L, 1004L, 1004L, 1004L, 1004L, 1004L, 1004L, 1004L, 1005L, 1005L, 1005L, 1005L, 1005L, 1005L, 1005L, 1005L, 1005L, 1005L, 1007L, 1007L, 1007L, 1007L, 1007L, 1007L, 1007L, 1007L, 1007L, 1007L, 1009L, 1009L, 1009L, 1009L, 1009L, 1009L, 1009L, 1009L, 1009L, 1009L, 1010L, 1010L, 1010L, 1010L, 1010L, 1010L, 1010L, 1010L, 1010L, 1010L, 1011L, 1011L, 1011L, 1011L, 1011L, 1011L, 1011L, 1011L, 1011L, 1011L, 1013L, 1013L, 1013L, 1013L, 1013L, 1013L, 1013L, 1013L, 1013L, 1013L), date = structure(c(1488884400, 1490612474, 1507460497, 1502276146, 1514372627, 1512644789, 1500980863, 1503572707, 1513940711, 1496660730, 1495796861, 1512644789, 1488884400, 1504436115, 1502276146, 1495796118, 1494068453, 1504868786, 1513940711, 1511780307, 1511348810, 1503572707, 1497524848, 1507028336, 1491476744, 1503572707, 1492340161, 1501844755, 1505300762, 1503140790, 1509620381, 1488884400, 1487156167, 1510052273, 1491476744, 1494068453, 1513940711, 1489748810, 1498388749, 1509620381, 1500980120, 1511780307, 1502708860, 1489748810, 1501412778, 1504436115, 1495796861, 1493204748, 1510484382, 1487156167, 1508324436, 1500548201, 1513940711, 1505732183, 1490612474, 1496660730, 1511348810, 1514372627, 1494068453, 1510052273, 1500548201, 1513076347, 1508756553, 1510484382, 1504436858, 1504004193, 1494932749, 1508324436, 1512644789, 1504868786, 1507460497, 1504004193, 1503140790, 1500980120, 1512212632, 1491476744, 1513940711, 1508756553, 1504436115, 1490612474, 1495796861, 1509188631, 1508756553, 1486292805, 1504004193, 1498388749, 1495796861, 1486292805, 1513940711, 1499684790), class = c("POSIXct", "POSIXt"), tzone = "UTC"), band = c("fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5", "fit1", "fit1", "fit1", "fit1", "fit1", "fit5", "fit5", "fit5", "fit5", "fit5"), value = c(0.496538754230172, 0.503271496428091, 0.97387311299285, 0.580658673638122, 0.55924511798107, 0.832069876834949, 0.669456383223215, 1.12835570514478, 0.650077806710299, 0.380956367547047, 0.315803532869213, 0.792491389890908, 0.542150595815071, 1.03016500582205, 0.761751198659722, 0.367933240661702, 0.478285303617102, 1.68901870452092, 0.740965064159661, 1.09028738312622, 0.822334909416119, 0.758342181009204, 0.404208383270466, 0.892795714415756, 0.452540219822814, 1.15220190981348, 0.522093412373678, 0.953592910857701, 1.27850667816495, 1.10756222303339, 0.722797148902218, 0.465842402588039, 0.524130056243481, 0.724757971315511, 0.401849347220063, 0.455169211763473, 0.736683498842155, 0.530595901306756, 0.598435246507131, 0.855911625573028, 0.459872179640563, 0.851473466057886, 0.600348304937791, 0.484896112230185, 0.491357621589034, 1.21884821937325, 0.408355867626313, 0.541537217668289, 1.20173675518489, 0.61126928681528, 1.02122136799224, 0.489289990779144, 0.829092258901136, 0.88152853467569, 0.528559966420024, 0.544164467022259, 1.15093592993106, 0.876559089290843, 0.582149928218707, 1.26592404446571, 0.479960992971744, 0.840894959543198, 1.00459298341354, 0.98285777345435, 0.754965044767638, 1.14971147250154, 0.678568628236206, 1.38981008816777, 0.989354634818581, 1.25116433808614, 1.2142398253614, 1.03201975237089, 0.928602154928637, 0.642961745200205, 0.842888403466734, 0.649606669375906, 0.724490820076092, 1.68294181717141, 1.83216850101507, 0.69741924948021, 0.268972923828825, 1.16584414990533, 1.20604228862346, 0.586060027904748, 1.16356144256577, 0.52670838257608, 0.382147314320451, 0.668308513834733, 0.78509264848017, 0.733357618207109)), row.names = c(NA, -90L), class = c("grouped_df", "tbl_df", "tbl", "data.frame"), vars = c("id", "band"), drop = TRUE, indices = list(0:4, 5:9, 10:14, 15:19, 20:24, 25:29, 30:34, 35:39, 40:44, 45:49, 50:54, 55:59, 60:64, 65:69, 70:74, 75:79, 80:84, 85:89), group_sizes = c(5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L), biggest_group_size = 5L, labels = structure(list(id = c(1001L, 1001L, 1002L, 1002L, 1004L, 1004L, 1005L, 1005L, 1007L, 1007L, 1009L, 1009L, 1010L, 1010L, 1011L, 1011L, 1013L, 1013L), band = c("fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5")), row.names = c(NA, -18L), class = "data.frame", vars = c("id", "band"), drop = TRUE, indices = list(0:4, 5:9, 10:14, 15:19, 20:24, 25:29, 30:34, 35:39, 40:44, 45:49, 50:54, 55:59, 60:64, 65:69, 70:74, 75:79, 80:84, 85:89), group_sizes = c(5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L), biggest_group_size = 5L, labels = structure(list(merge_id = c(1001L, 1001L, 1002L, 1002L, 1004L, 1004L, 1005L, 1005L, 1007L, 1007L, 1009L, 1009L, 1010L, 1010L, 1011L, 1011L, 1013L, 1013L), band = c("fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5", "fit1", "fit5")), row.names = c(NA, -18L), class = "data.frame", vars = c("merge_id", "band"), drop = TRUE)))
ggplot(df, aes(x = date, y = value, colour = band)) +
geom_point() +
geom_line() +
facet_wrap(~id)
Однако это становится громоздким, и графики становятся слишком маленькими, когда их слишком много id
, поэтому я хочу визуально проверить случайное подмножество. Я ожидаю, что следующее вернет только три из id
, но вместо этого мы получаем четыре id
, и мы даже не получаем все band
для каждого {{X4} }. Здесь я выбираю seed 1234, но вы получите другой результат, если продолжите повторный запуск с разными начальными числами, с разными компоновками комбинаций идентификаторов диапазонов.
set.seed(1234)
ggplot(
data = df %>% filter(id %in% sample(unique(df$id), 3)), # filtering to subset of 3 ids
mapping = aes(x = date, y = value, colour = band)
) +
geom_point() +
geom_line() +
facet_wrap(~id)
Обратите внимание, что это работает, если я делаю выборку вне вызова ggplot()
. (Это был бы желаемый результат)
set.seed(1234)
some_ids <- sample(unique(df$id), 3) # moved sample() outside of ggplot()
ggplot(
data = df %>% filter(id %in% some_ids),
mapping = aes(x = date, y = value, colour = band)
) +
geom_point() +
geom_line() +
facet_wrap(~id)
Почему это происходит? Я не вижу логической разницы между двумя вариантами. Это определенно связано с sample
, а не с частью unique(df$id)
, так как вы можете заменить это на c(1001, 1002, 1004, 1005, 1007, 1009, 1010, 1011, 1013)
и все равно получить проблему. Я также осознаю, что это может быть как-то связано с моими конкретными данными, поскольку я попытался создать аналогичное представление со встроенным набором данных, но я не могу представить, что это будет, поскольку это уже довольно ограниченное подмножество.
РЕДАКТИРОВАТЬ: например, я не могу воспроизвести эту ошибку, если использую этот еще более урезанный набор данных. Я сбит с толку, потому что я не могу сказать никакой разницы между этим набором данных и тем, что находится в моем dput
, за исключением фактических значений.
df2 <- tibble(
id = rep(1:9, each = 5, times = 2),
date = rep(seq(as.POSIXct("2018-01-01 00:00:00"), by = "month", length.out = 5), times = 18),
band = rep(c("b1", "b2"), each = 45),
value = c(rnorm(45, 0), rnorm(45, 1))
)
1 ответ
TL; DR: выражение фильтра оценивается несколько раз, поэтому вы не должны использовать недетерминированное выражение.
Не уверен, что этого достаточно для ответа, но если вы попытаетесь запустить свой пример с разными начальными числами, вы заметите, что количество диаграмм меняется с каждым начальным числом. Это говорит о том, что количество идентификаторов, которые мы фильтруем для фрейма данных, изменяется с каждым начальным значением, предполагая, что sample
на самом деле вызывается несколько раз. Мы можем подтвердить это, создав функцию, которая заменяет sample
:
sample_out <- function(data, n) {
print("running sample_out ")
return (sample(data, n))
}
А затем использовать его вместо sample
:
ggplot(
data = df %>% filter(id %in% sample_out(unique(df$id), 3)),
mapping = aes(x = date, y = value, colour = band)
)
Вы увидите, что sample_out
на самом деле вызывается несколько раз. В моем сеансе это вызывается 18 раз с данными выше, независимо от начального числа. Поэкспериментируя с разными размерами фреймов данных, кажется, что sample
будет вызываться (row_count / 5) раз. Это означает, что filter
каким-то образом оценивает свой аргумент несколько раз. Полный ответ объяснил бы, почему это происходит с filter
, но здесь я немного заблудился. Я считаю, что соответствующий источник находится здесь:
https://github.com/tidyverse/dplyr/blob/master/R/tbl-df.r#L55
filter.tbl_df <- function(.data, ..., .preserve = TRUE) {
// elided
out <- filter_impl(.data, quo)
filter_impl
в основном вызывает реализацию C ++, и я думаю, что ключевой строкой является:
https://github.com/tidyverse/dplyr/blob/master/src/filter.cpp#L408
template <typename SlicedTibble>
SEXP filter_template(const SlicedTibble& gdf, const NamedQuosure& quo) {
// elided
Proxy call_proxy(quo.expr(), gdf, quo.env()) ;
// elided
int ngroups = gdf.ngroups() ;
// elided
for (int i = 0; i < ngroups; i++, ++git) {
// elided
LogicalVector g_test = check_result_lgl_type(call_proxy.get(indices));
// elided
}
// elided
}
Обратите внимание, что для каждой группы таблицы выполняется call_proxy.get
. Я предполагаю, что мы видим, как sample_out
вызывается 18 раз, потому что в соответствующей таблице 18 групп.
В любом случае, на этот вопрос, вероятно, можно было бы быстро и авторитетно ответить, разместив сообщение в соответствующем контакте сообщества dplyr. В моем приключении, когда я узнал о dyplr, я не смог найти предупреждения об этом, так что, возможно, я что-то упустил. В документации dplyr
говорится, что его оценка немного отличается от той, которую можно использовать для: https://dplyr.tidyverse.org/articles/programming.html.
Большинство функций dplyr используют нестандартную оценку (NSE). Это универсальный термин, означающий, что они не следуют обычным правилам оценки R. Вместо этого они фиксируют введенное вами выражение и оценивают его произвольным образом.
Похожие вопросы
Связанные вопросы
Новые вопросы
r
R - это бесплатный язык программирования с открытым исходным кодом и программная среда для статистических вычислений, биоинформатики, визуализации и общих вычислений. Пожалуйста, предоставьте минимальные и воспроизводимые примеры вместе с желаемым результатом. Используйте dput () для данных и укажите все небазовые пакеты с вызовами library (). Не вставляйте изображения для данных или кода, вместо этого используйте блоки кода с отступом. По вопросам, связанным со статистикой, используйте https://stats.stackexchange.com.