model {

for (i in 1:N) {
	Y[i] ~ dbern(p[i])
	logit(p[i]) <- w[i]
	mu[i] <- beta[1]+beta[2]*LivingArea[i]/1000+beta[3]*Age[i]
	useless[i] <- HalfBaths[i]+x[i]+y[i]+Age[i]+OtherArea[i]+Beds[i]+Baths[i]+HalfBaths[i]
}

for (i in 1:3) {beta[i] ~ dnorm(0.0,0.001)}
w[1:N] ~ spatial.exp(mu[], x[], y[], spat.prec, phi, 1) 	
#for (i in 1:N) {w[i] ~ dnorm(0.0,spat.prec)}
phi ~ dunif(0.1,10)	
	
spat.prec ~ dgamma(0.1, 0.1)
sigmasq <- 1/spat.prec
	
}

#Initial values
list(beta=c(0,0,0), phi=5.0, spat.prec=1.0,w = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0))

list(beta=c(0,0,0,0))
list(beta=c(5,5,5,5))

#Data values
list(N = 50, LivingArea = c(1606, 
979, 1144, 2154, 2992, 1174, 1600, 1020, 1864, 1526, 2408, 1545, 
2205, 2583, 1705, 2998, 1900, 2200, 1789, 1787, 2121, 1535, 1010, 
1546, 1413, 2014, 1358, 3418, 1836, 1200, 3052, 1963, 1709, 2477, 
2642, 2147, 1500, 1700, 1885, 1918, 2066, 2064, 1762, 1328, 1376, 
2830, 1485, 1593, 1624, 1533), OtherArea = c(794, 452, 584, 786, 
1522, 410, 550, 529, 1114, 706, 613, 599, 635, 634, 572, 754, 
850, 768, 596, 275, 256, 478, 386, 791, 660, 695, 567, 1016, 
556, 1042, 1161, 819, 570, 994, 1008, 512, 410, 832, 920, 559, 
521, 750, 865, 646, 655, 934, 694, 698, 700, 512), Age = c(11, 
2, 0, 13, 3, 0, 7, 13, 19, 10, 0, 6, 12, 12, 15, 23, 10, 14, 
20, 42, 17, 14, 14, 14, 10, 25, 5, 10, 13, 20, 5, 8, 20, 11, 
11, 24, 24, 50, 12, 15, 18, 5, 16, 15, 17, 11, 0, 4, 8, 10), 
Beds = c(4, 2, 3, 4, 5, 3, 3, 3, 3, 3, 4, 3, 4, 3, 3, 6, 3, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 2, 4, 3, 3, 3, 4, 3, 3, 3, 3, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3), Baths = c(2, 1, 2, 2, 3, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2), HalfBaths = c(0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0), Y = c(0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0), x = c(3.95789322711243, 
-1.65953903186291, -4.65039033698281, 3.40190163834073, 5.50124918904885, 
-4.38198060447266, 4.67684786776564, 4.28381933087544, -4.76542307948716, 
3.29645495771072, 7.3226009453725, 3.21018040083212, -9.5105237078035, 
-9.84603587344221, -0.653002534948152, 0.65070188010523, -0.892654081832554, 
1.08207466449688, -1.06520319558976, -8.27392172587868, 3.96747928898814, 
2.17488571829025, 4.29340539275114, 1.96399235703295, 1.04373041699543, 
-1.50616204185711, 4.54264300151125, 5.0698764046572, 1.88730386202869, 
-8.9353599952804, -8.08220048837143, -6.72097970206655, -1.19940806184415, 
5.79841710718611, 3.01845916332487, -2.05256756875448, -7.92882349836427, 
-6.12664386579204, -0.0778388224250416, 2.27074633704455, 2.21322996579170, 
4.44678238275694, -0.662588596822492, 5.41497463217025, 0.852009179488183, 
-6.5771887739351, 1.80102930515145, -1.40071536122711, 4.67684786776564, 
4.14961446461968), y = c(-6.86199459609406, 4.26750561169679, 
7.05961021972557, 0.116599000078903, 0.205554941393838, -7.86719673295503, 
-1.02870874435364, 0.789328306274477, 10.5166604890855, -1.61692990629987, 
-5.77784406131572, -0.242560612980816, -5.38087817319692, -5.04284559619954, 
3.90612210010388, 0.312302070972312, 12.1801365916785, -4.50243825270997, 
9.44918919330357, -3.5973115498283, 0.363451737228557, 12.3758396625719, 
0.137726036141240, 5.51511268864143, 4.1963408586446, -1.13323197539893, 
0.471310816073034, 0.723723299554674, -4.24335407362955, 5.05031789527003, 
-6.29601241947656, -6.08029426178721, -0.719586848283376, 0.204442992127441, 
-2.83562630231740, 3.71486682627645, 5.18486375650922, -0.643974298165603, 
4.00063778775119, -0.142485179001515, 1.16850300613053, 0.436840388813536, 
-4.87494125696685, -7.23116175255174, 2.20372777318541, -7.14109386197041, 
-1.27778538003609, -6.93649519694543, -1.00869365755770, -2.06726935920771
))