# Generate an "mbnma.network" object that stores data in the correct format
network <- mbnma.network(triptans)
# Generate a network plot at the dose/treatment level
plot(network, level="treatment")
# Generate a network plot at the agent level
plot(network, level="agent", remove.loops=TRUE)
# \donttest{
# Perform "split" NMA to examine dose-response relationship
nma <- nma.run(network)
plot(nma)
# Analyse data using mbnma.run() with an Emax dose-response function
# and common treatment effects
result <- mbnma.run(network, fun=demax(),
method="common")
# Generate forest plots for model results
plot(result)
# Rank results and plot rankograms
ranks <- rank(result)
plot(ranks, params="emax")
# Predict responses
pred <- predict(result, E0=0.2)
# Plot predicted response with "split" NMA results displayed
plot(pred, overlay.split=TRUE)
# }
Run the code above in your browser using DataLab