当我们facet在ggplot中使用选项时,我们会在标题周围找到一个漂亮的灰色框(3,4,5).
library(ggplot2)
data(mtcars)
ggplot(mtcars, aes(cyl)) + geom_bar() + facet_wrap(~gear) + theme_bw()
当我们不使用该facet选项时,我们如何在图表标题周围放置一个类似的灰色框？
ggplot(mtcars, aes(cyl)) + geom_bar() + theme_bw() +
ggtitle("How do you put a grey box around me??")
这个问题是从这里的问题而来的 ，因为我们想要绘制一个多变量模型。
假设我们要绘制包发行说明survreg中示例中的生存函数
模型在哪里
survreg(Surv(time, status) ~ ph.ecog + age + strata(sex), lung)
据我了解，以下问题似乎解决了该问题，但没有提供数据集并且仅使用一个变量。在此示例中，我们有两个连续变量和一个因子
根据以下问题，我们创建一些虚拟数据。然后将其转换为可以理解的格式ggplot2，并生成一个简单的图表，显示var随时间的变化。
test_data <-
data.frame(
var0 = 100 + c(0, cumsum(runif(49, -20, 20))),
var1 = 150 + c(0, cumsum(runif(49, -10, 10))),
var2 = 120 + c(0, cumsum(runif(49, -5, 10))),
date = seq(as.Date("2002-01-01"), by="1 month", length.out=100)
)
#
library("reshape2")
library("ggplot2")
#
test_data_long <- melt(test_data, id="date") # convert to long format
ggplot(data=test_data_long,
aes(x=date, y=value, colour=variable)) +
geom_line() + theme_bw()
我想在同一张图中绘制三个的平均值var，并显示平均值的置信区间。可能有 +-1SD。为此，我认为stat_summary()可以使用该函数，如此处和此处所述。
通过添加以下任一命令，我没有获得平均值，也没有获得置信区间。任何建议将不胜感激。
stat_summary(fun.data=mean_cl_normal)
#stat_summary(fun.data ="mean_sdl", mult=1, geom …我正在尝试仅提取每个行b是唯一值时的行a。
这是一些示例数据
a <- c(1,1,2,2,3,3,4,4,5,5,5,6,6,7,7,8,8,9,9,9,9,9,10,10,10)
b <- c(1,2,1,1,5,5,6,1,1,1,3,2,2,1,1,2,3,1,2,3,4,4,1,2,2)
df1 <- data.frame(a, b)
并使用dplyr包
library(dplyr)
Unique <- df1 %>%
group_by(a) %>%
filter(n_distinct(b))
所需的输出应该是数据帧长度 18