169 похожих чатов

Всем доброго понедельника! Есть у меня data.table из неотрицательных чисел такого

вида:
data.table(x = c(0, 0, 0, 0, 5:1), y = c(0, 0, 1:5, 0, 0), z = c(0, 1:8))

Нужно убрать все "ведущие" нули в каждом векторе (их каждый раз разное количество), образующиеся хвосты можно заменить нулями:
data.table(x = c(5:1, 0, 0, 0, 0), y = c(1:5, 0, 0, 0, 0), z = c(1:8, 0))

Решение у меня пока нашлось такое:
dt2 <-
dt1[, lapply(.SD, function(x)
shift(
x,
n = min(which(x > 0)),
fill = 0,
type = "lead"
))]

Вроде работает, но может я что-то упустил... Кто сможет подсказать?

5 ответов

17 просмотров

0. В целом, идея неплоха. 1. Немного режет which и min. 2. К чему здесь data.table? Если это вектора разной длины, каждый со своим жизненным циклом, то укладка их в прямоугольное преставление (а это тоже список векторов) создает ненужные ассоциации о связях в строках. Лучше честно писать, что это список векторов. Я бы делал так: library(data.table) dt1 <- data.table(x = c(0, 0, 0, 0, 5:1), y = c(0, 0, 1:5, 0, 0), z = c(0, 1:8)) dt1 %>% purrr::map(~.[cumsum(.) > 0])

Эдик-Амперян Автор вопроса
Ilya Shutov
0. В целом, идея неплоха. 1. Немного режет which и...

Илья, спасибо Чтение набора данных было из файла (бинарного, специфического формата "фортрановского типа") Поэтому образовалось это "прямоугольное представление", т.е. там уже нули были в каждом векторе предопределенной длины UPDATE: мне потом еще этот "список векторов" в файл отправить, поэтому я решил, что "прямоугольный" вариант data.table отправить через fwrite лучше, чем по списку бежать

Эдик Амперян
Илья, спасибо Чтение набора данных было из файла (...

вот и проявляется специфика задачи. для разового преобразования не суть важно решение, если оно успешное. само по себе это было похоже на разреженные матрицы. Если дальше предполагались вычисления с ними — возникли бы доп. вопросы.

Эдик-Амперян Автор вопроса
Ilya Shutov
вот и проявляется специфика задачи. для разового п...

Да, спасибо еще раз Специфика задачи есть, и цель еще иногда определяет средтва. Кстати, cumsum позволяет решить в стиле base R - если, например, только data.frame оперировать... и хвосты нулями забивать df1 <- data.frame( x = c(0, 0, 0, 1:5, 0, 0), y = c(0, 1:5, 0, 0, 0, 0), z = c(0, 0, 0, 0, 0, 1:5) ) df2 <- sapply( df1, function(x) { cs <- cumsum(x) c(x[cs > 0], x[cs <= 0]) } ) df3 <- as.data.frame(df2)

Эдик-Амперян Автор вопроса
Ilya Shutov
вот и проявляется специфика задачи. для разового п...

И вот еще один нюанс - мне попался вектор из нулей! min(which(0 > 0)) - это фейл! а вот sum(cumsum(0 <= 0)) - уже без ошибок

Похожие вопросы

Обсуждают сегодня

30500 за редактор? )
Владимир
47
а через ESC-код ?
Alexey Kulakov
29
Чёт не понял, я ж правильной функцией воспользовался чтобы вывести отладочную информацию? но что-то она не ловится
notme
18
У меня есть функция где происходит это: write_bit(buffer, 1); write_bit(buffer, 0); write_bit(buffer, 1); write_bit(buffer, 1); write_bit(buffer, 1); w...
~
13
Недавно Google Project Zero нашёл багу в SQLite с помощью LLM, о чём достаточно было шумно в определённых интернетах, которые сопровождались рассказами, что скоро всех "ибешни...
Alex Sherbakov
5
Ребят в СИ можно реализовать ООП?
Николай
33
Как передать управляющий символ в открытую через CreateProcess консоль? Собсна, есть процедура: procedure TRedirectThread.WriteData(Data: OEMString); var Written: Cardinal;...
Serjone
6
в JclConsole объявлено так: function CtrlHandler(CtrlType: DWORD): BOOL; stdcall; - где ваше объявление с stdcall? у вас на картинке нет stdcall
Karagy
8
https://github.com/erlang/otp/blob/OTP-27.1/lib/kernel/src/logger_h_common.erl#L174 https://github.com/erlang/otp/blob/OTP-27.1/lib/kernel/src/logger_olp.erl#L76 15 лет назад...
Maksim Lapshin
20
program test; {$mode delphi} procedure proc(v: int32); overload; begin end; procedure proc(v: int64); overload; begin end; var x: uint64; begin proc(x); end. Уж не знаю...
notme
6
Карта сайта