margint.rob {rmargint} | R Documentation |
This function computes robust marginal integration procedures for additive models.
margint.rob(Xp, yp, point = NULL, windows, prob = NULL, sigma.hat = NULL, win.sigma = NULL, epsilon = 1e-06, type = "0", degree = NULL, typePhi = "Huber", k.h = 1.345, k.t = 4.685, max.it = 20, qderivate = FALSE, orderkernel = 2, Qmeasure = NULL)
Xp |
Matrix (n by p) of explanatory variables. |
yp |
Vector of responses (missing values are allowed). |
point |
Matrix of points where predictions will be computed and returned. |
windows |
Vector or a squared matrix of bandwidths for the smoothing estimation procedure. |
prob |
Probabilities of observing each response (n). Defaults to |
sigma.hat |
Estimate of the residual standard error. If |
win.sigma |
Vector of bandwidths for estimating sigma.hat. If |
epsilon |
Convergence criterion. |
type |
Three different type of estimators can be selected: type |
degree |
Degree of the local polynomial smoother in the direction of interest when using the estimator of type |
typePhi |
One of either |
k.h |
Tuning constant for a Huber-type loss function. Defaults to |
k.t |
Tuning constant for a Tukey-type loss function. Defaults to |
max.it |
Maximum number of iterations for the algorithm. |
qderivate |
If TRUE, it calculates |
orderkernel |
Order of the kernel used in the nuisance directions when using the estimator of type |
Qmeasure |
A matrix of points where the integration procedure ocurrs. Defaults to |
This function computes three types of robust marginal integration procedures for additive models.
A list with the following components:
mu |
Estimate for the intercept. |
g.matrix |
Matrix of estimated additive components (n by p). |
sigma.hat |
Estimate of the residual standard error. |
prediction |
Matrix of estimated additive components for the points listed in the argument point. |
mul |
A vector of size p showing in each component the estimated intercept that considers only that direction of interest when using the type |
g.derivative |
Matrix of estimated derivatives of the additive components (only when qderivate is |
prediction.derivate |
Matrix of estimated derivatives of the additive components for the points listed in the argument point (only when qderivate is |
Xp |
Matrix of explanatory variables. |
yp |
Vector of responses. |
Alejandra Martinez, ale_m_martinez@hotmail.com, Matias Salibian-Barrera
Boente G. and Martinez A. (2017). Marginal integration M-estimators for additive models. TEST, 26(2), 231-260. https://doi.org/10.1007/s11749-016-0508-0
function.g1 <- function(x1) 24*(x1-1/2)^2-2 function.g2 <- function(x2) 2*pi*sin(pi*x2)-4 set.seed(140) n <- 150 x1 <- runif(n) x2 <- runif(n) X <- cbind(x1, x2) eps <- rnorm(n,0,sd=0.15) regresion <- function.g1(x1) + function.g2(x2) y <- regresion + eps bandw <- matrix(0.25,2,2) set.seed(8090) nQ <- 80 Qmeasure <- matrix(runif(nQ*2), nQ, 2) fit.rob <- margint.rob(Xp=X, yp=y, windows=bandw, type='alpha', degree=1, Qmeasure=Qmeasure)