title: Why are my predicted values constant?
body: <h1>mydata</h1>
  <pre><code>mydata_x &lt;- rep(seq(200,600,by=200),times=1,each=4)
  mydata_y &lt;- rep(seq(600,1200,by=200),times=4,each=1)
  mydata_z &lt;- c(0,0,0,0,0,0,0,529,0,0,0,0)
  mydata &lt;- data.frame(x=mydata_x,y=mydata_y,z=mydata_z)
  </code></pre>
  <p>I've taken samples from a virtual field and now want to predict neighboring values.  It's been done by the two methods described below... both predicting a constant value for each point.  I don't understand why.</p>
  <p>Keep in mind that I'm a new R user and love it.  I've been trying to follow an example I found online to perform my task.  Using the example dataset, I have no problem achieving their outcome.  However, I'm unable to apply it to my own data.  See the example here:  <a href="http://www.stat.ucla.edu/%7Enchristo/statistics_c173_c273/c173c273_lec11_w11.pdf" rel="nofollow noreferrer">http://www.stat.ucla.edu/~nchristo/statistics_c173_c273/c173c273_lec11_w11.pdf</a></p>
  <pre><code>require(geoR)
  
  #create x,y coordinates for locations to predict
  
  kx &lt;- rep(seq(250,700,by=100),times=1,each=length(seq(650,1250,by=100)))
  ky &lt;- rep(seq(650,1250,by=100),times=length(seq(250,700,by=100)),each=1)  
  k_df &lt;- data.frame(x=kx,y=ky) #coordinates as data.frame
  k_matrix &lt;- as.matrix(cbind(kx,ky)) #coordinates as.matrix
  
  #convert sample data (top of post) as.geodata
  
  b &lt;- as.geodata(mydata)
  
  #predict values
  
  prediction &lt;- ksline(b,cov.model = &quot;gaussian&quot;,cov.pars = c(10,3.33),locations=k_df)
  prediction$predict
  </code></pre>
  <p>All 35 outcomes return 44.0833 #I just don't get it!!!</p>
  <p>I don't understand how I am returning a single value across the field when there are so many zeros.  I would expect to see a single region of decreasing values.</p>
  <h1>Thoughts?</h1>