remove(list=ls()) library(BRugs) #Library needed #Example of simulated squared distances d2=c(24555, 22512, 100685, 11678, 49589, 29577, 3533, 46040, 10313, 121275, 23488, 29810, 10635, 58222, 55974, 3333, 6723, 20479, 29257, 29702, 38641, 22270, 7018, 5947, 58308, 4810, 25044, 21188, 22501, 187, 12149, 67319) #Example of simulated times-at-liberty t=c(838, 618, 798, 135, 481, 349, 646, 61, 244, 498, 166, 373, 298, 689, 491, 143, 982, 720, 942, 69, 968, 243, 743, 759, 189, 364, 717, 571, 945, 62, 449, 876) #INITIAL VALUES initials<- function(){list( a=runif(1,19000,21000), k=runif(1,0,1))} #(a is epsilon/k) #DATA data<-list( d2=d2, t=t, n=length(t)) #MODEL model <- function(){ for( i in 1 : n) { mu[i] <- 2*a*(1 - exp(-k*t[i])) b[i]<-1/mu[i] d2[i] ~ dexp(b[i])} a~dunif(0,1000000000) k~dunif(0,10000) epsilon<-a*k radius<-sqrt((log(1-0.95))/-(1/a))#in meters, distance embracing 95% positions } #RUNNING THE MODEL MCMC<-BRugsFit( modelFile=model, data=data, inits=initials, numChains = 5, parametersToSave=c("a","radius","k","epsilon"), nBurnin = 10000, nIter = 30000, nThin = 10, DIC = FALSE, working.directory=NULL, digits = 5, BRugsVerbose=FALSE) samplesHistory("*",beg=0,thin=0,mfrow = c(3, 1)) samplesDensity("radius",beg=0) samplesStats("radius",beg=0)