help
list(hyper_eta=1,hyper_tau=0.1,mu=c(100,80),c=c(2,2),rho=c(0.2,0.2),q=c(0.1,0.2),
p=structure(.Data=c(0.1,0.2,0.0.3,0.4),.
Dim=c(2,2))
,Q=structure(.Data=c(12,23,33,44),.Dim=c(2,2)))
how to set value for this Q ?
MODEL
model {
##### Hyperpriors #######
# One over the typical concentration parameter
# Supplementary Text, eqn. 10
hyper_eta ~ dexp(1)
# Gene-gene variability in isoform noise
# Supplementary Text, eqn. 12
hyper_tau ~ dexp(0.1)
# Supplementary Text, eqn. 11
lognorm_param <- log(1/hyper_eta) - 0.5 * 1/hyper_tau
for (g in 1:nGene) {
###### Priors #######
#Reparametrization according to eqn. 5 & 6
alpha[g] <- c[g] * rho[g] / (sqrt(rho[g]*sqrt(rho[g]) + (1-rho[g])*(1-rho[g])))
beta[g] <- c[g] * (1-rho[g]) / (sqrt(rho[g]*sqrt(rho[g]) + (1-rho[g])*(1-rho[g])))
#Eqn. 11
c[g] ~ dlnorm(lognorm_param , hyper_tau)
#Eqn. 8
rho[g] ~ dbeta(1,1)
#Eqn. 15
q[g] ~ dbeta(1,1)
#Eqn. 14
mu[g] ~ dunif(0.01,10000)
#Reparametrization for negative binomial distribution
size[g] <- q[g] * mu[g] / (1-q[g])
for (i in 1:nCell) {
##### Latent Variables ######
#Isoform preference
#eqn. 5 (truncation avoids numerical qlems)
p[i,g] ~ dbeta(alpha[g], beta[g]) T(0.0001,0.9999)
#mRNA expression level
#eqb. 13
Q[i,g] ~ dbin(q[g], size[g])
#Isoform number
#eqn 2 & 3
M[1,i,g] ~dbin(p[i,g],Q[i,g]) #number of isoform 1
M[2,i,g] <- Q[i,g] - M[1,i,g]
for (j in 1:nIso) {
###### observables #######
#eqn. 1
N[j,i,g] ~ dbin(rteff[i],M[j,i,g]) #number of molecules we can observe
}
}
}
}
-------------------------------------------------------------------
This list is for BUGS users to support each other, ask questions or share ideas.
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|