<html><head></head><body><div style="color:#000; background-color:#fff; font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px"><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span>Dear </span><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;" id="yui_3_16_0_ym19_1_1495120744730_7201">Dr. Benedikt Gräler,</span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;"><br></span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;" id="yui_3_16_0_ym19_1_1495120744730_7202">Thank you for your advice. Now, I have found some directs to move on. </span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;"><br></span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;" id="yui_3_16_0_ym19_1_1495120744730_7353">Actually, I do not have enough background to analyze by myself. Based on your suggestions, I have plotted the variogram. I think there are daily and weekly seasonality.</span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;"><br></span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;" id="yui_3_16_0_ym19_1_1495120744730_7380">According to your suggestions, I am trying to model daily seasonality. It would be great if you can give me more advice.</span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;"><br></span></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif">With Best Regards,</font></div><div id="yui_3_16_0_ym19_1_1495120744730_6868" dir="ltr"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif">Ankur</font></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 13px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> On Wednesday, May 17, 2017 6:10 AM, Dr. Benedikt Gräler <b.graeler@52north.org> wrote:<br></font></div> <br><br> <div class="y_msg_container"><div dir="ltr">Dear Ankur,<br clear="none"><br clear="none">I see two issues from your plot:<br clear="none">i) there seems to be a seasonality in your data (maybe daily?) as the <br clear="none">sample variogram grows up to half the length of your temporal axis and <br clear="none">then drops again<br clear="none">ii) the starting values of your fitting routine are far off from what <br clear="none">can be deduced from the sample variogram. Hence, the numerical routines <br clear="none">have a hard time to identify the correct values. Additionally, the <br clear="none">parameters are of different orders of magnitude further complicating the <br clear="none">numerical fit. Manual re-scaling might help.<br clear="none"><br clear="none">Try adding the parameter "scales=list(arrows=F)" for a visual assessment <br clear="none">of the empirical variogram; in your case<br clear="none"><br clear="none">plot(var, wireframe=T, scales=list(arrows=F))<br clear="none"><br clear="none">This reveals that the variogram indicates the same strength of <br clear="none">dependence for values 24 and 0 hours apart -> daily cycle? I'd suggest <br clear="none">to model the daily cycle first and then to model the spatio-temporal <br clear="none">variogram of the residuals.<br clear="none"><br clear="none">HTH,<br clear="none"><br clear="none"> Ben<br clear="none"><br clear="none"><br clear="none"><div class="yqt2363048959" id="yqtfd87299"><br clear="none">On 08/05/2017 06:33, Ankur Sarker via R-sig-Geo wrote:<br clear="none">> Hi,<br clear="none">> <br clear="none">> Can anyone suggest me the most appropriate variogram model to fit my data?<br clear="none">> <br clear="none">> I have tried three different models and results are too bad. Here is the <br clear="none">> comparison of different variogram models as attached.<br clear="none">> <br clear="none">> Thanks,<br clear="none">> Ankur<br clear="none">> <br clear="none">> <br clear="none">> _______________________________________________<br clear="none">> R-sig-Geo mailing list<br clear="none">> <a shape="rect" ymailto="mailto:R-sig-Geo@r-project.org" href="mailto:R-sig-Geo@r-project.org">R-sig-Geo@r-project.org</a><br clear="none">> <a shape="rect" href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a></div><br clear="none">> <br clear="none"><br clear="none">-- <br clear="none">Dr. Benedikt Gräler<br clear="none">52°North Initiative for Geospatial Open Source Software GmbH<br clear="none">Martin-Luther-King-Weg 24<br clear="none">48155 Muenster, Germany<br clear="none"><br clear="none">E-Mail: <a shape="rect" ymailto="mailto:b.graeler@52north.org" href="mailto:b.graeler@52north.org">b.graeler@52north.org</a><br clear="none">Fon: +49-(0)-251/396371-39<br clear="none">Fax: +49-(0)-251/396371-11<br clear="none"><br clear="none"><a shape="rect" href="http://52north.org/" target="_blank">http://52north.org/</a><br clear="none">Twitter: @FiveTwoN<br clear="none"><br clear="none">General Managers: Dr. Albert Remke, Dr. Andreas Wytzisk<br clear="none">Local Court Muenster HRB 10849<div class="yqt2363048959" id="yqtfd14767"><br clear="none"></div></div><div class="yqt2363048959" id="yqtfd92488">_______________________________________________<br clear="none">R-sig-Geo mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-sig-Geo@r-project.org" href="mailto:R-sig-Geo@r-project.org">R-sig-Geo@r-project.org</a><br clear="none"><a shape="rect" href="https://stat.ethz.ch/mailman/listinfo/r-sig-geo" target="_blank">https://stat.ethz.ch/mailman/listinfo/r-sig-geo</a></div><br><br></div> </div> </div> </div></div></body></html>