2015-12-20 28 views
5

我試圖創建代表重複措施的變量「波浪」。具體來說,我試圖創建代表變量1 - 10,11 - 20 ... 91 - 100的平均值的連續變量。請注意,「...」代表波浪3到9的變量,因爲避免輸入這些是我的目標!突變多個/連續列(使用dplyr或base R)

下面是一個例子的數據幀,df,用10行100列:

mat <- matrix(runif(1000, 1, 10), ncol = 100) 
df <- data.frame(mat) 
dim(df) 
> 10 100 

我用dplyr功能mutate其中工程一旦所有的變量類型,但時間密集型和容易出錯。我一直沒能找到一種途徑,無需手動鍵入列的名稱這樣做,因爲我開始做如下(注意,「...」到9象徵波3):

df <- df %>% 
     mutate(wave_1 = (X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10)/10, 
      wave_2 = (X11 + X12 + X13 + X14 + X15 + X16 + X17 + X18 + X19 + X20)/10, 
      ... 
      wave_10 = (X91 + X92 + X93 + X94 + X95 + X96 + X97 + X98 + X99 + X100)/10) 

你可以mutate突變多個/連續列'dplyr'?其他方法也是受歡迎的。

+1

它是否必須與dplyr? –

+0

不,謝謝 - 另一種解決方案也會很棒 –

回答

5

這裏是與封裝zoo的一種方法:

library(zoo) 
t(rollapply(t(df), width = 10, by = 10, function(x) sum(x)/10)) 

這裏是與基礎R做到這一點的一種方法:

splits <- 1:100 
dim(splits) <- c(10, 10) 
splits <- split(splits, col(splits)) 
results <- do.call("cbind", lapply(splits, function(x) data.frame(rowSums(df[,x]/10)))) 
names(results) <- paste0("wave_", 1:10) 
results 

另一個用鹼R(G的禮貌非常簡潔的方式格羅騰迪克):

t(apply(df, 1, tapply, gl(10, 10), mean)) 

這裏是dplyr和的解決方案:

library(dplyr) 
library(tidyr) 
df$row <- 1:nrow(df) 
df2 <- df %>% gather(column, value, -row) 
df2$column <- cut(as.numeric(gsub("X", "", df2$column)),breaks = c(0:10*10)) 
df2 <- df2 %>% group_by(row, column) %>% summarise(value = sum(value)/10) 
df2 %>% spread(column, value) %>% select(-row) 
+0

謝謝,我會等待看看有沒有人使用dplyr做出回答,如果沒什麼問題的話 –

+1

@JoshuaRosenberg肯定,不需要匆忙 –

+0

@ G.Grothendieck相當不錯 –