R/starts_uni_estimate_prepare_fitting_prior_model.R

In alexanderrobitzsch/STARTS: Functions for the STARTS Model

## File Name: starts_uni_estimate_prepare_fitting_prior_model.R
## File Version: 0.082


starts_uni_estimate_prepare_fitting_prior_model <- function( pars_est,
    prior_var_trait, prior_var_ar, prior_var_state, prior_a, sd0, estimator)
{
    prior_model <- NULL
    if ( estimator !='ML'){
        prior_model <- list()
        prior_model[['var_trait']] <- starts_uni_estimate_prepare_fitting_prior_variance(
                                            prior_var=prior_var_trait, sd0=sd0)
        prior_model[['var_ar']] <- starts_uni_estimate_prepare_fitting_prior_variance(
                                            prior_var=prior_var_ar, sd0=sd0)
        prior_model[['var_state']] <- starts_uni_estimate_prepare_fitting_prior_variance(
                                            prior_var=prior_var_state, sd0=sd0)
        vec <- round( c(NA, (prior_a[1]+1)*prior_a[2], (prior_a[1]+1)*(1-prior_a[2])), 4 )
        prior_model[[ 'a' ]] <- list( 'dbeta', as.list(vec) )
        pars_remove <- names(pars_est)[ ! pars_est ]
        if ( length(pars_remove) > 0 ){
            for (pp in pars_remove){
                prior_model[[pp]] <- NULL
            }
        }
    }
    return(prior_model)
}