profile likelihood

Ben Bolker · 2011-12-09T21:39:47Z

plocq hotmail.com> writes: > > Hi, > > I try to use the function profile() of the SpatialExtremes' package to > obtain the profile likelihood of parameters for an extreme values fit based > on Poisson process : > > fit > But when I call "profile(fit)", I obtain the following error (even if I > precise others arguments of the function) : > > [1] "profiling loc" > Erreur dans nlpot(p[1], p[2], ...) : > argument(s) inutilis?(s)

Ben Bolker

Fri, Dec 9, 2011 1:39 PM

plocq <less.is.mo_ <at> hotmail.com> writes:

We need a reproducible example, please.

  fpot() is in the evd package, not in the SpatialExtremes package
(you didn't say explicitly that it was, but this would be helpful
information).

  This minimal example works:

library(evd)
uvdata <- rgpd(100, loc = 0, scale = 1.1, shape = 0.2)
M1 <- fpot(uvdata, 1)
library(SpatialExtremes)
profile(M1)

  So you need to provide us with a minimal example that doesn't work.

profile likelihood

Thread (3 messages)