R语言中的检验方法总结

转自:http://blog.csdn.net/hawksoft/article/details/7679263


1、W检验(Shapiro–Wilk (夏皮罗–威克尔 ) W统计量检验)

     检验数据是否符合正态分布,R函数:shapiro.test().
     结果含义:当p值小于某个显著性水平α(比如0.05)时,则认为

    样本不是来自正态分布的总体,否则则承认样本来自正态分布的总体。


2、K检验(经验分布的Kolmogorov-Smirnov检验)

      R函数:ks.test(),如果P值很小,说明拒绝原假设,表明数据不符合F(n,m)分布。


3、相关性检验:
     R函数:cor.test()

cor.test(x, y,
alternative = c("two.sided", "less", "greater"),
method = c("pearson", "kendall", "spearman"),
exact = NULL, conf.level = 0.95, ...)

结果含义:如果p值很小,则拒绝原假设,认为x,y是相关的。否则认为是不相关的。


4、T检验
用于正态总体均值假设检验,单样本,双样本都可以。      

t.test()

t.test(x, y = NULL,
alternative = c("two.sided", "less", "greater"),
mu = 0, paired = FALSE, var.equal = FALSE,
conf.level = 0.95, ...)

结果意义:P值小于显著性水平时拒绝原假设,否则,接受原假设。具体的假设要看所选择的是双边假设还是单边假设(又分小于和大于)

 

5、正态总体方差检验

t.test(x, y = NULL,
alternative = c("two.sided", "less", "greater"),
mu = 0, paired = FALSE, var.equal = FALSE,
conf.level = 0.95, ...)

结果含义:P值小于显著性水平时拒绝原假设,否则,接受原假设。具体的假设要看所选择的是双边假设还是单边假设(又分小于和大于)

6、二项分布总体假设检验

binom.test(x, n, p = 0.5,
alternative = c("two.sided", "less", "greater"),
conf.level = 0.95)

原假设:p=p0,p<p0,p<p0 计算结果p-值很小,表示拒绝假设,否则为接受假设.


7、Pearson 拟合优度χ2检验
chisq.test(x, y = NULL, correct = TRUE,
p = rep(1/length(x), length(x)), rescale.p = FALSE,
simulate.p.value = FALSE, B = 2000)     

 原假设H0:X符合F分布。
  p-值小于某个显著性水平,则表示拒绝原假设,否则接受原假设。


8、Fisher精确的独立检验:
fisher.test(x, y = NULL, workspace = 200000, hybrid = FALSE,
control = list(), or = 1, alternative = "two.sided",
conf.int = TRUE, conf.level = 0.95)
原假设:X,Y相关。


9、McNemar检验:

mcnemar.test(x, y = NULL, correct = TRUE)

原假设:两组数据的频数没有区别。


10、秩相关检验
cor.test(x, y,
alternative = c("two.sided", "less", "greater"),
method = "spearman", conf.level = 0.95, ...)

原假设:x,y相关.


11、Wilcoxon秩检验
wilcox.test(x, y = NULL,
alternative = c("two.sided", "less", "greater"),
mu = 0, paired = FALSE, exact = NULL, correct = TRUE,
conf.int = FALSE, conf.level = 0.95, ...)
原假设:中位数大于,小于,不等于mu.

  • 5
    点赞
  • 65
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,我会为你编写一个自定义函数来计算符号秩和检验统计量W-plus的各个分位数,并实现例7.6.6的计算过程。 首先,我们需要导入一些必要的库: ```python import numpy as np from itertools import combinations ``` 然后,定义一个自编函数 `wilcoxon_quantiles()` 来计算W-plus的分位数。该函数接受两个参数:样本大小 `n` 和重复次数 `num_repeats`。其,样本大小 `n` 是茆诗松书的例子的样本大小。 ```python def wilcoxon_quantiles(n, num_repeats=10000): w_values = [] for _ in range(num_repeats): # 生成两组随机数(0和1),每组大小为n group1 = np.random.randint(2, size=n) group2 = np.random.randint(2, size=n) # 计算符号秩和检验统计量W-plus W_plus = sum(group1 > group2) w_values.append(W_plus) # 返回分位数 return np.percentile(w_values, [0.1, 1, 5, 10, 25, 50, 75, 90, 95, 99, 99.9]) ``` 接下来,我们可以调用 `wilcoxon_quantiles()` 函数来计算W-plus的各个分位数。例如,对于书例7.6.6的样本大小为10的情况,可以这样调用函数: ```python quantiles = wilcoxon_quantiles(n=10) print("W-plus的分位数:", quantiles) ``` 这样,你就可以得到W-plus的各个分位数。 请注意,由于W-plus的分布不易求出,我们使用了随机模拟的方法来估计其分位数。因此,每次运行代码得到的分位数可能会有轻微的差异。 希望这个编程实现能够帮助到你!如果有任何问题,请随时提问。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值