Kovarianssifunktio

**Kovarianssifunktio** on todennäköisyyslaskennassa, erityisesti stokastiikassa ja tilastotieteessä käytettävä riippuvuuden määrän ilmaisemiseen käytettävä mitta. Stokastisissa prosesseissa, jotka muistuttavat esimerkiksi aikasarjoja, käytetään autokorrelaatiota prosessin sisäisen riippuvuuden mittaamiseksi. Satunnaiskentissä, jotka muodostetaan moniston \mathbb{R}^n avulla, satunnaismuuttujien välillä olevan riippuvuuden määrää mitataan spatiaalisen autokorrelaation avulla. Riippuvuutta ilmaisevan kovarianssifunktion rinnalla käytetään myös korrelogrammia ja variogrammia.

## Määritelmä
### Kovarianssin määritelmä
Kahden satunnaismuuttujan X ja Y välinen riippuvuus määritellään kovarianssin avulla. Kovarianssi merkitään ja määritellään
:\sigma(X,Y)=\sigma_{XY}=cov(X,Y)=E[(X-E[X])(Y-E[X])].
Lausekkeet voidaan kertoa auki
:=E[(X-E[X])(Y-E[X])]=E[XY-XE[Y]-YE[X]+E[X]E[Y]].
Koska odotusarvo E on lineaarinen operaattori ja jos sen arvot merkitään E[X]=\mu_X ja E[Y]=\mu_Y, saadaan
:=E[XY-XE[Y]-YE[X]+E[X]E[Y]]=E[XY-Y\mu_X-X\mu_Y+\mu_X\mu_Y]
:=E[XY]-E[Y]\mu_X-E[X]\mu_Y+\mu_X\mu_Y=E[XY]-\mu_X\mu_Y.

### Satunnaiskentän kovarianssi
Satunnaiskenttien kohdassa x olevat tai indeksillä x merkityt satunnaismuuttujat voidaan esittää merkinnällä
:Z_x=Z(x). 
Silloin merkinnät Z(x) ja Z(y) tarkoittavat eri pisteissä x ja y olevia satunnaismuuttujia. Näiden välinen riippuvuus on ollut tapana merkitä
:\sigma(Z(x),Z(y))=\sigma(Z_x,Z_y)=\sigma(x,y)=\sigma_{xy}.
Käytännön sovelluksissa, jollainen on esimerkiksi geostatistiikka, satunnaiskenttiä Z(x) ovat esimerkiksi maaston alueita (x\in \mathbb{R}^2) tai maanalaisia tilavuuksia (x\in \mathbb{R}^3). Tutkittavaa ilmiötä vastaavan matemaattisen satunnaiskentän stationäärisyys voi olla niin vahvaa, että kovarianssi riippuu pelkästään joko pisteiden välisestä etäisyydestä tai etäisyyden suuntavektorista. Tällaisissa tapauksissa kovarianssin arvot ovat esimerkiksi etäisyyden ja suunnan funktioita \sigma(h,\alpha) tai pelkästään etäisyyden funktioita \sigma(h). Tämän tyyppisiä kovariansseja kutsutaan **kovarianssifunktioiksi** ja niitä merkitään \sigma:n sijasta tavallisesti C(h).

## Ominaisuuksia
### Peruspiirteet
Jos satunnaismuuttujan varianssin arvo merkitään \sigma^2(X)=\sigma^2, niin satunnaismuuttujan kovarianssi itsensä kanssa on \sigma(X,X)=E[XX]-\mu_X\mu_X=E[X^2]-\mu^2=\sigma^2 (varianssin määritelmän mukaan). Satunnaiskentissä, jossa satunnaismuuttujien riippuvuus heikkenee etäisyyden kasvaessa, on yleensä suurin riippuvuus lähekkäin sijaitsevilla satunnaismuuttujilla. Kovarianassifunktiossa tämä tapahtuu etäisyysarvoilla nolla eli C(0), jolloin kovarianssi saa varianssin arvon eli C(0)=\sigma^2.

Monesti halutaan, vaikka kokeellisesti eli tilastollisesti määritellyt kovarianssien arvot ovat tietyillä etäisyyksillä negatiivisia, kovarianssifunktion olevan positiividefiniitti eli C(h) \ge 0. Mikäli negatiivisia riippuvuuksia löytyy, tarkoittaa se, että toisen satunnaismuuttujan arvot vähenisi kun toisen kasvaisi, ja päinvastoin. Tämä on luonnonsovelluksissa harvinaista, mutta taloussovelluksissa yleisempää.

