Description Usage Arguments Value Author(s) References Examples
Perform MCMC iterations of the model, as described in the reference.
1 2 3 4 5 6 7 8 9 10 | iBAT(Y, X, distance, disfix, intercept=1, xi,
R=-1, tran, mu, sigma=((rgamma(4,1,1))^(-0.5)),
cmu=1/1000000, c=10, delta=3, d, e=0.001, f=0.999,
alpha=20, deltak=c(-1,0,0.58,1), tauk=c(1,1,1,2),
upp_bounds=c(-0.1, 0.1, 0.73, Inf),
low_bounds=c(-Inf, -0.1, 0.1, 0.73),
alpha_IG=c(1,1,1,1), beta_IG=c(1,1,1,1),
low_IG=c(0.41,0.41,0.41,1), a=c(1,1,1,1),
niter=500000, burnin=200000, Cout=1000,
phi=0.5, pR=0.4, selectioncgh=-1, pXI=0.6, indep=0)
|
Y |
Matrix of gene expression data |
X |
Matrix of CGH data |
distance |
Vector of distance between CGH probes |
disfix |
Length of the chromosome under investigation |
intercept |
If set to one an intercept is included in the regression model |
xi |
Initialized matrix of latent states |
R |
Initialized association matrix in a vector form. Default set to -1, that automatically creates a vector with all the positions set to zero |
tran |
Initialized transition matrix |
mu |
Initialized state specific mean vector |
sigma |
Initialized state specific standard deviation vector |
cmu |
Parameter that controls the variance of the prior on the intercept |
c |
Parameter that determines the shrinkage in the model |
delta |
Parameter of the Inverse-Gamma prior on the error variance |
d |
Parameter of the Inverse-Gamma prior on the error variance |
e |
Parameter of the Beta prior on the inclusion probability |
f |
Parameter of the Beta prior on the inclusion probability |
alpha |
Parameter that regulates the strength of the independent part of the mixture |
deltak |
Vector of mean of the prior on the state specific mean |
tauk |
Vector of sd of the prior on the state specific mean |
upp_bounds |
Vector of upper bounds of the prior on the state specific mean |
low_bounds |
Vector of lower bounds of the prior on the state specific mean |
alpha_IG |
Vector of parameters of the prior on the state specific standard deviation |
beta_IG |
Vector of parameters of the prior on the state specific standard deviation |
low_IG |
Truncation of the prior on the state specific standard deviation |
a |
Vector of parameters of the prior on the transition matrix |
niter |
Number of Monte Carlo Markov Chain iterations |
burnin |
Burn-in |
Cout |
Print the number of iterations ran every Cout iterations |
phi |
Probability of an A/D step |
pR |
Parameter of the distribution used to select the rows to be updated at every MCMC iteration |
selectioncgh |
Number of samples not in neutral state in order to consider a CGH as a potential candidate for association with gene expression. Default set to -1 that automatically set it to 10% of the samples |
pXI |
Parameter of the distribution used to select the rows to be updated at every MCMC iteration |
indep |
If set to an integer different from zero, run the analysis with an independent prior, see reference. |
The output consists of an R list composed by 4*niter+3 elements objects, where niter is the number of MCMC iterations. The first niter objects of the list are vectors, each containing the positions of the association matrix set to one, at the corresponding MCMC iteration. Each of the following niter objects of the list are the transition matrices at the corresponding MCMC iteration, while the third and the fourth set of niter objects are the vectors of state specific mean and state specific variance, respectively. The last three objects of the list consist of three matrices counting the number of times the corresponding latent state has been set to 1,3 and 4, respectively.
Alberto Cassese
Cassese A, Guindani M, Tadesse M, Falciani F, Vannucci M. A hierarchical Bayesian model for inference of copy number variants and their association to gene expression. Annals of Applied Statistics, 8(1), 148-175.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ## Not run:
data(NCI_60)
Y <- NCI_60$Affy
X <- NCI_60$aCGH
distance <- NCI_60$distance
disfix <- 146274826
xi <- InitXi(X)
tran <- Tran(xi)
mu <- InitMu()
d=0.2587288
Y <- Center(Y)
res <- iBAT(Y=Y,X=X,distance=distance,disfix=disfix,xi=xi,tran=tran,mu=mu,d=d)
summRes <- Inference(res,G=dim(Y)[[2]],M=dim(X)[[2]],niter=niter,burnin=bi,threshold=0.5)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.