<div dir="ltr">Dear Jon and Roman<div><br></div><div>Thanks for the reply.</div><div><br></div><div>Im attaching the code.</div><div><br></div><div>You may be right that the variogram may be problematic but it fitted well visually.</div><div><br></div><div>I think also that the prediction grid was not well created. I wiould be grateful to have to hints.code how to create a useful grid and also for the spplot or other alternatives</div><div><br></div><div>I can send more info if required.</div><div><br></div><div>Cheers</div><div><br></div><div>Anand</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 4, 2015 at 1:39 PM, Jon Skoien <span dir="ltr"><<a href="mailto:jon.skoien@jrc.ec.europa.eu" target="_blank">jon.skoien@jrc.ec.europa.eu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    Hi,<br>
    <br>
    You have attached some data, but we still dont know exactly what you
    do with the data (do you use the logarithm, how many bins, did you
    project the data set...?). However, having had a quick look using
    the LatLong coordinates (not really recommended), I can see that the
    variogram from this data set is not particularly well defined, and
    it could well be that you fitted a model with pure nugget effect. In
    that case you will get the same prediction everywhere and no color
    differences. <br>
    <br>
    Have a look at <br>
    > summary(map)<br>
    and particularly the summary of var1.pred. <br>
    You could also plot the variogram of the data together with your
    model, vmf, and see how they fit. <br>
    <br>
    I am not sure what you mean by not being able to create a prediction
    grid with data included. Do you mean that you want a grid which
    exactly overlays some of the data points? Matching one point should
    not be difficult, 2 and 3 should be possible, but you might get some
    problems with precision (FAQ 7.31). <br>
    <br>
    Best wishes,<br>
    Jon<div><div class="h5"><br>
    <br>
    <br>
    <br>
    <div>On 8/3/2015 5:30 PM, Anand Sookun
      wrote:<br>
    </div>
    <blockquote type="cite">
      <div dir="ltr">Hi
        <div><br>
        </div>
        <div>Thank you Roman for your kind response.</div>
        <div><br>
        </div>
        <div>The data is as attached which is part of the whole dataset.</div>
        <div><br>
        </div>
        <div>My other problem is that im not able to create an
          appropriate prediction grid with some data included.</div>
        <div><br>
        </div>
        <div>The codes which failed to show colour ramps are</div>
        <div><br>
        </div>
        <div>map <- krige(logcarbseq ~ 1, locations = carbon_seq,
          newdata = map,model = vmf)</div>
        <div>str(map)</div>
        <div>print(spplot(map, "var1.pred", asp=1,
          col.regions=bpy.colors(64),main="OK prediction, carbon")) </div>
        <div>print(spplot(map,
          "var1.var",col.regions=cm.colors(64),asp=1, main="OK
          prediction variance, carbon"))</div>
        <div><br>
        </div>
        <div>I look forward to hopefully get some solutions.</div>
        <div><br>
        </div>
        <div>Thanks again and warm regards</div>
        <div><br>
        </div>
        <div>Anand</div>
      </div>
      <div class="gmail_extra"><br>
        <div class="gmail_quote">On Mon, Aug 3, 2015 at 4:09 PM, Roman
          Luštrik <span dir="ltr"><<a href="mailto:roman.lustrik@gmail.com" target="_blank">roman.lustrik@gmail.com</a>></span>
          wrote:<br>
          <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
            <div dir="ltr">What is the variability of the plotted
              variable?
              <div><br>
              </div>
              <div>Cheers,</div>
              <div>Roman</div>
            </div>
            <div class="gmail_extra"><br>
              <div class="gmail_quote"><span>On Sun, Aug 2,
                  2015 at 11:03 AM, Anand Sookun <span dir="ltr"><<a href="mailto:asookun@gmail.com" target="_blank"></a><a href="mailto:asookun@gmail.com" target="_blank">asookun@gmail.com</a>></span>
                  wrote:<br>
                </span>
                <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span>Hi<br>
                    <br>
                    I have kriged my data but the map shows uniform
                    colours - no gradual<br>
                    contour colour plot on map<br>
                    <br>
                    code: print(spplot(k40, "var1.pred", asp=1,<br>
                    col.regions=bpy.colors(64),main="OK prediction,
                    log-ppm Zn"))<br>
                    <br>
                    thnaks for any help<br>
                    <br>
                    Anand Sookun<br>
                    --<br>
-------------------------------------------------------<br>
                  </span>
                  *Anand Sookun (**MPhil. **PhD Candidate, University of
                  Mauritius)*<br>
                  *Director cum Freelance Consultant*<br>
                  *Soft Sustainable Services (Mauritius)*<br>
                  *Research and Consultancy*<span><br>
                    <br>
                    Website: <a href="http://sustain-ability.webs.com/" rel="noreferrer" target="_blank">http://sustain-ability.webs.com/</a><br>
                    <br>
                    <br>
                    <a href="https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic" rel="noreferrer" target="_blank">https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic</a><br>
                    <br>
                  </span>
                          [[alternative HTML version deleted]]<br>
                  <br>
                  _______________________________________________<br>
                  R-sig-Geo mailing list<br>
                  <a href="mailto:R-sig-Geo@r-project.org" target="_blank">R-sig-Geo@r-project.org</a><br>
                  <a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" rel="noreferrer" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a><span><font color="#888888"><br>
                    </font></span></blockquote>
              </div>
              <span><font color="#888888"><br>
                  <br clear="all">
                  <div><br>
                  </div>
                  -- <br>
                  <div>In God we trust, all others bring data.</div>
                </font></span></div>
          </blockquote>
        </div>
        <br>
        <br clear="all">
        <div><br>
        </div>
        -- <br>
        <div>
          <div dir="ltr">
            <div>
              <div dir="ltr">
                <div>
                  <div dir="ltr">
                    <div>
                      <div dir="ltr">
                        <div>
                          <div dir="ltr">
                            <div>-------------------------------------------------------</div>
                            <div><font color="#6fa8dc" size="1"><i><b>Anand
                                    Sookun (</b></i><i><b>MPhil. </b></i><i><b>PhD
                                    Candidate, University of Mauritius)</b></i></font></div>
                            <div><i><b><font color="#6fa8dc" size="1">Director
                                    cum Freelance Consultant</font></b></i></div>
                            <div><i><b><font color="#6fa8dc" size="1">Soft
                                    Sustainable Services (Mauritius)</font></b></i></div>
                            <div><i><b><font color="#6fa8dc" size="1">Research
                                    and Consultancy</font></b></i></div>
                            <div><font color="#6fa8dc" size="1"><br>
                              </font></div>
                            <font color="#6fa8dc" size="1">Website: <a href="http://sustain-ability.webs.com/" target="_blank"></a><a href="http://sustain-ability.webs.com/" target="_blank">http://sustain-ability.webs.com/</a><span style="font-family:helvetica,arial;line-height:12px;text-align:center"> </span></font></div>
                          <div dir="ltr"><font color="#6fa8dc" size="1"><span style="font-family:helvetica,arial;line-height:12px;text-align:center">
                                <div dir="ltr"><br>
                                </div>
                              </span></font></div>
                          <div dir="ltr">
                            <div style="text-align:center"><font color="#6fa8dc" size="1" face="helvetica, arial"><span style="line-height:12px"><br>
                                </span></font></div>
                            <div><span style="text-align:center"><font color="#000000" face="helvetica,
                                  arial"><span style="font-size:12px;line-height:12px"><a href="https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic" target="_blank"></a><a href="https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic" target="_blank">https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic</a></span></font><br>
                              </span></div>
                            <div><span style="text-align:center"><font color="#000000" face="helvetica,
                                  arial"><span style="font-size:12px;line-height:12px"><br>
                                  </span></font></span></div>
                          </div>
                        </div>
                      </div>
                    </div>
                  </div>
                </div>
              </div>
            </div>
          </div>
        </div>
      </div>
      <br>
      <fieldset></fieldset>
      <br>
      <pre>_______________________________________________
R-sig-Geo mailing list
<a href="mailto:R-sig-Geo@r-project.org" target="_blank">R-sig-Geo@r-project.org</a>
<a href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a>
</pre>
    </blockquote>
    <br>
    </div></div><span class="HOEnZb"><font color="#888888"><pre cols="72">-- 
Jon Olav Skøien
Joint Research Centre - European Commission
Institute for Environment and Sustainability (IES)
Climate Risk Management Unit

Via Fermi 2749, TP 100-01,  I-21027 Ispra (VA), ITALY

<a href="mailto:jon.skoien@jrc.ec.europa.eu" target="_blank">jon.skoien@jrc.ec.europa.eu</a>
Tel:  <a href="tel:%2B39%200332%20789205" value="+390332789205" target="_blank">+39 0332 789205</a>

Disclaimer: Views expressed in this email are those of the individual and do not necessarily represent official views of the European Commission. </pre>
  </font></span></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>-------------------------------------------------------</div><div><font color="#6fa8dc" size="1"><i><b>Anand Sookun (</b></i><i><b>MPhil. </b></i><i><b>PhD Candidate, University of Mauritius)</b></i></font></div><div><i><b><font color="#6fa8dc" size="1">Director cum Freelance Consultant</font></b></i></div><div><i><b><font color="#6fa8dc" size="1">Soft Sustainable Services (Mauritius)</font></b></i></div><div><i><b><font color="#6fa8dc" size="1">Research and Consultancy</font></b></i></div><div><font color="#6fa8dc" size="1"><br></font></div><font color="#6fa8dc" size="1">Website: <a href="http://sustain-ability.webs.com/" style="margin:0px;padding:0px;border:0px;font-family:'Helvetica Neue',helvetica,sans-serif;font-weight:700;line-height:1.4em;vertical-align:baseline;text-align:center" target="_blank">http://sustain-ability.webs.com/</a><span style="font-family:helvetica,arial;line-height:12px;text-align:center"> </span></font></div><div dir="ltr"><font color="#6fa8dc" size="1"><span style="font-family:helvetica,arial;line-height:12px;text-align:center"><div dir="ltr"><br></div></span></font></div><div dir="ltr"><div style="text-align:center"><font color="#6fa8dc" face="helvetica, arial" size="1"><span style="line-height:12px"><br></span></font></div><div><span style="text-align:center"><font color="#000000" face="helvetica, arial"><span style="font-size:12px;line-height:12px"><a href="https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic" target="_blank">https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic</a></span></font><br></span></div><div><span style="text-align:center"><font color="#000000" face="helvetica, arial"><span style="font-size:12px;line-height:12px"><br></span></font></span></div></div></div></div></div></div></div></div></div></div></div>
</div>