<div dir="ltr"><div>Dear All,</div><div>I used gstat package of R to fit variogram model to experimental variogram having normalized the raw data as</div><div><font color="#000000" size="3" face="Times New Roman">

</font></div><p style="margin:0cm 0cm 10pt" class="MsoNormal"><font color="#000000"><span style="line-height:115%;font-size:12pt"><font face="Calibri">VN = (VP -VMean)/VSd</font></span><span style="line-height:115%;font-family:"Calibri","sans-serif";font-size:11pt">
 
 
  
  
  
  
  
  
  
  
  
  
  
  
 
 
 

 
</span></font></p><div><font color="#000000" size="3" face="Times New Roman">

</font></div><p style="margin:0cm 0cm 10pt" class="MsoNormal"><span style="line-height:115%;font-size:12pt"><font color="#000000"><font face="Calibri">Where VN is the normalized yield value, VP is the
original yield value at sampling point, VMean is average yield<span>  </span>value, and VSd is the standard deviation of
yield value.</font></font></span></p><p style="margin:0cm 0cm 10pt" class="MsoNormal"><span style="line-height:115%;font-size:12pt"><font color="#000000"><font face="Calibri">The reason for normalization is to be able to assess the stability of spatial pattern of yield over a number of years. Having normalized the data, I expected the sill i.e. the variance of the data to be one but to my surprise the semivariogram plot shows sill of 0.8. What is responsible for this?</font></font></span></p>
<p style="margin:0cm 0cm 10pt" class="MsoNormal"><span style="line-height:115%;font-size:12pt"><font color="#000000" face="Calibri">The result of fitting exponential function showed a nugget effect of zero as showed below</font></span></p>
<span style="line-height:115%;font-size:12pt"><span style="font:13px/15px "Lucida Console";color:rgb(0,0,0);text-transform:none;text-indent:0px;letter-spacing:normal;word-spacing:0px;white-space:pre-wrap;border-collapse:separate;font-size-adjust:none;font-stretch:normal;background-color:rgb(225,226,229)"><pre style="margin:0px;line-height:1.2;font-family:"Lucida Console";font-size:10pt!important;white-space:pre-wrap!important;outline-style:none" tabindex="0">
<span style="color:blue">## Exponential variogram model
</span><span style="color:blue;white-space:pre">> </span><span style="color:blue">exp.var<-vgm(psill = 0.8, model="Exp", range= 15, nugget= 0.1)
</span><span style="color:blue;white-space:pre">> </span><span style="color:blue">exp.mod<-fit.variogram(empvar,model=exp.var)
</span><span style="color:blue;white-space:pre">> </span><span style="color:blue">print(exp.mod)
</span>  model     psill    range
1   Nug 0.0000000  0.00000
2   Exp 0.8066099 11.36279</pre><div> </div><div>Does it mean the exponential structure explained 100 percent of total variation in a normalized yield variable? I am confused! Have a look at the attached file for the semivariogram plot. Do I need to detrend before fitting the model? it seems the plot increases monotonically without level off. Please advise.</div>
<div> </div><div>Thank you</div><div>Moshood</div></span><div></div></span><div><font color="#000000" size="3" face="Times New Roman">

</font></div><div><font color="#000000" size="3" face="Times New Roman"></font> </div></div>