R的百分比分布

时间:2016-09-06 14:34:11

标签: r distribution percentage

我具有以百分比表示的参数(天然气混合物组成)的分布。如何测试这些数据的分布参数(它应该是伽玛,正态分布或对数正态分布)并根据R中的参数生成随机组合?

1 个答案:

答案 0 :(得分:1)

对于CrossValidated,这可能是一个更好的问题,但是:

  • 根据合适的优点,从一系列可能的分布中进行选择通常是个好主意。相反,您应该根据数据的定性特征进行选择,如下所示:

distributions

令人沮丧的是,这张图表实际上并没有为您的数据提供最佳选择(合成,连续,在0和1之间[或0和100]之间),这是一个Beta distribution(虽然有如果您的样本中的值恰好为0或100,则会出现技术问题。

在R:

## some arbitrary data
z <- c(2,8,40,45,56,58,70,89)
## fit (beta values must be in (0,1), not (0,100), so divide by 100)
(m <- MASS::fitdistr(z/100,"beta",start=list(shape1=1,shape2=1)))
## sample 1000 new values
z_new <- 100*rbeta(n=1000,shape1=m$estimate["shape1"],
                   shape2=m$estimate["shape2"])