Learn R Programming

PReMiuM (version 3.2.13)

summariseVarSelectRho: summariseVarSelectRho

Description

This function summarises the posterior distribution of rho, a parameter for variable selection only.

Usage

summariseVarSelectRho(runInfoObj)

Value

A list with the following elements.

rho

A matrix that has as many columns as the number of covariates and as many rows as the number of sweeps. This matrix records the samples from the posterior distribution of rho for each covariate at each sweep.

rhoMean

Vector with the column means of the matrix rho above. Each value corresponds to the posterior mean of rho for each covariate.

rhoMedian

Vector with the column medians of the matrix rho above. Each value corresponds to the posterior median of rho for each covariate.

rhoLowerCI

Vector with the column lower confidence intervals of the matrix rho above. Each value corresponds to the lower confidence interval of the posterior distribution of rho for each covariate.

rhoUpperCI

Vector with the column upper confidence intervals of the matrix rho above. Each value corresponds to the upper confidence interval of the posterior distribution of rho for each covariate.

Arguments

runInfoObj

Object of type runInfoObj

Authors

David Hastie, Department of Epidemiology and Biostatistics, Imperial College London, UK

Silvia Liverani, Department of Epidemiology and Biostatistics, Imperial College London and MRC Biostatistics Unit, Cambridge, UK

Maintainer: Silvia Liverani <liveranis@gmail.com>

References

Silvia Liverani, David I. Hastie, Lamiae Azizi, Michail Papathomas, Sylvia Richardson (2015). PReMiuM: An R Package for Profile Regression Mixture Models Using Dirichlet Processes. Journal of Statistical Software, 64(7), 1-30. tools:::Rd_expr_doi("10.18637/jss.v064.i07").

Examples

Run this code
if (FALSE) {
inputs <- generateSampleDataFile(clusSummaryVarSelectBernoulliDiscrete())

runInfoObj<-profRegr(yModel=inputs$yModel, 
    xModel=inputs$xModel, nSweeps=10, nClusInit=15, 
    nBurn=20, data=inputs$inputData, output="output", 
    covNames = inputs$covNames, varSelect="Continuous")

rho<-summariseVarSelectRho(runInfoObj)
}

Run the code above in your browser using DataLab