Message-ID: <4CC089B0.8010208@cornell.edu>
Date: 2010-10-21T18:42:56Z
From: Zia Ahmed
Subject: Help: CDF plot
In-Reply-To: <4C75ECB4.60503@csiro.au>
I want to write a R- function to make a CDF plot of all realizations
(sim1, sim2, sim3........) with CDF of observed data (suppose zinc of
muse data). Help will be appreciated
Thanks
Zia
library(gstat)
data(meuse)
data(meuse.grid)
coordinates(meuse) = ~x + y
coordinates(meuse.grid) = ~x + y
lzn.vgm <- variogram(log(zinc) ~ 1, meuse)
lzn.fit <- fit.variogram(lzn.vgm, model = vgm(1, "Sph", 900, 1))
lzn.sim <- krige(log(zinc) ~ 1, meuse, meuse.grid, v.fit, nsim = 6, nmax
= 40)
/lzn.sim at data <- exp(lzn.sim at data)/
/ /
/ /
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20101021/9857b9c7/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: zua3.vcf
Type: text/x-vcard
Size: 281 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20101021/9857b9c7/attachment.vcf>