dsi*_*cha 6 algorithm statistics optimization montecarlo
生成随机列联表的有效方法是什么？列联表被定义为矩形矩阵,使得每行的总和是固定的,并且每列的总和是固定的,但是只要每行和列的总和是正确的,各个元素可以是任何一个.
请注意,生成随机列联表非常容易,但我正在寻找比天真算法更有效的东西.
小智 6
查看R 的networksis包的代码可能会有所帮助.我认为有效的计算需要花哨的马尔可夫链序列重要性重采样技术,所以如果你能避免它,你可能想避免重新实现它.
编辑:相关论文是Chen,Diaconis,Holmes和Liu(2005).用作者的话来说,"[o] ur方法与其他现有的基于Monte Carlo的算法相比,有时会比其他几个数量级更有效."