alqrfe (version 1.3)

plot_alpha: plot multiple penalized quantile regression - alpha

Description

plot QR intercepts for several taus

Usage

plot_alpha(
  Beta,
  tau = 1:9/10,
  D,
  ylab = expression(alpha[1]),
  col = 2,
  lwd = 1,
  lty = 2,
  pch = 1,
  cex.axis = 1,
  cex.lab = 1,
  main = ""
)

Arguments

Beta

Numeric array, with three dimmensions: 1) tau, 2) coef., lower bound, upper bound, 3) exploratory variables.

tau

Numeric vector, identifies the percentiles.

D

intercept's number.

ylab

y legend

col

color.

lwd

line width.

lty

line type.

pch

point character.

cex.axis

cex axis length.

cex.lab

cex axis length.

main

title.

Examples

n = 10
m = 5
d = 4
N = n*m
L = N*d
x = matrix(rnorm(L), ncol=d, nrow=N)
subj = rep(1:n, each=m)
alpha = rnorm(n)
beta = rnorm(d)
eps = rnorm(N)
y = x %*% beta  + matrix(rep(alpha, each=m) + eps)
y = as.vector(y)

Beta = mqr_alpha(x,y,subj,tau=1:9/10, method="qr", ngrid = 10)
plot_alpha(Beta,tau=1:9/10,D=1)