We use cookies to give you the best experience possible. By continuing we’ll assume you’re on board with our cookie policy

Check Writers' Offers

Hire a Professional Writer Now

The input space is limited by 250 symbols

Choose 3 Hours or More.
Back
2/4 steps

How Many Pages?

Back
3/4 steps

Back
Get Offer

# Traumatic Systolic Blood Pressure ModelConsider observed Essay

Paper type: Essay Pages: 5 (1102 words)

Views: 234

Traumatic Systolic Blood Pressure Model

Consider observed systolic blood pressure monitored during traumatic events on a representative sample of size n from a given population. We assume the response is associated with a time where, We consider a model

where is an unknown smooth function and an error vector. For Bayesian approach to inference, we must put a prior on the unknown functions as well as all parameters.

We consider Gaussian process prior (Rasmussen and Williams 2006) on the functional term

where represents Gaussian process with mean function and covariance function .

For simplicity, we use a mean zero Gaussian process. That is, we set The covariance function will be chosen to have a stationary parametric form. In particular, is modelled with the Gaussian covariance function,

where is the process variance and is a spatial dependence parameter. It is well known that the derivative of a Gaussian process with covariance function (2) exists in mean square and it is also a Gaussian process with covariance function as the second derivative given by The process derivative is characterized with variance .

This variance indicates how the spatial dependence parameter relates to the dynamic properties of the trend. Let be a set of individual specific covariates which is informative about the dynamics of traumatic events. Since each subject in the representative sample is characterized with a unique covariate set, we propose to model the spatial dependence parameter based on some important covariate statistics instead of the covariates themselves.

In particular, we consider normal model for the spatial dependence parameter priorly,

where, is a vector of coefficients modelled with a normal prior where and are known. We consider inverse gamma models for where are known. The remaining standard deviation parameters are modelled with

distributions. where all hyperparameters are known and denotes the half-t distribution with scale and degrees of freedom . A half-t model can be modelled as a scale mixture of inverse gamma distributions (Wand et al. 2011). Let be auxiliary variables. Then the half-t priors in terms of scale mixtures of inverse gamma models are;

These alternative forms of the Half-t priors will be employed to simplify computations.

Variational Bayes Approach for Gaussian processes

Consider likelihood function for sample from a population with parameter Let be the prior model for the parameter Then inference for the parameter in the Bayesian paradigm is made via the posterior distribution

where Usually, the posterior distribution in most applications is intractable due the model evidence except very simple conjugate models.

Variational Bayes methods (Attias 2000), approximate the true intractable posterior distribution with a distribution belonging to some tractable class The distribution, can be allowed to follow some parametric family such as multivariate normal, or can be divided into blocks with the assumption of posterior independence between the blocks. Then a choice for is made to approximate as well as possible, via minimizing the Kullback-Leibler divergence between the two distributions,

Based on (3), the logarithm of the model can be expressed as

Due to the non-negativity of the Kullback-Leibler divergence, that is,

serves as a lower bound on where denotes expectation with respect to . From (5), minimizing the Kullback-Leibler divergence between and is equivalent to maximizing the lower bound . Note that when Otherwise, will be close to and a good approximation is attained for purposes such as Bayesian model choice if is close to (Bishop, (2006) and Ormerod and Wand, (2010)).

Now we apply the variational approximation technique to our traumatic systolic blood pressure model (1). Let denotes the full set of parameters in the model, we consider the following variational approximation to the posterior distribution

## where

and

Let denotes the set of all variational parameters in the model. We will write for suppressing the dependence on when needed. The lower bound can be derived in closed form based on the above variational posterior. We optimize this lower bound with respect to the variational parameters via variational Bayes iterative algorithm, where the following convention for notation has been adopted. and The expectations involved in the above notations are evaluated based on Lemma 1 with details provided in the appendix. The full set of the variational updates are given in Algorithm 1 below.

Algorithm 1 Variational Bayes Algorithm for TSBP model in (1)

## Initialize:

Set

Do until the change in the lower bound is less than a specified tolerance:

















Posterior inference with Hybrid Adaptive Metropolis Sampler

Our MCMC algorithm uses Gibbs steps for all parameters except the parameters and for which adaptive Metropolis steps are employed. The derivation of the full conditional distributions for the Gibbs updates is given in the next subsection.

1. Set tuning parameter and acceptance rate; Select a batching schedule

2. Initialization:

### Set and select starting values

3. Metropolis-within-Gibbs update: Update using Metropolis as follows;

## where

4. Update using Gibbs as follows.

a)

b)

c)

d)

e)

f)

5. Adaptation step: Update the proposal distribution:

a)

Traumatic Systolic Blood Pressure ModelConsider observed. (2019, Dec 17). Retrieved from https://studymoose.com/traumatic-systolic-blood-pressure-modelconsider-observed-best-essay

How to Avoid Plagiarism
• Use multiple resourses when assembling your essay
• Use Plagiarism Checker to double check your essay
• Get help from professional writers when not sure you can do it yourself 