data(rats)
fit <- network(rats,specifygraph=FALSE,inspect=FALSE)
fit.prior <- jointprior(fit,12)
W1 <- fit$nodes$W1
W1 <- cond.node(W1,fit,fit.prior)
W1.post <- postc0c(W1$condprior[[1]]$mu,
W1$condprior[[1]]$tau,
W1$condprior[[1]]$rho,
W1$condprior[[1]]$phi,
rats[,W1$idx])
Run the code above in your browser using DataLab