edgep {BiDAG} | R Documentation |
This function estimates the posterior probabilities of edges by averaging over a sample of DAGs obtained via an MCMC scheme.
edgep(MCMCchain, pdag = FALSE, burnin = 0.2, endstep = 1)
MCMCchain |
an object of class |
pdag |
logical, if TRUE (FALSE by default) all DAGs in the MCMCchain are first converted to equivalence class (CPDAG) before the averaging |
burnin |
number between |
endstep |
number between |
a square matrix with dimensions equal to the number of variables; each entry [i,j]
is an estimate of the posterior probability of the edge from node i
to node j
Polina Suter
Bostonscore<-scoreparameters("bge", Boston) ## Not run: samplefit<-orderMCMC(Bostonscore, iterations=25000,chainout=TRUE) edgesposterior<-edgep(samplefit, pdag=TRUE, burnin=0.2) ## End(Not run)