rinter | R Documentation |
Generate one (or several) realisation(s) of the inhibition or contagious process in a region S x T.
rinter(npoints,s.region,t.region,hs="step",gs="min",thetas=0, deltas,ht="step",gt="min",thetat=1,deltat,recent="all",nsim=1, discrete.time=FALSE,replace=FALSE,inhibition=TRUE,...)
npoints |
Number of points to simulate. |
s.region |
Two-column matrix specifying polygonal region containing
all data locations.
If |
t.region |
Vector containing the minimum and maximum values of
the time interval.
If |
hs, ht |
Function which depends on the distance between points
and |
thetas, thetat |
Parameters of |
deltas, deltat |
Spatial and temporal distance of inhibition. |
gs, gt |
Compute the probability of acceptance of a new point from
|
recent |
If “ |
nsim |
Number of simulations to generate. Default is 1. |
discrete.time |
If TRUE, times belong to N, otherwise belong to R^+. |
replace |
Logical. If TRUE allows times repeat. |
inhibition |
Logical. If TRUE, an inhibition process is generated. Otherwise, it is a contagious process. |
... |
Additional parameters if |
A list containing:
xyt |
Matrix (or list of matrices if |
s.region, t.region |
Parameters passed in argument. |
Edith Gabriel <edith.gabriel@inrae.fr>, Peter J Diggle.
plot.stpp
, animation
and stan
for plotting space-time point patterns.
# simple inhibition process inh1 = rinter(npoints=200,thetas=0,deltas=0.05,thetat=0,deltat=0.001, inhibition=TRUE) plot(inh1$xyt,style="elegant") # inhibition process using hs and ht defined by the user hs = function(d,theta,delta,mus=0.1) { res=NULL a=(1-theta)/mus b=theta-a*delta for(i in 1:length(d)) { if (d[i]<=delta) res=c(res,theta) if (d[i]>(delta+mus)) res=c(res,1) if (d[i]>delta & d[i]<=(delta+mus)) res=c(res,a*d[i]+b) } return(res) } ht = function(d,theta,delta,mut=0.3) { res=NULL a=(1-theta)/mut b=theta-a*delta for(i in 1:length(d)) { if (d[i]<=delta) res=c(res,theta) if (d[i]>(delta+mut)) res=c(res,1) if (d[i]>delta & d[i]<=(delta+mut)) res=c(res,a*d[i]+b) } return(res) } d=seq(0,1,length=100) plot(d,hs(d,theta=0.2,delta=0.1,mus=0.1),xlab="",ylab="",type="l", ylim=c(0,1),lwd=2,las=1) lines(d,ht(d,theta=0.1,delta=0.05,mut=0.3),col=2,lwd=2) legend("bottomright",col=1:2,lty=1,lwd=2,legend=c(expression(h[s]), expression(h[t])),bty="n",cex=2) inh2 = rinter(npoints=100, hs=hs, gs="min", thetas=0.2, deltas=0.1, ht=ht, gt="min", thetat=0.1, deltat=0.05, inhibition=TRUE) animation(inh2$xyt, runtime=15, cex=0.8) # simple contagious process for given spatial and temporal regions data(northcumbria) cont1 = rinter(npoints=100, s.region=northcumbria, t.region=c(1,200), thetas=0, deltas=10000, thetat=0, deltat=10, recent=1, inhibition=FALSE) plot(cont1$xyt,pch=19,s.region=cont1$s.region,mark=TRUE,mark.col=4)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.