在R中查找描述数据的分布(qqPlot包冲突)

时间:2014-02-02 18:00:37

标签: r package distribution

我通过算法传递随机值生成数据,我想了解结果的分布。我已经阅读了关于Anderson.Darling / Shapiro-Wilk和QQPlots的讨论  here 因此会考虑自己警告不要作出毫无根据的结论。但是,我使用qqPlot函数在两个不同的库library(car)library(qualityTools)

中得到了不同的结果

事实证明,我没有足够的声誉来发布图片,因此我无法显示我的数据(我从一个更大的集合中抽取了10,000个点,并且不知道上传这样的集合是否合适或可能),或相关的情节。正常情况下,对数正态和伽马(其他一些工作,但明显不适合数据,其他人给出错误,根本没有绘图),它跟随伽玛最接近的一些偏差开始向分位数的上端开始) 。但作为一个例子,请说我检查以下内容:

library(MASS)
library(qualityTools)
library(gplots)
set.seed(123)
x <- rgamma(1000,200,200)
par(mfrow <- c(2, 2))
qqPlot(x,"normal")
qqPlot(x,"log-normal")
qqPlot(x,"gamma")
plot(density(x))

产生一件事,同时清除加载的包并用car包重新启动R

library(car)
library(gplots)
set.seed(123)
x <- rgamma(1000,200,200)
par(mfrow <- c(2, 2))
qqPlot(x,"norm")
qqPlot(x,"lnorm")
qqPlot(x,"gamma")
plot(density(x))

产生别的东西。由于需要形状参数,因此缺少伽玛分布。我假设car版本中没有设置默认值或假设qualityTools中的参数,我可能还会猜测car版本更强大(黑色更少)框)。我下载了qualityToolsMASS的源代码(在运行qqPlot时自动加载)但在任何文件中找不到qqPlot(Windows文件内容搜索tar'd文件夹)。我的最终目标是找到一个模拟我的随机数据(或尽可能接近)的分布。这种方法(尝试分发)是否有效?如果是这样,我如何获得检查其他分布所需的参数(例如形状,dfs等)的初始估计值?任何帮助将不胜感激,谢谢。

0 个答案:

没有答案