计算概率密度函数（PDF）的最大值：代码错误排查求助

让我们一步步来梳理你的问题和代码里的问题：

核心问题拆解 #

• 函数定义完全偏离需求
  你要计算的是概率密度函数f(x) = 3(1-x)^2的最大值，但你定义的integral函数却是3*x*(1-x)^2——这和你给出的PDF完全不是同一个函数，相当于你在计算一个无关函数的最大值，结果自然不对。

• 调用max()的方式错误
  就算函数定义正确，你写的max(integral(x))里的x根本没有被定义，R无法识别这个变量。而且max()是用来求一组离散数值的最大值，如果你想找连续函数的极值，直接这么用是行不通的，得用专门的优化工具，或者先做数学分析。

正确解决方法 #

方法1：数学分析（最直接高效）

你的PDF是f(x) = 3(1-x)^2，定义域x∈[0,1]：

1. 对函数求导：f’(x) = -6(1-x)
2. 在x∈[0,1)区间内，f’(x) ≤ 0，说明函数在定义域内单调递减
3. 因此最大值出现在定义域的左端点x=0处，代入得f(0)=3*(1-0)^2=3

方法2：用R代码验证

如果你想用代码确认结果，正确的步骤如下：

# 正确定义目标PDF函数
pdf_func <- function(x) {
  3 * (1 - x)^2
}

# 使用optimize工具找最大值：因为optimize默认找最小值，所以我们优化函数的负值
max_result <- optimize(function(x) -pdf_func(x), interval = c(0, 1))

# 输出结果
cat("最大值点：", max_result$minimum, "\n")
cat("最大值：", -max_result$objective, "\n")

运行这段代码会得到最大值点为0，最大值为3，和数学分析的结果一致。

另外补充：如果你之前混淆了PDF和CDF（累积分布函数），原PDF的CDF是F(x) = 1 - (1-x)^3，它在x∈[0,1]的最大值是1（当x=1时），但这和你最初的需求不符，所以要明确自己的计算目标。

内容的提问来源于stack exchange，提问作者The Statistician Magician