Split xts oggetto per eventi

Question

Ho un oggetto xts con una manciata di indicatori di eventi. Devo dividerlo per determinati eventi, in modo che tutte le voci da un determinato evento fino all'evento successivo vengano salvate nello stesso xts, creando infine un elenco di oggetti xts, con ciascuno contenente l'evento come ultima voce senza altri eventi dello stesso tipo.

un esempio:

ts = as.Date(Sys.Date()-99:0) 
e1 = numeric(100);e1[10*1:10]=1 
e2 = numeric(100);e2[15*1:6]=1 
y = 1:100 # just a sample content 
xs = as.xts(cbind(e1,e2,y),order.by=ts) 

ee = e1*e2==1 # the event in which both e1 and e2 are 1, should happen at 30,60,90 

# here should be splitting function that gets xs and ee as parameters 
# and should return a list of 4 xts: the first with the entries 1 through 30, 
# the second with entries 31 to 60, the third with entries 61 to 90, and the last 
# with entries 91 to 100 

I vostri suggerimenti saranno molto apprezzati.

Answer 1

Utilizzare cumsum(ee) per creare la variabile di raggruppamento, quindi chiamare semplicemente split. È necessario apportare una piccola modifica all'uscita di cumsum, poiché si desidera che i valori TRUE siano l'ultima osservazione nel gruppo (anziché il primo).

split(xs, c(0,head(cumsum(ee),-1))) 
split(xs,rev(cumsum(rev(ee)))) # grouping factors reversed