ggplot小提琴情节不绘图

时间:2019-01-21 09:09:20

标签: r ggplot2

我尝试使用带有代码的R包ggplot2绘制小提琴图

norm2 = function(v) return(sqrt(sum(v*v)))
myfct = function(d) {
  vec_length = Inf
  while (vec_length > 1){
    vec_length = norm2(runif(n=d,min=-1,max=1))
  }
  return(vec_length)
}

df = data.frame(x = rep.int(1:5, 2))
df$vec_length = sapply(df$x, myfct)
ggplot(df, aes(factor(x),vec_length)) + geom_violin(trim=FALSE)

但我明白了

Warning:
In max(data$density) :
  no non-missing argument for max; return -Inf

我的情节是

enter image description here

我做错了什么?

1 个答案:

答案 0 :(得分:2)

每个vec_length的数据中只有两个x(y)。这是一个“特殊情况”,小提琴会缩成一条线。在这种情况下,人们可能也将geom_violin()实施为geom_line(),但这并不是这样实现的:

library(ggplot2)
ggplot(df1, aes(factor(x), vec_length)) + geom_line()

enter image description here

要拉小提琴,您至少需要三个y值:

df2 <- data.frame(x=rep.int(1:5, 3))
df2$vec_length <- sapply(df2$x, myfct)
ggplot(df2, aes(factor(x), vec_length)) + geom_violin(trim=FALSE)

enter image description here


数据

library("SpatioTemporal")
set.seed(42)
myfct <- function(d) {
  vec_length <- Inf
  while (vec_length > 1){
    vec_length <- norm2(runif(n=d, min=- 1, max=1))
  }
  return(vec_length)
}

df1 <- data.frame(x=rep.int(1:5, 2))
df1$vec_length <- sapply(df1$x, myfct)