TITEL              : res.vor, res.vor2, res.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)

res.vor <- interp(ko.vor2$x, ko.vor2$y, ko.vor2$okest1)
persp(res.vor, box=F, zlim=c(-1,-0.65))
contour(res.vor)
image(res.vor)

res.vor2 <- res.vor
res.vor2$z[werte.vor$z==1] <- NA
persp(res.vor2)

## Vorhersagefehler für die Residuen
ko.vor2$ukvar <- sqrt(ko.vor2$ukvar)
res.sd <- interp(ko.vor2$x, ko.vor2$y, ko.vor2$ukvar)
persp(res.sd, box=F)
contour(res.sd)
image(res.sd)
 
...                    
--------------------------------------------------------------------------------
* / 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








