RegMax | R Documentation |
Conduct the simultaneous inference for a set of regression coefficients in the null hypothesises H01 that assume the set of regression coefficients components are all zeroes.
RegMax(X, Y, G1, Nsplit = 5, sub.frac=0.5, alpha=0.05, seed=1, standardized=FALSE)
X |
a |
Y |
a |
G1 |
a positive vector with values between 1 and p, the set of H01. |
Nsplit |
a positive integer, the random split times used, default as 5. |
sub.frac |
a positive number between 0 and 1, the proportion of the sample used in the stage I. |
alpha |
a positive real, the significance level. |
seed |
a non-negative integer, the random seed. |
standardized |
a logical value, whether standerdize the covariates matrix in the stage I. |
return a vector with names 'CriticalValue', 'TestStatistic', 'reject_status', 'p-value' if Nsplit=1, and 'reject_status' and 'adjusted_p-value' if Nsplit>1.
nothing
Liu Wei
Liu, W., Lin, H., Liu, J., & Zheng, S. (2020). Two-directional simultaneous inference for high-dimensional models. arXiv preprint arXiv:2012.11100.
gendata_Reg
### Example n <- 50; p <- 20; i <- 1 s0 <- 5 # First five components are nonzeros rho <- 1; dat1 <- gendata_Reg(n, p, s0, seed=i, rho) # ex1: H01 is false RegMax(dat1$X, dat1$Y, 1:p) # ex1: H01 is true RegMax(dat1$X, dat1$Y, p)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.