切片采样算法在光释光年龄模型参数估计中的应用
彭俊, 董治宝, 韩凤清

Application of slice sampling method for optimizing OSL age models used for equivalent dose determination
Jun PENG, Zhibao DONG, Fengqing HAN
表1 多变量切片采样算法流程
Tab.1 A general procedure for applying the Slice sampling method to sampling from a multivariable model
输入:概率模型,迭代总次数nsim,初始参数状态X0,模型参数总 数N;
输出:随机采样样本及参数估计结果
1. 分配用于随机样本储存的数据矩阵(chain):chain=matrix(nrow=nsim, ncol=N);
2. 通过如下迭代方式对参数轮流更新:
do i=1, nsim
do j=1, N
根据切片采样算法以条件概率P(Xj/X0[-j])生成随机 变量Xj
X0[j] = Xj
chain[i,j] = Xj
end do
end do
3. 根据随机样本数据chain 进行参数估计