R语言
r2dtable
位于 stats
包(package)。 说明
使用 Patefield 算法生成具有给定边际的随机 2 向表。
用法
r2dtable(n, r, c)
参数
n |
一个非负数字,给出要绘制的表格的数量。 |
r |
长度至少为 2 的非负向量,给出行总数,强制为 |
c |
长度至少为 2 的非负向量,给出列总计,强制为 |
值
长度为 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 rWishart 随机 Wishart 分布式矩阵
- R relevel 因子水平重新排序
- R reorder.default 因子水平重新排序
- R rect.hclust 在层次集群周围绘制矩形
- R replications 条款的重复次数
- R reshape 重塑分组数据
- R read.ftable 操作平面列联表
- R reorder.dendrogram 重新排序树状图
- R residuals 提取模型残差
- R runmed 运行中位数 – 稳健散点图平滑
- R stlmethods STL 对象的方法
- R medpolish 矩阵的中值波兰(稳健双向分解)
- R naprint 调整缺失值
- R summary.nls 总结非线性最小二乘模型拟合
- R summary.manova 多元方差分析的汇总方法
- R formula 模型公式
- R nls.control 控制 nls 中的迭代
- R aggregate 计算数据子集的汇总统计
- R deriv 简单表达式的符号和算法导数
- R kruskal.test Kruskal-Wallis 秩和检验
- R quade.test 四方测试
- R decompose 移动平均线的经典季节性分解
- R plot.stepfun 绘制阶跃函数
- R alias 查找模型中的别名(依赖项)
- R qqnorm 分位数-分位数图
注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Random 2-way Tables with Given Marginals。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。