predCurves | R Documentation |
This function draws predicted curves (lines) against an explanatory variable for each category of a categorical variable.
predCurves(
object,
indata = NULL,
xname,
group = NULL,
legend = TRUE,
legend.space = "top",
legend.ncol = 2,
...
)
object |
Either an |
indata |
A data.frame object containing the data. If not specified, data is extracted from
the |
xname |
The name of variable to be plotted. |
group |
A character string or a sequence of length equivalent to rows of data to plot.
|
legend |
A logical value indicating whether a legend for |
legend.space |
A character string specifies one of the four sides,
which can be one of |
legend.ncol |
An integer specifies a number of columns, possibly
divided into blocks, each containing some rows. Default,
|
... |
Other arguments to be passed to |
Zhang, Z., Charlton, C.M.J., Parker, R.M.A., Leckie, G., and Browne, W.J. (2016) Centre for Multilevel Modelling, University of Bristol.
predLines
## Not run:
library(R2MLwiN)
# NOTE: if MLwiN not saved in location R2MLwiN defaults to, specify path via:
# options(MLwiN_path = 'path/to/MLwiN vX.XX/')
# If using R2MLwiN via WINE, the path may look like this:
# options(MLwiN_path = '/home/USERNAME/.wine/drive_c/Program Files (x86)/MLwiN vX.XX/')
## Read alevchem data
data(alevchem, package = "R2MLwiN")
alevchem$gcseav <- alevchem$gcse_tot/alevchem$gcse_no - 6
# Avoids warning when fitting factor as continuous response:
alevchem$a_point_num <- as.numeric(alevchem$a_point)
## Example: A-level Chemistry
(mymodel <- runMLwiN(a_point_num ~ 1 + gcseav + I(gcseav^2) + I(gcseav^3)
+ gender + (1 | pupil), estoptions = list(EstM = 1, resi.store = TRUE),
data = alevchem))
predCurves(mymodel, xname = "gcseav", group = "genderfemale")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.