在R中，以绘图方式在密度中添加痕迹

Question

使用R中的图包，我想做一个图。实际上，我需要在我的图中再添加一条密度线。我有一个按地理区域划分的一些上市公司的收入信息数据。像这样的东西

head(data)
id    income region 
  1     4556     1
  2     6545     1
  3    65465     2
  4    54555     1
  5    71442     2
  6     5645     6

首先，我用下面的密度图分析了5个和6个地区的收入

reg56<- data[data$region %in% c(5,6) , ]
dens <- with(reg56, tapply(income, INDEX = region, density))
df <- data.frame(
x = unlist(lapply(dens, "[[", "x")),
y = unlist(lapply(dens, "[[", "y")),
cut = rep(names(dens), each = length(dens[[1]]$x))
)

# plot the density 
p<- plot_ly(df, x = x, y = y, color = cut) 

但是，我想要的不止这些。我想加上总收入，即所有地区的收入。我今天晚上试过了

data$aux<- 1
dens2 <- with(data, tapply(income, INDEX = 1, density)) 
df2 <- data.frame(
 x = unlist(lapply(dens2, "[[", "x")),
 y = unlist(lapply(dens2, "[[", "y")),
 cut = rep(names(dens2), each = length(dens2[[1]]$x)) )

p<- plot_ly(df, x = x, y = y, color = cut) 
p<-  add_trace(p, df2, x = x, y = y, color = cut)  
p
Error in FUN(X[[i]], ...) : 
'options' must be a fully named list, or have no names (NULL)

有什么解决办法吗？

Answer 1

因为您没有命名传递给add_trace的参数，所以它将它们解释为对应于默认参数顺序。add_trace的用法是

添加_跟踪（p=last_plot（），…，组，颜色，颜色，符号，符号，大小，数据=NULL，求值=FALSE）

因此，在您提供data.framedf2作为第二个参数的函数调用中，假设它对应于...参数，该参数必须是一个命名列表。您需要指定data=df2，以便add_trace理解这个参数是什么。

让我们生成一些虚拟数据来演示

library(plotly)
set.seed(999)
data <- data.frame(id=1:500, income = round(rnorm(500,50000,15000)), region=sample(6,500,replace=T) )

现在，（在计算了df和df2之后，如您的示例所示）：

p <- plot_ly(df, x = x, y = y, color = cut) %>%
  add_trace(data=df2, x = x, y = y, color = cut)  
p