2015-04-21 11 views
5

Ho un frame di dati impostato in questo modo:Calcolando una media in una cornice di dati sulla base di posizioni da colonne separate

N1 <- c(1,2,4,3,2,3,4,5,4,3,4,5,4,5,6,8,9) 
Start <- c("","Start","","","","","","","Start","","","","Start","","","","") 
Stop <- c("","","","","Stop","","","","","","Stop","","","","Stop","","") 

Con N1 essendo i miei dati di interesse. Vorrei calcolare la media di una stringa di numeri basata sulle posizioni "Start" e "Stop" nelle prossime due colonne.

Le stringhe come definiti da "Start" e "Stop" apparirebbe in questo modo:

2,4,3,2 
4,3,4 
4,5,6 

Quindi il mio risultato finale dovrebbe essere di 3 mezzi:

2.75,3.6,5 

risposta

5

si può provare:

mapply(function(start, stop){ 
      mean(N1[start:stop]) 
     }, 
     start=which(Start!=""), 
     stop=which(Stop!="")) 

#[1] 2.750000 3.666667 5.000000 
2

puoi provare anche rollapply

library(zoo) 
x <- sort(c(which(Stop != ""), which(Start != ""))) # indices of Start and Stop 
rollapply(x, 2, FUN = function(y) mean(N1[y[1]:y[2]]), by=2) 
[1] 2.750000 3.666667 5.000000 
4
library(data.table) # need latest 1.9.5+ 

# set up data to have all 1's column for the period we're interested in and 0 otherwise 
d = data.table(N1, event = cumsum((Start != "") - c(0, head(Stop != "", -1)))) 

d[, mean(N1), by = .(event, rleid(event))][event == 1, V1] 
#[1] 2.750000 3.666667 5.000000 

# or equivalently 
d[, .(event[1], mean(N1)), by = rleid(event)][V1 == 1, V2] 
Problemi correlati