[R] Weibull point process

[Torbjørn Ergon]

 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 
 <vicvoncastle at 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