fracpoly_mvshape_plot overlays the fractional polynomial plot with the mvshape plot.
Usage
fracpoly_mvshape_plot(
fracpoly,
degree = "both",
mvshape,
xref = NULL,
logx = FALSE,
logy = FALSE,
pref_x = "x",
pref_y = "y",
xbreaks = NULL,
ybreaks = NULL,
cicolour = "grey",
betacolour = "black",
intcolour = "grey",
mvbetacolour = "red",
mvcicolour = "red",
xlim = NULL,
ylim = NULL
)Arguments
- fracpoly
a fracpoly object.
- degree
the fractonal polyniomal degree to be plotted (options: 1, 2 and both).
- mvshape
a mvshape object.
- xref
the reference point for the figures of binary outcomes.
- logx
plot the x-axis on the log-scale.
- logy
plot the y-axis on the log-scale.
- pref_x
the prefix for the x-axis.
- pref_y
the prefix for the y-axis.
- xbreaks
breaks in the x-axis.
- ybreaks
breaks in the y-axis.
- cicolour
colour of the 95% confidence intervals.
- betacolour
colour of the regression line.
- intcolour
colour of the intercept line.
- xlim
x-axis limits.
Author
James Staley jrstaley95@gmail.com
Examples
if (FALSE) {
# Data
y <- rnorm(5000)
x <- rnorm(5000, 10, 1)
c1 <- rbinom(5000, 1, 0.5)
c2 <- rnorm(5000)
study <- c(rep("study1", 1000), rep("study2", 1000), rep("study3", 1000), rep("study4", 1000), rep("study5", 1000))
covar <- data.frame(c1 = c1, c2 = c2, study = study)
# Analyses
fp <- fracpoly(y = y, x = x, covar = covar, family = "gaussian")
mvma <- mvshape(y = y, x = x, covar = covar[, c("c1", "c2")], study = study, family = "gaussian")
fracpoly_mvshape_plot(fracpoly = fp, mvshape = mvma)
}