TITEL              : werte.vor, werte.vor2, werte.sd
KEYWORDS (statist) : Geostatistik; Universal-Kriging
ANWENDUNGSGEBIET   : Serie 2, Aufgabe 3
VERW.(VORL./J.-S.) : NDK 5/1997-99
Ueb.aufgabe 	   : "/u/sfs/kurse/ndk/ndk5/Ueb/Zeitreihen2/ueb-99-zeit-2c-2.tex"
--------------------------------------------------------------------------------
KURZE BESCHREIBUNG : Matrix fuer Contour-Plot und Perspektiv-Plot fuer den
Universal-Kriging Vorhersage (werte.vor, werte.vor2) und Standardfehler der
Vorhersage (werte.sd)
                     
                     
                     
SPEZ. EIG.         : 
                     
--------------------------------------------------------------------------------

Einlesen mit S-plus:


ko.vor2 <- read.table(pred.dat, header="T")

ko.vor2$okest[ko.vor$okest==-99] <- NA
ko.vor2 <- na.omit(ko.vor2)
werte.vor <- interp(ko.vor2$x, ko.vor2$y, ko.vor2$okest)

werte.vor2 <- werte.vor
werte.vor2$z[werte.vor$z==1] <- NA

ko.vor2$okvar <- sqrt(ko.vor2$okvar)
werte.sd <- interp(ko.vor2$x, ko.vor2$y, ko.vor2$okvar)

persp(werte.vor, box=F, zlim=c(-1,-0.65))
contour(werte.vor)
image(werte.vor)

...
                     
--------------------------------------------------------------------------------
* / LITERATUR      :  /
                     
                     
                   * D=Daten, vA=volle Analyse, tA=teilweise Analyse,
                     K=Kommentar

--------------------------------------------------------------------------------
BEMERKUNGEN        : Siehe kobalt.doc und pred.doc
                     
                     
--------------------------------------------------------------------------------
USERS              : R. Invernizzi
DATUM/UNTERSCHRIFT : Feb. 1999