Description Usage Arguments Value Examples

Generic function that can be used to create and fetch the reference model
structure for all those objects that have this method. All these
implementations are wrappers to the `init_refmodel`

-function so
the returned object has the same type.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 | ```
get_refmodel(object, ...)
## S3 method for class 'refmodel'
get_refmodel(object, ...)
## S3 method for class 'vsel'
get_refmodel(object, ...)
## Default S3 method:
get_refmodel(
object,
data,
y,
formula,
ref_predfun,
proj_predfun,
div_minimizer,
fetch_data,
family = NULL,
wobs = NULL,
folds = NULL,
cvfits = NULL,
offset = NULL,
cvfun = NULL,
dis = NULL,
...
)
## S3 method for class 'stanreg'
get_refmodel(
object,
data = NULL,
ref_predfun = NULL,
proj_predfun = NULL,
div_minimizer = NULL,
folds = NULL,
...
)
init_refmodel(
object,
data,
formula,
family,
ref_predfun = NULL,
div_minimizer = NULL,
proj_predfun = NULL,
folds = NULL,
extract_model_data = NULL,
cvfun = NULL,
cvfits = NULL,
dis = NULL,
...
)
``` |

`object` |
Object on which the reference model is created. See possible types below. |

`...` |
Arguments passed to the methods. |

`data` |
Data on which the reference model was fitted. |

`y` |
Target response. |

`formula` |
Reference model's lme4-like formula. |

`ref_predfun` |
Prediction function for the linear predictor of the reference model. |

`proj_predfun` |
Prediction function for the linear predictor of the projections. |

`div_minimizer` |
Maximum likelihood estimator for the underlying projection. |

`fetch_data` |
Wrapper function for fetching the data without directly accessing it. It should have a prototype fetch_data(data, data_points, newdata = NULL), where data_points is a vector of data indices and newdata, if not NULL, is a data frame with new data for testing. |

`family` |
A family object that represents the observation model for the reference model. |

`wobs` |
A weights vector for the observations in the data. The default is a vector of ones. |

`folds` |
Only used for K-fold variable selection. It is a vector of fold indices for each data point in data. |

`cvfits` |
Only used for K-fold variable selection. A list of K-fold fitted objects on which reference models are created. |

`offset` |
A vector of offsets per observation to add to the linear predictor. |

`cvfun` |
Only used for K-fold variable selection. A function that, given a folds vector, fits a reference model per fold and returns the fitted object. |

`dis` |
A dispersion vector for each observation. |

`extract_model_data` |
A function with prototype extract_model_data(object, newdata, wrhs, orhs), where object is a reference model fit, newdata is either NULL or a data frame with new observations, wrhs is a right hand side formula to recover the weights from the data frame and orhs is a right hand side formula to recover the offset from the data frame. |

An object of type `refmodel`

(the same type as returned by
init_refmodel) that can be passed to all the functions that take the
reference fit as the first argument, such as varsel,
cv_varsel, project, proj_predict and
proj_linpred.

1 2 3 4 5 6 7 8 9 10 11 12 13 | ```
if (requireNamespace('rstanarm', quietly=TRUE)) {
### Usage with stanreg objects
dat <- data.frame(y = rnorm(100), x = rnorm(100))
fit <- rstanarm::stan_glm(y ~ x, family = gaussian(), data = dat)
ref <- get_refmodel(fit)
print(class(ref))
# variable selection, use the already constructed reference model
vs <- varsel(ref)
# this will first construct the reference model and then execute
# exactly the same way as the previous command (the result is identical)
vs <- varsel(fit)
}
``` |

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.