各位老师，您好！我在运行中，遇到与jdd12345同样的问题。是哪里有问题？敬请帮助分析解答，非常感谢！
model {for (i in 1:K) {
for (j in 1:n) {
Y[i, j] ~ dnorm(eta[i, j], tauC)
eta[i, j] <- phi[i,1] *pow( (1 - phi[i,2]*exp(-phi[i,3]*x[j])),(1/(1-phi[i,4])))
}
phi[i, 1] <- exp(theta[i, 1]) # theta
phi[i, 2] <- exp(theta[i, 2])
phi[i, 3] <- exp(theta[i, 3])
phi[i, 4] <- exp(theta[i, 4])
theta[i, 1:4] ~ dmnorm(mu[1:4], tau[1:4, 1:4])
}
mu[1:4] ~ dmnorm(mean[1:4], prec[1:4, 1:4])
tau[1:4, 1:4] ~ dwish(R[1:4, 1:4], 4)
sigma2[1:4, 1:4] <- inverse(tau[1:4, 1:4])
for (i in 1:4) {sigma[i] <- sqrt(sigma2[i, i])}
tauC ~ dgamma(1.0E-2, 1.0E-2)
sigmaC <- 1 / sqrt(tauC)
}
list(n=60,K=13, x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60),
Y = structure(
.Data = c(800.536825,1525.59445,2633.110575,3566.913525,4498.27325,6004.741225,8157.4308,9557.8418, 10899.79175,12625.83518,14081.38665,15545.78068,17322.0922,18665.69535,20032.2003,21110.92155,22436.96173, 23686.02318,24541.17493,25060.80845,26097.68073,27731.68565,29250.06345,30919.91538,32630.24343,34365.56075, 35794.0928,37436.1842,38694.5139,40055.60015,41164.2769,42260.48438,43369.74715,44191.62455,45348.8596, 46364.43835,47445.96525,49037.38713,50331.62473,51722.33915,53399.84303,55163.91988,57171.27898,58327.3395, 59889.58755,61434.72335,63573.06908,66444.2944,69175.76825,71793.35148,74539.68338,77315.97395,80537.1838, 84242.45823,88504.49108,93122.18475,96883.6548,NA,NA,NA,
482.395675,1079.468075,1694.264575,2302.594825,4007.591175,5637.091025,7467.6437,9599.88015,11454.5557, 13125.08133,14497.03468,16532.00403,18368.06478,20392.5255,22535.36158,24193.71098,25000.0735,26041.62333, 26996.94358,28209.32235,29732.87215,30950.89638,32099.02475,33684.18393,35164.62385,37020.50448,38802.30095, 41257.31685,43682.64615,45548.3721,47246.9724,49398.5842,51017.64818,53158.76688,54887.815,56939.82185, 59469.04663,60941.64405,62657.24423,64586.40773,66079.96553,68012.46528,69325.32108,70868.68155,73229.15503, 75738.39898,77594.92355,79045.46405,81077.2353,82713.186,84772.2391,87220.47228,89766.88898,93319.5261,96033.62188, 98934.1758,NA,NA,NA,NA,
160.623475,552.569425,1216.120225,1808.4028,2313.87215,3120.8359,4218.6755,4906.415725,5773.1126,6463.750325, 7023.367125,7369.043525,7727.4324,8021.9512,8409.906125,8736.711625,9201.82325,9669.249175,10344.70285,
10924.83623,11603.63393,12615.1722,13362.8952,14139.81685,14916.6887,15830.9034,16727.85115,17407.36963, 18363.87105,19281.38798,20083.4464,21147.56025,21919.03768,22813.5424,23868.12205,25159.91385,26328.7602, 27623.67875,28691.48483,29659.74415,30889.38618,32004.73288,33609.92625,35026.77513,36328.14468,38084.4545, 39800.67823,42109.89958,44727.49305,47758.79593,50096.18235,53448.0804,56916.2315,60301.57705,63667.7271, 66516.6564,69548.3245,NA,NA,NA,
119.939325,506.393125,1313.323475, 2241.21935, 3401.8085, 4474.89325, 5483.73475, 6255.96945, 7042.247475, 7852.919, 8459.56575, 8976.406775, 9634.75375, 10260.01058, 10679.08375, 11261.42708, 11886.31583, 12841.79915, 13765.72478, 14423.31413, 15225.62065, 16362.8887, 17209.4482, 18139.10915, 19252.18033, 20720.9117, 22133.77903, 23498.99543, 24892.95345, 26481.17725, 27734.67555, 28893.4829, 30003.732, 31327.823, 32823.494, 34208.39233, 36015.95865, 37531.3536, 39391.03545, 42083.24733, 45163.2416, 48485.78645, 52747.08335, 56065.68985, 60425.8641, 64294.86348, 67159.93175, 70282.48415, 72695.71848, 75401.98808, 78017.1394, 80351.7753, 82885.01018, 86167.38845, 88848.0566, 91561.92063, 94040.63325, NA, NA, NA,
23.609525, 160.779, 328.57835, 592.29885, 1470.721275, 2505.83915, 3465.48495, 4562.878625, 5645.712, 6565.4725, 7423.082825, 8222.5859, 8807.754575, 9326.7274, 10016.6901, 10407.13563, 10844.16465, 11266.92293, 11613.8569, 11941.87488, 12313.97085, 12572.6302, 12737.44365, 12885.5228, 13018.33888, 13205.399, 13453.25423, 13790.66615, 14102.18853, 14400.34323, 14680.45903, 14942.51958, 15178.65765, 15420.44375, 15640.3549, 15804.81235, 16024.84518, 16232.97105, 16403.29838, 16703.00575, 16858.48085, 17018.69795, 17308.91075, 17524.57125, 17822.45143, 18073.6797, 18466.42745, 18828.8378, 18972.77485, 19114.4829, 19255.90173, 19406.30185, 19583.52728, 19901.85225, 20168.0341, 20346.66143, 20517.9394, NA, NA, NA,
148.0028, 590.0434, 1403.84175, 2201.31495, 3396.185375, 4648.859925, 6030.58975, 7150.30715, 8146.51395, 8834.6353, 9525.793725, 10176.42323, 10678.3648, 11308.79225, 11824.16668, 12538.87093, 13332.46768, 14021.11945, 14498.43065, 14976.70198, 15775.01838, 16416.53303, 17041.75705, 17902.82053, 18836.97868, 19725.3275, 20481.78915, 21299.31368, 22165.67285, 22912.8064, 23522.3051, 24101.11293, 24633.98495, 25280.55873, 25831.3144, 26574.51125, 27277.85383, 27985.4916, 28794.26485, 29751.838, 30456.4676, 31410.57388, 32278.65233, 32962.78485, 33874.72633, 34828.3402, 36008.32998, 37456.22625, 39139.6124, 40704.03693, 42604.61503, 44437.55343, 46423.33885, 48828.49813, 51208.56158, 53961.83515, NA, NA, NA, NA,
121.39095, 475.61765, 1315.6825, 2251.339725, 3193.256575, 4156.637675, 5250.3781, 6432.332125, 7270.970775, 8011.379575, 8611.8002, 9060.87005, 9530.0311, 9990.5896, 10492.8371, 10918.1986, 11462.0812, 12193.89353, 12957.96465, 13455.5432, 14063.68855, 14883.44985, 15535.28478, 16338.96648, 17181.77918, 18301.19385, 19409.39095, 20655.35943, 21868.14, 23300.53568, 24458.73115, 25491.28368, 26639.15728, 27873.4465, 29583.04253, 30897.02438, 32615.04975, 34187.64103, 35861.7292, 37992.34035, 40564.77128, 42243.60728, 44889.1885, 46185.98843, 47937.35525, 50554.46473, 52234.77563, 54666.58795, 57247.22428, 59459.71768, 61225.15903, 63570.297, 66185.99048, 69037.97805, 73251.80478, 76384.02973, 79766.4511, NA, NA, NA,
15.0545, 156.733175, 532.0181, 1163.536725, 2000.03315, 2424.7499, 3988.4589, 5181.039775, 6679.05655, 8016.961425, 9260.232725, 9972.745775, 10492.40353, 11296.28143, 12263.29383, 12948.07508, 13484.9526, 13854.22988, 14274.90718, 14740.19572, 15379.0258, 15880.764, 16278.80568, 16746.7459, 17242.32263, 17720.3809, 18280.35333, 18869.38495, 19478.71313, 20047.72938, 20680.71115, 21450.87558, 22029.92758, 22362.04055, 22640.93263, 22889.3092, 23220.59053, 23459.625, 23705.6484, 23937.43843, 24164.3869, 24340.76853, 24545.9718, 24812.64638, 25074.92468, 25334.10895, 25603.97173, 25897.43288, 26205.58095, 26435.844, 26669.86845, 26922.53125, 27141.7909, 27409.72843, NA, NA, NA, NA, NA, NA,
64.797475, 353.762975, 1093.39955, 2360.37665, 3747.1143, 5943.56015, 7968.373275, 10274.47293, 11562.37693, 12747.39018, 13692.6225, 14388.4377, 14762.67725, 15133.0987,
15627.63693, 16393.63245, 17149.04688, 17834.87298, 18305.86545, 18611.5706, 19254.29593, 20539.31175, 21699.50228, 22980.9599, 23968.24925, 24744.2012, 25454.1283, 26026.60675, 26917.09958, 27879.25825, 28601.5909, 29238.36848, 29959.45565, 31018.18945, 32527.18623, 34268.92933, 36312.30165, 37974.31178, 39294.06273, 40705.836, 42135.97573, 43228.86325, 44546.45185, 45716.7825, 46761.16105, 48101.12195, 49311.48595, 50656.94583, 51679.63175, 52437.89068, 53151.67935, 54005.05368, 54844.4638, 55964.82163, 57189.53148, 58235.90645, 59126.7894, NA, NA, NA,
10.2574, 42.338675, 72.68375, 193.7074, 401.0172, 1134.712475, 2149.34545, 3454.42585, 4773.88765, 6202.820975, 7407.9966, 8302.15535, 9387.1098, 10094.60078, 10905.26573, 11809.27223, 12449.07345, 12994.68875, 13464.14448, 13918.02293, 14545.61555, 15092.5654, 15488.4468, 15844.17198,16298.51388, 16694.47175, 17133.71363, 17703.33133, 18274.0536, 18837.49105, 19401.16885, 20117.1802, 20974.766, 21455.46193, 21985.56195, 22592.5785, 22949.7722, 23428.22048, 23859.37883, 24256.2315, 24674.13673, 25024.76305, 25340.9898, 25629.60095, 25906.1707, 26158.21393, 26376.93183, 26606.80488, 26911.84848, 27209.48773, 27490.73465, 27811.60103, 28060.74898, 28292.4402, 28562.68713, 28838.30238, 29040.2797, 29235.91295, 29475.42288, 29694.48263,
48.673, 284.25815, 863.240975, 1926.145725, 3165.958975, 5256.085825, 7725.554, 9936.3292, 12645.93238, 15010.22265, 17061.85253, 19574.8577, 22424.65218, 24729.1242, 27138.63818, 28615.46603, 30998.85255, 33252.1544, 35855.9953, 37372.08238, 39116.54823, 40781.75893, 42502.36775, 44413.3961, 46215.35758, 47533.2626, 49093.55428, 51032.1148, 52762.5402, 54978.24153, 56692.99433, 58116.88193, 59698.85585, 61384.62348, 63417.9137, 65777.96808, 68094.83115, 70889.04468, 73220.38198, 76379.63413, 80973.59043, 83681.35548, 86825.66935, 89187.23255, 91961.6621, 95403.939, 98885.79055, 102942.7377, 106809.7003, 110623.4302, 114078.6281, 118498.0183, 123321.7504, 128896.758, 134796.139, 140399.3125, 144849.0123, NA, NA, NA,
38.42925, 217.66235, 704.406875, 1370.219825, 2451.3492, 3823.370025, 4762.134975, 5620.4607, 6293.478825, 7102.91395, 7855.58525, 8673.30735, 9566.556075, 10311.95543, 10897.07218, 11511.027, 12196.34928, 13001.49955, 13736.5999, 14482.03045, 15918.49468, 17939.19575, 19294.23685, 20783.66155, 22578.50118, 24628.64363, 26598.79058, 28090.60363, 29773.8428, 30999.66418, 32107.36588, 33207.5038, 34239.9671, 35418.23145, 36902.17893, 38285.00693, 39745.03905, 41013.77495, 41993.26948, 42916.96765, 44185.08045, 45218.96105, 46729.33662, 47972.90348, 49357.76575, 50675.27613, 51887.2592, 52995.87153, 54256.96078, 55700.43883, 56830.45485, 58151.22815, 59656.14153, 61196.0301, 62956.60333, 64797.37763, 66415.85438, NA, NA, NA,
220.7985, 787.11695, 1381.431925, 2188.786725, 3647.072425, 5178.129525, 6639.63025, 8144.5204, 9604.853975, 11144.3953, 12019.8523, 12681.31355, 13453.28045, 14022.8392, 14612.95083, 15197.96748, 15794.00153, 16524.97223, 17634.3075, 18432.59385, 19890.46618, 21944.61818, 23208.77598, 24768.17425, 26273.3395, 27642.01563, 29167.58168, 30421.91905, 31658.09938, 32925.51768, 34102.89265, 35157.43853, 36267.9057, 37175.30618, 38302.05465, 39388.70765, 40463.56563, 41701.77558, 42594.16095, 43679.05325, 45429.61763, 46694.1779, 48302.36248, 49422.7665, 50538.18963, 52247.71483, 54097.94993, 55651.7662, 57382.26653, 59213.6578, 60636.9664, 62350.0335, 63887.6461, 65738.02775, 67865.54415, 69794.51858, 72333.96133, NA, NA, NA),
.Dim = c(13,60)),
mean=c(0,0,0,0),
R=structure(.Data=c(0.1,0,0,0,
0,0.1,0,0,
0,0,0.1,0,
0,0,0,0.1), .Dim=c(4,4)),
prec=structure(.Data= c(1.0E-3, 0, 0,0,
0,1.0E-3, 0,0,
0,0,1.0E-3,0,
0,0,0,1.0E-3), .Dim = c(4,4)))
Inits
list( theta=structure(
.Data=c(12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01,
12,0.001,2,0.01),
.Dim=c(13,4)),
mu=c(12,0.001,2,0.01),
tau=structure(.Data=c(0.1,0,0,0,
0,0.1,0,0,
0,0,0.1,0,
0,0,0,0.1),.Dim=c(4,4)),
tauC=2000)
TRAP 20 (precondition violated)
Math.Power [000001D5H]
.x REAL -0.3339797871997909
.y REAL 2.439812672847587
MathSpecfunc.Power [000008C1H]
.absY REAL 2.439812672847587
.n INTEGER 2
.power REAL 0.1320756233643248
.x REAL -0.3339797871997909
.y REAL 2.439812672847587
GraphStack.Node.Value [00000F97H]
.i INTEGER 13
.j INTEGER 4
.k INTEGER 3
.node GraphStack.Node [02AFD6D0H]
.op INTEGER 21
.p GraphNodes.Node [02B5B0F0H]
.top INTEGER 1
.work ARRAY 10 OF REAL Elements
GraphNormal.Node.Likelihood [00000284H]
.logTau REAL -19.91873285926873
.mu REAL 4703.107024087158
.node GraphNormal.Node [02B42C50H]
.tau REAL 2.235652155511597E-9
.x REAL 800.536825
GraphStochastic.Node.Conditional [000004BCH]
.likelihood GraphStochastic.List [02BFC070H]
.log REAL -614.7253757546347
.node GraphStochastic.Node [02B41D60H]
UpdaterMetnormal.Updater.MCMC [000001B8H]
.alpha REAL -0.0747334161769686
.newDen REAL -608.533817200884
.newVal REAL -0.149468385663475
.oldDen REAL -608.0509908114547
.oldVal REAL 0.7544179870690695
.overRelax BOOLEAN TRUE
.prior GraphStochastic.Node [02B41D60H]
.rate REAL 369.51
.res INTEGER 0
.updater UpdaterMetnormal.Updater [02C5A800H]
BugsUpdaters.Updater.MCMC [000005B8H]
.chain INTEGER 0
.depth INTEGER 1
.i INTEGER 2
.name BugsNames.Name [02BC5D80H]
.ok BOOLEAN TRUE
.overRelax BOOLEAN TRUE
.res INTEGER 0
.size INTEGER 52
.string ARRAY 120 OF CHAR "" ...
.updater BugsUpdaters.Updater [02C371E0H]
BugsSampler.MCMC [0000036DH]
.chain INTEGER 0
.cursor BugsSampler.List [02C5B780H]
.depth INTEGER 1
.depth0 INTEGER 1
.depth1 INTEGER 3
.offset INTEGER 6416052
.ok BOOLEAN TRUE
.overRelax BOOLEAN TRUE
.updater BugsUpdaters.Updater [02C371E0H]
BugsCmds.Action.Step [00000B1FH]
.a BugsCmds.Action [02B710B0H]
.chain INTEGER 0
.i INTEGER 1
.j INTEGER 406
.numChains INTEGER 1
BugsCmds.Action.Do [00000C58H]
.a BugsCmds.Action [02B710B0H]
.elapsedTime LONGINT 634921F163480384H
.profileList ARRAY 2024 OF CHAR "" ...
.res INTEGER 0
.s ARRAY 120 OF CHAR 17CX, 1AX ...
Services.Exec [00000136H]
.a Services.Action [02B710B0H]
.t POINTER [63480384H]
扫码加好友，拉您进群
全部版块
我的主页
收藏
