model { # offspring for (i in 1:300) { weight[i] ~ dnorm(mu[i], tau.within) mu[i] <- nu + alpha[sire[i]] + beta[dam[i]] } # sires for (i in 1:6) { alpha[i] ~ dnorm(0.0, tau.sire) } alpha[7] <- -sum(alpha[1:6]) # dams for (i in 1:52) { beta[i] ~ dnorm(0.0, tau.dam) } beta[53] <- -sum(beta[1:52]) # priors nu ~ dnorm(0, 0.001) sd.sire ~ dunif(0, 10) sd.dam ~ dunif(0, 10) sd.within ~ dunif(0, 10) # observational components v.sire <- sd.sire*sd.sire v.dam <- sd.dam*sd.dam v.within <- sd.within*sd.within # precisions tau.sire <- 1.0/v.sire tau.dam <- 1.0/v.dam tau.within <- 1.0/v.within } sire[] dam[] weight[] 1 1 99.2789343018915 1 1 92.8290811609914 1 1 92.1977611557311 1 1 96.2203721102781 1 1 97.435349690979 1 1 99.5683756647001 1 1 92.2818423839868 1 1 99.5110665900414 1 2 103.230303160319 1 2 112.099498381456 1 2 100.050051149295 1 2 109.136852501619 1 2 108.070188557951 1 2 107.3213363051 1 2 104.369222048311 1 2 107.287163679546 1 2 108.116646908932 1 2 101.277201319335 1 3 97.0380991400381 1 3 98.9362261287801 1 4 95.3038053010358 1 4 99.9383310752465 1 4 95.058437743349 1 4 95.1358805398451 1 4 99.6866322402375 1 5 103.418074372728 1 5 99.2689027878344 1 5 106.152993333834 1 5 95.7637834454326 1 5 97.6875907623822 1 6 102.398105247554 1 6 104.010146654669 1 6 107.254733184839 1 6 104.218780012531 1 6 95.1780637145864 1 7 104.744351967969 1 7 111.426480907402 1 7 100.278024116546 1 7 101.113146204984 1 7 100.174552488987 1 8 96.3018244826266 1 8 96.7786453150454 1 8 102.040224207027 1 8 101.906272140626 1 8 104.79517923878 2 9 103.84432153948 2 9 102.52311847906 2 9 107.625406267505 2 9 96.1341774551378 2 9 105.619228120084 2 9 104.362308917399 2 9 106.065105278472 2 9 102.92139354849 2 10 102.739350909711 2 10 94.8255641175595 2 10 100.659017898143 2 10 99.4375570029191 2 10 102.370162026357 2 10 101.908121938339 2 10 104.902784232212 2 10 98.984289977038 2 10 100.625107830063 2 10 97.8564461178563 2 11 99.4310294946206 2 11 99.0715624600573 2 12 99.3896916519929 2 12 99.7290173948314 2 12 101.964387224715 2 12 98.3564189100951 2 12 101.317050349436 2 13 97.1091582742227 2 13 104.014674455748 2 13 101.75423163567 2 13 96.4691397547025 2 13 98.1351329759152 2 14 101.716342426621 2 14 99.0731015591762 2 14 94.6563927805436 2 14 100.443941040031 2 14 100.607978676336 2 15 99.3799443070043 2 15 104.797609950732 2 15 98.4754851382706 2 15 103.000843268357 2 15 105.350983935119 2 16 98.0937881477792 2 16 99.1265453664283 2 16 98.9874210313283 2 16 102.444617348049 2 16 102.956939364977 3 17 108.125324944251 3 17 96.8826539766912 3 17 102.341947564029 3 17 98.6899643421196 3 17 108.70412076864 3 17 104.291580205875 3 17 106.833413413647 3 17 107.4495987968 3 18 101.478496852206 3 18 108.528743257664 3 18 104.047194849274 3 18 105.453676585458 3 18 108.484435736869 3 18 104.737687660416 3 18 99.2345321936847 3 18 100.189460607996 3 18 106.519019447624 3 18 104.272713435934 3 19 101.073700651745 3 19 107.624652453828 3 20 107.506675626372 3 20 101.088910708414 3 20 103.921678974506 3 20 102.059539395336 3 20 102.444598501587 3 21 94.8315941832359 3 21 99.4873189299915 3 21 103.319300835455 3 21 102.578671869247 3 21 97.9551676467937 3 22 107.570938483335 3 22 100.567802257186 3 22 107.893101738335 3 22 110.513725257452 3 22 101.15033538097 3 23 98.2792537924256 3 23 99.1629587001827 3 23 97.040554281022 3 23 90.4384802270158 3 23 97.3491381798554 3 24 105.217895708239 3 24 98.4599701404825 3 24 101.863478931962 3 24 102.982728544086 3 24 93.8526351095429 4 25 96.3550093225691 4 25 93.509472832981 4 25 91.1207932579228 4 25 100.625590240756 4 25 95.3949238862748 4 25 92.9949158847995 4 25 95.5848646767758 4 25 97.0222377268572 4 26 95.3117676330789 4 26 91.7099314454402 4 26 97.7829225732262 4 26 96.167465051509 4 26 92.7227609711859 4 26 93.1140869477763 4 26 97.784468065118 4 26 97.7906052983482 4 26 90.0520923107166 4 26 92.1327438996688 4 27 93.3727745622639 4 27 99.4362180242547 4 28 98.976060008652 4 28 93.6522251325325 4 28 100.630012118654 4 28 93.4083347644252 4 28 95.8333799338561 4 29 95.0008516359806 4 29 96.5880147284064 4 29 91.2350406111584 4 29 89.8887067192343 4 29 93.4339244042573 4 30 98.6816835305515 4 30 91.421946537241 4 30 95.9979954034806 4 30 97.1966345927588 4 30 102.831181645 4 31 99.6920160584322 4 31 98.4908908426844 4 31 97.6716512257098 4 31 98.3828149703803 4 31 94.0827580340892 4 32 90.6715055792914 4 32 95.5849666697151 4 32 93.2403152901777 4 32 85.5469239320921 4 32 89.1232046729652 5 33 96.5467445125667 5 33 90.0868364292052 5 33 101.940346970098 5 33 99.9283685054313 5 33 102.679726350312 5 33 100.051254293833 5 33 95.4901079568733 5 33 95.2291100356104 5 34 95.6601264017848 5 34 99.8201405671877 5 34 103.16484664031 5 34 100.771976401298 5 34 100.204902256575 5 34 100.438808348453 5 34 97.8867120217271 5 34 100.56237386476 5 34 98.3381883705712 5 34 97.7489042279231 5 35 100.563439175521 5 35 101.471618758563 5 36 98.606319505759 5 36 95.5545374764659 5 36 101.809042802703 5 36 97.6539315118755 5 36 105.312767549519 5 37 100.791282585641 5 37 101.20250838184 5 37 100.323599598345 5 37 95.3631130018725 5 37 97.8194272391681 5 38 99.2558188428742 5 38 96.4212677223435 5 38 92.8946701969386 5 38 98.3876024767126 5 38 102.515328205616 5 39 97.7187510930811 5 39 105.730865760526 5 39 94.2439881427354 5 39 97.570852183436 5 39 97.7237052527154 5 40 103.244321923243 5 40 101.06065240344 5 40 99.6303474718485 5 40 103.566473786173 5 40 102.382229666422 6 41 102.783571179046 6 41 102.057525186727 6 41 100.733607807739 6 41 97.0390772625521 6 41 97.0149094820694 6 41 95.6012148782019 6 41 103.222217291171 6 41 100.617230601629 6 42 102.758652246596 6 42 94.7260771819642 6 42 97.5286123898256 6 42 102.152059585886 6 42 94.9606684835349 6 42 100.443610777302 6 42 100.139920019299 6 42 100.504304429444 6 42 102.19485854991 6 42 98.7359019184512 6 43 99.3444773728923 6 43 103.004314024007 6 44 101.594139395699 6 44 101.483301319968 6 44 105.67317409516 6 44 105.101638698918 6 44 98.4934993272372 6 45 94.4119022458766 6 45 105.258066593973 6 45 104.713106811365 6 45 105.620231310106 6 45 100.643008495523 6 46 101.818896132192 6 46 101.901826748305 6 46 103.833041477753 6 46 99.5194491514286 6 46 103.458879254412 6 47 94.5697129922525 6 47 96.2540995949465 6 47 95.1196907908487 6 47 95.0548081187196 6 47 97.3532999361522 7 48 101.931050433174 7 48 99.4102116229249 7 48 97.9044267666817 7 48 99.1060726133379 7 48 99.2159400899229 7 48 103.315963334111 7 48 102.231632021853 7 48 102.688955385621 7 49 100.532398906546 7 49 100.538339222786 7 49 103.47633392613 7 49 101.017014401513 7 49 102.139371114488 7 49 97.7563283358633 7 49 102.202720631733 7 49 102.900505664288 7 49 101.631437508119 7 49 106.246369334446 7 50 97.3363613712318 7 50 104.46740233381 7 51 92.7627489755597 7 51 88.2311303368408 7 51 107.116583153958 7 51 102.410098419439 7 51 96.378785050323 7 52 98.1762730625568 7 52 96.203860940023 7 52 99.1156334881672 7 52 96.5419115014224 7 52 102.100956942292 7 53 117.842216117441 7 53 124.174111459722 7 53 115.02591859104 7 53 111.885901553834 7 53 115.320308219732 END