data(ozone)
## Use half of the observations for demonstration
set.seed(27695)
train.id <- sort(sample(1:nrow(ozone),ceiling(nrow(ozone)/2)))
cossoObj <- cosso(x=ozone[train.id,2:5],y=ozone[train.id,1],nbasis=50)
## Compute model parameters
predict.cosso(cossoObj,M=2,type="c")
## Make prediction on full data
predict.cosso(cossoObj,ozone[,2:5],M=2,type="c")
Run the code above in your browser using DataLab