get_data {dynamite} | R Documentation |
Extract the Model Data of the Dynamite Model
Description
Returns the input data to the Stan model. Mostly useful for debugging.
Usage
get_data(x, ...)
## S3 method for class 'dynamiteformula'
get_data(x, data, time, group = NULL, ...)
## S3 method for class 'dynamitefit'
get_data(x, ...)
Arguments
x |
[dynamiteformula or dynamitefit ] The model formula or an
existing dynamitefit object. See dynamiteformula() and dynamite() .
|
... |
Ignored.
|
data |
[data.frame , tibble::tibble , or data.table::data.table ]
The data that contains the variables in the model in long format.
Supported column types are integer , logical , double , and
factor . Columns of type character will be converted to factors.
Unused factor levels will be dropped. The data can contain missing
values which will simply be ignored in the estimation in a case-wise
fashion (per time-point and per channel). Input data is converted to
channel specific matrix representations via stats::model.matrix.lm() .
|
time |
[character(1) ] A column name of data that denotes the
time index of observations. If this variable is a factor, the integer
representation of its levels are used internally for defining the time
indexing.
|
group |
[character(1) ] A column name of data that denotes the
unique groups or NULL corresponding to a scenario without any groups.
If group is NULL , a new column .group is created with constant
value 1L is created indicating that all observations belong to the same
group. In case of name conflicts with data , see the group_var element
of the return object to get the column name of the new variable.
|
Value
A list
containing the input data to Stan.
See Also
Model outputs
as.data.frame.dynamitefit()
,
as.data.table.dynamitefit()
,
as_draws_df.dynamitefit()
,
coef.dynamitefit()
,
confint.dynamitefit()
,
dynamite()
,
get_code()
,
get_parameter_dims()
,
get_parameter_names()
,
get_parameter_types()
,
ndraws.dynamitefit()
,
nobs.dynamitefit()
Examples
data.table::setDTthreads(1) # For CRAN
d <- data.frame(y = rnorm(10), x = 1:10, time = 1:10, id = 1)
str(get_data(obs(y ~ x, family = "gaussian"),
data = d, time = "time", group = "id"
))
[Package
dynamite version 1.5.2
Index]