# Xavier Fernández i Marín # dg abr 21 15:44:13 CEST 2013 # Explain percentage of regulatory agencies in countries model { for (i in 1:d.I) { d.y.perc.betareg[i] ~ dbeta(a[i], b[i]) a[i] <- mu[i] * phi[d.ic[i]] b[i] <- (1 - mu[i]) * phi[d.ic[i]] logit(mu[i]) <- (theta[d.time[i]] ) + (beta[1] * log.gdp.cap.sc[i]) + (beta[2] * igo.number.sc[i]) + (beta[3] * d.W.dist.sc[i]) + (beta[4] * d.W.imp.sc[i]) + (beta[5] * d.W.igo.sector.economic.sc[i]) + (beta[6] * d.W.hist.sc[i]) } # phi factor for (c in 1:d.C) { phi[c] ~ dgamma(0.001, 0.001) } # theta by time theta[1] ~ dnorm(0, 0.001) for (yr in 2:d.Y) { theta[yr] ~ dnorm(theta[yr-1], tau.theta) } tau.theta <- pow(sigma.theta, -2) sigma.theta ~ dunif(0, 10) # control variables beta[1:6] ~ dmnorm(b0[1:6], B0[1:6,1:6]) # missing values for (i in 1:d.I) { log.gdp.cap.sc[i] ~ dnorm(0, 1) } }