说我的数据如下:
set.seed(123)
a <- sample(1:10, 50, replace = TRUE)
b <- sample(1:10, 50, replace = TRUE)
c <- sample(3:10, 50, replace = TRUE)
data <- data.frame(a,b,c)
data$AA <- cut(data$a, breaks = c(0,3,6,10), right = TRUE,
labels = c ("typ 1", "typ 2", "typ 3"))
data$BB <- cut(data$b, breaks = c(0,3,6,10), right = TRUE,
labels = c ("typ 1", "typ 2", "typ 3"))
data$CC <- cut(data$c, breaks = c(0,3,6,10), right = TRUE,
labels = c ("typ 1", "typ 2", "typ 3"))
我想绘制a,b
和c
(第1列:第3列)的typ1,2,3的堆积密度分布
所需的输出类似于
根据@MrFlick的建议,我需要转换为长格式。
我的新表应如下所示:
我尝试过使用
data$a <- as.factor(data$a)
data$b <- as.factor(data$b)
data$c <- as.factor(data$c)
data.2 <- melt(data[,-7], idvar=c("a", "b", "c"))
但是这会输出以下消息
Using a, b, c, AA, BB, CC as id variables