data(rats)
fit <- network(rats)
fit.prior <- jointprior(fit,12)
fit <- learn(fit,rats,fit.prior)$nw
rats.extra <- addarrow(fit,rats,fit.prior)$nw
plot(rats.extra)
rats.extra2<- addarrow(rats.extra[[2]],rats,fit.prior)$nw
plot(rats.extra2)
rats.extra3<- turnarrow(rats.extra2[[2]],rats,fit.prior)$nw
plot(rats.extra3)
rats.extra4<- removearrow(rats.extra2[[2]],rats,fit.prior)$nw
plot(rats.extra4)
Run the code above in your browser using DataLab