[Gstat-info] Non-Interactive curve fitting and prediction
Mach Nife
machnife at yahoo.com
Fri Feb 24 16:33:40 CET 2006
Hi Edzer,
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
different...
machnife
--- "Edzer J. Pebesma" <e.pebesma at geo.uu.nl> wrote:
> Mach Nife wrote:
>
> >Hi,
> >
> >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) +
> 300
> >Exp(10);
> >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
> it
> >takes this fitted formula to make the predictions?
> It
> >would be nice if it could all be done in one
> operation
> >(Fit the formula and use this formula right away to
> >make the predictions). I've tried adding this at
> the
> >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
> file
> with predictions.
>
> In R this would all be fairly easy.
> --
> Edzer
>
> >machnife
> >
> >__________________________________________________
> >Do You Yahoo!?
> >Tired of spam? Yahoo! Mail has the best spam
> protection around
> >http://mail.yahoo.com
> >_______________________________________________
> >Gstat-info mailing list
> >Gstat-info at geo.uu.nl
>
>http://mailman.geo.uu.nl/mailman/listinfo/gstat-info
> >
> >
>
>
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com
More information about the Gstat-info
mailing list