## Not run:
# #See ?'crisp-package' for a full example of how to use this package
#
# #generate data (using a very small 'n' for illustration purposes)
# set.seed(1)
# data <- sim.data(n = 15, scenario = 2)
#
# #fit model for a range of tuning parameters, i.e., lambda values
# #lambda sequence is chosen automatically if not specified
# crisp.out <- crisp(X = data$X, y = data$y)
# #or fit model and select lambda using 2-fold cross-validation
# #note: use larger 'n.fold' (e.g., 10) in practice
# crispCV.out <- crispCV(X = data$X, y = data$y, n.fold = 2)
#
# #plot the estimated relationships between two predictors and outcome
# #do this for a specific fit
# plot(crisp.out, lambda.index = 25)
# #or for the fit chosen using cross-validation
# plot(crispCV.out)
# ## End(Not run)
Run the code above in your browser using DataLab