#
#  MDS Model for Nations Data
#
model{
#
# Fix USSR
#
        x[10,1] <- 0.000
        x[10,2] <- 0.000
#
#  Fix USA 2nd Dimension
#
         x[11,2] <- 0.000
#
for (i in 1:N-1){
          llh[i,i] <- 0.0
	  for (j in i+1:N){
#
#
dstar[i,j] ~ dlnorm(mu[i,j],tau)
                 mu[i,j] <- log(sqrt((x[i,1]-x[j,1])*(x[i,1]-x[j,1])+(x[i,2]-x[j,2])*(x[i,2]-x[j,2])))
                 llh[i,j] <- (log(dstar[i,j])-mu[i,j])*(log(dstar[i,j])-mu[i,j])
                 llh[j,i] <- (log(dstar[i,j])-mu[i,j])*(log(dstar[i,j])-mu[i,j])
	     }
	  }
#
# Borrowed from Simon Jackman
#
llh[N,N] <- 0.0
sumllh2 <- sum(llh[,])
#   
## priors
#
tau ~ dgamma(.1,.1)
#
  for(l in 1:2){x[1,l] ~ dnorm(0.0, 0.01)}
  x[2,1] ~ dnorm(0.0,0.01)
  for(l in 1:2){
      for(k in 3:6) {x[k,l] ~ dnorm(0.0, 0.01)}
  }
  x[7,2] ~ dnorm(0.0,0.01)
  x[8,1] ~ dnorm(0.0,0.01)
  x[9,2] ~ dnorm(0.0,0.01)
  for(l in 1:2){x[12,l] ~ dnorm(0.0, 0.01)}
#
#  Kludge to fix rotation -- set 1st and 2nd dimension coordinates of
#          4 Nations to fix the sign flips

   x[2,2] ~ dnorm(0.0,0.01)I(0,)
   x[7,1] ~ dnorm(0.0,0.01)I(0,)
   x[8,2] ~ dnorm(0.0,0.01)I(,0)
   x[9,1] ~ dnorm(0.0,0.01)I(,0)
   x[11,1] ~ dnorm(0.0,0.01)I(0,)
}