# \donttest{
# Define network
network <- mb.network(obesityBW_CFB, reference="plac")
# Run an MBNMA with a quadratic time-course function
quad <- mb.run(network,
fun=tpoly(degree=2, pool.1="rel", method.1="common",
pool.2="rel", method.2="common"),
intercept=TRUE)
# Predict responses
pred <- predict(quad, times=c(0:50), treats=c(1:5),
ref.resp = network$data.ab[network$data.ab$treatment==1,],
E0=10)
# Generate summary of predictions
summary(pred)
# }
Run the code above in your browser using DataLab