首页 > 解决方案 > 在 R 中为我的偏度测试设置可接受的区域

问题描述

我正在编写下面的函数,让我对带有复制test of skewness的样本向量进行 a。(10, 20, 50, 100)1000

library(moments)
out <- t(sapply(c(10, 20, 50, 100), function(x)
  table(replicate(1000, skewness(rgamma(n = x, shape = 3, rate = 0.5))) > 2)))
row.names(out) <- c(10, 20, 50, 100)
out

我的条件

我拒绝零假设的条件是统计数据必须满足两 (2) 个条件:

  1. 小于 -2
  2. 或大于+2。

我有的

但在我的R功能中,我只能描述第二个条件。

我想要的是

如何在我的函数中同时包含第一个和第二个条件?

标签: r

解决方案


也许添加abs将是满足这两个条件的最简单方法

out <- t(sapply(c(10, 20, 50, 100), function(x)
  table(abs(unlist(replicate(1000, skewness(rgamma(n = x, shape = 3, rate = 0.5))))) > 2)))
row.names(out) <- c(10, 20, 50, 100)
out

推荐阅读