Variogrammi

Variogrammi eli variogrammifunktio on 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äliset riippuvuus ilmaistaan spatiaalisen autokorrelaation avulla. Korrelaatiota ilmaisevan variogrammin rinnalla käytetään myös korrelogrammia ja kovarianssifunktiota. Kovarianssia käytetään sellaisten satunnaiskenttien yhteydessä, joilla on olemassa äärelliset odotusarvot ja toiset momentit. Jos toiset momentit ovat äärettömät, voidaan vielä käyttää variogrammia riippuvuuden määrän ilmaisemiseen. Variogrammeja käytettiin aluksi kaivosteollisuuden parissa geostatistiikassa, jonka puitteissa havaittiin, kuinka satunnaiskentän toinen momentti kasvoi kaivosalueen laajetessa liikaa. Variogrammien käyttö mahdollisti erityisesti kriging-estimoinnin, jolla suunniteltiin kaivoksen tuotantoa, mutta sen käyttö on levinnyt muillekin ympäristön suureita estimoiville aloille. Semivariogrammi on variogrammin puolikas, vaikka joissakin teksteissä nämä kaksi sekoitetaan toisiinsa. Tässä artikkelissa variogrammia merkitään 2\gamma ja semivariogrammia \gamma (lue: "gamma"). Vaikka variogrammilla on tärkeä asema geostatistiikan teoriassa, käytetään kuitenkin sovelluksissa lähinnä semivariogrammia. == Määritelmä == Satunnaiskentässä Z(x), missä x \in \mathbb{R}^n, eri kohdissa sijaitsevien satunnaismuuttujien erotuksen Z(y)-Z(x) varianssi (lyhennetään Var) antaa satunnaismuuttujien variogrammin arvon :2\gamma(x,y) = Var(Z(y)-Z(x)). Koska satunnaismuuttujan X varianssi määritellään toisen keskusmomentin avulla :\sigma_X^2=Var(X)=E[(X-\mu_X)^2], voidaan erotuksen varianssi myös ilmaista keskusmomentin avulla :2\gamma(x,y) =E[\\{(Z(y)-Z(x))-(\mu_y-\mu_x)\\}^2] = E[\\{(Z(y)-\mu_y)-(Z(x)-\mu_x)\\}^2]. Tätä kutsutaan yleisesti variogrammiksi. Sovelluksissa käytetään yleisemmin semivariogrammia, joka on variogrammin puolikas \gamma(x,y). Tätä variogrammin yleistä määritelmää voidaan edelleen tarkentaa, mikäli saadaan tarkempia tietoja satunnaiskentän stationäärisyydestä ja laadusta. Mikäli variogrammin arvot riippuvat vain satunnaismuuttujien välisestä etäisyysvektorista \bar h \in \mathbb{R},voidaan edelliset lausekkeet kirjoittaa :2\gamma(h) = Var(Z(x+\bar h)-Z(x)), tai :2\gamma(h) =E[\\{(Z(x+\bar h)-Z(x))-(\mu_{x+\bar h}-\mu_x)\\}^2] = E[\\{(Z(x+\bar h)-\mu_{x+\bar h})-(Z(x)-\mu_x)\\}^2]. == Stationääriset satunnaiskentät == Näin määriteltynä huomioidaan se mahdollisuus, että variogrammilla voi olla eri arvoja jokaisen lähiympäristön pisteparin välillä. Käytännön sovellutuksissa tällainen teoreettinen variogrammi ei ole käyttökelpoinen, sillä sen selville saaminen näytteitä tutkimalla on epätodennäköistä. Optimaalinen toive on, että satunnaiskenttä on jossakin määrin stationäärinen. On olemassa tilastollisia testejä riippuvuuksien paljastamiseksi. Eräs sellainen on Mantelin testi. === Stationäärisyyksiä === Vahva stationäärisyys () tarkoittaa sitä, että jokainen satunnaiskentän satunnaismuuttuja on jakaumaltaan identtinen toisen satunnaismuuttujan kanssa ja että mikä tahansa satunnaismuuttujien kokoelman yhteisjakauma on identtinen missä tahansa muualla avaruudellisesti samoin järjestäytyneen satunnaismuuttujaparven yhteisjakaumaan kanssa. Yksinkertainen seuraus on, että kaikki satunnaismuuttujien momentit ovat samat ja niistä erityisesti odotusarvo ja kovarianssi. Jos odotusarvo on äärellisinä olemassa ja se on E[Z(x)]=m, voidaan odotusarvo vähentää jokaisen satunnaismuuttujan arvosta pois Y(x)=Z(x)-m. Näin saadun modifioidun satunnaiskentän satunnaismuuttujien Y(x) odotusarvoksi tulee nolla E[Y(x)]=0. Kovarianssi pisteiden x ja y=x+h välille lasketaan silloin :C(h)=E[(Z(x)-m)(Z(x+h)-m)]=E[Y(x)Y(x+h)], jolloin varianssi tulee olla arvoltaan sama kuin kovarianssi etäisyydellä h=0 eli :C(0)=\sigma^2. Viimeinen ehto on myös kovarianssifunktion olemassaoloehto, sillä sitä ei ole olemassa ilman äärellistä varianssia. Vahvasti stationäärisen satunnaiskentän riippuvuusilmiö levittäytyy samanlaisena läpi avaruuden sen eri suunnissa ja on helposti määritettävissä tilastollisin keinoin. Käytännön sovellutuksissa tällaisia kenttiä ei kuitenkaan juuri tavata ja siksi työhypoteesiksi otetaan heikommin stationäärisiä satunnaiskenttiä. Heikosti stationäärinen () satunnaiskenttä ei ole yhtenäisesti jakautunut, koska satunnaiskentän satunnaismuuttujilla on vain samat odotusarvot ja varianssit (kaksi ensimmäistä momenttia) muttei samoja jakaumia. Tämä vastaa monien kenttätöiden tuloksia, jossa tilastoanalyysin tuloksena onnistutaan saadaan vain kaksi ensimmäistä momenttia. Muita satunnaiskentän piirteitä ei aina ole käytettävissä. Tarkemmissa geostatistisissa analyyseissä pyritään selvittämään vielä satunnaismuuttujien riippuvuuden määrä ja laatu. Lisäselvityksen tuloksena esitetään usein korrelogrammi tai kovarianssifunktio, jotka ilmaisevat pitoisuuden tai muun suureen jatkuvuuden määrää ja laatua. Tämäkään stationäärisyys ei kaikissa kenttätöissä toteudu. Esimerkiksi odotusarvo voi riippua paikasta eikä ole kaikkialla vakio ja eräillä kentillä satunnaismuuttujien varianssi kasvaa tutkimuskentän laajentuessa. Varianssi voi siis olla ääretön äärettömän laajalla tutkimuskentällä. Seuraava satunnaiskentän stationäärisyysehto onkin tällaisia tapauksia varten. Sisäisesti stationäärinen () satunnaiskenttä määritellään modifioidun satunnaiskentän avulla. Jos satunnaiskentän satunnaismuuttujien erotukset muodostavat heikosti stationäärisen satunnaiskentän, voivat erotusten varianssit olla äärellisinä olemassa. Näitä tapauksia varten lasketaan satunnaismuuttujien välisten riippuvuuksien määrä variogrammin avulla. Variogrammissa tietyn erotusvektorin erottamien satunnaismuuttujien varianssi muodostaa variogrammifunktion yhden arvon. Laskemalla kaikkien eri etäisyydellä olevien satunnaismuuttujien erotuksen varianssin saadaan variogrammin kaikki arvot. Koska variogrammi kehitettiin nimeen omaan sisäisesti stationäärisen satunnaiskentän riippuvuusfunktioksi, käytetään sitä siinä kovarianssifunktion sijasta. Se on kuitenkin myös käyttökelpoinen stationäärisemmissä satunnaiskentissä, joten seuraavassa käsitellään sen ominaisuuksia myös niissä. === Vahvasti stationääriset variogrammit === Vahvan stationäärisyyden satunnaiskentässä Z(x), jossa on olemassa äärellinen odotusarvo E[Z(x)]=m, voidaan odotusarvo vähentää jokaisen satunnaismuuttujan arvosta pois Y(x)=Z(x)-m. Tämän satunnaismuuttujien odotusarvo on silloin nolla E[Y(x)]=0. Jos satunnaiskentän Z(x) toiset momentit ovat äärellisinä olemassa, saadaan :E[(Z(x)-m)^2]=E[Y(x)^2]=\sigma^2 ja :C(h)=E[(Z(x)-m)(Z(x+h)-m)]=E[Y(x)Y(x+h)]. Variogrammi on tällaisessa satunnaiskentässä :2\gamma(h)= Var(Y(x+h)-Y(x))=E[((Y(x+h)-Y(x))-(m-m))^2]=E[(Y(x+h)-Y(x))^2] :::=E[Y(x+h)^2]-2E[Y(x+h)Y(x)]+E[Y(x)^2]=\sigma^2-2C(h)+\sigma^2 eli variogrammi on kovarianssifunktion yhdistelmäfunktio :2\gamma(h)=2[\sigma^2-C(h)]=2[C(0)-C(h)]. Semivariogrammi on silloin :\gamma(h)=\sigma^2-C(h)=C(0)-C(h). Yleensä kovarianssit saavat suurimman arvonsa arvolla h=0 ja muilla etäisyyden harvoilla pienempiä arvoja. Tämän vuoksi variogrammi on vahvasti stationäärisessä satunnaiskentässä ylhäältä rajoitettu funktio. Jos kovarianssi saa vain positiivisia arvoja, on semivariogrammin yläraja :\gamma(h)\le C(0)=\sigma^2, ja jos riippuvuus voi olla myös käänteistä ja |C(h)| \le C(0) niin yläraja on kaksinkertainen :\gamma(h)\le C(0)-(-C(0))=2C(0)=2\sigma^2. === Heikosti stationääriset variogrammit === Heikon stationäärisyyden satunnaiskentässä on vahvaan stationäärisyyteen verrattuna suuri ehtojen huojennus se, että vain satunnaiskentän odotusarvosta ja kovarianssifunktiosta on selvä käsitys, mutta jakaumien laadusta ei tiedetä mitään. Odotusarvot ovat tietyn rajatun ympäristön A sisällä vakio eli E[Z(x_0)]=m, kun x_0\subseteq A. Kovarianssin arvot ovat saman ympäristön sisällä :C(h)=E[Z(x_0)Z(x_0+h)]-m^2 ja sen vuoksi myös varianssi on äärellisenä olemassa. Koska vahvasti stationäärisessä satunnaiskentässä johdetaan odotusarvon, kovarianssifunktion ja variogrammin lausekkeet samanhenkisten olettamusten varassa, pätevät ne tässäkin stationäärisyyden hypoteesissä. Jos odotusarvo ei olekaan vakio, vaan muuttuu paikan mukaan muodostaen esimerkiksi pitoisuuksien trendin, eivät edellä kuvatut olettamukset pidä enää paikkaansa. Toisaalta varianssi voi olla arvoltaan jopa ääretön, jolloin kovarianssifunktio ei esitä riippuvuutta luotettavalla tavalla. Silloin tulee riippuvuusfunktiona käyttää variogrammia. === Sisäisen stationäärisyyden variogrammi === Mikäli lähekkäisissä pisteissä olevat satunnaismuuttujat muistuttavat laadultaan ja arvoiltaan toisiaan ja jos esimerkiksi samat odotusarvot, äärelliset toiset momentit tai varianssit eivät vaihtele satunnaiskentän eri osissa merkittävästi, voidaan silloin hyödyntää tilastollisesti koko satunnaiskentän mitattuja arvoja. Odotusarvot voivat muodostaa trendin eli odotusarvot ovat tietyssä suunnassa paikan funktio, joka näkyy jopa erotuksissa :E[Z(x+h)-Z(x)]=m(x). Tällaisissa tapauksissa vain variogrammi olemassa :Var((Z(x+h)-Z(x))^2)=2\gamma(h) ja se voidaan estimoida näyteotoksen perusteella koko satunnaiskentälle yhteisesti. Jos satunnaiskenttä on riittävän stationäärinen, riippuu variogrammi vain suuntavektorista \bar{h} (anisotropia) :2\gamma(\bar{h})=Var(Z(x+\bar{h})-Z(x)) tai, kun satunnaismuuttujien erotus on sama kaikissa suunnissa, variogrammi riippuu vain pisteiden välisestä etäisyydestä h = |x-y| (isotropia) :2\gamma(h)=Var(Z(y)-Z(x)). Sekä variogrammista että semivariogrammista käytetäänkin sovelluksissa kahta viimeiseksi esitettyä muotoa :\gamma(\bar{h})=\frac{1}{2}Var(Z(x)-Z(x+\bar{h})) ja :\gamma(h)=\frac{1}{2}Var(Z(y)-Z(x)). Yleensä satunnaiskenttä ei toteuta heikon stationäärisyysden ehtoja ja silloin variogrammi ei ole enää ylhäältä rajoitettu. Semivariogrammi muistuttaa silloin enemmän logaritmista \gamma(h)=3\alpha log|h| (de Wijsin malli) tai lineaarista \gamma(h)=A|h| mallia. == Variogrammin tulkinta == Variogrammi antaa tavallisesti pienimmän arvon satunnaismuuttujille, jotka sijaitsevat lähellä toisiaan. Niiden riippuvuus on suurta ja sen vuoksi niiden saamat arvot ovat lähellä toisiaan. Silloin satunnaismuuttujien erotuksen varianssi on myös pieni. Riippuvuus vähenee, kun satunnaismuuttujat etääntyvät toisistaan, jolloin variogrammin arvot suurenevat. Variogrammin arvo origossa onkin nolla ja se kasvaakin varsin yleisesti etäisyyden kasvaessa. == Variogrammin estimaattori == Variogrammin arvojen estimointi suoritetaan laskemalla sen arvoja näytteistä suoritetulla otoksella. Otos muodostetaan valikoimalla sopivan etäisyyden päässä toisistaan olevia näytteitä. Jos mitattava suure ei toista itseään (temperospatiaaliset satunnaismuuttujat) vaan se esiintyy pitkän aikaa samanlaisena, sanotaan sillä olevan vain yksi käytettävissä oleva reaalisaatio eli ulostulo. Silloin tulee satunnaiskentän olla stationäärinen ja kaikki näytteet edustavat tulkinnan mukaan samalla tavalla jakautuvia satunnaismuuttujia. Silloin näistä muodostettu variogrammi on luotettava riippuvuuden mitta. Satunnaismuuttujien Z(x) reaalisaatioita merkitään tässä z(x). === Määritelmä === Variogrammi estimoidaan mitatuista näytteistä laskemalla samalla etäisyydellä hsijaitsevien näytteiden arvojen erotukset z(x+h)-z(x). Sama etäisyys voidaan korvata samalla etäisyysvälillä, jolloin aineisto on etäisyyden mukaan luokiteltua. Variogrammin estimaatti 2\gamma^*(h) on näiden erotusten neliöiden keskiarvo. Semivarigrammin estimaatti \gamma^*(h) (merkitään myös \hat{\gamma}(h) ) on variogrammin puolikas: :\gamma^*(h)\approx \frac{1}{2N(h)}\sum (z(x+h)-z(x))^2, missä N(h) on etäisyydellä h olevien näytteiden lukumäärä. Jos satunnaiskenttä on anisotropinen, huomioidaan myös etäisyysvektorin suunta semivariogrammia laskettaessa. Silloin semivariogrammi on sekä etäisyyden että suunnan funktio. Estimoituihin semivariogrammin pisteisiin sovitetaan jokin ennalta valittu käyrä, joka antaa jatkuvan funktion ominaisuudet riippuvuuslaskuihin. === Esimerkki === Eräs varhaisia krigingin sovelluskohteita oli geologinen malminarviointi. Siinä kallioon kairataan syviä reikiä, joiden kiviytimet nostetaan maanpinnalle analysoitavaksi. Tutkimuksissa kairauksia on voitu tehdä satoja ja niiden analysoitujen ytimien osat muodostavat tuhansien näytteiden tietokannan. Koska kunkin näytteen sijainti kalliossa tunnetaan, voidaan niiden arvot hyödyntää riippuvuusanalyysissä ja krigingestimoinnissa. Näytteistä lasketaankin estimointia varten semivariogrammit, joista tutkitaan riippuvuuden laatu ja vaikutusala. == Variogrammien approksimaatiofunktioita == === Yleisiä ominaisuuksia === Huntu () tarkoittaa variogrammin ja semivariogrammin origossa olevaa epäjatkuvuuskohtaa ja se on seurausta satunnaiskentässä olevasta satunnaisuudesta, kun eri pisteiden arvojen välillä ei ole keskinäistä riippuvuutta, tai kun näytteiden suureiden arvoissa on mittausvirhettä. Kun lasketaan lähekkäisillä satunnaismuuttujilla semivariogrammin arvoja, saadaan raja-arvoksi huntua merkitsevä luku b :\lim_{h \to 0}\gamma(h)=b e0 vaikka \gamma(0)=0. Kynnys () tarkoittaa semivariogrammin saamaa maksimiarvoa niissä tapauksissa, kun satunnaiskenttä on vahvasti- tai heikosti stationäärinen, odotusarvo on sekä äärellinen että vakio ja varianssi on äärellisenä olemassa. Kynnysarvo C saavutetaan yleensä vaikutussäteen jälkeen tai asymptoottisesti hyvin kaukana :\lim_{h \to \infty}\gamma(h)=C. Vaikutussäde () kertoo etäisyyden, jonka jälkeen variogrammi tai semivariogrammi saavuttaa suurimman arvonsa. Semivariogrammilla kyseinen arvo on satunnaiskentän varianssi. Edelliset kolme ominaisuutta on mahdollista sisällyttää mihin tahansa sovituskäyrään, jolla approksimoidaan näytteistä saatua variogrammiaineistoa. Jos merkitään vaikutussäteen a sisällä variogrammifunktiota f(x), sen ulkopuolella kynnystä C ja origossa huntua b,saadaan paloittaisella määrittelyllä semivariogrammiksi :\gamma(h) = \begin{cases} 0, & \mbox{jos }h=0 \\\ b+f(h), & \mbox{jos }0 Semivariogrammilla voi edelleen olla vielä monta sisäkkäistä vaikutussädettä ja mutkikas lausekerakenne. === Paljon käytettyjä semivariogrammikäyriä === Alla on luettelo variogrammianalyysissä käytettävistä sovituskäyristä. Lausekkeissa vakio C tarkoittaa kynnyksen suuruutta, h etäisyyttä ja avaikutussädettä. Muut suureet, kuten esimerkiksi \alpha ja \beta, ovat käyrän sovittamisessa käytettäviä parametreja. semivariogrammien käyriä Käyrä Nimi Normi \delta Yhtälö \gamma(h) kynnys on C, ei vaikutussädettä puhdas huntu — \gamma(h)=\begin{cases} C, & \text{kun }h>0 \\\ 0, & \text{kun }h=0 \end{cases} käyrät, jotka saavauttavat kynnyksen C vaikutussäteen a etäisyydellä gaussin malli 1,731 \gamma(h)=C \left( 1-e^{-\left( \frac{h}{a} \right)^2} \right) kuutiollinen 1 \gamma\left(h\right)=C\left(7\left(\frac h a\right)^2-\frac{35}4\left(\frac ha\right)^3+\frac72\left(\frac ha\right)^5-\frac34\left(\frac ha\right)^7\right) eksponenttifunktion malli ≃ 2,996 \gamma(h)=C \left( 1-e^{-\frac{h}{a}} \right) pallomainen 1 \gamma(h)=\begin{cases} C \left( \frac 32 \frac{h}{a} - \frac 12 \left( \frac{h}{a} \right)^3 \right), & \text{kun }0 \leqslant h \leqslant a \\\ C, & \text{kun }h \geqslant a \end{cases} sinc-käyrä ≃ 20,371 \gamma\left(h\right)=C\left(1-\frac{sin\left(\frac ha\right)}{\frac ha}\right) epästationääristen satunnaiskenttien käyrät lineaarinen 1 \gamma(h)=C \frac ha potenssimalli 1 \gamma(h)=C \left(\frac ha\right)^b\ \text{, missä }\ 0 havemmin käytettävät käyrät eksponentiaalinen (yleinen) \sqrt[\alpha]{3} \gamma\left(h\right)=C\left(1-e^{-\left(\frac ha\right)^\alpha}\right) gamma \sqrt[\alpha]{20}-1 \gamma\left(h\right)=C\left(1-\frac1{\left(1+\frac ha\right)^\alpha}\right), \mbox{ missä }\alpha > 0 ensimmäisen lajin Besselin funktio 1 \gamma\left(h\right)=C\left(1-2^\alpha \Gamma\left(\alpha+1\right)\frac{J_\alpha\left(\frac ha\right)}{\left(\frac ha\right)^\alpha}\right),\mbox{ missä } \alpha > \frac{d}{2}-1 Matérnin kovarianssifunktiosta johdettu käyrä 1 \gamma\left(h\right)=C\left(1-\frac{\left(\frac ha\right)^\alpha}{2^{\alpha-1}\Gamma\left(\alpha\right)}K_{-\alpha}\left(\frac ha\right)\right), \mbox{ missä }\alpha > 0 Cauchyn yleistys \sqrt{\sqrt[\alpha]{20}-1} \gamma\left(h\right)=C\left(1-\left(1+\left(\frac ha\right)^2\right)^{-\alpha}\right), \mbox{ missä }\alpha > 0 == Lähteet == * * Heikkinen, Juha: Geostatistiikka , luentomoniste, Helsingin Yliopisto, 2006 === Viitteet === Luokka:Geostatistiikka