Am 24.02.2011 23:19, schrieb Edzer Pebesma:

On 02/24/2011 09:40 PM, Anna Gretschel wrote:
Dear List,

The following formula is given for a spherical model in variofit (geoR).

/rho(h) = 1 - 1.5 * (h/phi) + 0.5(h/phi)^3 if h<  phi , 0 otherwise

/Could anyone explain what the term "if h<phi, 0 otherwise" means? Which
part of the equation becomes 0 if h>phi?
rho(h)

And shouldn't there be a * between 0.5 and (h/phi)?
That would have been consistent, but as this is documentatation, I don't
think it is wrong.

Thanks in advance,
Anna

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Thank you! If you don't mind I would like to ask you another question:
I want to use this function to calculate a rmse (root mean square error) for my fitted models.

obs beeing the observed data (in my case variog$v) and pred beeing the predicted data I use the following formula:

rmse <- function(obs, pred) sqrt(mean((obs-pred)^2))

pred I calculated like this:

pre<- 1 - 1.5 * (h/phi) + 0.5(h/phi)^3 if h<  phi , 0 otherwise
with phi beeing the range and h beeing the distances (variog$u)



The output however is not reliable so I think I made a mistake but cannot figure out what is wrong. Do you have an idea what I where the error is in reasoning?

Thank you so much for allways replying to my beginners questions :)
Anna

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to