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


R SignRank Wilcoxon 有符号秩统计量的分布


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

说明

从大小为 n 的样本获得的 Wilcoxon Signed Rank 统计分布的密度、分布函数、分位数函数和随机生成。

用法

dsignrank(x, n, log = FALSE)
psignrank(q, n, lower.tail = TRUE, log.p = FALSE)
qsignrank(p, n, lower.tail = TRUE, log.p = FALSE)
rsignrank(nn, n)

参数

x, q

分位数向量。

p

概率向量。

nn

观察次数。如果是 length(nn) > 1 ,则长度被视为所需的数量。

n

样本中的观测值数量。正整数或此类整数的向量。

log, log.p

逻辑性;如果为 TRUE,则概率 p 以 log(p) 形式给出。

lower.tail

逻辑性;如果为 TRUE(默认值),则概率为 ,否则为

细节

该分布的获得如下。令 x 为来自关于原点对称的连续分布的大小为 n 的样本。那么 Wilcoxon 符号秩统计量就是 x[i] 为正的绝对值 x[i] 的秩之和。该统计量取 之间的值,其均值和方差分别为

如果前两个参数中的任何一个是向量,则使用回收规则对两者的所有组合进行计算,直到较长向量的长度。

dsignrank 给出密度,psignrank 给出分布函数,qsignrank 给出分位数函数,rsignrank 生成随机偏差。

结果的长度由 rsignranknn 确定,并且是其他函数的数值参数长度的最大值。

nn 之外的数字参数将被回收到结果的长度。仅使用逻辑参数的第一个元素。

例子

require(graphics)

par(mfrow = c(2,2))
for(n in c(4:5,10,40)) {
  x <- seq(0, n*(n+1)/2, length.out = 501)
  plot(x, dsignrank(x, n = n), type = "l",
       main = paste0("dsignrank(x, n = ", n, ")"))
}

作者

Kurt Hornik; efficiency improvement by Ivo Ugrina.

也可以看看

wilcox.test 根据数据计算统计量,查找 p 值等。

Distributions 用于标准分布,包括 dwilcox 用于两样本 Wilcoxon 秩和统计的分布。

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Distribution of the Wilcoxon Signed Rank Statistic。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。