Tablodaki her satırı, satır zaman damgaları

Question

arasındaki farka göre başka bir tabloda bir satıra eşleştirin Her biri aynı sistemin ayrı özniteliklerini ölçen iki eşit olmayan aralıklı zaman dizim var. İki dizinin veri noktaları aynı zamanda örneklenmez ve seri aynı uzunlukta değildir. Her satırı, A serisinden, zamana en yakın B satırına eşleştirmek istiyorum. Aklımda bulunan şey, B'ye en yakın satırdaki indeksleri içeren bir A sütunu eklemektir. Her iki seride de Unix zamanında ölçülen bir zaman sütunu vardır (ör. 1459719755). örneğin

, verilen iki veri seti

a time 
2 1459719755 
4 1459719772 
3 1459719773 

b time 
45 1459719756 
2 1459719763 
13 1459719766 
22 1459719774 

birinci veri kümesi B[4,]$timeA[2,]$time ve A[3,]$time yakın bir değere sahiptir, A[1,]$time yakın bir değere sahiptir B[1,]$time yana

a time  index 
2 1459719755 1 
4 1459719772 4 
3 1459719773 4 

için güncellenmelidir.

Bunu yapmak için uygun bir yol var mı? Böyle

Answer 1

deneyin şey:

(1+ecdf(bdat$time)(adat$time)*nrow(bdat)) 
[1] 1 4 4 

Neden bu iş ki? Ecdf işlevi, 0 ile 1 arasında bir değere sahip olan başka bir işlevi döndürür. Bu, ecdf'a ilk bağımsız değişken tarafından tanımlanan değer dağılımındaki yeni bir değerin "olasılık aralığı" [0,1] 'de "konumunu" döndürür. İfade gerçekten sadece bu işlevin sonucunu [1, nrow (bdat)] aralığına yeniden ölçeklendirmektir. (Bence zarif sayılır.)

Başka bir yaklaşım, bdat$time'un sıralanmış değerleri üzerine approxfun kullanmak olacaktır; Bunların yuvarlatılması gerekebilir. Onları indis olarak kullanmak yerine tamsayıya keser.

apf <- approxfun(x=sort(bdat$time), y=seq(length(bdat$time)) ,rule=2) 
apf(adat$time) 
#[1] 1.000 3.750 3.875 
round(apf(adat$time)) 
#[1] 1 4 4 

Her iki durumda da "sipariş istatistiğinden" sıralı bir değer öngörüyorsunuz. İkinci durumda, bağların arzu ettiğiniz şekilde işlendiğini kontrol etmelisiniz.