Description Usage Arguments Details Value Source Examples
This function fits a Bayesian circular regression model based on the projected normal distribution.
1 2 3 4 5 6 7 8 9 |
pred.I |
model equation for effects of component 1. |
data |
the dataframe used for analysis. |
pred.II |
model equation for effects of component 2. |
its |
output iterations of the MCMC sampler. |
burn |
number of burn-in iterations. |
n.lag |
amount of lag for the iterations and burn-in. |
seed |
user-specified random seed. |
Because the model is based on the projected normal distribution, a
model equation has to be given for two components. By default the equation
of the second component pred.II
is set to be equal to that of the
first component. Note that the circular outcome needs to be measured in
radians on a scale from 0 to 2π. For more information about the
projected normal distribution see Presnell, Morrisson & Littell (1998).
A tutorial on how to use this function can be found in Cremers & Klugkist (2018). More details on the sampling algorithm and interpretation of the coefficients from the model can be found in Cremers, Mulder & Klugkist (2018) and Cremers, Mainhard & Klugkist (2018). The uninformative priors for the regression coefficients of the two components are set to N(0, 10000).
A bpnr
object, which can be further analyzed using the
associated functions traceplot.bpnr
, BFc.bpnr
,
coef_lin.bpnr
, coef_circ.bpnr
,
fit.bpnr
and print.bpnr
.
A bpnr
object contains the following elements (some elements are not
returned if not applicable)
beta1
A matrix of posterior samples for the
coefficients beta1
of the first component.
beta2
A matrix
of posterior samples for the coefficients beta2
for the second
component.
Likelihood
A matrix containing the posterior density values for all individuals in the dataset for all iterations. The rowsums of this matrix are the likelihood values for all iterations
its
Number of output iterations.
n.lag
One in
n.lag
iterations will be saved as output iteration. Set lag to 1 to
save all iterations (default).
burn-in
Burn-in time for the MCMC sampler.
p1
Number of parameters predicting the first component.
p2
Number of parameters predicting the second component.
theta
The circular outcome vector measured in radians.
a.x
A matrix of posterior samples for a.x
which
describes the location of the inflection point of the regression curve on
the axis of the predictor.
a.c
A matrix of posterior samples
for a.c
which describes the location of the inflection point of the
regression curve on the axis of the circular outcome.
b.c
A
matrix of posterior samples for b.c
which describes the slope of the
tangent line at the inflection point.
SAM
A matrix of posterior samples for the circular regression slopes at the mean.
AS
A matrix of posterior samples for the average slopes of the circular regression.
SSDO
A matrix of posterior samples for the signed shortest distance to the origin.
circ.diff
A matrix of posterior samples for the circular difference between levels of categorical variables and the intercept.
Call
The matched call.
lin.coef.I
The mean, mode, standard deviation and 95
confidence interval of the highest posterior density of the linear
coefficients for beta1
.
lin.coef.II
The mean, mode,
standard deviation and 95
density of the linear coefficients for beta2
.
circ.coef
The mean, mode, standard deviation and 95
confidence interval of the highest posterior density for the a.x
,
a.c
, b.c
, AS
, SAM
and SSDO
of the
circular coefficients.
circ.coef.cat
The mean, mode, standard deviation and 95 density the circular difference between levels of categorical variables and the intercept.
circ.coef.means
The mean, mode, standard deviation and 95 circular means of the categorical variables.
model.fit
A list of information criteria for assessment of model fit.
mm
A list of information, model matrices, sample size, etc. on the specified model.
Cremers, J., Mulder, K.T. & Klugkist, I. (2018). Circular interpretation of regression coefficients. British Journal of Mathematical and Statistical Psychology, 71(1), 75-95.
Cremers, J., Mainhard, M.T. & Klugkist, I. (2018). Assessing a Bayesian Embedding Approach to Circular Regression Models. Methodology, 14, 69-81.
Cremers, J. & Klugkist, I. (2018). One direction? A tutorial for circular data with examples in cognitive psychology. Frontiers in Psychology: Cognitive Science.
Presnell, B., Morrison, S.P. & Littell, R.C. (1998). Projected multivariate linear models for directional data. Journal of the American Statistical Association, 93 (443), 1068 - 1077.
1 2 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.