Thanks Ken!
This was not exactly what I was after - see my reply to Rolf Turner.
Sorry for not explaining this well enough.
Cheers,
Torbjørn
On Wed, 14 Sep 2011 16:37:48 -0400, Ken Hutchison
<vicvoncas...@gmail.com> wrote:
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
______________________________________________
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.
