0

我正在使用 ggplot2 绘制回归线。我还想用置信区间绘制平均点，以显示它是否显着。例如：

library("ggplot2")
library("gridExtra")
library("epicalc")
options(digits=2)

subset1 <- subset(na.omit(iris), Species == "setosa")
subset2 <- subset(na.omit(iris), Species == "versicolor")
subset3 <- subset(na.omit(iris), Species == "virginica")

meanx <- c(ci(subset1$Sepal.Length)$mean,
               ci(subset2$Sepal.Length)$mean,
               ci(subset3$Sepal.Length)$mean)

meany <- c(ci(subset1$Sepal.Width)$mean,
           ci(subset2$Sepal.Width)$mean,
           ci(subset3$Sepal.Width)$mean)

Species <- factor(c("setosa", "versicolor", "virginica"))
meanmatrix <- as.data.frame(cbind(Species, meanx, meany))

lowerx <- c(ci(subset1$Sepal.Length)$lower95ci,
            ci(subset2$Sepal.Length)$lower95ci,
            ci(subset3$Sepal.Length)$lower95ci)

upperx <- c(ci(subset1$Sepal.Length)$upper95ci,
            ci(subset2$Sepal.Length)$upper95ci,
            ci(subset3$Sepal.Length)$upper95ci)

lowery <- c(ci(subset1$Sepal.Width)$lower95ci,
            ci(subset2$Sepal.Width)$lower95ci,
            ci(subset3$Sepal.Width)$lower95ci)

uppery <- c(ci(subset1$Sepal.Width)$upper95ci,
            ci(subset2$Sepal.Width)$upper95ci,
            ci(subset3$Sepal.Width)$upper95ci)

px <- ggplot(data = meanmatrix, geom = 'blank',
             aes(y = meanx, x = meany,color = factor(Species)))
pbx <- px + 
  geom_point(size = 5) +
  geom_errorbar(aes(ymin=lowerx, ymax=upperx), colour="black", width=.1) +
  scale_color_manual(values = c("#00FFFF", "#FFFF00", "#00FF00")) +
  theme(panel.background = element_rect(fill='white', colour='red'),
        axis.title.x = element_blank(),
        axis.title.y = element_blank(),
        legend.position = "none") +
  coord_flip()


py <- ggplot(data = meanmatrix, geom = 'blank',
             aes(y = meany, x = meany,color = factor(Species)))
pby <- py + 
  geom_point(size = 5) +
  geom_errorbar(aes(ymin=lowery, ymax=uppery), colour="black", width=.1) +
  scale_color_manual(values = c("#00FFFF", "#FFFF00", "#00FF00")) +
  theme(panel.background = element_rect(fill='white', colour='red'),
        axis.title.x = element_blank(),
        axis.title.y = element_blank(),
        legend.position = "none")

p <- ggplot(iris, aes(x = Sepal.Length, y = Sepal.Width,
                    color = factor(Species)))

p0 <- p + 
  scale_color_manual(values = c("#00FFFF", "#FFFF00", "#00FF00")) +
  scale_linetype_manual(breaks = c("0","1"), values = c(1,2), labels = c("male", "female")) +
  geom_smooth(method = "lm",se = FALSE, size = 1.2) +
  theme(panel.background = element_rect(fill='white', colour='red'),
        axis.title.x = element_blank(),
        axis.title.y = element_blank(),
        legend.position = "none")

grid.newpage() 
pushViewport(viewport(layout = grid.layout(nrow=3, ncol=3)))
print(p0,vp = viewport(layout.pos.row = 1:2, layout.pos.col = 2:3))
print(pby,vp = viewport(layout.pos.row = 1:2, layout.pos.col = 1))    
print(pbx,vp = viewport(layout.pos.row = 3, layout.pos.col = 2:3))

三个地块的尺度是不同的。我怎样才能使它们通用，以便我可以比较它们？谢谢。在此处输入图像描述

4

1 回答 1

1

scale_x_continuous就像 Ernest A. 评论的那样，您可以使用和手动更改 x 和 y 轴的比例scale_y_continuous。只需将breaks参数设置为相同的值。

在绘图本身中绘制不确定性条，或绘制包括 95% 置信区间的回归线也可能更容易。

于 2012-12-03T18:40:56.750 回答