LM: Lindsey's Method

In SunJinglan/StatComp21080: Density estimations using Lindsey's method and penalized Gaussian mixtures

Description

A density estimation using Lindsey's method with Gibbs sampling

Usage

LM(df, n = 20, a = min(df), b = max(df), iter = 5000, burnin = 1000)

Arguments

df	the sample to be estimated
n	the number of equal width bins
a	the minimum of the range to estimate the density
b	the maximum of the range to estimate the density
iter	the number of Gibbs samples
burnin	the number of between-sample random numbers

Value

a list where t represents the possilble values and lm represents the density of that value

Examples

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

## End(Not run)

SunJinglan/StatComp21080 documentation built on Dec. 24, 2021, 1:24 a.m.