InnProd | R Documentation |
Inference for weighted inner product of the regression vectors in high dimensional generalized linear regressions
InnProd(
X1,
y1,
X2,
y2,
G,
A = NULL,
model = c("linear", "logistic", "logistic_alter"),
intercept = TRUE,
beta.init1 = NULL,
beta.init2 = NULL,
split = TRUE,
lambda = NULL,
mu = NULL,
prob.filter = 0.05,
rescale = 1.1,
tau = c(0.25, 0.5, 1),
alpha = 0.05,
verbose = FALSE
)
X1 |
Design matrix for the first sample, of dimension |
y1 |
Outcome vector for the first sample, of length |
X2 |
Design matrix for the second sample, of dimension |
y2 |
Outcome vector for the second sample, of length |
G |
The set of indices, |
A |
The matrix A in the quadratic form, of dimension |
model |
The high dimensional regression model, either |
intercept |
Should intercept(s) be fitted for the initial estimators (default = |
beta.init1 |
The initial estimator of the regression vector for the 1st data (default = |
beta.init2 |
The initial estimator of the regression vector for the 2nd data (default = |
split |
Sampling splitting or not for computing the initial estimators.
It take effects only when |
lambda |
The tuning parameter in fitting initial model. If |
mu |
The dual tuning parameter used in the construction of the projection direction. If |
prob.filter |
The threshold of estimated probabilities for filtering observations in logistic regression. (default = 0.05) |
rescale |
The factor to enlarge the standard error to account for the finite sample bias. (default = 1.1) |
tau |
The enlargement factor for asymptotic variance of the
bias-corrected estimator to handle super-efficiency. It allows for a scalar
or vector. (default = |
alpha |
Level of significance to construct two-sided confidence interval (default = 0.05) |
verbose |
Should intermediate message(s) be printed. (default = |
est.plugin |
The plugin(biased) estimator for the inner product form of the regression vectors restricted to |
est.debias |
The bias-corrected estimator of the inner product form of the regression vectors |
se |
Standard errors of the bias-corrected estimator,
length of |
ci.mat |
The matrix of two.sided confidence interval for the inner product
form of the regression vector; row corresponds to different values of
|
X1 = matrix(rnorm(100*5), nrow=100, ncol=5)
y1 = -0.5 + X1[,1] * 0.5 + X1[,2] * 1 + rnorm(100)
X2 = matrix(rnorm(90*5), nrow=90, ncol=5)
y2 = -0.4 + X2[,1] * 0.48 + X2[,2] * 1.1 + rnorm(90)
G = c(1,2)
A = matrix(c(1.5, 0.8, 0.8, 1.5), nrow=2, ncol=2)
Est = InnProd(X1, y1, X2, y2, G, A, model="linear")
Est$ci
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.