Как в R итерироваться по разницам между элементами вектора

Как вариант

Prelude> let sample = take 10 $ scanl (+) 0 [1,2..]
Prelude> sample
[0,1,3,6,10,15,21,28,36,45]
Prelude> let mysticFunction = id
Prelude> zipWith (-) (tail sample) sample
[1,2,3,4,5,6,7,8,9]
Prelude> map mysticFunction $ zipWith (-) (tail sample) sample
[1,2,3,4,5,6,7,8,9]

На R сам переведешь :)

надо прощее писать :)

sapply(1:(length(ivec)-1), function(x) ivec[x+1] - ivec[x])

vmyfun<-Vectorize(myfun)
vmyfun(ivec[-1], ivec[-length(ivec)])

Принудительная векторизация понятно дело необходимое только когда исходная функция не понимает векторные аргументы

можно еще вот так

mapply(myfun, ivec[-1], ivec[-length(ivec)])

ovec <- mapply('-', ivec[-1], ivec[-length(ivec)])

Именно это я и хотел, спасибо :)

надо прощее писать :)

Вот кстати научи делать слайсы...

Когда надо порезать data.frame с заведомо известными столбцами — волосы мягкие и шелковистые:

> df <- data.frame(time = c(1:10), data = c(10:19))
> df$data[df$time >= 4 & df$time <= 6]

А когда надо порезать матрицу, то даже с заведомо известными столбцами код превращается в мешанину:

> mx <- cbind(time = c(1:10), data = c(10:19))
> mx[mx[,"time"] >= 4 & mx[,"time"] <= 6, "data"]

На R сам переведешь :)

zipWith — вот это не осилил :(

но смотрится норм :)

используй пакет dplyr

Prelude> :t zipWith
zipWith :: (a -> b -> c) -> [a] -> [b] -> [c]

Берем функцию от двух аргументов и два списка, результат - список из поэлементного применения функции.

Вот кстати научи делать слайсы...

уже сказали про plyr

но есть еще полезный dsl для агрегации, это «сладкая парочка» melt и cast (лучше сразу reshape2)

то есть zipWith это mapply() ? :)

Именно это я и хотел, спасибо :)

для такой функции быстрее всё таки вот так

> "-"(1:9,2:10)
[1] -1 -1 -1 -1 -1 -1 -1 -1 -1

вопрос должен звучать как, «т.е. mapply это zipWith, который не проверяет типы»? Да, mapply это почти, что zipWith, во всяком случае так же применяет типы, однако далее есть разница, mapply не фьюзится (промежуточные структуры создаются), mapply не ленивый (ну, если честно, в R есть подобие ленивости, но серьёзно его рассматривать нельзя).

Можно же функцию diff() использовать.

1) mapply() не ленивый, поскольку вызывается сишная реализация без ленивости

2) а ленивость в R почти никто не учитывает когда пишет код :(

Можно же функцию diff() использовать.

для "-" да, но вопрос был для произвольной пользовательской функции?

вот например

http://cran.r-project.org/web/packages/lazyeval/vignettes/lazyeval.html

но вопрос был для произвольной пользовательской функции?

Да

ну и ещё один способ это library(iterators)

http://cran.r-project.org/web/packages/iterators/vignettes/writing.pdf

если не хватает памяти, и функции тяжелые можно параллелить вычисления

Для сравнения корректнее использовать microbenchmark всё же - повторения и случайность попыток нивелируют многие побочные эффекты.

случайность и попытки я обеспечиваю вручную :)

Так и запишем: не критичен к собственному состоянию :)