Bir zaman serisinin bir bitişiklik matrisinin özfonksiyonları?


15

Basit bir zaman serisini düşünün:

> tp <- seq_len(10)
> tp
 [1]  1  2  3  4  5  6  7  8  9 10

bu zaman serileri için örnekler arasındaki zamansal bağlantıları temsil eden bir bitişiklik matrisini hesaplayabiliriz. Bu matrisi hesaplarken, 0 zamanında hayali bir alan ekliyoruz ve bu gözlem ile 1. zamandaki ilk gerçek gözlem arasındaki bağlantı, bağlantı 0 olarak bilinir. Zaman 1 ile zaman 2 arasında, bağlantı 1'dir ve bu şekilde devam eder. Zaman yönlü bir işlem olduğundan, siteler sitenin "yukarı akışındaki" bağlantılara (etkilenen) bağlanır. Bu nedenle her site 0 bağlantısına bağlıdır, ancak 9 bağlantı yalnızca 10 numaralı siteye bağlanmıştır; site 10 dışındaki her siteden sonra geçici olarak oluşur. Bu şekilde tanımlanan bitişiklik matrisi aşağıdaki gibi oluşturulur:

> adjmat <- matrix(0, ncol = length(tp), nrow = length(tp))
> adjmat[lower.tri(adjmat, diag = TRUE)] <- 1
> rownames(adjmat) <- paste("Site", seq_along(tp))
> colnames(adjmat) <- paste("Link", seq_along(tp)-1)
> adjmat
        Link 0 Link 1 Link 2 Link 3 Link 4 Link 5 Link 6 Link 7
Site 1       1      0      0      0      0      0      0      0
Site 2       1      1      0      0      0      0      0      0
Site 3       1      1      1      0      0      0      0      0
Site 4       1      1      1      1      0      0      0      0
Site 5       1      1      1      1      1      0      0      0
Site 6       1      1      1      1      1      1      0      0
Site 7       1      1      1      1      1      1      1      0
Site 8       1      1      1      1      1      1      1      1
Site 9       1      1      1      1      1      1      1      1
Site 10      1      1      1      1      1      1      1      1
        Link 8 Link 9
Site 1       0      0
Site 2       0      0
Site 3       0      0
Site 4       0      0
Site 5       0      0
Site 6       0      0
Site 7       0      0
Site 8       0      0
Site 9       1      0
Site 10      1      1

SVD, bu matrisin farklı zamansal ölçekler olarak varyasyonun Özfonksiyonlarına ayrışmasını sağlar. Gösterir ekstre işlevleri aşağıdaki şekil (dan SVD$u)

> SVD <- svd(adjmat, nu = length(tp), nv = 0)

Özfonksiyonlar

Özfonksiyonlar çeşitli zamansal ölçeklerde periyodik bileşenlerdir. Denemek tp <- seq_len(25)(veya daha uzun), bunu yukarıda gösterdiğim daha kısa örnekten daha iyi gösterir.

Bu tür analizlerin istatistiklerde uygun bir adı var mı? Tekil Spektrum Analizine benziyor, ancak bu gömülü bir zaman serisinin (sütunları zaman serisinin gecikmeli versiyonları olan bir matris) ayrışmasıdır.

Arka plan: Bu fikri, asimetrik Özvektör Haritaları (AEM) adlı uzamsal ekolojiden bilinen bir yöne sahip bir uzamsal süreci göz önünde bulunduran ve bir örneğin bağlanabileceği 1'ler içeren bir uzamsal örnek dizisi arasında bir bitişik matris oluşturan bir fikri değiştirerek ortaya çıktım. bir bağlantıya ve yapamayacağı bir 0'a, bağlantıların sadece "aşağı akım" 'a bağlanabilmesi - dolayısıyla analizin asimetrik doğasına. Yukarıda tarif ettiğim, AEM yönteminin tek boyutlu bir versiyonudur. İlgileniyorsanız , AEM yönteminin bir kopyasını burada bulabilirsiniz .

Şekil şu şekilde üretildi:

layout(matrix(1:12, ncol = 3, nrow = 4))
op <- par(mar = c(3,4,1,1))
apply(SVD$u, 2, function(x, t) plot(t, x, type = "l", xlab = "", ylab = ""),
      t = tp)
par(op)
layout(1)

Analiz nedir? Şunu mu demek istediniz: "Bu özfonksiyonlar için bir kullanım var mı?" "En" matrislerin özfonksiyonları, daha yüksek sayıların daha hızlı salınması ile salınır; ama size ne anlatıyorlar?
petrelharp

Hayır bu yaklaşım hakkında daha fazla bilgi edinin.
Monica'yı eski durumuna getirin. G. Simpson

Yanıtlar:


Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.