R 为什么我的图中省略了负值的幂

发布于11月20日

我想描绘一个图表，它是一些值的幂，以一种连续的方式.

f <- function(x) {
  return (x^(1/3))
}  

ggplot(data.frame(x=seq(-5,5,length.out=10)), aes(x)) +
  stat_function(fun=f)

这只显示x > 0的值，但如果我使用在线图形工具Desmos，它也会显示x < 0的值.在交互式会话中，-3^(1/3)返回-1.44225.

为什么R省略它，我怎么也能描绘负的x值呢？

推荐答案

负数的分数次方是复数.如果要在绘图中包含这些组件，可以分别显示虚部和实部:

real_part <- function(x) {
  Re(complex(length(x), x, 0)^(1/3))
}

imaginary_part <- function(x) {
  Im(complex(length(x), x, 0)^(1/3))
}

为了清楚起见，我们可以使用标签进行绘制:

library(geomtextpath)
#> Loading required package: ggplot2

ggplot(data.frame(x = seq(-5, 5, length.out = 10)), aes(x)) +
  stat_function(fun = real_part, aes(label = "Real part"), geom = "textpath",
                hjust = 0.8, vjust = -0.2, size = 6) +
  stat_function(fun = imaginary_part, color = "red", linetype = 2, hjust = 0.2,
                aes(label = "Imaginary part"), geom = "textpath", vjust = -0.2,
                size = 6) +
  theme_minimal(base_size = 20)