如何在R语言绘图中根据指定y值(0.5)求解对应x值

一、你的代码为什么没成功 #

用x[which(exp(-(lambda_wei*x)^kappa_wei) == 0.5)]找不到结果，核心原因是你生成的x是离散序列，生存函数的计算值精确等于0.5的概率几乎为0，which()找不到匹配的索引，自然返回空值。

二、两种可行解法 #

1. 解析解法（推荐，精确计算） #

对于指数和威布尔分布的生存函数，我们可以直接通过数学公式推导x的精确值，完全不需要依赖离散序列：

指数分布（Wykladniczy）

指数分布生存函数为S(x) = exp(-λx)，令S(x)=0.5，解方程可得：

# 计算指数分布的中位生存期（对应y=0.5的x值）
x_exp <- log(2)/lambda_wyk
x_exp

威布尔分布（Weibull）

威布尔分布生存函数为S(x) = exp(-(λx)^k)，令S(x)=0.5，解方程可得：

# 计算威布尔分布的中位生存期（对应y=0.5的x值）
x_weib <- (log(2)^(1/kappa_wei))/lambda_wei
x_weib

2. 数值插值解法（适用于无解析解的情况） #

如果遇到无法推导解析公式的函数，可以用插值方法找近似x值，步骤如下：

# 计算威布尔生存函数在x序列上的所有y值
y_weib <- exp(-(lambda_wei*x)^kappa_wei)
# 用approx函数直接根据目标y值反求x，method指定线性插值
x_interp <- approx(y_weib, x, xout=0.5, method="linear")$y
x_interp

注：approx()函数中，前两个参数是已知的y-x对应序列，xout是你要找的目标y值，返回的$y就是对应的x值。

三、验证结果（可选） #

可以把计算出的x值画到原图上，验证是否正确：

# 添加中位生存期的垂直线
abline(v=x_exp, col='red', lty=2)
abline(v=x_weib, col='blue', lty=2)
# 添加文字标注
text(x_exp, 0.5, paste("指数中位:", round(x_exp,1)), col='red', pos=4)
text(x_weib, 0.5, paste("威布尔中位:", round(x_weib,1)), col='blue', pos=4)

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