Satunnaiskentissä riippuvuus yleensä heikkenee etäisyyden kasvaessa. Jos riippuvuuden heikkeneminen ja etäisyyden kasvaminen on koko ajan samansuuntaista, on kovarianssi etäisyyden suhteen vähenevä funktio. Kovarianssifunktion maksimi eli varianssi saavutetaan silloin origossa C(0)=\sigma^2.

### Stationäärisyysoletuksien vaikutuksia
Satunnaiskentän vahva stationäärisyys takaa kaikille satunnaismuuttujille saman jakauman, samat momentit ja satunnaismuuttujien moniulotteiset jakaumat ovat invariantit paikan valinnalle. Nämä ehdot takaavat myös kovarianssifunktiolle saman arvon samalla etäisyydellä olevalle satunnaismuuttujaparille. Jos niiden toinen momentti on lisäksi äärellinen, on niillä äärellinen odotusarvo, joka on vakio, ja kovarianssifunktion arvo riippuu vain erotusvektorista. Tällaisia ideaalisia satunnaiskenttiä ei luonnosta kuitenkaan löydy.

Satunnaiskentän heikko stationäärisyys ei takaa samaa moniulotteista jakaumaa, joka on invariantti paikan valinnalle. Satunnaiskentän kaikilla satunnaismuuttujilla on kuitenkin äärellinen odotusarvo E[Z(x)]=\mu ja toinen momentti E[Z(x)^2]. Kovarianssifunktio riippu vain etäisyysvektorista. Muita satunnaiskentän ominaisuuksia ei taata. Niiden seurauksena satunnaiskentän varianssi on vakio \sigma^2=C(0) ja korrelogrammi on rajoitettuna olemassa
:\rho(h)=\frac{C(h)}{C(0)}. 
Kovarianssifunktio on näiden lisäksi etäisyysvektorin suhteen symmetrinen eli C(h)=C(-h) ja sen arvot ovat |C(h)|\le C(0). Kaikki edellä kerrotut vahvan ja heikon stationäärisyyden ominaisuudet ovat voimassa Gaussisella satunnaiskentällä.

Heikkokin stationäärisyys on käytännän kannalta erittäin vahva oletus ja tulee harvoin täytyttyä. Siksi on kehitetty sovelluksia varten sellainen stationäärisyysehto, jossa odotusarvo voi vaihdella ja toinen momentti voi olla ääretön.

Satunnaiskentän sisäinen stationäärisyys on määritelty niin, ettei kovarianssia voi käyttää, sillä se muuttuu negatiiviseksi. Kovarianssifunktion sijasta käytetään näissä tapauksissa variogrammia.

## Kovarianssin estimaattori
Yleensä kovarianssifunktio määritetään tutkimuskentältä z(x) määritetyistä näytteistä z_i. Aluksi tulee kovarianssifunktion määrittelyssä huomioida sekä etäisyys että suunta, joista muodostetaan omat luokat ja joiden kovarianssit lasketaan luokittain. Jos kovarianssi on kuitenkin samanlainen eri suunnissa, voidaan eri suuntien etäisyydet yhdistää samaan etäisyysluokkaan. Kullekin etäisyysluokalle h_klasketaan sitten kovarianssien arvot
:C(h_k)=\frac{\sum_{\forall i, j} (z_i-\mu)(z_j-\mu)}{n-1}.
Kun kaikki etäisyysluokat h_k on määritetty, sovitetaan niihin kovarianssifunktion lauseke niin, että käyrä kulkee mahdollisimman lähellä kaikkia estimoituja etäisyysluokkien kovarianssiarvoja. Tätä käyrää kutsutaan kokeelliseksi eli empiiriseksi kovarianssifunktioksi.

## Lähteet
* 
* Koistinen, Esko: Geomatematiikan menetelmiä ja sovelluksia malmivaratutkimuksissa, Tutkimusraportti nro 52, Geologinen tutkimuslaitos, 1981
* Heikkinen, Juha: Geostatistiikka , luentomoniste, Helsingin yliopisto, 2006
### Viitteet