help: krigged map do not show colours

Hi,

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.

Have a look at

summary(map)

and particularly the summary of var1.pred.
You could also plot the variogram of the data together with your model,
vmf, and see how they fit.

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).

Best wishes,
Jon




On 8/3/2015 5:30 PM, Anand Sookun wrote:

Hi

Thank you Roman for your kind response.

The data is as attached which is part of the whole dataset.

My other problem is that im not able to create an appropriate prediction
grid with some data included.

The codes which failed to show colour ramps are

map <- krige(logcarbseq ~ 1, locations = carbon_seq, newdata = map,model =
vmf)
str(map)
print(spplot(map, "var1.pred", asp=1, col.regions=bpy.colors(64),main="OK
prediction, carbon"))
print(spplot(map, "var1.var",col.regions=cm.colors(64),asp=1, main="OK
prediction variance, carbon"))

I look forward to hopefully get some solutions.

Thanks again and warm regards

Anand

On Mon, Aug 3, 2015 at 4:09 PM, Roman Lu?trik <roman.lustrik at gmail.com>
wrote:

What is the variability of the plotted variable?

Cheers,
Roman

On Sun, Aug 2, 2015 at 11:03 AM, Anand Sookun < <asookun at gmail.com>
asookun at gmail.com> wrote:

Hi

I have kriged my data but the map shows uniform colours - no gradual
contour colour plot on map

code: print(spplot(k40, "var1.pred", asp=1,
col.regions=bpy.colors(64),main="OK prediction, log-ppm Zn"))

thnaks for any help

Anand Sookun
--
-------------------------------------------------------
*Anand Sookun (**MPhil. **PhD Candidate, University of Mauritius)*
*Director cum Freelance Consultant*
*Soft Sustainable Services (Mauritius)*
*Research and Consultancy*

Website: http://sustain-ability.webs.com/



https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic

        [[alternative HTML version deleted]]

--
In God we trust, all others bring data.

--
-------------------------------------------------------
*Anand Sookun (**MPhil. **PhD Candidate, University of Mauritius)*
*Director cum Freelance Consultant*
*Soft Sustainable Services (Mauritius)*
*Research and Consultancy*

Website:  <http://sustain-ability.webs.com/>
http://sustain-ability.webs.com/



<https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic>
https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic