[Gstat-info] Non-Interactive curve fitting and prediction
machnife at yahoo.com
Fri Feb 24 16:33:40 CET 2006
I would really appreciate if you would have the
kindness to write the equivalent script for R. If I
want to be converted to R this would help a lot.
Also, I'm wondering, if I want to use the Exponential
model, in a totally non-interactive way, what would be
a good initial formula? I noticed that with different
initial formulas, the result of the fitting is
--- "Edzer J. Pebesma" <e.pebesma at geo.uu.nl> wrote:
> Mach Nife wrote:
> >I've used this command file:
> >data(density): 'density.txt', x=2, y=3, v=1, min=4,
> >max=40, radius=1000;
> >variogram(density): 'density.sem': 0.14 Nug(0) +
> >method: semivariogram;
> >set fit=7;
> >And this command:
> >gstat -l cmd_c.log cmd_c.gst
> >The fitted formula is shown in the cmd_c.log file
> as :
> >424437 Nug(0) - 87.4479 Exp(-12.784)
> >Is there a way to modify the above command file, so
> >takes this fitted formula to make the predictions?
> >would be nice if it could all be done in one
> >(Fit the formula and use this formula right away to
> >make the predictions). I've tried adding this at
> >bottom of the command file, without success:
> >predictions(density): 'Dens_Krig';
> >Is it possible?
> No; you'd need to do some scripting to pull the
> fitted model (which is in a syntax gstat can read)
> out of the log file and put it into the next command
> with predictions.
> In R this would all be fairly easy.
> >Do You Yahoo!?
> >Tired of spam? Yahoo! Mail has the best spam
> protection around
> >Gstat-info mailing list
> >Gstat-info at geo.uu.nl
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
More information about the Gstat-info