Description Usage Arguments Value References Examples

Implements a screening procedure proposed by Liu, Li and Wu(2014) for varying coefficient models with ultra-high dimensional predictors.

The function code is adapted from the relevant authors' code. Special thanks are due to Jingyuan Liu for providing some of the code upon which this function is based.

1 | ```
screenVCM(X, Y, U)
``` |

`X` |
Matrix of predictors to be screened. There should be one row for each observation. |

`Y` |
Vector of responses. It should have the same length as the number of rows of X. |

`U` |
Covariate, with which coefficient functions vary. |

A list with following components:

`CORR_sq:` |
A vector of the unconditioned squared correlation with length equal to the number of columns in the input matrix X. The hgh the unconditioned squared correlation is, the more desirable it is to retain the corresponding X covariate in a later predictive model. |

`rank:` |
Vector for the rank of the predictors in terms of the conditional correlation
( |

Liu, J., Li, R., & Wu, R. (2014). Feature selection for varying coefficient models with ultrahigh-dimensional covariates. Journal of the American Statistical Association, 109: 266-274. <DOI:10.1080/01621459.2013.850086>

1 2 3 4 5 6 7 8 9 | ```
set.seed(12345678)
data1 <- simulateVCM(p=250,trueIdx = c(2, 20, 80, 120, 200))
screenResults<- screenVCM(X = data1$X, Y = data1$Y, U = data1$U)
rank <- screenResults$rank
print(which(rank <= 10)) # This prints the numbers of the columns having the best rank.
trueIdx <- c(2, 20, 80, 120, 200) # These were the column numbers (indices) of
# the nonnull predictors in the simulated data;
print(rank[trueIdx]) # It can be seen that these predictors were all given high rank,
# showing that the function worked well.;
``` |

```
[1] 2 19 20 21 80 116 117 120 130 200
[1] 3 1 5 4 2
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.