gibbsR: A Gibbs sampler using R

In Lrrzz/StatComp21064: Example functions for 'Statistical Computing' course

Description

A Gibbs sampler using R

Usage

gibbsR(N, thin, n, aa, bb)

Arguments

N	the number of samples
thin	the number of between-sample random numbers
n	the self-selected fixed parameter 'n' in ex9.8
aa	the self-selected fixed parameter 'a' in ex9.8
bb	the self-selected fixed parameter 'b' in ex9.8

Value

a random sample of size n

Examples

## Not run: 
rnR <- gibbsR(100,10,22,2,2)
par(mfrow=c(2,1));
plot(rnR[,1],type='l')
plot(rnR[,2],type='l')

## End(Not run)

Lrrzz/StatComp21064 documentation built on Dec. 23, 2021, 10:21 p.m.