<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from rtf -->
<style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<font face="Calibri, sans-serif" size="2">
<div>Hello</div>
<div> </div>
<div>The variograms function in gstat is returning semi-variance estimates which appear to be two orders of magnitude lower than what I expect from calculating the variance of the data in the attached space time object created with the STFDF function:</div>
<div> </div>
<div>library(gstat)</div>
<div>library(spacetime)</div>
<div>breaks<-seq(0.5,12.5,1)</div>
<div>var_all <- variogram(VMC~1,sp_time_var, bounds=breaks, cutoff = 12.5) # VMC is the numeric data to calculate a sp-time varigram</div>
<div>plot(var_all) # note the low values ranges on the scale <0.15</div>
<div>var(sp_time_var@data$VMC) # variance is 11.5</div>
<div> </div>
<div>So I would expect the variance in the upper part of the graph to approach a value around 11.</div>
<div> </div>
<div>Have I mis-specified the variogram function? I know from gstat that one had to specify a gstat object, but here I thought specifying the spacetime object was the right way to proceed and that VMC~1 was the correct syntax?</div>
<div> </div>
<div>Any advice much appreciated.</div>
<div>Barry</div>
<div> </div>
<div> </div>
<div> </div>
</font>
<br />--
<br />This message (and any attachments) is for the recipient only. NERC
<br />is subject to the Freedom of Information Act 2000 and the contents
<br />of this email and any reply you make may be disclosed by NERC unless
<br />it is exempt from release under the Act. Any material supplied to
<br />NERC may be stored in an electronic records management system.</body>
</html>