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
<mailto: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
<mailto: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]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org <mailto:R-sig-Geo at r-project.org>
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
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/
https://www.linkedin.com/profile/view?id=161468467&trk=nav_responsive_tab_profile_pic
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
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 jon.skoien at jrc.ec.europa.eu Tel: +39 0332 789205 Disclaimer: Views expressed in this email are those of the individual and do not necessarily represent official views of the European Commission. [[alternative HTML version deleted]]