# NOT RUN {
# }
# NOT RUN {
set.seed(1973)
## generate an array with two constant blocks
data(MajorAlly)
Y <- MajorAlly
fit <- NetworkChange(newY, R=2, m=2, mcmc=G, initial.s = initial.s,
burnin=G, verbose=0, v0=v0, v1=v1)
drawRegimeRaw(fit, newY)
# }
Run the code above in your browser using DataLab