M：加速EM、MM及其他固定点算法收敛.pdfVIP

SQUAREM：加速EM、MM及其他不动点算法的收敛

RaviVaradhan

1SQUAREM概述

“

SQUAREM”是一个旨在加速缓慢收敛收缩映射的软件包。它可以用于加速如EM（期

望最大化）算法、MM（重大化最小化）算法以及其他非线性固定点迭代（例如用于寻找

主特征向量的幂法）等慢速、线性收敛的收缩映射。它采用了一种称为平方外推法（

SQUAREM）的新方法，该方法由Varadhan和Roland（《斯堪的纳维亚统计学》，

35:335‑353）以及Roland、Vardhan和Frangakis（《数值算法》，44:159‑172）提

出。

本软件包中的函数可通过以下方式使用：

library(SQUAREM)

您可以查看该包的基本信息，包括所有可用的函数：

help(package=SQUAREM)

虽然不是必需的，但如果您希望完全重现本指南中的示例，请执行以下操作：

require(setRNG)setRNG(list(kind=Wichmann‑Hill,normal.kind=Box‑Muller,

seed=123))

，示例需要按照这里的顺序运行（或者至少是生成随机数的部分）。对于某些示例