如何使用 ggplot plot 改变 R 中一个变量的值的方程

Question

我有一个方程，它给出的值介于 1 和 0 之间，并且想在 R 到 plot 中使用 ggplot 如何改变方程中的一个变量（年龄）会改变 Z78E6221F6393D1356681DB398CZF1 值。 我希望 X 轴是从 0 到 100 的变量范围，并且 plot 是 X 轴上方程的 output 值。

我像这样设置传递给方程的变量：

age = 100
d <- (365/4)*age
k <- 5
N <-8
m <- 15000000
Nm <- N*m
p <- 1 - (1 - (1 - (1 - u)^d)^k)^Nm

我确信在 R 中必须有一种简单的方法来给出一系列年龄值以传递给方程，并且 plot 将 P 的结果值作为一条线。 但我不确定如何最好地了解它。 我是否需要事先填写范围内的值表，还是可以将范围本身放在 ggplot 命令中？

提前致谢

Answer 1

您可以使用stat_function 。

p <- function(age, u, k, N, m)
{
  1 - (1 - (1 - (1 - u)^((365/4) * age))^k)^(N * m)
}

ggplot(data.frame(age = 1:100), aes(x = age)) + 
  stat_function(fun = p, args = list(u = 0.000005, k = 5, N = 8, m = 15000000)) +
  labs(y = "p")

问题中没有定义变量u ，但这是一条逻辑曲线，将u设置为 0.000005 会将范围的中间值置于 p = 0.5，因此u可能接近该值。