Description Usage Arguments Details Value Author(s) See Also Examples
Similar to usual predict methods and
predict.cv.grpreg
in grpreg
package.
1 2 3 |
object |
A fitted " |
newdata |
Optionally, a matrix or data frame where to predict. If omits, the fitted predictors are used. |
type |
The type of prediction: " |
... |
Not used. |
This function gives the predictions at newdata
or all predictors if the
argument newdata
is not supplied. Typically, type = "response"
is
used for linear or poisson regression, and type = "class"
or
type = "probability"
is used for logistic regression.
The predicted values depending on the type.
Debin Qiu, Jeongyoun Ahn
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | library(MASS)
set.seed(23)
n <- 30 # sample size
p <- 3 # number of predictors in each group
J <- 50 # group size
group <- rep(1:J,each = 3) # group indices
X <- mvrnorm(n,seq(0,5,length.out = p*J),diag(p*J))
beta <- runif(12,-2,5)
mu <- X%*%matrix(c(beta,rep(0,p*J-12)),ncol = 1)
# linear regression with family = "gaussian"
y <- mu + rnorm(n)
## with cross-validation
gss11 <- grpss(X,y,group,select = TRUE,cross.validation = TRUE)
predict(gss11) # fitted values
predict(gss11, newdata = X[1,]) # predicted values at given 'newdata'
# logistic regression with family = "binomial"
set.seed(23)
y1 <- rbinom(n,1,1/(1 + exp(-mu)))
gss21 <- grpss(X,y1,group, criterion = "gDC",select = TRUE,
family = "binomial")
predict(gss21)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.