extractsamples {BAMBI} | R Documentation |

## Extract MCMC samples for parameters from an angmcmc object

### Description

Extract MCMC samples for parameters from an angmcmc object

### Usage

```
extractsamples(object, par.name, comp.label, chain.no, drop = TRUE, ...)
```

### Arguments

`object` |
angular MCMC object |

`par.name` |
vector of names of parameters for which point estimates are to be computed. If |

`comp.label` |
vector of component labels (positive integers, e.g., |

`chain.no` |
vector of chain numbers whose samples are to be be used. in the estimation. By default all chains are used. |

`drop` |
logical. Should the dimension of the output be dropped, if |

`...` |
additional arguments to be passed to the function. |

### Details

The default for both `par.name`

and `comp.label`

are the all possible choices
available in `object`

.

### Value

Returns a four dimensional array with

dimension 1 - model parameters and mixing proportions dimention 2 - components dimension 3 - MCMC iterations dimension 4 - chain number

### Examples

```
# first fit a vmsin mixture model
# illustration only - more iterations needed for convergence
fit.vmsin.20 <- fit_vmsinmix(tim8, ncomp = 3, n.iter = 20,
n.chains = 1)
# extract Markov chain realizations for kappa1 from component 1
extr_kappa1_1 <- extractsamples(fit.vmsin.20, "kappa1", 1)
# for kappa1 from component from all components
extr_kappa1 <- extractsamples(fit.vmsin.20, "kappa1")
# for all parameters in component 1
extr_1 <- extractsamples(fit.vmsin.20, comp.label = 1)
```

*BAMBI*version 2.3.5 Index]