[R-sig-Geo] geoR, Kriging, non-standard models

[Gildas Mazo]

Dear P.J. Ribeiro,

thank you for your answer. Actually I have realized I have misswritten my question. I indeed fit the two power models using variofit(), but I don't know how to perform the kriging map based on the combination (4.1.2) ??

Thank you,
Gildas

----- Mail original -----
> De: "Paulo Justiniano Ribeiro Jr" <paulojus at c3sl.ufpr.br>
> À: "Gildas Mazo" <gildas.mazo at inria.fr>
> Cc: r-sig-geo at r-project.org
> Envoyé: Jeudi 8 Décembre 2011 11:59:53
> Objet: Re: [R-sig-Geo] geoR, Kriging, non-standard models
> Gildas
> 
> regarding with geoR you can only "fit" such models visually using
> lines.variomodel() or alike.
> The automatic fitting functions (variofit(), likfit())
> do not fit "nested" models.
> 
> 
> 
> 
> Paulo Justiniano Ribeiro Jr
> LEG (Laboratorio de Estatistica e Geoinformacao)
> Universidade Federal do Parana
> Caixa Postal 19.081
> CEP 81.531-990
> Curitiba, PR - Brasil
> Tel: (+55) 41 3361 3573
> VOIP: (+55) (41) (3361 3600) 1053 1066
> Fax: (+55) 41 3361 3141
> e-mail: paulojus AT ufpr br
> http://www.leg.ufpr.br/~paulojus
> 
> On Tue, 6 Dec 2011, Gildas Mazo wrote:
> 
> > Dear R users,
> >
> > I am learning geostatistics and was trying to analyse the
> > Wolfcamp-Aquifer data from Cressie's famous book (1994) using the
> > geoR
> > package. However I don't know how to fit the semivariogram model
> > found
> > by Cressie (4.1.2) page 217. It is a smooth combination of two power
> > models. How can I do that with geoR ? Or gstat ?
> >
> > Thanks,
> > Gildas
> >
> > --
> >
> > Gildas Mazo
> > Doctorant Equipe-Projet MISTIS
> > INRIA Grenoble Rhone-Alpes
> > Inovall?e, 655, av. de l'Europe, Montbonnot
> > 38334 Saint-Ismier Cedex France
> > Bureau: G-206
> > Tel: +33 (0) 4 56527104
> >
> > 	[[alternative HTML version deleted]]
> >
> >

-- 
Gildas Mazo
Doctorant Equipe-Projet MISTIS
INRIA Grenoble Rhone-Alpes
Inovallée, 655, av. de l'Europe, Montbonnot
38334 Saint-Ismier Cedex France
Bureau: G-206
Tel: +33 (0) 4 56527104