# WinBUGS 1.4 code for joint model (Model XI in Guo and Carlin paper) # Written by Xu Guo, 3/1/03 model{ for (i in 1:N) { for (j in 1:M) { Y[i, j] ~ dnorm(muy[i, j], tauz) muy[i, j]<-beta1[1]+beta1[2]*t[j]+beta1[3]*t[j]*randgrp1[i]+ beta1[4]*gender1[i]+beta1[5]*prevoi1[i]+beta1[6]*stratum1[i]+U[i,1]+U[i,2]*t[j] } surt[i] ~ dweib(p,mut[i]) I(surt.cen[i],) log(mut[i])<-beta2[1]+beta2[2]*randgrp1[i]+beta2[3]*gender1[i]+ beta2[4]*prevoi1[i]+beta2[5]*stratum1[i]+r1*U[i, 1]+r2*U[i, 2] U[i,1:2] ~ dmnorm(U0[],tau[,]) } p <- 1 # Use this for Exponential model #p ~ dgamma(1,1) # Use this for full Weibull model sigmaz<-1/tauz sigma[1:2,1:2]<-inverse(tau[,]) sigma1<-sigma[1,1] sigma2<-sigma[2,2] sigma12<-sigma[1,2] cor<-sigma12/(sqrt(sigma1*sigma2)) #priors tau[1:2,1:2] ~ dwish(R[,], 23) beta1[1:6]~dmnorm(betamu1[],Sigma1[,]) tauz~dgamma(0.1, 0.1) beta2[1:5]~dmnorm(betamu2[],Sigma2[,]) r1~dnorm(0, 0.01) r2~dnorm(0, 0.01) } #data list(N=467, M=5, betamu1=c(0,0,0,0,0,0), betamu2=c(0,0,0,0,0), Sigma1=structure(.Data=c(0.01,0,0,0,0,0, 0,0.01,0,0,0,0, 0,0,0.01,0,0,0, 0,0,0,0.01,0,0, 0,0,0,0,0.01,0, 0,0,0,0,0,0.01),.Dim=c(6,6)), Sigma2=structure(.Data=c(0.01,0,0,0,0, 0,0.01,0,0,0, 0,0,0.01,0,0, 0,0,0,0.01,0, 0,0,0,0,0.01),.Dim=c(5,5)), U0=c(0,0), R=structure(.Data=c(100, 0, 0, 100), .Dim=c(2,2)), t=c(0,2,6,12,18), Y=structure( .Data=c(10.67707825, NA, 8.426149773, 9.433981132, NA, 6.32455532, NA, 8.124038405, 4.582575695, 5, 3.464101615, 3.605551276, 6.164414003, NA ,NA, 3.872983346, 4.582575695, 2.645751311, 1.732050808, NA, 7.280109889, 8.602325267, 8.602325267, 6.708203933, NA, 4.582575695, NA, NA, NA, NA, 6.782329983, 5.385164807, 4.472135955, 3.16227766, NA, 3.464101615, 1 ,0 ,NA, NA, 0, 5.744562647, NA, NA, NA, 12.12435565, 13.41640787, 10.53565375, 7.483314774, NA, 6.8556546, NA ,NA ,NA ,NA, 13.60147051 ,16.18641406 ,13.3041347, 9.539392014 ,11.74734012, 2.449489743 ,4, 2.449489743,NA ,NA, 2 ,2.828427125, NA ,NA, NA, 11.74734012, 9.848857802, 9.110433579 ,8.888194417 ,NA, 1.732050808, 1, 5.099019514, NA, NA, 6.557438524, 4.898979486, 5, 5 ,NA, 5.65685425, 5.385164807 ,3, 2.645751311 ,NA, 3, NA, NA, NA, NA, 19.05255888, 19.07878403 ,16.43167673, 16.70329309, NA, 7.211102551, 5.385164807, 3.31662479, NA, NA, 1.732050808, 1.732050808, 1.414213562, 2, NA, 3.872983346, NA, NA, NA, NA, 14.03566885, 15.96871942, NA, 15.32970972, 14.35270009, 2.236067978, NA, NA, NA, NA, 17.17556404, 19.84943324, 20.27313493, 17.05872211, 13.60147051, 4.123105626, NA, 2.236067978, 1.414213562, 1.732050808, 16.88194302, NA, NA, NA, NA, 4.582575695, NA, NA, NA, NA, 5.744562647, NA, 4.358898944, NA, NA, 9.110433579, 10.58300524, 5.567764363, NA, NA, 1.732050808, NA ,1, NA, NA, 7.745966692, 6.557438524, NA, NA, 3.464101615, 1 ,NA ,1.414213562 ,NA ,2.645751311, 12.52996409, 10.04987562 ,7.141428429, 6.8556546, 6.32455532, 0, NA, 1, NA, NA, 3.872983346, 3, 2.236067978, NA, NA, 16.73320053, 15.26433752, NA, NA, NA, 10, NA, 7.810249676, 4.472135955, NA, 16.21727474, 14.28285686, NA, NA, NA, 4.69041576, 3.31662479, NA, 1.414213562, NA, 8.306623863, NA, NA, NA, NA, 16.4924225, NA, NA, NA, NA, 4.582575695, 5.385164807, 3.464101615, NA, NA, 1.732050808, NA, NA, NA, NA, 6.244997998, 4.898979486, 2.645751311, NA, NA, 4, 3.605551276, NA, NA, NA, 1, 1.414213562, NA, NA, NA, 1, NA, NA, NA, NA, 13.78404875, 14.38749457, 12.12435565, 4, NA, 3, 4.472135955, NA, NA, NA, 1.414213562, 0, NA, NA, NA, 3.605551276, NA, NA, NA, NA, 10.53565375, 13.67479433, 11.40175425, 9, NA, 2.449489743, 3.16227766, 4.795831523, 2, NA, 6.92820323, 7.280109889, 3.605551276, 4.472135955, NA, 3.31662479, 5.65685425, 3.16227766, NA, NA, 2.449489743, 1.732050808, 1.732050808, NA, NA, 11.40175425, 5.196152423, 5.744562647, NA, NA, 2.828427125, NA, 2.236067978, 1, NA, 12.40967365, 13.3041347, NA, NA, NA, 16.46207763, 15.8113883, 11.18033989, 13, NA, 13.92838828, 17.52141547, 11.91637529, 9.949874371, NA, 5.65685425, 5.567764363, 2.449489743, NA, NA, 3.872983346, NA, NA, NA, NA, 7.071067812, NA, 4.582575695, NA, NA, 4.123105626, NA, 2.449489743, NA, NA, 13.85640646, 8.831760866, 5.744562647, 3.31662479, NA, 2.236067978, NA, NA, NA, NA, 6.782329983, 7.810249676, 8.888194417, NA, NA, 2, 2.236067978, NA, NA ,NA, 15.03329638 ,13.67479433, 13.78404875, 13.67479433 ,NA, 1, NA, NA ,NA ,NA, 4.582575695, 6.08276253, 2, 2.828427125, NA, 10.48808848, 8.124038405, 7.615773106, NA, NA, 4.358898944, 4.358898944, 4, 3.16227766, NA, 3.741657387, 1, 1.414213562, NA, NA, 7.211102551, 11.70469991, 10.19803903, 11.40175425, NA, 4.582575695, 5.567764363, 5.65685425, 2.828427125, NA, 10.95445115, 17.32050808, 14.83239697, 16.4924225, NA, 2.236067978, 2.236067978, NA, NA, NA, 6.480740698, 5.830951895, 3.872983346, NA, 2.645751311, 7.071067812, 7.071067812, NA, NA, NA, 8, 3.872983346, NA, NA, NA, 3.872983346, 7.348469228, 6.708203933, 5, NA, 15.55634919, 19.74841766, 16.24807681, 18.22086716, NA, 6.244997998, 1.414213562, 4.898979486, NA, NA, 2.828427125, 5.477225575, 1.732050808, NA, NA, 2, 3.31662479, 1.414213562, 4.582575695, 4.358898944, 14.10673598, 14.24780685, 16.91153453, 15.96871942, 13.52774926, 8.062257748, 7.141428429, 6.782329983, 2.236067978, 6.08276253, 6.708203933, 6.92820323, 7.745966692, 7.810249676, 6.480740698, 17.05872211, 17.72004515, 21.28379665, 19.54482029, 16.52271164, 6.244997998, 5.477225575, 4.582575695, 4, NA, 10.7238053, 9.486832981, 11.78982612, 11.87434209, NA, 14.96662955, NA ,19.05255888, 16.18641406, NA, 6.32455532, 4, NA, NA, NA, 6.32455532, 5.099019514, 6, NA, NA, 4.898979486, NA, 1.732050808, NA ,NA, 4.582575695, 5.916079783, 5.291502622, NA, NA, 13.26649916, 14.07124728, 9.899494937, 4.69041576, NA, 2.449489743, 1.414213562, NA, NA, NA, 3.605551276, 12.72792206, NA, 6.480740698, NA, 3.872983346, 4.472135955, 4.472135955, NA, NA, 3.31662479 ,2.236067978 ,NA ,NA, NA, 16.58312395, 20.37154879, 16.85229955, 17.54992878, NA, 7.071067812, 8.485281374, 6, NA, NA, 1.414213562, NA, NA, NA, NA, 7.810249676, NA, NA, NA, NA, 16.55294536, 16.2788206, NA, NA, NA, 5.291502622, 6.480740698, 4 ,NA, NA, 6.403124237, 5.477225575, NA ,NA ,NA, 7.348469228, NA ,NA ,NA ,NA, 2 ,2, 1.414213562, NA, NA, 7.348469228, 4.582575695, NA ,NA ,NA, 5.385164807, 4.582575695, NA ,NA ,NA, 10, NA, 3.741657387 ,5.385164807, NA, 3.16227766, NA, NA ,NA ,NA, 12.489996, NA, 6.244997998 ,NA ,NA, 0 ,NA, NA, NA, NA, 7.681145748, 6.557438524 ,7.810249676, 5 ,NA, 2.449489743, 2 ,NA, NA, NA, 8.426149773, NA, NA, 2.449489743 ,NA, 5.196152423, 2.828427125 ,3.605551276, 4.242640687, NA, 3.872983346, 1.732050808, 1.732050808, 2, NA, 7.937253933, 5.477225575, 4.69041576, 2, NA, 5, 12.76714534, 10.63014581, 7.681145748, NA, 3.741657387, NA, NA, NA, NA, 4.582575695, 3.872983346, 2.236067978, NA ,NA, 14.03566885 ,12.40967365 ,11.61895004 ,10.48808848 ,5.567764363, 11.40175425, 24.12467616, NA,NA ,NA, 4 ,6.08276253, 3.16227766, 2.828427125 ,NA, 6.08276253 ,5.099019514 ,0, 2.828427125, NA, 3.16227766 ,3.464101615 ,3, NA, NA, 7.071067812, 6.32455532, 4.472135955 ,3.16227766, NA, 10.44030651, 6.480740698, 4.242640687, NA ,NA, 9.433981132, 4.358898944, 6.557438524, 3 ,2.449489743, 7.211102551, 3.872983346, 5.099019514, NA, NA, 6.480740698 ,4.795831523, 5.477225575, 3.31662479 ,3.872983346, 2.645751311, NA, 1, NA ,NA, 3.16227766, 0, 2.236067978, NA, NA, 15.23154621 ,14.86606875, 12.88409873, 11.87434209, NA, 10.81665383 ,4.123105626, NA, 7.416198487, NA, 16.97056275 ,16.37070554, 12.24744871, 17.2626765, NA, 12.9614814, 15.03329638 ,11.35781669, 12,NA, 13 ,12.20655562, 9.643650761, 6.708203933 ,NA, 3.464101615 ,4 ,0, NA ,NA, 11.22497216 ,15.68438714, 12.36931688 ,11.83215957 ,NA, 11.66190379 ,10.58300524, NA ,9.273618496, NA, 1 ,1.732050808 ,0 ,0 ,NA, 9.539392014, 6.08276253, NA ,NA, NA, 16.73320053 ,NA, 17.8605711 ,NA, NA, 12.64911064 ,10.7238053 ,7.549834435 ,8.306623863 ,NA, 6.782329983 ,5.567764363, 2.828427125, 3.16227766, NA, 2.236067978 ,1.732050808 ,1.732050808, 2.645751311, NA, 0 ,NA, NA ,NA ,NA, 15.68438714 ,18.24828759, NA, 15.23154621 ,13.3041347, 14.96662955 ,NA ,14.28285686, 11.87434209 ,NA, 15.49193339 ,18.97366596, NA ,17.60681686 ,NA, 14.07124728 ,13.11487705, NA ,16.91153453, NA, 8.124038405 ,9.433981132, 9 ,4.242640687 ,NA, 6.32455532, 3.872983346 ,5.830951895, 2.236067978, NA, 3.16227766 ,5.477225575, NA ,NA, NA, 5.477225575 ,NA, NA, NA ,NA, 10.95445115 ,12.64911064,10.95445115 ,11.40175425, 8.94427191, 7.071067812 ,NA ,4.472135955, 3.16227766, NA, 6.8556546 ,10.48808848, 7.937253933, 4.472135955 ,NA, 4.472135955, 3.16227766, 4.472135955 ,NA ,NA, 5.477225575 ,6.32455532 ,9.746794345, 7.071067812 ,NA, 8.94427191 ,7.071067812, 7.745966692, 3.16227766, NA, 5.477225575, 2, 6.32455532, 3.16227766, NA, 8.717797887, 10.48808848, NA, NA, NA, 15.8113883, NA, 14.49137675, 14.83239697, NA, 5.477225575, 6.32455532, 6.32455532, 17.88854382, NA, 10, 6.32455532, NA ,NA, NA, 11.83215957, 10.24695077, NA ,9.433981132 ,NA, 4.472135955 ,4.472135955, 3.16227766, 4.472135955, NA, 16.43167673, 17.02938637, 14.83239697, 15.49193339 ,NA, 6.32455532, 4.472135955 ,5.477225575 ,NA ,NA, 13.78404875 ,9.486832981, 11.40175425, 8.94427191, NA, 6.32455532, NA ,NA ,NA, NA, 6.32455532, 8.94427191, 4.472135955, 4.472135955, NA, 8.366600265, 6.32455532, 5.477225575, 0, NA, 13.41640787, 16.1245155 ,12.64911064, 13.41640787 ,NA, 0 ,NA, 7.745966692, 8.366600265, NA, 13.78404875, 16.73320053, 11.40175425, 14.14213562, NA, 16.43167673, 19.23538406, 18.97366596, 18.70828693, NA, 2.645751311, 2, 2.449489743, 16.94107435, NA, 19.23538406, NA, NA ,NA, NA, 4.898979486, 1.732050808, 3.31662479, 2.645751311 ,NA, 12 ,15.39480432, 15.55634919, 10.90871212, NA, 3.605551276, 5.196152423, NA ,4.358898944, NA, 10.19803903, 13.19090596, 13.92838828 ,NA ,NA, 5 ,0, 3 ,4.795831523 ,NA, 10.48808848 ,15.8113883 ,13.03840481, 8.94427191, NA, 4.472135955, NA, NA ,NA, NA, 11.40175425, NA, NA ,NA, NA, 12.64911064 ,15.16575089 ,17.32050808, 14.14213562, NA, 13.78404875, 17.32050808, 12.24744871, 15.49193339 ,NA, 4.472135955 ,3.16227766, 0 ,2 ,NA, 2.236067978, 0, NA, NA, NA, 7.681145748, NA, NA, NA, NA, 2.236067978, 3.741657387, 2.236067978, NA, NA, 5.291502622, NA, 6.782329983, NA, NA, 2.449489743, 2.236067978, NA, NA, NA, 10.34408043, 11.22497216, 8.485281374, 3.741657387, NA, 4.795831523, NA, NA, NA, NA, 2, 2.236067978, 2, 3 ,NA, 5.65685425, 6.244997998, 6.782329983 ,5.477225575, 5.830951895, 11.95826074, 12.40967365, 16.24807681, 13.26649916, NA, 3.741657387, 2.828427125, NA ,NA, NA, 16.52271164 ,19.57038579 ,16.21727474 ,19.20937271, NA, 8.485281374 ,4.472135955, 3 ,NA, NA, 8.774964387, 4.358898944, NA, NA, NA, 9.746794345 ,15.29705854, 14.69693846, 9.899494937, NA, 9.695359715, 7.280109889 ,5.385164807, 2.236067978, NA, 12.84523258, NA, 5.916079783, NA, NA, 2.449489743, 4.242640687, NA, NA, NA, 16.76305461, 14.38749457, 12.64911064 ,8.426149773,7.874007874, 4.123105626, 3.16227766 ,NA ,NA, NA, 2.645751311 ,3.872983346 ,3.31662479, NA, NA, 10.44030651, 9.591663047 ,9.219544457, NA ,8.246211251, 14.89966443 ,11.66190379 ,NA, 13.37908816, NA, 2.449489743 ,3.16227766 ,2.828427125 ,2.236067978, 3.31662479, 3.16227766 ,2.828427125, 1.732050808 ,NA, NA, 4.582575695, 4.123105626, 2.828427125, NA, NA, 3.16227766, NA, NA, NA, NA, 0, 4.582575695, NA, NA, NA, 3.605551276, NA, NA, NA, NA, 3.741657387, NA, NA ,NA, NA, 3 ,NA ,3.16227766, 5.477225575, NA, 3.31662479, 2.645751311, 3.741657387, 4 ,NA, 7.280109889 ,5.196152423, NA, NA,NA, 11.44552314 ,11.35781669, 11.04536102, 8.124038405, NA, 1.732050808 ,2.645751311, NA ,NA, NA, 14.69693846, 15.77973384 ,3.741657387 ,11.78982612 ,NA, 2.828427125 ,2.645751311, NA ,NA ,NA, 6.08276253 ,5.477225575 ,4.582575695, 4 ,NA, 14.14213562 ,NA ,14.24780685 ,14.4568323, NA, 5.65685425 ,7.549834435 ,3.605551276, 2.645751311 ,NA, 8.185352772 ,6.32455532 ,5.65685425 ,7.348469228 ,NA, 15.32970972 ,NA ,11.5758369 ,7.280109889, NA, 6.557438524 ,7.416198487,NA, NA, NA, 5.916079783, 4.358898944 ,3.31662479, NA, NA, 15.09966887 ,15.32970972 ,5.385164807, 17.1464282, NA, 10.44030651 ,11.13552873, 9.38083152 ,NA ,NA, 3.31662479 ,NA ,NA ,NA,NA, 9.055385138 ,NA ,NA ,NA ,NA, 10.58300524 ,3, 2.236067978, NA,NA, 13.45362405 ,10.48808848 ,NA ,NA, NA, 8.717797887 ,4.242640687, 3.31662479, 2.645751311, NA, 12.32882801 ,12.9614814, 12, 8.246211251, NA, 2.645751311 ,2.828427125, 2.236067978 ,1.414213562, NA, 15.16575089 ,15.90597372 ,4.123105626, 13.03840481, NA, 1.414213562 ,1.732050808 ,1.414213562 ,2 ,NA, 1 ,1.414213562 ,1, 1 ,NA, 3, 4.69041576, 2.828427125, 3.31662479, NA, 16.76305461, 17.32050808, 19.18332609, 14.2126704, 11.40175425, 2.236067978 ,NA, NA ,NA, NA, 3.464101615, 5.916079783 ,5.830951895 ,3.741657387, 3.741657387, 16.64331698, 9.16515139, NA, NA ,NA, 2 ,1.732050808 ,1 ,1, NA, 1.414213562 ,1.414213562 ,1.414213562, 1.414213562 ,NA, 15.93737745 ,16.85229955, 17.05872211 ,14.28285686 ,NA, 7.141428429, 8.660254038, 8.831760866, NA, NA, 2.449489743 ,1.414213562 ,1 ,NA ,NA, 4.358898944 ,5.196152423 ,5.830951895 ,3.464101615, NA, 4 ,3.16227766 ,1.732050808 ,2, NA, 1.414213562 ,2.449489743, 1 ,1, NA, 2.828427125 ,3.464101615, 2 ,NA, NA, 4.123105626 ,2 ,2.645751311, NA, NA, 3.31662479 ,2.828427125, NA, 2.236067978, NA, 6.403124237 ,7.211102551, 5.477225575, 1, NA, 7.348469228 ,6.244997998, 4.795831523, 3.741657387, 4.358898944, 4 ,NA ,3.16227766 ,1.414213562 ,2, 8.366600265 ,5.916079783 ,7, 5.099019514, NA, 3.16227766 ,NA ,NA,NA ,NA, 3.16227766 ,3.16227766, 2 ,NA ,NA, 2.828427125 ,3 ,2.449489743 ,3 ,NA, 3.464101615 ,NA ,2.449489743 ,NA, NA, 3.16227766 ,2.828427125 ,1.732050808, 2.645751311, NA, 4.123105626, 1.732050808 ,1 ,NA ,NA, 3.872983346 ,2.236067978, 1 ,1, NA, 5.567764363 ,5.099019514 ,5.196152423, 3.464101615, NA, 6.32455532 ,3 ,NA ,NA ,NA, 6.633249581 ,5.477225575 ,NA, NA, NA, 6.32455532 ,2.449489743, NA ,NA ,NA, 3.464101615 ,NA, NA, NA, NA, 8.94427191 ,5.916079783, 6.8556546, 1.732050808, NA, 9.949874371, 14.14213562 ,12.32882801 ,13.82027496, NA, 2.828427125 ,NA, NA, 2.236067978, NA, 1.732050808, NA ,NA ,NA, NA, 5.830951895 ,5.830951895, 6 ,3.605551276, NA, 7.745966692 ,5.744562647 ,6.08276253, 3.741657387, NA, 6.164414003 ,4.123105626 ,3.464101615, 7.416198487, NA, 5.567764363 ,4.69041576, 3.16227766 ,NA ,NA, 5.830951895 ,3.605551276, 3 ,NA,NA, 3.31662479, 3.31662479 ,1.732050808 ,NA ,NA, 8.774964387, 7.141428429 ,6.08276253, NA, NA, 3 ,4.582575695 ,NA ,NA ,NA, 11.04536102 ,10.58300524 ,8.774964387, 8.485281374, NA, 16.52271164 ,11.35781669 ,12.76714534 ,12.40967365, NA, 9 ,7.211102551 ,5.099019514, NA ,NA, 3.16227766, 2.645751311, 3.741657387, 3.16227766, NA, 2 ,NA ,2.449489743, NA, NA, 5.916079783, NA ,3.31662479 ,2, NA, 8.062257748 ,7.615773106, 4.69041576, 3.872983346, NA, 2.645751311, 2 ,2.449489743, NA, NA, 13.41640787 ,15.42724862 ,13.07669683 ,8.062257748 ,NA, 2, 2.828427125 ,1.732050808 ,NA ,NA, 9.643650761 ,11.91637529 ,NA ,3.872983346, NA, 17.23368794 ,11.04536102 ,9.433981132, 6.32455532, NA, 17.32050808 ,18.16590213, 17.02938637, 16.43167673 ,18.43908892, 3.605551276 ,3.16227766 ,NA ,NA,NA, 2.828427125 ,3.872983346, 4.123105626, 3.464101615, 2.828427125, 1.732050808 ,1.732050808, 2.828427125, 2.449489743 ,1.732050808, 2.645751311 ,3 ,1.414213562 ,NA, NA, 16.52271164 ,12.32882801 ,8.062257748 ,6.164414003 ,4.582575695, 11.70469991 ,11, 12.08304597, 11.74734012 ,NA, 6.480740698 ,NA ,NA, NA, NA, 2.828427125 ,NA ,NA, NA ,NA, 7.874007874 ,NA ,NA ,NA ,NA, 1.732050808, 2, 1.732050808 ,2 ,NA, 3.605551276 ,2.645751311 ,1.732050808 ,NA ,NA, 6.08276253 ,4.123105626 ,4.472135955 ,1 ,NA, 2.449489743 ,1.732050808 ,1, NA ,NA, 11.3137085 ,NA, 4.898979486, 3.872983346 ,NA, 8 ,5.291502622 ,NA , NA, NA, 1.414213562, NA, NA, NA, NA, 11.18033989 ,12.12435565 ,14.24780685 ,10.63014581, NA, 4.242640687 ,5.477225575 ,4.69041576, NA ,NA, 4.582575695 ,NA, 2.828427125 ,2.645751311, NA, 6.8556546 ,7.483314774,NA ,NA ,NA, 2.236067978 ,2.236067978 ,2 ,NA ,NA, 3.31662479 ,3.464101615 ,2.828427125, NA ,NA, 13.56465997 ,11.74734012, 6.08276253, NA, NA, 3.464101615 ,NA, NA ,NA, NA, 2.645751311 ,3.16227766, 3, 2, NA, 6.32455532 ,5.385164807 ,5.291502622, 4.472135955,NA, 10 ,12.04159458 ,9.38083152, 6.708203933, NA, 10.63014581 ,10.77032961, 13.15294644, 11.22497216, NA, 8.602325267 ,7.483314774 ,5.385164807 ,5.196152423, NA, 12.68857754 ,13.26649916 ,NA ,5.477225575 ,NA, 2.236067978 ,NA, NA ,NA ,NA, 5.567764363 ,4.898979486 ,5.567764363, NA,NA, 15.26433752,16.40121947 ,13.49073756 ,8.426149773, NA, 7.937253933 ,8.426149773, 4 ,4 ,NA, 15.71623365 ,17 ,14.79864859 ,NA ,NA, 3.16227766 ,4.472135955 ,NA ,NA ,NA, 17.05872211 ,19.15724406, 17.97220076, 16.34013464 ,NA, 5.477225575, 6.403124237 ,4.69041576 ,4, NA, 2.236067978, NA,NA ,NA, NA, 10.39230485,21.79449472 ,16.85229955 ,12.28820573, NA, 3.605551276 ,3.872983346 ,NA,NA, NA, 18.46618531,NA , NA ,NA ,NA, 2.449489743 ,2.449489743, NA, NA, NA, 2.645751311 ,3.464101615 ,3, 0 ,NA, 4.795831523 ,7.348469228, NA ,NA ,NA, 5.196152423 ,5.744562647 ,4.242640687, 2.236067978, 3.464101615, 7.211102551 ,6.92820323, 9.16515139, 7.937253933 ,4.472135955, 0 ,NA ,3.872983346 ,NA ,NA, 0 ,NA, NA,NA ,NA, 2 ,3.16227766, 0, NA, NA, 3.741657387, NA, NA, 0 ,NA, 6.92820323 ,11.95826074 ,NA, NA, NA, 7 ,8.774964387 ,NA ,9.486832981 ,NA, 2.236067978 ,1.414213562 ,NA ,NA ,NA, 7.937253933 ,10.48808848 ,7.416198487 ,5.099019514 ,NA, 2.449489743 ,0 ,NA ,2.236067978,NA, 2.828427125 ,0 ,NA, NA ,NA, 3.16227766, NA, 2.645751311 ,NA ,NA, 15.39480432 ,17.49285569 ,16.91153453 ,16.61324773, NA, 8.124038405 ,8.366600265 ,NA ,NA ,NA, 12.4498996 ,10.19803903 ,NA, 13.56465997 ,NA, 16.76305461 ,NA ,NA ,NA, NA, 0 ,NA, NA, NA ,NA, 14.17744688 ,15.42724862 ,15.49193339, 20.4450483, NA, 2.449489743, 0 ,NA,NA ,NA, 4.472135955 ,3.464101615 ,3.464101615,NA ,NA, 5.196152423, 4.472135955 ,3.31662479 ,NA ,NA, 10.48808848 ,7.549834435, 6.403124237,NA ,NA, 6.08276253 ,6.32455532 ,NA ,NA ,NA, 1.732050808 ,2.236067978 ,1 ,NA, NA, 13.26649916 ,14.38749457 ,19.57038579 ,12.68857754, NA, 14.73091986 ,9.486832981,9.848857802, 3 ,NA, 9.16515139 ,6.244997998 ,NA,NA, NA, 4.69041576 ,3.741657387, 2.449489743, NA, NA, 2.828427125, 5.385164807 ,4 ,2.645751311 ,NA, 8.94427191,7.141428429 ,NA, NA, NA, 14.96662955 ,13.92838828 ,11.18033989 ,NA, NA, 3.605551276 ,4.358898944,NA ,1.732050808 ,1.414213562, 9.797958971 ,7.280109889, NA, NA ,NA, 3.16227766,4.123105626 ,3 ,NA ,NA, 5.477225575 ,NA, NA, NA, NA, 7.549834435 ,7.416198487, NA ,3.31662479, NA, 1.414213562 ,1.732050808 ,3.31662479, NA, NA, 4.358898944 ,NA ,NA ,NA ,NA, 6.244997998 ,3.741657387, 4 ,NA ,NA, 16.97056275, 17.97220076 ,18.89444363 ,19.20937271 ,NA, 3.872983346 ,3.31662479 ,1.732050808 ,NA, NA, 6.557438524, 3.605551276, NA, NA, NA, 3.741657387,3.741657387, 3.16227766, NA, NA, 11.18033989 ,NA ,NA ,NA ,NA, 2.645751311 ,2.645751311, 1.732050808 ,3.31662479 ,NA, 3.741657387, 3, 3.16227766, 2.645751311, NA, 1.414213562, 2.828427125 ,NA ,NA, NA, 3.31662479, 2.645751311, 2.828427125, 0, NA, 2.236067978, 2.645751311, NA ,NA, NA, 6.782329983, 6.8556546, 5.916079783, 5, NA, 1.732050808, 1.414213562, NA, NA ,NA, 11.48912529, 9.486832981, 7.745966692, 4.795831523, NA, 9.16515139, 7.549834435, 7.348469228, NA ,NA, 5.099019514, NA, NA, NA, NA, 15.45962483, 14.14213562, 9.433981132, 7.141428429, NA, 2, NA ,NA ,NA, NA, 5.385164807, 4.795831523, 5.916079783, NA, NA, 4.582575695, 10.63014581, 2.236067978, NA, NA, 12.88409873, NA, NA, NA, NA, 16.52271164, 9.327379053, 18.94729532, 16.88194302, NA, 13.45362405, NA, NA, NA ,NA, 1.414213562, 2 ,1.414213562, 0, NA, 5.291502622, 6.557438524 ,5.291502622 ,5.291502622 ,NA, 8.366600265, 7.745966692, 8.124038405, 7, NA, 8.717797887 ,10.24695077, 6.32455532, NA ,NA, 8.831760866 ,NA, NA, NA ,NA, 8.660254038 ,11.48912529, NA ,NA, NA, 6 ,NA ,NA ,NA ,NA, 3.464101615 ,2.449489743 ,4, 3.605551276, NA, 6.32455532 ,NA, 5.291502622, NA ,NA, 4, 7.071067812, 4.123105626, 2.236067978, NA, 2.236067978, NA ,NA, NA ,NA, 1.732050808, 3.464101615, NA ,NA ,NA, 7.483314774, 6.708203933 ,3, 2.828427125 ,NA, 8.124038405, 8.831760866 ,7.071067812 ,5.65685425 ,NA, 8.062257748, 9.797958971 ,7.745966692 ,9.539392014,NA, 16.4924225 ,15.49193339 ,15.29705854, 9.899494937 ,NA, 7.745966692 ,4.898979486 ,9.486832981, 8.717797887, NA, 10.24695077 ,16.4924225 ,19.36491673 ,12.24744871 ,NA, 6 ,5.65685425 ,6.480740698,NA ,NA, 4.898979486 ,2 ,NA ,NA, NA, 17.32050808 ,16.43167673 ,11.95826074, 12.64911064 ,NA, 4.472135955,4.242640687 ,3.16227766 ,NA,NA, 3.605551276,3.741657387 ,2.449489743 ,2.645751311,2.645751311, 6.480740698, 4.472135955,5.477225575 ,4.242640687, NA, 5.196152423 ,6.480740698, 2.449489743 ,2.449489743,NA, 2.645751311 ,3 ,4.582575695 ,NA ,NA, 8.366600265, 7.937253933, 5.477225575 ,NA ,NA, 12.60952021, 8.544003745 ,NA ,NA ,NA, 3.741657387,4 ,2.236067978, NA ,NA, 11.74734012 ,11.83215957 ,12.40967365 ,10.95445115 ,11.22497216, 13 ,12, 9.16515139, NA ,NA, 6.32455532 ,6, NA ,3.16227766, NA, 6.32455532 ,NA ,NA ,NA ,NA, 4.898979486, 3.872983346 ,2.236067978, 2.236067978 ,1.732050808, 4.582575695,3.464101615 ,4.898979486 ,10.09950494 ,NA, 2.449489743, 4.898979486 ,NA,NA ,NA, 3.16227766 ,4.472135955 ,3.741657387 ,1.732050808, NA, 13.78404875 ,11.66190379 ,NA, NA, NA, 14.83239697 ,20, 15.49193339 ,NA ,NA, 15.71623365,17.32050808 ,14.83239697, 15.55634919, NA, 3 ,5.477225575,3.16227766, NA, NA, 1.732050808 ,4.582575695 ,NA,NA ,NA, 2.645751311 ,3.16227766 ,2.449489743,2.645751311 ,NA, 4.898979486, NA ,NA ,NA ,NA, 3.16227766 ,2.449489743, 3.464101615,NA ,NA, 3 ,3.16227766, NA,NA ,NA, 5.477225575 ,2.236067978,NA ,NA ,NA), .Dim=c(467,5)), surt=c(NA, NA, 18.53, NA, NA, 1.9, NA, 9.57, 11.57, NA, 2.4, NA, 12.23, 2.6, NA, 7.1, NA, NA, 2.73, NA, 8.07, NA, 12.53, NA, 4.57, NA, NA, NA, 6.13, NA, 16.27, 7.17, NA, NA, NA, 7.67, 14.87, NA, NA, NA, NA, NA, NA, 8.47, 8.43, 11.07, 11.07, 5.9, 2, NA, 7.57, 5.57, 4.93, NA, NA, NA, 12.2, NA, 10.57, 15.7, NA, NA, NA, 7.17, 2.6, 15.57, 11.13, NA, 0.77, NA, 6, NA, 1.3, NA, 12.67, NA, 10.33, NA, NA, NA, 3.13, NA, 4.6, NA, NA, NA, 10.4, 9.77, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 10.3, NA, 5.1, NA, 11.03, 11.17, NA, NA, 1.07, 4.67,NA, NA,NA,NA,10.23,15.1,NA,NA,7.97,NA, 3.83,NA,NA,NA,NA,NA,NA,15.47,3.2,NA, NA,NA,NA,NA,13.87,NA,NA,NA,11.3,NA, 13.2,10.93,NA,16.4,NA,NA,NA,14.63,NA,NA, NA,10.23,NA,NA,NA,NA,6.63,NA,NA,NA, NA,NA,NA,5.43,10.3,NA,17.8,NA,16.63,NA, NA,NA,9.37,NA,NA,9.37,NA,NA,NA,11.93, NA,3.13,NA,NA,NA,NA,NA,NA,NA,NA, NA,NA,NA,NA,NA,NA,1.7,NA,NA,NA, NA,7.03,3.3,11.6,NA,4.33,NA,4.53,NA,NA, NA,2.47,NA,12.3,11.9,NA,NA,9.93,12.27,NA, NA,11.2,NA,NA,NA,NA,NA,11.47,10.7,2.17, 2.7,NA,NA,17.03,14.7,13.3,NA,10.87,NA,NA, NA,NA,NA,7.93,NA,NA,NA,14.43,NA,NA, NA,NA,NA,NA,NA,NA,NA,NA,NA,4.53, NA,3.37,NA,NA,NA,NA,9.8,NA,NA,NA, 12.63,10.03,NA,NA,NA,NA,NA,5.03,NA,NA, 11,NA,8.13,NA,NA,NA,5.53,3.57,2.57,NA, NA,NA,1.43,NA,NA,NA,12.37,8.8,NA,9.57, 10.63,NA,NA,NA,NA,NA,12.07,NA,9.6,NA, 11.3,NA,13.9,NA,4.97,NA,18.07,NA,NA,NA, NA,3.63,3.97,NA,10.07,NA,11.67,NA,6.13,2.03, NA,NA,NA,NA,7.4,10.8,NA,1.13,NA,NA,NA, NA,NA,NA,1.07,7.17,NA,NA,NA,11.5,NA, NA,NA,NA,6.3,0.47,4.1,NA,8.17,NA,NA,6.9, 3.43,7.97,15.33,5.37,NA,3.2,NA,NA,18.63, 8.47,NA,11.93,NA,7.83,1.63,NA,6.07,8.3,NA, 13.6,NA,NA,NA,NA,8.33,10.7,NA,9.67,NA, NA,7.7,8.73,2.13,NA,11.63,3,NA,NA,9.63, 15.07,9.93,NA,NA,NA,6.6,NA,5.23,NA,3.4, NA,12.2,6.87,NA,1.03,NA,10.6,4.2,NA,NA, NA,NA,NA,9.73,NA,NA,1.43,NA,8,NA, 2.47,9.07,14.13,15.17,NA,16.63,NA,NA,7.3,5.57, NA,7.73,NA,NA,NA,11.5,12.47,5.07,12.7,NA, NA,19.07,0.9,NA,13.1,NA,17.27,NA,9.33,NA, 10.23,3.3,NA,2,13.3,NA,3.9), surt.cen=c(16.97,19,0,12.7,15.13,0,14.33,0,0,16.43, 0,18.1,0,0,12.23,0,16.67,13.7,0,14.4, 0,17.4,0,20.27,0,20.27,18.13,19.57,0,18.87, 0,0,18.57,18.5,18.4,0,0,17.87,17.7,17.4, 17.37,17.37,17.23,0,0,0,0,0,0,14.4, 0,0,0,16.03,12.27,15.03,0,12.27,0,0, 12.3,14.17,12.7,0,0,0,0,14.17,0,14.07, 0,13.93,0,13.93,0,13.23,0,12.47,12.47,17.93, 0,18.03,0,15.97,14.9,15.6,0,0,18.13,17.93, 18.13,18.03,17.87,15.57,17.43,17.67,15.1,16.3,13.03,0, 15.03,0,12.9,0,0,14.07,18.4,0,0,12.2, 17.83,16.23,12.5,0,0,12.5,12.93,0,14.33, 0,12.5,12.2,12.27,12.77,16.47,20.23,0,0,13.5,19.2,12.23, 17.8,16.7,0,16.3,16.23,20.4,0,18.6,0,0,17.17,0,16.23,16.03, 16.03,0,15.27,14.83,14.6,0,13.83,14.13,13.13,13.03,0,18.37, 16.2,13.47,13.2,18.17,15.57,0,0,18.37,0,17.8,0,17.37,15.9, 17.33,0,16.97,17.1,0,15.8,15.5,15.1,0,14.63,0,14.57,14.33, 13.93,13.37,13.23,12.97,14.87,1.83,13.73,16.93,16.17,12.27, 12.53,12.77,0,2.17,12.3,12.23,12.23,0,0,0,12.2,0,17.57,0, 13.73,18.8,18.57,0,12.43,0,0,13.37,17.67,0,0,18.33,17.67,0, 18.13,18.57,17.93,16.23,17.2,0,0,0,0,12.67,12.7,0,0,0,13.87, 0,16.73,13.2,13.9,15.37,15.97,0,16,15.6,14.63,0,14.1,14.07, 12.47,12.67,12.33,12.2,16.17,15.53,14.4,12.7,19.27,0,20.67, 0,12.23,19.27,14.2,12.23,0,14.2,12.73,12.73,0,0,13.27,19.77, 19.53,18.37,15.93,0,20.27,12.93,0,19.57,0,21.4,21.13,21.13, 0,0,0,12.27,19.57,13.9,0,12.2,15.97,17.83,0,0,16.2,0,0,13.87, 13.43,13.03,12.33,12.3,0,12.33,0,15.97,0,20.87,0,18.97,0, 18.17,0,15.93,18.37,15.97,17.7,0,0,17.37,0,17,0,17.2,0,0, 16.47,16.23,16.03,15.77,0,0,12.57,0,13.2,13.47,13.9,13.93, 14.6,14.43,0,0,15.33,15.57,12.23,0,12.23,12.27,12.23,12.53, 0,0,0,12.23,0,20.5,19.8,0,0,0,0,0,13.9,0,17.83,19.93,0,0,16.87,0,18.77,0,0,13.5,0,0,19.8,0,17.23,16.5,13.4,17.63,0,0,12.53, 0,12.2,18.17,0,0,0,13,0,0,12.7,16.93,0,0,0,16.5,16.97,14.4, 0,12.53,0,13,0,16.03,0,0,15.03,0,14.1,0,0,13.6,13,13.5, 13.2,13.2,0,13.6,20.63,0,12.43,0,12.47,0,0,0,0,17.7,0,14.3, 12.43,0,0,17.2,0,19.8,16.47,13.5,0,0,0,0,19.03,7.93,0,0, 19.7,0,20.87,0,13.37,0,12.3,0,0,14.33,0,0,15.07,0), randgrp1=c(0,1,1,0,1,0,0,1,0,1,0,1,0,1,1,0,0,0,1,0,1,1,1,0,0,1,0,1,1,0,0,1,1,0,1,0,0,1,0,1,0,0,0,1, 0,1,0,0,1,1,1,0,0,0,0,0,1,1,0,0,0,0,0,1,1,1,1,1,0,1,0,0,1,1,1,0,0,1,0,0,1,1,1,0,1,1,0,0,0,1,0,1,0,1, 0,0,1,0,0,0,1,1,0,0,0,1,1,0,0,1,1,0,1,1,0,0,1,0,1,1,1,0,0,1,0,0,1,1,0,0,1,0,0,1,0,0,1,1,0,1,1,1,0,1, 0,0,0,1,1,0,1,0,0,1,1,1,1,1,1,0,0,1,1,1,0,0,1,0,0,0,0,1,0,0,0,0,0,0,1,0,1,1,0,1,1,0,0,1,1,0,1,1,1, 0,1,1,0,1,1,0,1,0,1,1,1,0,1,0,1,1,0,0,1,0,1,0,1,1,0,0,1,1,0,1,0,1,1,0,0,0,1,1,0,0,1,1,1,1,0,0,0,0, 0,1,1,1,1,0,0,1,0,0,0,0,0,0,0,0,1,0,0,1,0,1,0,1,1,0,1,1,0,1,1,1,1,1,1,0,1,0,0,1,1,1,1,1,0,1,1,1,0,0, 0,1,0,0,0,1,1,1,1,0,0,1,0,0,0,0,0,0,0,0,0,0,1,0,1,1,1,0,1,0,1,0,0,0,1,1,1,1,0,0,0,1,0,0,1,0,1,0, 1,1,0,0,0,1,1,1,1,0,0,0,1,1,1,1,0,1,1,1,1,0,0,1,0,0,0,1,0,0,0,0,1,0,1,1,0,1,1,0,0,0,0,1,0,0,0,0,1, 0,1,0,0,1,0,1,0,1,1,0,1,1,1,1,1,1,0,0,1,0,0,1,1,1,0,0,1,1,1,0,1,0,0,1,0,1,0,0,0,0,0,1,1,0,0,1,0,1,1, 0,0,0,0,1,0,1,1,1,1,1,1,0,1,1,0,0,0,1,1,1,1,1,1,0,1,0,0), gender1=c(1,1,-1,1,1,-1,1,-1,1,1, 1,1,1,1,1,1,1,1,-1,1, 1,1,-1,1,-1,1,1,1,1,-1, -1,1,1,1,1,1,1,-1,1,1, 1,-1,1,1,1,1,-1,1,1,-1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,-1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,-1, 1,1,1,1,1,1,1,1,1,1, -1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,-1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, -1,1,1,1,1,1,1,1,1,1, 1,-1,-1,-1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,-1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,-1, 1,1,1,1,-1,1,1,-1,-1,1, 1,1,-1,1,1,1,1,-1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, -1,1,1,1,1,1,1,1,1,1, 1,1,1,-1,1,1,1,1,1,1, -1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,-1,-1,1,-1,-1,-1, 1,-1,1,1,1,1,1,1,1,-1, 1,1,1,-1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,-1,1,1,1,1, 1,1,-1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,1,-1,1, 1,1,1,1,1,1,1,1,1,1, 1,1,1,1,1,1,1,-1,1,1, 1,1,1,1,1,1,-1,1,1,1, 1,1,1,1,1,1,1,1,1,1, 1,-1,1,1,1,1,1,1,1,1, 1,1,1,1,-1,1), stratum1=c(-1,-1,-1,1,1,1,-1,-1,1,-1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,-1,-1,1,1,-1,-1,-1,1,-1,-1,-1,-1,-1,-1,-1, -1,-1,-1,-1,1,-1,1,1,-1,-1,-1,-1,1,1,-1,1,-1,-1,-1,1,-1,-1,1,-1,1,1,1,1,1,-1,-1,-1,-1,-1,1,-1,1,-1,-1,-1,-1,1,1,-1,-1,-1, 1,1,1,-1,-1,1,-1,-1,-1,-1,1,1,-1,1,1,1,-1,1,-1,-1,1,1,1,1,1,-1,1,-1,1,1,1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,-1, -1,-1,-1,-1,-1,1,-1,-1,1,-1,-1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,1,1,1,-1,1,-1,1,-1,-1,1,1,-1,-1,1,1,1,-1,-1,-1,1, -1,-1,-1,1,-1,-1,-1,-1,-1,-1,1,-1,1,-1,-1,-1,-1,-1,-1,1,1,1,-1,1,-1,1,1,1,-1,1,-1,1,-1,-1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1, 1,-1,-1,-1,-1,1,-1,1,-1,1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,-1,1,-1,1,1,1,-1,1,1,-1,-1,-1,-1,-1,1,1,-1,-1,1,-1,-1,-1,-1,1,-1,-1,1, 1,-1,1,-1,1,-1,1,-1,1,1,-1,-1,1,1,1,-1,1,-1,-1,-1,-1,1,1,1,1,-1,-1,-1,1,1,1,1,-1,1,-1,1,-1,-1,-1,-1,1,1,1,1,-1,-1,1,-1,-1, 1,1,1,-1,-1,1,-1,1,-1,-1,-1,1,-1,1,1,-1,-1,-1,-1,-1,-1,-1,1,-1,1,1,-1,-1,-1,-1,-1,-1,1,1,-1,1,-1,1,1,-1,-1,-1,-1,-1,-1,1,1, 1,-1,-1,1,-1,1,-1,1,-1,-1,-1,1,-1,-1,1,1,-1,1,-1,1,-1,-1,1,-1,-1,1,-1,-1,-1,1,-1,1,-1,1,-1,1,-1,-1,1,-1,-1,1,1,-1,-1,-1, -1,-1,-1,-1,-1,1,-1,-1,-1,-1,1,-1,1,-1,1,1,1,-1,-1,-1,1,-1,-1,1,1,-1,-1,-1,-1,1,-1,-1,1,1,-1,1,1,-1,-1,-1,1,-1,-1,1, 1,1,-1,-1,1,-1,1), prevoi1=c(1,-1,1,1,1,1,1,-1,1,1,1,-1,-1,-1,1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1,-1,1,1,1,1,-1,1,1,-1,-1,-1,1,-1,-1,-1,1, -1,1,1,-1,-1,1,1,1,1,1,1,1,-1,1,1,1,-1,1,1,1,1,1,1,1,1,1,-1,1,1,1,-1,1,-1,-1,-1,1,1,1,1,1,1,1,1,-1,-1,-1,1,-1,1,-1,-1,1, 1,-1,1,1,1,-1,1,1,-1,1,1,1,1,1,-1,1,1,1,1,1,1,-1,1,-1,1,1,1,1,-1,1,1,1,-1,-1,-1,1,1,-1,-1,-1,-1,-1,1,1,-1,1,-1,-1,1,1,-1, -1,-1,-1,-1,-1,-1,1,1,-1,-1,-1,1,1,1,1,1,-1,1,-1,1,1,-1,1,1,-1,-1,1,1,1,-1,1,-1,1,-1,1,-1,1,-1,-1,1,1,1,1,1,-1,1,-1,1,-1,-1, -1,1,1,1,1,1,1,-1,1,1,1,-1,1,-1,1,1,-1,1,1,1,1,1,1,-1,-1,-1,-1,1,1,1,1,1,1,1,1,-1,1,-1,1,1,-1,-1,-1,-1,1,1,-1,-1,1,-1,1, -1,1,1,1,-1,1,1,1,-1,1,1,1,1,1,-1,1,1,1,1,-1,1,1,-1,-1,1,1,-1,1,-1,1,1,1,1,1,1,1,1,1,1,1,-1,1,1,1,-1,1,1,1,1,1,1,-1, -1,1,1,1,1,-1,1,-1,1,-1,-1,-1,1,1,1,1,1,1,-1,1,-1,1,1,1,1,-1,-1,1,-1,1,-1,1,1,1,-1,1,1,-1,1,-1,1,-1,1,1,1,1,1,1,-1,-1,1,1, -1,-1,1,1,1,1,1,1,1,1,1,-1,-1,1,-1,1,1,1,-1,-1,1,-1,1,-1,1,1,1,1,1,-1,-1,1,1,1,1,-1,1,1,-1,1,1,1,1,1,1,-1,1,1,1, -1,1,1,1,1,1,1,1,-1,1,1,-1,1,1,-1,-1,-1,1,1,1,1,-1,-1,-1,1,1,1,1,1,1,1,1,-1,-1,1,-1,1,1,1,-1,-1,1,1,1,-1,1,1,1,1,1,1,-1, -1,-1,1,1,-1,1,1,1,-1,1,1,1,1)) #initial values list(beta1=c(0,0,0,0,0,0),tauz=1,beta2=c(0,0,0,0,0),r1=0,r2=0, U=structure(.Data=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, 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, 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, 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, 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, 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, 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, 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, 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, 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,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,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,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,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,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,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,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,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,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),.Dim=c(467,2))) list(beta1=c(2,2,2,2,2,2),tauz=2,beta2=c(2,2,2,2,2),r1=2,r2=2, U=structure(.Data=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, 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, 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, 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, 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, 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, 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, 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, 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, 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,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,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,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,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,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,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,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,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,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),.Dim=c(467,2))) list(beta1=c(-2,-2,-2,-2,-2,-2),tauz=3,beta2=c(-2,-2,-2,-2,-2),r1=-2,r2=-2, U=structure(.Data=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, 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, 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, 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, 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, 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, 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, 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, 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, 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,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,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,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,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,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,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,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,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,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),.Dim=c(467,2)))