Unconstrained quantile regression with multiple compositional predictors | R Documentation |
Unconstrained quantile regression with multiple compositional predictors.
ulc.rq2(y, x, z = NULL, tau = 0.5, xnew = NULL, znew = NULL)
y |
A numerical vector containing the response variable values. |
x |
A list with multiple matrices with the predictor variables, the compositional data. No zero values are allowed. |
z |
A matrix, data.frame, factor or a vector with some other covariate(s). |
tau |
The quantile to be estimated, a number between 0 and 1. |
xnew |
A matrix containing a list with multiple matrices with compositional data whose response is to be predicted. If you have no new data, leave this NULL as is by default. |
znew |
A matrix, data.frame, factor or a vector with the values of some other covariate(s). If you have no new data, leave this NULL as is by default. |
The function performs the unconstrained log-contrast quantile regression model.
The logarithm of the compositional predictor variables is used (hence no zero
values are allowed). The response variable is linked to the log-transformed data
without the constraint that the sum of the regression coefficients
equals 0. If you want the regression without the zum-to-zero contraints see
lc.rq2
. Extra predictors variables are allowed as well, for
instance categorical or continuous.
A list including:
mod |
The object as returned by the function quantreg::rq(). This is useful for hypothesis testing purposes. |
be |
The unconstrained regression coefficients. Their sum does not equal 0. |
est |
If the arguments "xnew" and znew were given these are the predicted or estimated values, otherwise it is NULL. |
Michail Tsagris.
R implementation and documentation: Michail Tsagris mtsagris@uoc.gr.
Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.
Koenker R. W. and Bassett G. W. (1978). Regression Quantiles, Econometrica, 46(1): 33–50.
Koenker R. W. and d'Orey V. (1987). Algorithm AS 229: Computing Regression Quantiles. Applied Statistics, 36(3): 383–393.
ulc.rq, lc.rq
y <- rnorm(150)
x <- list()
x1 <- as.matrix(iris[, 2:4])
x1 <- x1 / rowSums(x1)
x[[ 1 ]] <- x1
x[[ 2 ]] <- rdiri(150, runif(4) )
x[[ 3 ]] <- rdiri(150, runif(5) )
mod <- ulc.rq2(y, x)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.