Learn R Programming

MixStable (version 0.1.0)

metropolis_hastings: Metropolis-Hastings MCMC for stable mixture clustering

Description

Performs Metropolis-Hastings sampling to estimate parameters of a two-component alpha-stable mixture model. Proposals are generated for each parameter and accepted based on likelihood ratios. Cluster assignments are updated at each iteration.

Usage

metropolis_hastings(fct, iterations, lok, aa = c(1, 1), proposal_std = 0.1)

Value

List of sampled parameters and weights across iterations.

Arguments

fct

Function to estimate initial parameters.

iterations

Number of MCMC iterations.

lok

Numeric vector of observations.

aa

Prior parameters for Dirichlet distribution.

proposal_std

Standard deviation for proposal distributions.