20210928_W1_Exercises

hist(rnorm(1000))

install.packages("mlmRev", repos='https://cloud.r-project.org')

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

## package 'mlmRev' successfully unpacked and MD5 sums checked
## 
## The downloaded binary packages are in
##  C:\Users\USER\AppData\Local\Temp\Rtmp6vywJp\downloaded_packages

library(mlmRev)

## 載入需要的套件：lme4

## 載入需要的套件：Matrix

install.packages("lme4")

## Warning: 正在使用 'lme4' 這個程式套件，因此不會被安裝

install.packages("Matrix")

## Warning: 正在使用 'Matrix' 這個程式套件，因此不會被安裝

library(lattice)
dta <- Chem97
dta$age = dta$age + 18.5
View(dta)

histogram(~gcsescore | factor(score), data=dta)

tp1 <- histogram(~gcsescore | factor(score), data=dta)

# (col.position, row.position, n.col, n.row)

plot(tp1, split=c(1,1,1,1))

tp2 <- densityplot(~ gcsescore, 
                   groups=score, 
                   data=dta, 
                   plot.points=F, 
                   auto.key=list(space="right",title="score"))

# put the second trellis obj at row 2
plot(tp2, split=c(1,1,1,1), newpage=F)

tp3 <- bwplot(~gcsescore |factor(score), data = dta)

#try boxplot
plot( tp3, split = c(1, 1, 1, 1), newpage = F)

plot(tp1, split=c(1,1,1,3))
plot(tp2, split=c(1,2,1,3), newpage=F)
plot( tp3, split = c(1, 3, 1, 3), newpage = F)