I don't know if this settles the matter, but if you are modeling Weibull as
the Log-Intensity I.E. in a non-homegenous poisson process you can:
#Where T is the length of the series desired
Weibull.process<-function(T,shape,scale)
{
logLambda=rweibull(T,shape,scale)
Lambda=exp(logLambda)
Point.Process.Counts=rpois(T,Lambda)
return(Point.Process.Counts)
}
I haven't actually tried this code (may contain clerical errors) but I hope
it gets you on the right track.
Good luck,
Ken Hutchison
2011/9/14 Torbjørn Ergon <t.h.er...@bio.uio.no>
> Dear list,
>
> I'm looking for a function to generate (simulate) a random Weibull point
> process. Can anyone help?
>
> Cheers,
>
> Torbjørn Ergon, University of Oslo
>
> ______________________________**________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
> PLEASE do read the posting guide http://www.R-project.org/**
> posting-guide.html <http://www.R-project.org/posting-guide.html>
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
