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))

3つのプロットのスケールは異なります。それらを比較できるように、どうすればそれらをユニバーサルにすることができますか？ありがとう。ここに画像の説明を入力してください

4

1 に答える 1

1

scale_x_continuousErnest A. がコメントしたように、 と を使用して x 軸と y 軸のスケールを手動で変更できますscale_y_continuous。breaks引数を同じ値に設定するだけです。

プロット自体に不確実性バーをプロットするか、95% 信頼区間を含む回帰直線をプロットする方が簡単な場合もあります。

于 2012-12-03T18:40:56.750 に答える