求助：R中Weibull分布参数计算输入mu=0、sd=1返回NA值

首先，你遇到的NA结果是因为输入的mu=0不符合两参数Weibull分布（loc=0）的基本性质：
当loc=0时，Weibull分布的所有取值都大于0，它的均值公式为 scale * Gamma(1 + 1/shape)，其中Gamma函数在自变量大于0时结果恒为正，同时shape和scale参数本身也必须大于0，因此Weibull分布的均值必然是正数，你输入mu=0时，函数找不到符合条件的shape和scale，所以返回NA。

mixdist::weibullpar函数的输入参数要求 #

• 当设置loc=0时：
  • mu必须是正数，这是由Weibull分布的支撑集（取值范围）决定的
  • sd必须是正数，且要满足Weibull分布的方差约束：方差 = scale²[Gamma(1+2/shape) - (Gamma(1+1/shape))²]，不是所有正数mu和sd的组合都能对应有效的shape和scale，但大部分合理的正数组合可以正常计算
• 如果使用三参数Weibull（loc≠0），loc需要小于mu（因为分布的取值范围是(loc, +∞)，均值必须大于loc）

正确使用示例 #

# 合法的正数均值和标准差
mu = 1
sd = 2
sp = mixdist::weibullpar(mu, sd, loc = 0)$shape
sc = mixdist::weibullpar(mu, sd, loc = 0)$scale

# 查看结果
sp
sc

给初学者的提示 #

• 先熟悉Weibull分布的基本定义：两参数Weibull（loc=0）是定义在正实数域上的分布，所有统计量（均值、中位数等）都必须为正
• 可以手动推导均值和方差公式，结合Gamma函数的性质，就能明白为什么mu不能为0
• 用R内置函数验证：比如生成Weibull随机数，计算其均值，结果肯定是正数

# 生成1000个Weibull随机数，计算均值
set.seed(123)
wb_data = rweibull(1000, shape=2, scale=1)
mean(wb_data) # 结果约为0.886，是正数

内容的提问来源于stack exchange，提问作者Timon