#
#  MDS Model for Color Circle Data
#
model{

#  Fix 490 at 0.0, 0.0 and 600 (Yellow) 2nd dimension 0.0
#
         x[5,1] <- .000
	 x[5,2] <- .000
         x[10,2] <- 0.0

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])
	     }
#             sumllh[i] <- sum(llh[])
	  }
#
#  Borrowed From Simon Jackman
#
                 llh[N,N] <- 0.0
   sumllh2 <- sum(llh[,])
#
   
  ## priors
   tau ~ dgamma(.1,.1)
   sigma <- sqrt(1/tau)
#
   x[1,1] ~ dnorm(0.0,0.01)
   x[1,2] ~ dnorm(0.0,0.01)
   x[2,1] ~ dnorm(0.0,0.01)
   x[3,1] ~ dnorm(0.0,0.01)
   x[4,1] ~ dnorm(0.0,0.01)
   x[4,2] ~ dnorm(0.0,0.01)
   x[6,1] ~ dnorm(0.0,0.01)
   x[7,1] ~ dnorm(0.0,0.01)
   x[7,2] ~ dnorm(0.0,0.01)
   x[8,1] ~ dnorm(0.0,0.01)
   x[9,1] ~ dnorm(0.0,0.01)
   x[9,2] ~ dnorm(0.0,0.01)
   x[11,1] ~ dnorm(0.0,0.01)
   x[11,2] ~ dnorm(0.0,0.01)
   x[12,1] ~ dnorm(0.0,0.01)
   x[12,2] ~ dnorm(0.0,0.01)
   x[13,1] ~ dnorm(0.0,0.01)
   x[13,2] ~ dnorm(0.0,0.01)
   x[14,2] ~ dnorm(0.0,0.01)
#
#  Kludge to fix rotation -- set 1st and 2nd dimension coordinate of
#          14th point to be positive
#
   x[2,2] ~ dnorm(0.0,0.01)I(0,)
   x[3,2] ~ dnorm(0.0,0.01)I(0,)
   x[6,2] ~ dnorm(0.0,0.01)I(,0)
   x[8,2] ~ dnorm(0.0,0.01)I(,0)
   x[10,1] ~ dnorm(0.0,0.01)I(0,)
   x[14,1] ~ dnorm(0.0,0.01)I(0,)
}