当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R r2dtable 具有给定边际的随机 2 向表


R语言 r2dtable 位于 stats 包(package)。

说明

使用 Patefield 算法生成具有给定边际的随机 2 向表。

用法

r2dtable(n, r, c)

参数

n

一个非负数字,给出要绘制的表格的数量。

r

长度至少为 2 的非负向量,给出行总数,强制为 integer 。总和必须与 c 相同。

c

长度至少为 2 的非负向量,给出列总计,强制为 integer

长度为 n 的列表,包含生成的表作为其组件。

例子

## Fisher's Tea Drinker data.
TeaTasting <-
matrix(c(3, 1, 1, 3),
       nrow = 2,
       dimnames = list(Guess = c("Milk", "Tea"),
                       Truth = c("Milk", "Tea")))
## Simulate permutation test for independence based on the maximum
## Pearson residuals (rather than their sum).
rowTotals <- rowSums(TeaTasting)
colTotals <- colSums(TeaTasting)
nOfCases <- sum(rowTotals)
expected <- outer(rowTotals, colTotals) / nOfCases
maxSqResid <- function(x) max((x - expected) ^ 2 / expected)
simMaxSqResid <-
    sapply(r2dtable(1000, rowTotals, colTotals), maxSqResid)
sum(simMaxSqResid >= maxSqResid(TeaTasting)) / 1000
## Fisher's exact test gives p = 0.4857 ...

参考

Patefield, W. M. (1981). Algorithm AS 159: An efficient method of generating r x c tables with given row and column totals. Applied Statistics, 30, 91-97. doi:10.2307/2346669.

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Random 2-way Tables with Given Marginals。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。