0924_HOMEWORK1

data management and graphics package

install.packages("mlmRev",repos = "http://cran.us.r-project.org")

## 將程式套件安載入 'C:/Users/a0905/AppData/Local/R/win-library/4.2'
## (因為 'lib' 沒有被指定)

## 程式套件 'mlmRev' 開啟成功，MD5 和檢查也透過
## 
## 下載的二進位程式套件在
##  C:\Users\a0905\AppData\Local\Temp\RtmpMpQO7J\downloaded_packages 裡

library(mlmRev)

## 載入需要的套件：lme4

## 載入需要的套件：Matrix

library(tidyverse)

## ── Attaching packages
## ───────────────────────────────────────
## tidyverse 1.3.2 ──

## ✔ ggplot2 3.3.6      ✔ purrr   0.3.4 
## ✔ tibble  3.1.8      ✔ dplyr   1.0.10
## ✔ tidyr   1.2.1      ✔ stringr 1.4.1 
## ✔ readr   2.1.2      ✔ forcats 0.5.2 
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ tidyr::expand() masks Matrix::expand()
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ✖ tidyr::pack()   masks Matrix::pack()
## ✖ tidyr::unpack() masks Matrix::unpack()

data(Gcsemv, package="mlmRev")
dta <- Gcsemv

compute averages by school

dta_a <- dta %>%
        group_by(school) %>%
        summarize(ave_written = mean(written, na.rm=TRUE),
                  ave_course = mean(course, na.rm=TRUE))

superimpose two plots

ggplot(data=dta, aes(x=written, y=course)) +
 geom_point(color="skyblue") +
 stat_smooth(method="lm", formula=y ~ x, se=F, col="skyblue") +
 geom_point(data=dta_a, aes(ave_written, ave_course), color="steelblue") +
 stat_smooth(data=dta_a, aes(ave_written, ave_course),
             method="lm", formula= y ~ x, se=F, color="steelblue") +
 labs(x="written", 
      y="course") +
 theme_bw()

## Warning: Removed 382 rows containing non-finite values (stat_smooth).

## Warning: Removed 382 rows containing missing values (geom_point).

##