visW {DrBats} R Documentation

## Plot the estimates for the latent factors

### Description

Plot the estimates for the latent factors

### Usage

```visW(mcmc.output, Y, D, chain = 1, factors = c(1, 2))
```

### Arguments

 `mcmc.output` an mcmc list as produced by clean.mcmc `Y` the matrix of data `D` the number of latent factors `chain` the chain to plot (default = 1) `factors` a vector indicating the factors to plot (default = c(1, 2))

### Value

res.W a data frame containing the estimates for the factors, and their lower and upper bounds

Inertia the percentage of total inertia captured by each of the factors

### Author(s)

Gabrielle Weinrott

### Examples

```data("toydata")
data("stanfit")
codafit <- coda.obj(stanfit) ## convert to mcmc.list
W.res <- visW(codafit, Y = toydata\$Y.simul\$Y, D = toydata\$wlu\$D,
chain = 1, factors = c(1, 2))

## plot the results

data <- data.frame(time = rep(1:9, 2), W.res\$res.W)
ggplot2::ggplot() +
ggplot2::geom_step(data = data, ggplot2::aes(x = time, y = Estimation, colour = Factor)) +
ggplot2::geom_step(data = data, ggplot2::aes(x = time, y = Lower.est, colour = Factor),
linetype = "longdash") +
ggplot2::geom_step(data = data, ggplot2::aes(x = time, y = Upper.est, colour = Factor),
linetype = "longdash")

```

[Package DrBats version 0.1.5 Index]