我尝试使用带有代码的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
我的情节是
我做错了什么?
答案 0 :(得分:2)
每个vec_length
的数据中只有两个x
(y)。这是一个“特殊情况”,小提琴会缩成一条线。在这种情况下,人们可能也将geom_violin()
实施为geom_line()
,但这并不是这样实现的:
library(ggplot2)
ggplot(df1, aes(factor(x), vec_length)) + geom_line()
要拉小提琴,您至少需要三个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)
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)