r语言plot画图叠加

在R语言中，plot函数是最常用的画图函数之一。它可以画出散点图、线图、柱状图等各种类型的图形。但是，有时候我们需要在同一张图上叠加多个图形，比如在一张散点图上叠加线性回归线，或者在一张柱状图上叠加置信区间。这时候，我们就需要用到plot函数的一些高级技巧，来实现多个图形的叠加。

叠加散点图和线性回归线

首先，我们来看一个例子。假设我们有一组数据，其中x和y是两个变量，我们想要画出它们的散点图，并在同一张图上叠加一条线性回归线。我们可以按照以下步骤来实现：

1. 生成数据

我们先生成一个随机的数据集，包含100个x和y的值：

``` r

set.seed(123)

x <- rnorm(100)

y <- 2*x + rnorm(100)

```

2. 画散点图

我们使用plot函数来画出散点图，其中x轴是x变量，y轴是y变量：

``` r

plot(x, y, main="Scatterplot Example", xlab="x", ylab="y")

```

这会生成一个简单的散点图，如下图所示：


3. 添加线性回归线

我们可以使用abline函数来添加一条线性回归线。这个函数可以接受两个参数，第一个是斜率，第二个是截距。我们可以使用lm函数来拟合一条线性回归线，并从中获得斜率和截距：

``` r

lm.fit <- lm(y ~ x)

abline(lm.fit, col="red")

```

这会在原来的散点图上添加一条红色的线性回归线，如下图所示：


叠加柱状图和置信区间

除了在散点图上叠加线性回归线，我们还可以在柱状图上叠加置信区间。假设我们有一个数据集，其中x是一个因子变量，y是一个数值变量，我们想要画出它们的柱状图，并在同一张图上叠加y的置信区间。我们可以按照以下步骤来实现：

1. 生成数据

我们先生成一个随机的数据集，包含5个因子水平和它们对应的y的值：

``` r

set.seed(123)

x <- factor(c("A", "B", "C", "D", "E"))

y <- c(3.2, 4.1, 3.5, 5.7, 4.8)

```

2. 画柱状图

我们使用barplot函数来画出柱状图，其中x轴是x变量，y轴是y变量：

``` r

bp <- barplot(y, main="Barplot Example", xlab="x", ylab="y")

```

这会生成一个简单的柱状图，如下图所示：


3. 添加置信区间

我们可以使用segments函数来添加置信区间。这个函数可以接受四个参数，分别是x1、y1、x2、y2，表示一条线段的起点和终点。我们可以使用qt函数来计算置信区间的上下界，并从中获得y的置信区间：

``` r

alpha <- 0.05

n <- length(y)

se <- sd(y) / sqrt(n)

t <- qt(1-alpha/2, n-1)

ci <- t * se

segments(bp, y + ci, bp, y - ci, col="red", lwd=2)

```

这会在原来的柱状图上添加红色的置信区间，如下图所示：


叠加多个图形

除了在同一张图上叠加两个图形外，我们还可以在同一张图上叠加多个图形。比如，在一张散点图上叠加两条线性回归线，或者在一张柱状图上叠加多组数据的置信区间。实现起来也很简单，只需要按照上面的步骤，依次添加每个图形即可。