我具有以百分比表示的参数(天然气混合物组成)的分布。如何测试这些数据的分布参数(它应该是伽玛,正态分布或对数正态分布)并根据R中的参数生成随机组合?
答案 0 :(得分:1)
对于CrossValidated,这可能是一个更好的问题,但是:
令人沮丧的是,这张图表实际上并没有为您的数据提供最佳选择(合成,连续,在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"])