INITS: # we give uninformative values of 0 to each of the model coefficients list(alpha=c(0,0,0,0,0,0,0,0,0,0,0,0), #assigning initial values of 1 to each of the variance terms for spatial (tau.phi) and non-spatial (tau.theta) # random effects helps convergence and is only vaguely informative tau.phi = c(1,1,1,1,1,1,1,1,1,1), tau.theta = c(1,1,1,1,1,1,1,1,1,1), # this first Markov chain assigns 0 to the spatial residual random effect terms in control arm (phiC) phiC1=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), ############################################# # phi2 inits phiC2=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), ############################################## # phi3 inits phiC3=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), ############################################## # phi4 inits phiC4=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), ############################################## # phi5 inits phiC5=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), ############################################# # phi6 inits phiC6=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), ############################################## # phi7 inits phiC7=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), ############################################## # phi8 inits phiC8=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), ############################################## # phi9 inits phiC9=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), ############################################# # phi10 inits phiC10=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), ############################################## # this first Markov chain assigns 0 to the spatial residual random effect terms in treatment arm (phiT) phiT1=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), ############################################# # phi2 inits phiT2=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), ############################################## # phi3 inits phiT3=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), ############################################## # phi4 inits phiT4=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), ############################################## # phiT5 inits phiT5=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), ############################################# # phiT6 inits phiT6=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), ############################################## # phiT7 inits phiT7=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), ############################################## # phiT8 inits phiT8=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), ############################################## # phiT9 inits phiT9=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), ############################################# # phiT10 inits phiT10=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), ############################################# ############################################################################ # now start all the theta inits ############################################# # theta1 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta1 corresponds to the first time period theta1=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) , ############################################# # theta2 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta2 corresponds to the second time period theta2=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) , ############################################## # theta3 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta3 corresponds to the third time period theta3=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) , ############################################## # theta4 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta4 corresponds to the fourth time period theta4=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) , ############################################## # in this chain all non-spatial random effects are assigned 0 as initial value # where theta5 corresponds to the fifth time period theta5=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) , ############################################# # theta6 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta6 corresponds to the sixth time period theta6=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) , ############################################## # theta7 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta7 corresponds to the seventh time period theta7=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) , ############################################## # theta8 inits # in this chain all non-spatial random effects are assigned 0 as initial value # where theta8 corresponds to the eighth time period theta8=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) , ############################################## # in this chain all non-spatial random effects are assigned 0 as initial value # where theta9 corresponds to the ninth time period theta9=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) , ############################################# # in this chain all non-spatial random effects are assigned 0 as initial value # where theta10 corresponds to the tenth time period theta10=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) ) ############################################# ############################################## ####################################################################### # inits for second Markov chain # the variables have same meanings as documented in the first Markov chain # but different starting values to provide disparate chains that better # demonstrate convergence as indicated in the output document list(alpha=c(0.5,0.5,0.5,0.5,0.5,0.5, 0.5,0.5,0.5,0.5, 0.5, 0.5), tau.phi = c(0.25,0.25,0.25,0.25, 0.25,0.25,0.25,0.25,0.25,0.25), tau.theta = c(0.25,0.25,0.25,0.25, 0.25,0.25,0.25,0.25,0.25,0.25), phiC1=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# # phi2 inits phiC2=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi3 inits phiC3=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi4 inits phiC4=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi5 inits phiC5=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# # phi6 inits phiC6=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi7 inits phiC7=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi8 inits phiC8=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi9 inits phiC9=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# # phi10inits phiC10=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# phiT1=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# # phi2 inits phiT2=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi3 inits phiT3=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi4 inits phiT4=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi5 inits phiT5=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# # phi6 inits phiT6=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi7 inits phiT7=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi8 inits phiT8=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################## # phi9 inits phiT9=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# # phi10 inits phiT10=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005), ############################################# ############################################################################ # now start all the theta inits ############################################# # theta1 inits theta1=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################# # theta2 inits theta2=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################## # theta3 inits theta3=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################## # theta4 inits theta4=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################## theta5=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################# # theta6 inits theta6=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################## # theta7 inits theta7=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################## # theta8 inits theta8=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################## theta9=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) , ############################################# theta10=c( 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005) ) ############################################# ####################################################################### # inits for third Markov chain # the variables have same meanings as documented in the first Markov chain # but different starting values to provide disparate chains that better # demonstrate convergence as indicated in the output document list(alpha=c(-0.5,-0.5,-0.5,-0.5,-0.5,-0.5, -0.5,-0.5,-0.5,-0.5, -0.5, -0.5), tau.phi = c(0.5,0.5,0.5,0.5, 0.5,0.5,0.5,0.5,0.5,0.5), tau.theta = c(0.5,0.5,0.5,0.5, 0.5,0.5,0.5,0.5,0.5,0.5), phiC1=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################# # phi2 inits phiC2=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi3 inits phiC3=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi4 inits phiC4=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## phiC5=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################# # phi6inits phiC6=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi7 inits phiC7=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi8inits phiC8=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi9 inits phiC9=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi10 inits phiC10=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## phiT1=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################# # phi2 inits phiT2=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi3 inits phiT3=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi4 inits phiT4=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################ phiT5=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################# # phi6 inits phiT6=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi7 inits phiT7=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi8 inits phiT8=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################ # phi9 inits phiT9=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## # phi10inits phiT10=c( -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04, -1.00E-04), ############################################## ############################################# # theta1 inits theta1=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################# # theta2 inits theta2=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################## # theta3 inits theta3=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################## # theta4 inits theta4=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################## theta5=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################# # theta6 inits theta6=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################## # theta7 inits theta7=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################## # theta8 inits theta8=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################## theta9=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) , ############################################# theta10=c( -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03, -5.00E-03) ) ############################################# DAS ENDE