model{ for (i in 1:3875){ x1[i]<-beta1[1]+transparency[wdinumber[i],year[i]]*beta1[2] x2[i]<-beta2[1]+transparency[wdinumber[i],year[i]]*beta2[2] x3[i]<-beta3[1]+transparency[wdinumber[i],year[i]]*beta3[2] x4[i]<-beta4[1]+transparency[wdinumber[i],year[i]]*beta4[2] x5[i]<-beta5[1]+transparency[wdinumber[i],year[i]]*beta5[2] x6[i]<-beta6[1]+transparency[wdinumber[i],year[i]]*beta6[2] x7[i]<-beta7[1]+transparency[wdinumber[i],year[i]]*beta7[2] x8[i]<-beta8[1]+transparency[wdinumber[i],year[i]]*beta8[2] x9[i]<-beta9[1]+transparency[wdinumber[i],year[i]]*beta9[2] x10[i]<-beta10[1]+transparency[wdinumber[i],year[i]]*beta10[2] x11[i]<-beta11[1]+transparency[wdinumber[i],year[i]]*beta11[2] x12[i]<-beta12[1]+transparency[wdinumber[i],year[i]]*beta12[2] x13[i]<-beta13[1]+transparency[wdinumber[i],year[i]]*beta13[2] x14[i]<-beta14[1]+transparency[wdinumber[i],year[i]]*beta14[2] x15[i]<-beta15[1]+transparency[wdinumber[i],year[i]]*beta15[2] x16[i]<-beta16[1]+transparency[wdinumber[i],year[i]]*beta16[2] x17[i]<-beta17[1]+transparency[wdinumber[i],year[i]]*beta17[2] x18[i]<-beta18[1]+transparency[wdinumber[i],year[i]]*beta18[2] x19[i]<-beta19[1]+transparency[wdinumber[i],year[i]]*beta19[2] x20[i]<-beta20[1]+transparency[wdinumber[i],year[i]]*beta20[2] x21[i]<-beta21[1]+transparency[wdinumber[i],year[i]]*beta21[2] x22[i]<-beta22[1]+transparency[wdinumber[i],year[i]]*beta22[2] x23[i]<-beta23[1]+transparency[wdinumber[i],year[i]]*beta23[2] x24[i]<-beta24[1]+transparency[wdinumber[i],year[i]]*beta24[2] x25[i]<-beta25[1]+transparency[wdinumber[i],year[i]]*beta25[2] x26[i]<-beta26[1]+transparency[wdinumber[i],year[i]]*beta26[2] x27[i]<-beta27[1]+transparency[wdinumber[i],year[i]]*beta27[2] x28[i]<-beta28[1]+transparency[wdinumber[i],year[i]]*beta28[2] x29[i]<-beta29[1]+transparency[wdinumber[i],year[i]]*beta29[2] x30[i]<-beta30[1]+transparency[wdinumber[i],year[i]]*beta30[2] x31[i]<-beta31[1]+transparency[wdinumber[i],year[i]]*beta31[2] x32[i]<-beta32[1]+transparency[wdinumber[i],year[i]]*beta32[2] x33[i]<-beta33[1]+transparency[wdinumber[i],year[i]]*beta33[2] x34[i]<-beta34[1]+transparency[wdinumber[i],year[i]]*beta34[2] x35[i]<-beta35[1]+transparency[wdinumber[i],year[i]]*beta35[2] x36[i]<-beta36[1]+transparency[wdinumber[i],year[i]]*beta36[2] x37[i]<-beta37[1]+transparency[wdinumber[i],year[i]]*beta37[2] x38[i]<-beta38[1]+transparency[wdinumber[i],year[i]]*beta38[2] x39[i]<-beta39[1]+transparency[wdinumber[i],year[i]]*beta39[2] x40[i]<-beta40[1]+transparency[wdinumber[i],year[i]]*beta40[2] x41[i]<-beta41[1]+transparency[wdinumber[i],year[i]]*beta41[2] x42[i]<-beta42[1]+transparency[wdinumber[i],year[i]]*beta42[2] x43[i]<-beta43[1]+transparency[wdinumber[i],year[i]]*beta43[2] x44[i]<-beta44[1]+transparency[wdinumber[i],year[i]]*beta44[2] x45[i]<-beta45[1]+transparency[wdinumber[i],year[i]]*beta45[2] x46[i]<-beta46[1]+transparency[wdinumber[i],year[i]]*beta46[2] x47[i]<-beta47[1]+transparency[wdinumber[i],year[i]]*beta47[2] x48[i]<-beta48[1]+transparency[wdinumber[i],year[i]]*beta48[2] x49[i]<-beta49[1]+transparency[wdinumber[i],year[i]]*beta49[2] x50[i]<-beta50[1]+transparency[wdinumber[i],year[i]]*beta50[2] x51[i]<-beta51[1]+transparency[wdinumber[i],year[i]]*beta51[2] x52[i]<-beta52[1]+transparency[wdinumber[i],year[i]]*beta52[2] x53[i]<-beta53[1]+transparency[wdinumber[i],year[i]]*beta53[2] x54[i]<-beta54[1]+transparency[wdinumber[i],year[i]]*beta54[2] x55[i]<-beta55[1]+transparency[wdinumber[i],year[i]]*beta55[2] x56[i]<-beta56[1]+transparency[wdinumber[i],year[i]]*beta56[2] x57[i]<-beta57[1]+transparency[wdinumber[i],year[i]]*beta57[2] x58[i]<-beta58[1]+transparency[wdinumber[i],year[i]]*beta58[2] x59[i]<-beta59[1]+transparency[wdinumber[i],year[i]]*beta59[2] x60[i]<-beta60[1]+transparency[wdinumber[i],year[i]]*beta60[2] x61[i]<-beta61[1]+transparency[wdinumber[i],year[i]]*beta61[2] x62[i]<-beta62[1]+transparency[wdinumber[i],year[i]]*beta62[2] x63[i]<-beta63[1]+transparency[wdinumber[i],year[i]]*beta63[2] x64[i]<-beta64[1]+transparency[wdinumber[i],year[i]]*beta64[2] x65[i]<-beta65[1]+transparency[wdinumber[i],year[i]]*beta65[2] x66[i]<-beta66[1]+transparency[wdinumber[i],year[i]]*beta66[2] x67[i]<-beta67[1]+transparency[wdinumber[i],year[i]]*beta67[2] x68[i]<-beta68[1]+transparency[wdinumber[i],year[i]]*beta68[2] x69[i]<-beta69[1]+transparency[wdinumber[i],year[i]]*beta69[2] x70[i]<-beta70[1]+transparency[wdinumber[i],year[i]]*beta70[2] x71[i]<-beta71[1]+transparency[wdinumber[i],year[i]]*beta71[2] x72[i]<-beta72[1]+transparency[wdinumber[i],year[i]]*beta72[2] x73[i]<-beta73[1]+transparency[wdinumber[i],year[i]]*beta73[2] x74[i]<-beta74[1]+transparency[wdinumber[i],year[i]]*beta74[2] x75[i]<-beta75[1]+transparency[wdinumber[i],year[i]]*beta75[2] x76[i]<-beta76[1]+transparency[wdinumber[i],year[i]]*beta76[2] x77[i]<-beta77[1]+transparency[wdinumber[i],year[i]]*beta77[2] x78[i]<-beta78[1]+transparency[wdinumber[i],year[i]]*beta78[2] x79[i]<-beta79[1]+transparency[wdinumber[i],year[i]]*beta79[2] x80[i]<-beta80[1]+transparency[wdinumber[i],year[i]]*beta80[2] x81[i]<-beta81[1]+transparency[wdinumber[i],year[i]]*beta81[2] x82[i]<-beta82[1]+transparency[wdinumber[i],year[i]]*beta82[2] x83[i]<-beta83[1]+transparency[wdinumber[i],year[i]]*beta83[2] x84[i]<-beta84[1]+transparency[wdinumber[i],year[i]]*beta84[2] x85[i]<-beta85[1]+transparency[wdinumber[i],year[i]]*beta85[2] x86[i]<-beta86[1]+transparency[wdinumber[i],year[i]]*beta86[2] x87[i]<-beta87[1]+transparency[wdinumber[i],year[i]]*beta87[2] x88[i]<-beta88[1]+transparency[wdinumber[i],year[i]]*beta88[2] x89[i]<-beta89[1]+transparency[wdinumber[i],year[i]]*beta89[2] x90[i]<-beta90[1]+transparency[wdinumber[i],year[i]]*beta90[2] x91[i]<-beta91[1]+transparency[wdinumber[i],year[i]]*beta91[2] x92[i]<-beta92[1]+transparency[wdinumber[i],year[i]]*beta92[2] x93[i]<-beta93[1]+transparency[wdinumber[i],year[i]]*beta93[2] x94[i]<-beta94[1]+transparency[wdinumber[i],year[i]]*beta94[2] x95[i]<-beta95[1]+transparency[wdinumber[i],year[i]]*beta95[2] x96[i]<-beta96[1]+transparency[wdinumber[i],year[i]]*beta96[2] x97[i]<-beta97[1]+transparency[wdinumber[i],year[i]]*beta97[2] x98[i]<-beta98[1]+transparency[wdinumber[i],year[i]]*beta98[2] x99[i]<-beta99[1]+transparency[wdinumber[i],year[i]]*beta99[2] x100[i]<-beta100[1]+transparency[wdinumber[i],year[i]]*beta100[2] x101[i]<-beta101[1]+transparency[wdinumber[i],year[i]]*beta101[2] x102[i]<-beta102[1]+transparency[wdinumber[i],year[i]]*beta102[2] x103[i]<-beta103[1]+transparency[wdinumber[i],year[i]]*beta103[2] x104[i]<-beta104[1]+transparency[wdinumber[i],year[i]]*beta104[2] x105[i]<-beta105[1]+transparency[wdinumber[i],year[i]]*beta105[2] x106[i]<-beta106[1]+transparency[wdinumber[i],year[i]]*beta106[2] x107[i]<-beta107[1]+transparency[wdinumber[i],year[i]]*beta107[2] x108[i]<-beta108[1]+transparency[wdinumber[i],year[i]]*beta108[2] x109[i]<-beta109[1]+transparency[wdinumber[i],year[i]]*beta109[2] x110[i]<-beta110[1]+transparency[wdinumber[i],year[i]]*beta110[2] x111[i]<-beta111[1]+transparency[wdinumber[i],year[i]]*beta111[2] x112[i]<-beta112[1]+transparency[wdinumber[i],year[i]]*beta112[2] x113[i]<-beta113[1]+transparency[wdinumber[i],year[i]]*beta113[2] x114[i]<-beta114[1]+transparency[wdinumber[i],year[i]]*beta114[2] x115[i]<-beta115[1]+transparency[wdinumber[i],year[i]]*beta115[2] x116[i]<-beta116[1]+transparency[wdinumber[i],year[i]]*beta116[2] x117[i]<-beta117[1]+transparency[wdinumber[i],year[i]]*beta117[2] x118[i]<-beta118[1]+transparency[wdinumber[i],year[i]]*beta118[2] x119[i]<-beta119[1]+transparency[wdinumber[i],year[i]]*beta119[2] x120[i]<-beta120[1]+transparency[wdinumber[i],year[i]]*beta120[2] x121[i]<-beta121[1]+transparency[wdinumber[i],year[i]]*beta121[2] x122[i]<-beta122[1]+transparency[wdinumber[i],year[i]]*beta122[2] x123[i]<-beta123[1]+transparency[wdinumber[i],year[i]]*beta123[2] x124[i]<-beta124[1]+transparency[wdinumber[i],year[i]]*beta124[2] x125[i]<-beta125[1]+transparency[wdinumber[i],year[i]]*beta125[2] x126[i]<-beta126[1]+transparency[wdinumber[i],year[i]]*beta126[2] x127[i]<-beta127[1]+transparency[wdinumber[i],year[i]]*beta127[2] x128[i]<-beta128[1]+transparency[wdinumber[i],year[i]]*beta128[2] x129[i]<-beta129[1]+transparency[wdinumber[i],year[i]]*beta129[2] x130[i]<-beta130[1]+transparency[wdinumber[i],year[i]]*beta130[2] x131[i]<-beta131[1]+transparency[wdinumber[i],year[i]]*beta131[2] x132[i]<-beta132[1]+transparency[wdinumber[i],year[i]]*beta132[2] x133[i]<-beta133[1]+transparency[wdinumber[i],year[i]]*beta133[2] x134[i]<-beta134[1]+transparency[wdinumber[i],year[i]]*beta134[2] x135[i]<-beta135[1]+transparency[wdinumber[i],year[i]]*beta135[2] x136[i]<-beta136[1]+transparency[wdinumber[i],year[i]]*beta136[2] x137[i]<-beta137[1]+transparency[wdinumber[i],year[i]]*beta137[2] x138[i]<-beta138[1]+transparency[wdinumber[i],year[i]]*beta138[2] x139[i]<-beta139[1]+transparency[wdinumber[i],year[i]]*beta139[2] x140[i]<-beta140[1]+transparency[wdinumber[i],year[i]]*beta140[2] x141[i]<-beta141[1]+transparency[wdinumber[i],year[i]]*beta141[2] x142[i]<-beta142[1]+transparency[wdinumber[i],year[i]]*beta142[2] x143[i]<-beta143[1]+transparency[wdinumber[i],year[i]]*beta143[2] x144[i]<-beta144[1]+transparency[wdinumber[i],year[i]]*beta144[2] x145[i]<-beta145[1]+transparency[wdinumber[i],year[i]]*beta145[2] x146[i]<-beta146[1]+transparency[wdinumber[i],year[i]]*beta146[2] x147[i]<-beta147[1]+transparency[wdinumber[i],year[i]]*beta147[2] x148[i]<-beta148[1]+transparency[wdinumber[i],year[i]]*beta148[2] x149[i]<-beta149[1]+transparency[wdinumber[i],year[i]]*beta149[2] x150[i]<-beta150[1]+transparency[wdinumber[i],year[i]]*beta150[2] x151[i]<-beta151[1]+transparency[wdinumber[i],year[i]]*beta151[2] x152[i]<-beta152[1]+transparency[wdinumber[i],year[i]]*beta152[2] x153[i]<-beta153[1]+transparency[wdinumber[i],year[i]]*beta153[2] x154[i]<-beta154[1]+transparency[wdinumber[i],year[i]]*beta154[2] x155[i]<-beta155[1]+transparency[wdinumber[i],year[i]]*beta155[2] x156[i]<-beta156[1]+transparency[wdinumber[i],year[i]]*beta156[2] x157[i]<-beta157[1]+transparency[wdinumber[i],year[i]]*beta157[2] x158[i]<-beta158[1]+transparency[wdinumber[i],year[i]]*beta158[2] x159[i]<-beta159[1]+transparency[wdinumber[i],year[i]]*beta159[2] x160[i]<-beta160[1]+transparency[wdinumber[i],year[i]]*beta160[2] x161[i]<-beta161[1]+transparency[wdinumber[i],year[i]]*beta161[2] x162[i]<-beta162[1]+transparency[wdinumber[i],year[i]]*beta162[2] x163[i]<-beta163[1]+transparency[wdinumber[i],year[i]]*beta163[2] x164[i]<-beta164[1]+transparency[wdinumber[i],year[i]]*beta164[2] x165[i]<-beta165[1]+transparency[wdinumber[i],year[i]]*beta165[2] x166[i]<-beta166[1]+transparency[wdinumber[i],year[i]]*beta166[2] x167[i]<-beta167[1]+transparency[wdinumber[i],year[i]]*beta167[2] x168[i]<-beta168[1]+transparency[wdinumber[i],year[i]]*beta168[2] x169[i]<-beta169[1]+transparency[wdinumber[i],year[i]]*beta169[2] x170[i]<-beta170[1]+transparency[wdinumber[i],year[i]]*beta170[2] x171[i]<-beta171[1]+transparency[wdinumber[i],year[i]]*beta171[2] x172[i]<-beta172[1]+transparency[wdinumber[i],year[i]]*beta172[2] x173[i]<-beta173[1]+transparency[wdinumber[i],year[i]]*beta173[2] x174[i]<-beta174[1]+transparency[wdinumber[i],year[i]]*beta174[2] x175[i]<-beta175[1]+transparency[wdinumber[i],year[i]]*beta175[2] x176[i]<-beta176[1]+transparency[wdinumber[i],year[i]]*beta176[2] x177[i]<-beta177[1]+transparency[wdinumber[i],year[i]]*beta177[2] x178[i]<-beta178[1]+transparency[wdinumber[i],year[i]]*beta178[2] x179[i]<-beta179[1]+transparency[wdinumber[i],year[i]]*beta179[2] x180[i]<-beta180[1]+transparency[wdinumber[i],year[i]]*beta180[2] x181[i]<-beta181[1]+transparency[wdinumber[i],year[i]]*beta181[2] x182[i]<-beta182[1]+transparency[wdinumber[i],year[i]]*beta182[2] x183[i]<-beta183[1]+transparency[wdinumber[i],year[i]]*beta183[2] x184[i]<-beta184[1]+transparency[wdinumber[i],year[i]]*beta184[2] x185[i]<-beta185[1]+transparency[wdinumber[i],year[i]]*beta185[2] x186[i]<-beta186[1]+transparency[wdinumber[i],year[i]]*beta186[2] x187[i]<-beta187[1]+transparency[wdinumber[i],year[i]]*beta187[2] x188[i]<-beta188[1]+transparency[wdinumber[i],year[i]]*beta188[2] x189[i]<-beta189[1]+transparency[wdinumber[i],year[i]]*beta189[2] x190[i]<-beta190[1]+transparency[wdinumber[i],year[i]]*beta190[2] x191[i]<-beta191[1]+transparency[wdinumber[i],year[i]]*beta191[2] x192[i]<-beta192[1]+transparency[wdinumber[i],year[i]]*beta192[2] x193[i]<-beta193[1]+transparency[wdinumber[i],year[i]]*beta193[2] x194[i]<-beta194[1]+transparency[wdinumber[i],year[i]]*beta194[2] x195[i]<-beta195[1]+transparency[wdinumber[i],year[i]]*beta195[2] x196[i]<-beta196[1]+transparency[wdinumber[i],year[i]]*beta196[2] x197[i]<-beta197[1]+transparency[wdinumber[i],year[i]]*beta197[2] x198[i]<-beta198[1]+transparency[wdinumber[i],year[i]]*beta198[2] x199[i]<-beta199[1]+transparency[wdinumber[i],year[i]]*beta199[2] x200[i]<-beta200[1]+transparency[wdinumber[i],year[i]]*beta200[2] x201[i]<-beta201[1]+transparency[wdinumber[i],year[i]]*beta201[2] x202[i]<-beta202[1]+transparency[wdinumber[i],year[i]]*beta202[2] x203[i]<-beta203[1]+transparency[wdinumber[i],year[i]]*beta203[2] x204[i]<-beta204[1]+transparency[wdinumber[i],year[i]]*beta204[2] x205[i]<-beta205[1]+transparency[wdinumber[i],year[i]]*beta205[2] x206[i]<-beta206[1]+transparency[wdinumber[i],year[i]]*beta206[2] x207[i]<-beta207[1]+transparency[wdinumber[i],year[i]]*beta207[2] x208[i]<-beta208[1]+transparency[wdinumber[i],year[i]]*beta208[2] x209[i]<-beta209[1]+transparency[wdinumber[i],year[i]]*beta209[2] x210[i]<-beta210[1]+transparency[wdinumber[i],year[i]]*beta210[2] x211[i]<-beta211[1]+transparency[wdinumber[i],year[i]]*beta211[2] x212[i]<-beta212[1]+transparency[wdinumber[i],year[i]]*beta212[2] x213[i]<-beta213[1]+transparency[wdinumber[i],year[i]]*beta213[2] x214[i]<-beta214[1]+transparency[wdinumber[i],year[i]]*beta214[2] x215[i]<-beta215[1]+transparency[wdinumber[i],year[i]]*beta215[2] x216[i]<-beta216[1]+transparency[wdinumber[i],year[i]]*beta216[2] x217[i]<-beta217[1]+transparency[wdinumber[i],year[i]]*beta217[2] x218[i]<-beta218[1]+transparency[wdinumber[i],year[i]]*beta218[2] x219[i]<-beta219[1]+transparency[wdinumber[i],year[i]]*beta219[2] x220[i]<-beta220[1]+transparency[wdinumber[i],year[i]]*beta220[2] x221[i]<-beta221[1]+transparency[wdinumber[i],year[i]]*beta221[2] x222[i]<-beta222[1]+transparency[wdinumber[i],year[i]]*beta222[2] x223[i]<-beta223[1]+transparency[wdinumber[i],year[i]]*beta223[2] x224[i]<-beta224[1]+transparency[wdinumber[i],year[i]]*beta224[2] x225[i]<-beta225[1]+transparency[wdinumber[i],year[i]]*beta225[2] x226[i]<-beta226[1]+transparency[wdinumber[i],year[i]]*beta226[2] x227[i]<-beta227[1]+transparency[wdinumber[i],year[i]]*beta227[2] x228[i]<-beta228[1]+transparency[wdinumber[i],year[i]]*beta228[2] x229[i]<-beta229[1]+transparency[wdinumber[i],year[i]]*beta229[2] x230[i]<-beta230[1]+transparency[wdinumber[i],year[i]]*beta230[2] x231[i]<-beta231[1]+transparency[wdinumber[i],year[i]]*beta231[2] x232[i]<-beta232[1]+transparency[wdinumber[i],year[i]]*beta232[2] x233[i]<-beta233[1]+transparency[wdinumber[i],year[i]]*beta233[2] x234[i]<-beta234[1]+transparency[wdinumber[i],year[i]]*beta234[2] x235[i]<-beta235[1]+transparency[wdinumber[i],year[i]]*beta235[2] x236[i]<-beta236[1]+transparency[wdinumber[i],year[i]]*beta236[2] x237[i]<-beta237[1]+transparency[wdinumber[i],year[i]]*beta237[2] x238[i]<-beta238[1]+transparency[wdinumber[i],year[i]]*beta238[2] x239[i]<-beta239[1]+transparency[wdinumber[i],year[i]]*beta239[2] x240[i]<-beta240[1]+transparency[wdinumber[i],year[i]]*beta240[2] p1[i]<-1/(1+exp(-x1[i])) p2[i]<-1/(1+exp(-x2[i])) p3[i]<-1/(1+exp(-x3[i])) p4[i]<-1/(1+exp(-x4[i])) p5[i]<-1/(1+exp(-x5[i])) p6[i]<-1/(1+exp(-x6[i])) p7[i]<-1/(1+exp(-x7[i])) p8[i]<-1/(1+exp(-x8[i])) p9[i]<-1/(1+exp(-x9[i])) p10[i]<-1/(1+exp(-x10[i])) p11[i]<-1/(1+exp(-x11[i])) p12[i]<-1/(1+exp(-x12[i])) p13[i]<-1/(1+exp(-x13[i])) p14[i]<-1/(1+exp(-x14[i])) p15[i]<-1/(1+exp(-x15[i])) p16[i]<-1/(1+exp(-x16[i])) p17[i]<-1/(1+exp(-x17[i])) p18[i]<-1/(1+exp(-x18[i])) p19[i]<-1/(1+exp(-x19[i])) p20[i]<-1/(1+exp(-x20[i])) p21[i]<-1/(1+exp(-x21[i])) p22[i]<-1/(1+exp(-x22[i])) p23[i]<-1/(1+exp(-x23[i])) p24[i]<-1/(1+exp(-x24[i])) p25[i]<-1/(1+exp(-x25[i])) p26[i]<-1/(1+exp(-x26[i])) p27[i]<-1/(1+exp(-x27[i])) p28[i]<-1/(1+exp(-x28[i])) p29[i]<-1/(1+exp(-x29[i])) p30[i]<-1/(1+exp(-x30[i])) p31[i]<-1/(1+exp(-x31[i])) p32[i]<-1/(1+exp(-x32[i])) p33[i]<-1/(1+exp(-x33[i])) p34[i]<-1/(1+exp(-x34[i])) p35[i]<-1/(1+exp(-x35[i])) p36[i]<-1/(1+exp(-x36[i])) p37[i]<-1/(1+exp(-x37[i])) p38[i]<-1/(1+exp(-x38[i])) p39[i]<-1/(1+exp(-x39[i])) p40[i]<-1/(1+exp(-x40[i])) p41[i]<-1/(1+exp(-x41[i])) p42[i]<-1/(1+exp(-x42[i])) p43[i]<-1/(1+exp(-x43[i])) p44[i]<-1/(1+exp(-x44[i])) p45[i]<-1/(1+exp(-x45[i])) p46[i]<-1/(1+exp(-x46[i])) p47[i]<-1/(1+exp(-x47[i])) p48[i]<-1/(1+exp(-x48[i])) p49[i]<-1/(1+exp(-x49[i])) p50[i]<-1/(1+exp(-x50[i])) p51[i]<-1/(1+exp(-x51[i])) p52[i]<-1/(1+exp(-x52[i])) p53[i]<-1/(1+exp(-x53[i])) p54[i]<-1/(1+exp(-x54[i])) p55[i]<-1/(1+exp(-x55[i])) p56[i]<-1/(1+exp(-x56[i])) p57[i]<-1/(1+exp(-x57[i])) p58[i]<-1/(1+exp(-x58[i])) p59[i]<-1/(1+exp(-x59[i])) p60[i]<-1/(1+exp(-x60[i])) p61[i]<-1/(1+exp(-x61[i])) p62[i]<-1/(1+exp(-x62[i])) p63[i]<-1/(1+exp(-x63[i])) p64[i]<-1/(1+exp(-x64[i])) p65[i]<-1/(1+exp(-x65[i])) p66[i]<-1/(1+exp(-x66[i])) p67[i]<-1/(1+exp(-x67[i])) p68[i]<-1/(1+exp(-x68[i])) p69[i]<-1/(1+exp(-x69[i])) p70[i]<-1/(1+exp(-x70[i])) p71[i]<-1/(1+exp(-x71[i])) p72[i]<-1/(1+exp(-x72[i])) p73[i]<-1/(1+exp(-x73[i])) p74[i]<-1/(1+exp(-x74[i])) p75[i]<-1/(1+exp(-x75[i])) p76[i]<-1/(1+exp(-x76[i])) p77[i]<-1/(1+exp(-x77[i])) p78[i]<-1/(1+exp(-x78[i])) p79[i]<-1/(1+exp(-x79[i])) p80[i]<-1/(1+exp(-x80[i])) p81[i]<-1/(1+exp(-x81[i])) p82[i]<-1/(1+exp(-x82[i])) p83[i]<-1/(1+exp(-x83[i])) p84[i]<-1/(1+exp(-x84[i])) p85[i]<-1/(1+exp(-x85[i])) p86[i]<-1/(1+exp(-x86[i])) p87[i]<-1/(1+exp(-x87[i])) p88[i]<-1/(1+exp(-x88[i])) p89[i]<-1/(1+exp(-x89[i])) p90[i]<-1/(1+exp(-x90[i])) p91[i]<-1/(1+exp(-x91[i])) p92[i]<-1/(1+exp(-x92[i])) p93[i]<-1/(1+exp(-x93[i])) p94[i]<-1/(1+exp(-x94[i])) p95[i]<-1/(1+exp(-x95[i])) p96[i]<-1/(1+exp(-x96[i])) p97[i]<-1/(1+exp(-x97[i])) p98[i]<-1/(1+exp(-x98[i])) p99[i]<-1/(1+exp(-x99[i])) p100[i]<-1/(1+exp(-x100[i])) p101[i]<-1/(1+exp(-x101[i])) p102[i]<-1/(1+exp(-x102[i])) p103[i]<-1/(1+exp(-x103[i])) p104[i]<-1/(1+exp(-x104[i])) p105[i]<-1/(1+exp(-x105[i])) p106[i]<-1/(1+exp(-x106[i])) p107[i]<-1/(1+exp(-x107[i])) p108[i]<-1/(1+exp(-x108[i])) p109[i]<-1/(1+exp(-x109[i])) p110[i]<-1/(1+exp(-x110[i])) p111[i]<-1/(1+exp(-x111[i])) p112[i]<-1/(1+exp(-x112[i])) p113[i]<-1/(1+exp(-x113[i])) p114[i]<-1/(1+exp(-x114[i])) p115[i]<-1/(1+exp(-x115[i])) p116[i]<-1/(1+exp(-x116[i])) p117[i]<-1/(1+exp(-x117[i])) p118[i]<-1/(1+exp(-x118[i])) p119[i]<-1/(1+exp(-x119[i])) p120[i]<-1/(1+exp(-x120[i])) p121[i]<-1/(1+exp(-x121[i])) p122[i]<-1/(1+exp(-x122[i])) p123[i]<-1/(1+exp(-x123[i])) p124[i]<-1/(1+exp(-x124[i])) p125[i]<-1/(1+exp(-x125[i])) p126[i]<-1/(1+exp(-x126[i])) p127[i]<-1/(1+exp(-x127[i])) p128[i]<-1/(1+exp(-x128[i])) p129[i]<-1/(1+exp(-x129[i])) p130[i]<-1/(1+exp(-x130[i])) p131[i]<-1/(1+exp(-x131[i])) p132[i]<-1/(1+exp(-x132[i])) p133[i]<-1/(1+exp(-x133[i])) p134[i]<-1/(1+exp(-x134[i])) p135[i]<-1/(1+exp(-x135[i])) p136[i]<-1/(1+exp(-x136[i])) p137[i]<-1/(1+exp(-x137[i])) p138[i]<-1/(1+exp(-x138[i])) p139[i]<-1/(1+exp(-x139[i])) p140[i]<-1/(1+exp(-x140[i])) p141[i]<-1/(1+exp(-x141[i])) p142[i]<-1/(1+exp(-x142[i])) p143[i]<-1/(1+exp(-x143[i])) p144[i]<-1/(1+exp(-x144[i])) p145[i]<-1/(1+exp(-x145[i])) p146[i]<-1/(1+exp(-x146[i])) p147[i]<-1/(1+exp(-x147[i])) p148[i]<-1/(1+exp(-x148[i])) p149[i]<-1/(1+exp(-x149[i])) p150[i]<-1/(1+exp(-x150[i])) p151[i]<-1/(1+exp(-x151[i])) p152[i]<-1/(1+exp(-x152[i])) p153[i]<-1/(1+exp(-x153[i])) p154[i]<-1/(1+exp(-x154[i])) p155[i]<-1/(1+exp(-x155[i])) p156[i]<-1/(1+exp(-x156[i])) p157[i]<-1/(1+exp(-x157[i])) p158[i]<-1/(1+exp(-x158[i])) p159[i]<-1/(1+exp(-x159[i])) p160[i]<-1/(1+exp(-x160[i])) p161[i]<-1/(1+exp(-x161[i])) p162[i]<-1/(1+exp(-x162[i])) p163[i]<-1/(1+exp(-x163[i])) p164[i]<-1/(1+exp(-x164[i])) p165[i]<-1/(1+exp(-x165[i])) p166[i]<-1/(1+exp(-x166[i])) p167[i]<-1/(1+exp(-x167[i])) p168[i]<-1/(1+exp(-x168[i])) p169[i]<-1/(1+exp(-x169[i])) p170[i]<-1/(1+exp(-x170[i])) p171[i]<-1/(1+exp(-x171[i])) p172[i]<-1/(1+exp(-x172[i])) p173[i]<-1/(1+exp(-x173[i])) p174[i]<-1/(1+exp(-x174[i])) p175[i]<-1/(1+exp(-x175[i])) p176[i]<-1/(1+exp(-x176[i])) p177[i]<-1/(1+exp(-x177[i])) p178[i]<-1/(1+exp(-x178[i])) p179[i]<-1/(1+exp(-x179[i])) p180[i]<-1/(1+exp(-x180[i])) p181[i]<-1/(1+exp(-x181[i])) p182[i]<-1/(1+exp(-x182[i])) p183[i]<-1/(1+exp(-x183[i])) p184[i]<-1/(1+exp(-x184[i])) p185[i]<-1/(1+exp(-x185[i])) p186[i]<-1/(1+exp(-x186[i])) p187[i]<-1/(1+exp(-x187[i])) p188[i]<-1/(1+exp(-x188[i])) p189[i]<-1/(1+exp(-x189[i])) p190[i]<-1/(1+exp(-x190[i])) p191[i]<-1/(1+exp(-x191[i])) p192[i]<-1/(1+exp(-x192[i])) p193[i]<-1/(1+exp(-x193[i])) p194[i]<-1/(1+exp(-x194[i])) p195[i]<-1/(1+exp(-x195[i])) p196[i]<-1/(1+exp(-x196[i])) p197[i]<-1/(1+exp(-x197[i])) p198[i]<-1/(1+exp(-x198[i])) p199[i]<-1/(1+exp(-x199[i])) p200[i]<-1/(1+exp(-x200[i])) p201[i]<-1/(1+exp(-x201[i])) p202[i]<-1/(1+exp(-x202[i])) p203[i]<-1/(1+exp(-x203[i])) p204[i]<-1/(1+exp(-x204[i])) p205[i]<-1/(1+exp(-x205[i])) p206[i]<-1/(1+exp(-x206[i])) p207[i]<-1/(1+exp(-x207[i])) p208[i]<-1/(1+exp(-x208[i])) p209[i]<-1/(1+exp(-x209[i])) p210[i]<-1/(1+exp(-x210[i])) p211[i]<-1/(1+exp(-x211[i])) p212[i]<-1/(1+exp(-x212[i])) p213[i]<-1/(1+exp(-x213[i])) p214[i]<-1/(1+exp(-x214[i])) p215[i]<-1/(1+exp(-x215[i])) p216[i]<-1/(1+exp(-x216[i])) p217[i]<-1/(1+exp(-x217[i])) p218[i]<-1/(1+exp(-x218[i])) p219[i]<-1/(1+exp(-x219[i])) p220[i]<-1/(1+exp(-x220[i])) p221[i]<-1/(1+exp(-x221[i])) p222[i]<-1/(1+exp(-x222[i])) p223[i]<-1/(1+exp(-x223[i])) p224[i]<-1/(1+exp(-x224[i])) p225[i]<-1/(1+exp(-x225[i])) p226[i]<-1/(1+exp(-x226[i])) p227[i]<-1/(1+exp(-x227[i])) p228[i]<-1/(1+exp(-x228[i])) p229[i]<-1/(1+exp(-x229[i])) p230[i]<-1/(1+exp(-x230[i])) p231[i]<-1/(1+exp(-x231[i])) p232[i]<-1/(1+exp(-x232[i])) p233[i]<-1/(1+exp(-x233[i])) p234[i]<-1/(1+exp(-x234[i])) p235[i]<-1/(1+exp(-x235[i])) p236[i]<-1/(1+exp(-x236[i])) p237[i]<-1/(1+exp(-x237[i])) p238[i]<-1/(1+exp(-x238[i])) p239[i]<-1/(1+exp(-x239[i])) p240[i]<-1/(1+exp(-x240[i])) q1[i]<-max(min(p1[i],1),0) q2[i]<-max(min(p2[i],1),0) q3[i]<-max(min(p3[i],1),0) q4[i]<-max(min(p4[i],1),0) q5[i]<-max(min(p5[i],1),0) q6[i]<-max(min(p6[i],1),0) q7[i]<-max(min(p7[i],1),0) q8[i]<-max(min(p8[i],1),0) q9[i]<-max(min(p9[i],1),0) q10[i]<-max(min(p10[i],1),0) q11[i]<-max(min(p11[i],1),0) q12[i]<-max(min(p12[i],1),0) q13[i]<-max(min(p13[i],1),0) q14[i]<-max(min(p14[i],1),0) q15[i]<-max(min(p15[i],1),0) q16[i]<-max(min(p16[i],1),0) q17[i]<-max(min(p17[i],1),0) q18[i]<-max(min(p18[i],1),0) q19[i]<-max(min(p19[i],1),0) q20[i]<-max(min(p20[i],1),0) q21[i]<-max(min(p21[i],1),0) q22[i]<-max(min(p22[i],1),0) q23[i]<-max(min(p23[i],1),0) q24[i]<-max(min(p24[i],1),0) q25[i]<-max(min(p25[i],1),0) q26[i]<-max(min(p26[i],1),0) q27[i]<-max(min(p27[i],1),0) q28[i]<-max(min(p28[i],1),0) q29[i]<-max(min(p29[i],1),0) q30[i]<-max(min(p30[i],1),0) q31[i]<-max(min(p31[i],1),0) q32[i]<-max(min(p32[i],1),0) q33[i]<-max(min(p33[i],1),0) q34[i]<-max(min(p34[i],1),0) q35[i]<-max(min(p35[i],1),0) q36[i]<-max(min(p36[i],1),0) q37[i]<-max(min(p37[i],1),0) q38[i]<-max(min(p38[i],1),0) q39[i]<-max(min(p39[i],1),0) q40[i]<-max(min(p40[i],1),0) q41[i]<-max(min(p41[i],1),0) q42[i]<-max(min(p42[i],1),0) q43[i]<-max(min(p43[i],1),0) q44[i]<-max(min(p44[i],1),0) q45[i]<-max(min(p45[i],1),0) q46[i]<-max(min(p46[i],1),0) q47[i]<-max(min(p47[i],1),0) q48[i]<-max(min(p48[i],1),0) q49[i]<-max(min(p49[i],1),0) q50[i]<-max(min(p50[i],1),0) q51[i]<-max(min(p51[i],1),0) q52[i]<-max(min(p52[i],1),0) q53[i]<-max(min(p53[i],1),0) q54[i]<-max(min(p54[i],1),0) q55[i]<-max(min(p55[i],1),0) q56[i]<-max(min(p56[i],1),0) q57[i]<-max(min(p57[i],1),0) q58[i]<-max(min(p58[i],1),0) q59[i]<-max(min(p59[i],1),0) q60[i]<-max(min(p60[i],1),0) q61[i]<-max(min(p61[i],1),0) q62[i]<-max(min(p62[i],1),0) q63[i]<-max(min(p63[i],1),0) q64[i]<-max(min(p64[i],1),0) q65[i]<-max(min(p65[i],1),0) q66[i]<-max(min(p66[i],1),0) q67[i]<-max(min(p67[i],1),0) q68[i]<-max(min(p68[i],1),0) q69[i]<-max(min(p69[i],1),0) q70[i]<-max(min(p70[i],1),0) q71[i]<-max(min(p71[i],1),0) q72[i]<-max(min(p72[i],1),0) q73[i]<-max(min(p73[i],1),0) q74[i]<-max(min(p74[i],1),0) q75[i]<-max(min(p75[i],1),0) q76[i]<-max(min(p76[i],1),0) q77[i]<-max(min(p77[i],1),0) q78[i]<-max(min(p78[i],1),0) q79[i]<-max(min(p79[i],1),0) q80[i]<-max(min(p80[i],1),0) q81[i]<-max(min(p81[i],1),0) q82[i]<-max(min(p82[i],1),0) q83[i]<-max(min(p83[i],1),0) q84[i]<-max(min(p84[i],1),0) q85[i]<-max(min(p85[i],1),0) q86[i]<-max(min(p86[i],1),0) q87[i]<-max(min(p87[i],1),0) q88[i]<-max(min(p88[i],1),0) q89[i]<-max(min(p89[i],1),0) q90[i]<-max(min(p90[i],1),0) q91[i]<-max(min(p91[i],1),0) q92[i]<-max(min(p92[i],1),0) q93[i]<-max(min(p93[i],1),0) q94[i]<-max(min(p94[i],1),0) q95[i]<-max(min(p95[i],1),0) q96[i]<-max(min(p96[i],1),0) q97[i]<-max(min(p97[i],1),0) q98[i]<-max(min(p98[i],1),0) q99[i]<-max(min(p99[i],1),0) q100[i]<-max(min(p100[i],1),0) q101[i]<-max(min(p101[i],1),0) q102[i]<-max(min(p102[i],1),0) q103[i]<-max(min(p103[i],1),0) q104[i]<-max(min(p104[i],1),0) q105[i]<-max(min(p105[i],1),0) q106[i]<-max(min(p106[i],1),0) q107[i]<-max(min(p107[i],1),0) q108[i]<-max(min(p108[i],1),0) q109[i]<-max(min(p109[i],1),0) q110[i]<-max(min(p110[i],1),0) q111[i]<-max(min(p111[i],1),0) q112[i]<-max(min(p112[i],1),0) q113[i]<-max(min(p113[i],1),0) q114[i]<-max(min(p114[i],1),0) q115[i]<-max(min(p115[i],1),0) q116[i]<-max(min(p116[i],1),0) q117[i]<-max(min(p117[i],1),0) q118[i]<-max(min(p118[i],1),0) q119[i]<-max(min(p119[i],1),0) q120[i]<-max(min(p120[i],1),0) q121[i]<-max(min(p121[i],1),0) q122[i]<-max(min(p122[i],1),0) q123[i]<-max(min(p123[i],1),0) q124[i]<-max(min(p124[i],1),0) q125[i]<-max(min(p125[i],1),0) q126[i]<-max(min(p126[i],1),0) q127[i]<-max(min(p127[i],1),0) q128[i]<-max(min(p128[i],1),0) q129[i]<-max(min(p129[i],1),0) q130[i]<-max(min(p130[i],1),0) q131[i]<-max(min(p131[i],1),0) q132[i]<-max(min(p132[i],1),0) q133[i]<-max(min(p133[i],1),0) q134[i]<-max(min(p134[i],1),0) q135[i]<-max(min(p135[i],1),0) q136[i]<-max(min(p136[i],1),0) q137[i]<-max(min(p137[i],1),0) q138[i]<-max(min(p138[i],1),0) q139[i]<-max(min(p139[i],1),0) q140[i]<-max(min(p140[i],1),0) q141[i]<-max(min(p141[i],1),0) q142[i]<-max(min(p142[i],1),0) q143[i]<-max(min(p143[i],1),0) q144[i]<-max(min(p144[i],1),0) q145[i]<-max(min(p145[i],1),0) q146[i]<-max(min(p146[i],1),0) q147[i]<-max(min(p147[i],1),0) q148[i]<-max(min(p148[i],1),0) q149[i]<-max(min(p149[i],1),0) q150[i]<-max(min(p150[i],1),0) q151[i]<-max(min(p151[i],1),0) q152[i]<-max(min(p152[i],1),0) q153[i]<-max(min(p153[i],1),0) q154[i]<-max(min(p154[i],1),0) q155[i]<-max(min(p155[i],1),0) q156[i]<-max(min(p156[i],1),0) q157[i]<-max(min(p157[i],1),0) q158[i]<-max(min(p158[i],1),0) q159[i]<-max(min(p159[i],1),0) q160[i]<-max(min(p160[i],1),0) q161[i]<-max(min(p161[i],1),0) q162[i]<-max(min(p162[i],1),0) q163[i]<-max(min(p163[i],1),0) q164[i]<-max(min(p164[i],1),0) q165[i]<-max(min(p165[i],1),0) q166[i]<-max(min(p166[i],1),0) q167[i]<-max(min(p167[i],1),0) q168[i]<-max(min(p168[i],1),0) q169[i]<-max(min(p169[i],1),0) q170[i]<-max(min(p170[i],1),0) q171[i]<-max(min(p171[i],1),0) q172[i]<-max(min(p172[i],1),0) q173[i]<-max(min(p173[i],1),0) q174[i]<-max(min(p174[i],1),0) q175[i]<-max(min(p175[i],1),0) q176[i]<-max(min(p176[i],1),0) q177[i]<-max(min(p177[i],1),0) q178[i]<-max(min(p178[i],1),0) q179[i]<-max(min(p179[i],1),0) q180[i]<-max(min(p180[i],1),0) q181[i]<-max(min(p181[i],1),0) q182[i]<-max(min(p182[i],1),0) q183[i]<-max(min(p183[i],1),0) q184[i]<-max(min(p184[i],1),0) q185[i]<-max(min(p185[i],1),0) q186[i]<-max(min(p186[i],1),0) q187[i]<-max(min(p187[i],1),0) q188[i]<-max(min(p188[i],1),0) q189[i]<-max(min(p189[i],1),0) q190[i]<-max(min(p190[i],1),0) q191[i]<-max(min(p191[i],1),0) q192[i]<-max(min(p192[i],1),0) q193[i]<-max(min(p193[i],1),0) q194[i]<-max(min(p194[i],1),0) q195[i]<-max(min(p195[i],1),0) q196[i]<-max(min(p196[i],1),0) q197[i]<-max(min(p197[i],1),0) q198[i]<-max(min(p198[i],1),0) q199[i]<-max(min(p199[i],1),0) q200[i]<-max(min(p200[i],1),0) q201[i]<-max(min(p201[i],1),0) q202[i]<-max(min(p202[i],1),0) q203[i]<-max(min(p203[i],1),0) q204[i]<-max(min(p204[i],1),0) q205[i]<-max(min(p205[i],1),0) q206[i]<-max(min(p206[i],1),0) q207[i]<-max(min(p207[i],1),0) q208[i]<-max(min(p208[i],1),0) q209[i]<-max(min(p209[i],1),0) q210[i]<-max(min(p210[i],1),0) q211[i]<-max(min(p211[i],1),0) q212[i]<-max(min(p212[i],1),0) q213[i]<-max(min(p213[i],1),0) q214[i]<-max(min(p214[i],1),0) q215[i]<-max(min(p215[i],1),0) q216[i]<-max(min(p216[i],1),0) q217[i]<-max(min(p217[i],1),0) q218[i]<-max(min(p218[i],1),0) q219[i]<-max(min(p219[i],1),0) q220[i]<-max(min(p220[i],1),0) q221[i]<-max(min(p221[i],1),0) q222[i]<-max(min(p222[i],1),0) q223[i]<-max(min(p223[i],1),0) q224[i]<-max(min(p224[i],1),0) q225[i]<-max(min(p225[i],1),0) q226[i]<-max(min(p226[i],1),0) q227[i]<-max(min(p227[i],1),0) q228[i]<-max(min(p228[i],1),0) q229[i]<-max(min(p229[i],1),0) q230[i]<-max(min(p230[i],1),0) q231[i]<-max(min(p231[i],1),0) q232[i]<-max(min(p232[i],1),0) q233[i]<-max(min(p233[i],1),0) q234[i]<-max(min(p234[i],1),0) q235[i]<-max(min(p235[i],1),0) q236[i]<-max(min(p236[i],1),0) q237[i]<-max(min(p237[i],1),0) q238[i]<-max(min(p238[i],1),0) q239[i]<-max(min(p239[i],1),0) q240[i]<-max(min(p240[i],1),0) var_1[i]~dbern(q1[i]) var_2[i]~dbern(q2[i]) var_3[i]~dbern(q3[i]) var_4[i]~dbern(q4[i]) var_5[i]~dbern(q5[i]) var_6[i]~dbern(q6[i]) var_7[i]~dbern(q7[i]) var_8[i]~dbern(q8[i]) var_9[i]~dbern(q9[i]) var_10[i]~dbern(q10[i]) var_11[i]~dbern(q11[i]) var_12[i]~dbern(q12[i]) var_13[i]~dbern(q13[i]) var_14[i]~dbern(q14[i]) var_15[i]~dbern(q15[i]) var_16[i]~dbern(q16[i]) var_17[i]~dbern(q17[i]) var_18[i]~dbern(q18[i]) var_19[i]~dbern(q19[i]) var_20[i]~dbern(q20[i]) var_21[i]~dbern(q21[i]) var_22[i]~dbern(q22[i]) var_23[i]~dbern(q23[i]) var_24[i]~dbern(q24[i]) var_25[i]~dbern(q25[i]) var_26[i]~dbern(q26[i]) var_27[i]~dbern(q27[i]) var_28[i]~dbern(q28[i]) var_29[i]~dbern(q29[i]) var_30[i]~dbern(q30[i]) var_31[i]~dbern(q31[i]) var_32[i]~dbern(q32[i]) var_33[i]~dbern(q33[i]) var_34[i]~dbern(q34[i]) var_35[i]~dbern(q35[i]) var_36[i]~dbern(q36[i]) var_37[i]~dbern(q37[i]) var_38[i]~dbern(q38[i]) var_39[i]~dbern(q39[i]) var_40[i]~dbern(q40[i]) var_41[i]~dbern(q41[i]) var_42[i]~dbern(q42[i]) var_43[i]~dbern(q43[i]) var_44[i]~dbern(q44[i]) var_45[i]~dbern(q45[i]) var_46[i]~dbern(q46[i]) var_47[i]~dbern(q47[i]) var_48[i]~dbern(q48[i]) var_49[i]~dbern(q49[i]) var_50[i]~dbern(q50[i]) var_51[i]~dbern(q51[i]) var_52[i]~dbern(q52[i]) var_53[i]~dbern(q53[i]) var_54[i]~dbern(q54[i]) var_55[i]~dbern(q55[i]) var_56[i]~dbern(q56[i]) var_57[i]~dbern(q57[i]) var_58[i]~dbern(q58[i]) var_59[i]~dbern(q59[i]) var_60[i]~dbern(q60[i]) var_61[i]~dbern(q61[i]) var_62[i]~dbern(q62[i]) var_63[i]~dbern(q63[i]) var_64[i]~dbern(q64[i]) var_65[i]~dbern(q65[i]) var_66[i]~dbern(q66[i]) var_67[i]~dbern(q67[i]) var_68[i]~dbern(q68[i]) var_69[i]~dbern(q69[i]) var_70[i]~dbern(q70[i]) var_71[i]~dbern(q71[i]) var_72[i]~dbern(q72[i]) var_73[i]~dbern(q73[i]) var_74[i]~dbern(q74[i]) var_75[i]~dbern(q75[i]) var_76[i]~dbern(q76[i]) var_77[i]~dbern(q77[i]) var_78[i]~dbern(q78[i]) var_79[i]~dbern(q79[i]) var_80[i]~dbern(q80[i]) var_81[i]~dbern(q81[i]) var_82[i]~dbern(q82[i]) var_83[i]~dbern(q83[i]) var_84[i]~dbern(q84[i]) var_85[i]~dbern(q85[i]) var_86[i]~dbern(q86[i]) var_87[i]~dbern(q87[i]) var_88[i]~dbern(q88[i]) var_89[i]~dbern(q89[i]) var_90[i]~dbern(q90[i]) var_91[i]~dbern(q91[i]) var_92[i]~dbern(q92[i]) var_93[i]~dbern(q93[i]) var_94[i]~dbern(q94[i]) var_95[i]~dbern(q95[i]) var_96[i]~dbern(q96[i]) var_97[i]~dbern(q97[i]) var_98[i]~dbern(q98[i]) var_99[i]~dbern(q99[i]) var_100[i]~dbern(q100[i]) var_101[i]~dbern(q101[i]) var_102[i]~dbern(q102[i]) var_103[i]~dbern(q103[i]) var_104[i]~dbern(q104[i]) var_105[i]~dbern(q105[i]) var_106[i]~dbern(q106[i]) var_107[i]~dbern(q107[i]) var_108[i]~dbern(q108[i]) var_109[i]~dbern(q109[i]) var_110[i]~dbern(q110[i]) var_111[i]~dbern(q111[i]) var_112[i]~dbern(q112[i]) var_113[i]~dbern(q113[i]) var_114[i]~dbern(q114[i]) var_115[i]~dbern(q115[i]) var_116[i]~dbern(q116[i]) var_117[i]~dbern(q117[i]) var_118[i]~dbern(q118[i]) var_119[i]~dbern(q119[i]) var_120[i]~dbern(q120[i]) var_121[i]~dbern(q121[i]) var_122[i]~dbern(q122[i]) var_123[i]~dbern(q123[i]) var_124[i]~dbern(q124[i]) var_125[i]~dbern(q125[i]) var_126[i]~dbern(q126[i]) var_127[i]~dbern(q127[i]) var_128[i]~dbern(q128[i]) var_129[i]~dbern(q129[i]) var_130[i]~dbern(q130[i]) var_131[i]~dbern(q131[i]) var_132[i]~dbern(q132[i]) var_133[i]~dbern(q133[i]) var_134[i]~dbern(q134[i]) var_135[i]~dbern(q135[i]) var_136[i]~dbern(q136[i]) var_137[i]~dbern(q137[i]) var_138[i]~dbern(q138[i]) var_139[i]~dbern(q139[i]) var_140[i]~dbern(q140[i]) var_141[i]~dbern(q141[i]) var_142[i]~dbern(q142[i]) var_143[i]~dbern(q143[i]) var_144[i]~dbern(q144[i]) var_145[i]~dbern(q145[i]) var_146[i]~dbern(q146[i]) var_147[i]~dbern(q147[i]) var_148[i]~dbern(q148[i]) var_149[i]~dbern(q149[i]) var_150[i]~dbern(q150[i]) var_151[i]~dbern(q151[i]) var_152[i]~dbern(q152[i]) var_153[i]~dbern(q153[i]) var_154[i]~dbern(q154[i]) var_155[i]~dbern(q155[i]) var_156[i]~dbern(q156[i]) var_157[i]~dbern(q157[i]) var_158[i]~dbern(q158[i]) var_159[i]~dbern(q159[i]) var_160[i]~dbern(q160[i]) var_161[i]~dbern(q161[i]) var_162[i]~dbern(q162[i]) var_163[i]~dbern(q163[i]) var_164[i]~dbern(q164[i]) var_165[i]~dbern(q165[i]) var_166[i]~dbern(q166[i]) var_167[i]~dbern(q167[i]) var_168[i]~dbern(q168[i]) var_169[i]~dbern(q169[i]) var_170[i]~dbern(q170[i]) var_171[i]~dbern(q171[i]) var_172[i]~dbern(q172[i]) var_173[i]~dbern(q173[i]) var_174[i]~dbern(q174[i]) var_175[i]~dbern(q175[i]) var_176[i]~dbern(q176[i]) var_177[i]~dbern(q177[i]) var_178[i]~dbern(q178[i]) var_179[i]~dbern(q179[i]) var_180[i]~dbern(q180[i]) var_181[i]~dbern(q181[i]) var_182[i]~dbern(q182[i]) var_183[i]~dbern(q183[i]) var_184[i]~dbern(q184[i]) var_185[i]~dbern(q185[i]) var_186[i]~dbern(q186[i]) var_187[i]~dbern(q187[i]) var_188[i]~dbern(q188[i]) var_189[i]~dbern(q189[i]) var_190[i]~dbern(q190[i]) var_191[i]~dbern(q191[i]) var_192[i]~dbern(q192[i]) var_193[i]~dbern(q193[i]) var_194[i]~dbern(q194[i]) var_195[i]~dbern(q195[i]) var_196[i]~dbern(q196[i]) var_197[i]~dbern(q197[i]) var_198[i]~dbern(q198[i]) var_199[i]~dbern(q199[i]) var_200[i]~dbern(q200[i]) var_201[i]~dbern(q201[i]) var_202[i]~dbern(q202[i]) var_203[i]~dbern(q203[i]) var_204[i]~dbern(q204[i]) var_205[i]~dbern(q205[i]) var_206[i]~dbern(q206[i]) var_207[i]~dbern(q207[i]) var_208[i]~dbern(q208[i]) var_209[i]~dbern(q209[i]) var_210[i]~dbern(q210[i]) var_211[i]~dbern(q211[i]) var_212[i]~dbern(q212[i]) var_213[i]~dbern(q213[i]) var_214[i]~dbern(q214[i]) var_215[i]~dbern(q215[i]) var_216[i]~dbern(q216[i]) var_217[i]~dbern(q217[i]) var_218[i]~dbern(q218[i]) var_219[i]~dbern(q219[i]) var_220[i]~dbern(q220[i]) var_221[i]~dbern(q221[i]) var_222[i]~dbern(q222[i]) var_223[i]~dbern(q223[i]) var_224[i]~dbern(q224[i]) var_225[i]~dbern(q225[i]) var_226[i]~dbern(q226[i]) var_227[i]~dbern(q227[i]) var_228[i]~dbern(q228[i]) var_229[i]~dbern(q229[i]) var_230[i]~dbern(q230[i]) var_231[i]~dbern(q231[i]) var_232[i]~dbern(q232[i]) var_233[i]~dbern(q233[i]) var_234[i]~dbern(q234[i]) var_235[i]~dbern(q235[i]) var_236[i]~dbern(q236[i]) var_237[i]~dbern(q237[i]) var_238[i]~dbern(q238[i]) var_239[i]~dbern(q239[i]) var_240[i]~dbern(q240[i]) } #Model Priors mu[1]<-0 mu[2]<-0 alpha[1,1] <-.005 alpha[1,2]<-0 alpha[2,1]<-0 alpha[2,2] <- .005 beta1[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta2[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta3[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta4[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta5[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta6[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta7[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta8[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta9[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta10[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta11[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta12[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta13[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta14[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta15[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta16[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta17[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta18[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta19[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta20[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta21[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta22[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta23[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta24[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta25[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta26[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta27[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta28[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta29[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta30[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta31[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta32[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta33[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta34[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta35[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta36[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta37[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta38[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta39[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta40[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta41[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta42[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta43[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta44[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta45[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta46[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta47[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta48[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta49[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta50[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta51[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta52[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta53[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta54[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta55[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta56[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta57[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta58[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta59[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta60[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta61[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta62[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta63[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta64[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta65[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta66[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta67[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta68[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta69[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta70[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta71[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta72[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta73[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta74[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta75[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta76[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta77[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta78[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta79[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta80[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta81[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta82[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta83[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta84[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta85[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta86[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta87[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta88[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta89[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta90[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta91[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta92[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta93[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta94[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta95[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta96[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta97[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta98[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta99[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta100[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta101[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta102[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta103[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta104[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta105[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta106[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta107[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta108[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta109[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta110[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta111[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta112[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta113[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta114[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta115[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta116[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta117[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta118[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta119[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta120[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta121[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta122[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta123[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta124[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta125[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta126[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta127[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta128[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta129[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta130[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta131[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta132[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta133[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta134[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta135[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta136[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta137[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta138[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta139[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta140[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta141[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta142[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta143[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta144[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta145[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta146[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta147[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta148[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta149[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta150[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta151[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta152[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta153[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta154[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta155[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta156[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta157[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta158[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta159[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta160[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta161[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta162[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta163[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta164[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta165[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta166[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta167[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta168[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta169[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta170[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta171[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta172[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta173[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta174[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta175[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta176[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta177[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta178[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta179[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta180[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta181[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta182[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta183[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta184[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta185[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta186[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta187[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta188[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta189[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta190[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta191[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta192[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta193[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta194[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta195[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta196[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta197[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta198[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta199[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta200[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta201[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta202[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta203[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta204[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta205[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta206[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta207[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta208[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta209[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta210[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta211[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta212[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta213[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta214[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta215[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta216[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta217[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta218[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta219[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta220[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta221[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta222[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta223[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta224[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta225[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta226[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta227[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta228[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta229[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta230[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta231[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta232[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta233[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta234[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta235[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta236[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta237[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta238[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta239[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) beta240[1:2] ~ dmnorm(mu[1:2],alpha[1:2,1:2]) #Transparency Priors #To avoid issues of index flipping, we constrain Cuba to have negative transparency values, #and Sweden to have positive values. #Start with diffuse priors on the initial transparency term. We recenter/scale these below. for(i in 1:26){ transinit[i] ~ dnorm(0,.01) } transinit[27] ~ dnorm(0,.01)T(,0) for(i in 28:106){ transinit[i] ~ dnorm(0,.01) } transinit[107] ~ dnorm(0,.01)T(0,) for(i in 108:125){ transinit[i] ~ dnorm(0,.01) } mean.trans<- mean(transinit[1:125]) sd.trans <- sd(transinit[1:125]) for(i in 1:125){ transcentered[i] <- (transinit[i]-mean.trans)*pow(sd.trans,-1) } for (i in 1:26){ transparency[i,1] <- transcentered[i] tau[i] ~ dgamma(20, 4) for (j in 2:31){ transparency[i,j] ~ dnorm(transparency[i,(j-1)], tau[i]) } } #Cuba transparency[27,1] <- transcentered[27] tau[27] ~ dgamma(20, 4) for (j in 2:31){ transparency[27,j] ~ dnorm(transparency[27,(j-1)], tau[27])T(,0) } for (i in 28:106){ transparency[i,1] <- transcentered[i] tau[i] ~ dgamma(20, 4) for (j in 2:31){ transparency[i,j] ~ dnorm(transparency[i,(j-1)], tau[i]) } } #Sweden transparency[107,1] <- transcentered[107] tau[107] ~ dgamma(20, 4) for (j in 2:31){ transparency[107,j] ~ dnorm(transparency[107,(j-1)], tau[107])T(0,) } for (i in 108:125){ transparency[i,1] <-transcentered[i] tau[i] ~ dgamma(20, 4) for (j in 2:31){ transparency[i,j] ~ dnorm(transparency[i,(j-1)], tau[i]) } } }