Análisis geoestadístico uso espacial.exp en WinBugs

Question

Tengo un modelo de regresión logística en GeoBugs para la estimación de los predictores de la prevalencia de una enfermedad. ¿Alguien puede decirme si hay una forma sencilla de determinar los límites inferior y superior de su phi espacial.exp.

Esto es parte de la modelo:

Infected[i] ~ dbin(p[i], tested[i])
logit(p[i]) <- alpha + beta1*covar1[i] + beta2*covar2[i]  + u[i] 
mu[i] <- 0
u[1:N] ~ spatial.exp(mu[], x[], y[], tau, phi, 1)

Esta como me lo han especificado los priores:

phi ~ dunif(0.25, 20)

En el momento en que mi selección ha sido arbitraria. El GeoBugs manual sólo va tan lejos como decir que los límites de la distribución está determinada por los datos.

Lo que me gustaría saber es que hay una dummies versión de la mejor manera para seleccionar los valores.

Answer 1

He trabajado esto a mí. El límite inferior para el phi puede ser estiamted de

-ln(0.5)/(max separating distance between points)

Para encontrar el máximo de la separación de la distancia he utilizado el siguiente código en R. Mis datos están en un archivo plano con la x y la y, coords cambiado el nombre a long y lat respectivamente:

data <- read.csv(file="file.csv", header=T, sep=",")
coords <- data.frame(data$long,data$lat) 
library(sp)
pointDist <- apply(coords, 1, function(eachPoint) spDistsN1(as.matrix(coords), eachPoint, longlat=TRUE))
distances <- as.vector(pointDist)
max(distances)