library("perryExamples")
data("coleman")
## evaluate MM regression models tuned for 85% and 95% efficiency
tuning <- list(tuning.psi = c(3.443689, 4.685061))
## via model fitting function
# perform cross-validation
# note that the response is extracted from 'data' in
# this example and does not have to be supplied
perryTuning(lmrob, formula = Y ~ ., data = coleman,
tuning = tuning, splits = foldControl(K = 5, R = 10),
cost = rtmspe, costArgs = list(trim = 0.1), seed = 1234)
## via function call
# set up function call
call <- call("lmrob", formula = Y ~ .)
# perform cross-validation
perryTuning(call, data = coleman, y = coleman$Y,
tuning = tuning, splits = foldControl(K = 5, R = 10),
cost = rtmspe, costArgs = list(trim = 0.1), seed = 1234)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.