如何在R中使用自定义函数绘制指定PDF与CDF？

核心逻辑 #

只要你清楚目标PDF/CDF的数学表达式，完全可以在R里实现自定义绘制。核心就是把数学公式转成R可执行的函数，再结合绘图工具生成可视化结果。

步骤1：定义自定义PDF/CDF函数 #

无参数的简单分布示例 #

假设你看到的PDF是f(x)=2x（x∈[0,1]），对应的CDF是F(x)=x²（x∈[0,1]），可以这样定义函数：

# 自定义PDF：不在定义域内返回0
my_pdf <- function(x) {
  ifelse(x >= 0 & x <= 1, 2*x, 0)
}

# 自定义CDF：处理边界情况
my_cdf <- function(x) {
  ifelse(x < 0, 0, ifelse(x <=1, x^2, 1))
}

带参数的分布示例 #

如果目标分布有形状/尺度/位置参数（比如类似Beta分布的参数），直接把参数作为函数输入即可：

# 带形状参数a的PDF：f(x;a)=a*x^(a-1)，x∈[0,1]，a>0
my_pdf_param <- function(x, a) {
  ifelse(x >=0 & x <=1, a*x^(a-1), 0)
}

# 对应CDF
my_cdf_param <- function(x, a) {
  ifelse(x <0, 0, ifelse(x <=1, x^a, 1))
}

步骤2：生成绘图数据 #

先覆盖分布的有效范围生成x序列，再计算对应的PDF/CDF值：

# 生成x序列，包含定义域外的部分以展示边界行为
x_vals <- seq(-0.5, 1.5, by=0.01)

# 计算无参数的PDF/CDF值
pdf_vals <- sapply(x_vals, my_pdf)
cdf_vals <- sapply(x_vals, my_cdf)

# 计算带参数的PDF/CDF值（比如a=2）
pdf_vals_param <- sapply(x_vals, my_pdf_param, a=2)
cdf_vals_param <- sapply(x_vals, my_cdf_param, a=2)

步骤3：绘图（两种常用方式） #

Base R 绘图 #

# 绘制自定义PDF
plot(x_vals, pdf_vals, type="l", lwd=2, col="blue", 
     xlab="x", ylab="PDF", main="自定义PDF曲线")

# 绘制自定义CDF
plot(x_vals, cdf_vals, type="l", lwd=2, col="red", 
     xlab="x", ylab="CDF", main="自定义CDF曲线")

# 绘制带参数的PDF
plot(x_vals, pdf_vals_param, type="l", lwd=2, col="green", 
     xlab="x", ylab="PDF", main="带形状参数的自定义PDF (a=2)")

ggplot2 绘图 #

如果偏好ggplot2风格：

library(ggplot2)

# 整理成数据框
df <- data.frame(x = x_vals,
                 pdf = pdf_vals,
                 cdf = cdf_vals)

# 绘制PDF
ggplot(df, aes(x=x, y=pdf)) +
  geom_line(color="blue", linewidth=1) +
  labs(x="x", y="PDF", title="自定义PDF曲线") +
  theme_minimal()

# 绘制CDF
ggplot(df, aes(x=x, y=cdf)) +
  geom_line(color="red", linewidth=1) +
  labs(x="x", y="CDF", title="自定义CDF曲线") +
  theme_minimal()

参数转换的实用提示 #

如果对参数转换感到困惑，记住这两点：

• 先明确每个参数的数学含义（比如是控制曲线位置的位置参数、控制拉伸程度的尺度参数，还是控制形态的形状参数）
• 定义函数时把参数作为输入，测试极端参数值验证输出是否符合预期（比如参数取边界值时，PDF是否非负、CDF是否在0-1之间）

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