2020-04-20-In-class exercise-3

重新命名變項

library(datasets)
library(ggplot2)
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(tidyr)
dta <- datasets::USPersonalExpenditure
colnames(dta) <- c("year1940","year1945","year1950","year1955","year1960")
rownames(dta) <- c("Food_Tobacco","Household_Operation","Medical_Health","Personal_Care","Private_Education")
dta <- as.data.frame(dta)

取log，計算平均，嘗試畫圖

dat <- stack(log10(dta))

mean(dat$values)

## [1] 0.9183627

p0 <- ggplot(dat) + 
aes(x=values, y=ind) + 
  geom_point(size=rel(3)) + 
 geom_segment(aes(xend=dat$values, yend=dat$ind,y=0, x=0.9183627)) + 
  facet_grid(. ~ ind)

p0

## Warning: Use of `dat$values` is discouraged. Use `values` instead.

## Warning: Use of `dat$ind` is discouraged. Use `ind` instead.

重新一個一個慢慢畫，以year1945作為Y軸，以總體log平均為x軸 *目前想不到更好的方法

dta.1 <- (log(dta))

##畫出1945支出
p1 <- ggplot(dta.1) + 
aes(x=year1940, y=year1945) + 
  geom_point(size=rel(3)) + 
 geom_segment(aes(xend=year1940, yend=year1945,x=0.91)) + 
  labs(x="1940", 
       y="Personal Expenditure")+ 
scale_x_continuous(limits=c(-2, 5), 
                    breaks=seq(-1.1, 5, 1))

##畫出1945支出
p2 <- ggplot(dta.1) + 
aes(x=year1945, y=year1945) + 
  geom_point(size=rel(3)) + 
 geom_segment(aes(xend=year1945, yend=year1945,x=0.91)) + 
  labs(x="1945", 
       y="Personal Expenditure") + 
scale_x_continuous(limits=c(-2, 5), 
                    breaks=seq(-1.1, 5, 1))

##畫出1950支出
p3 <- ggplot(dta.1) + 
aes(x=year1950, y=year1945) + 
  geom_point(size=rel(3)) + 
 geom_segment(aes(xend=year1950, yend=year1945,x=0.91)) + 
  labs(x="1950", 
       y="Personal Expenditure")+ 
scale_x_continuous(limits=c(-2, 5), 
                    breaks=seq(-1.1, 5, 1))

##畫出1955支出
p4 <- ggplot(dta.1) + 
aes(x=year1955, y=year1945) + 
  geom_point(size=rel(3)) + 
 geom_segment(aes(xend=year1955, yend=year1945,x=0.91)) + 
  labs(x="1955", 
       y="Personal Expenditure")+ 
scale_x_continuous(limits=c(-2, 5), 
                    breaks=seq(-1.1, 5, 1))

##畫出1960支出
p5 <- ggplot(dta.1) + 
aes(x=year1960, y=year1945) + 
  geom_point(size=rel(3)) + 
 geom_segment(aes(xend=year1960, yend=year1945,x=0.91)) + 
  labs(x="1960", 
       y="Personal Expenditure")+ 
scale_x_continuous(limits=c(-2, 5), 
                    breaks=seq(-1.1, 5, 1))

##檔案合併
library(gridExtra)

## 
## Attaching package: 'gridExtra'

## The following object is masked from 'package:dplyr':
## 
##     combine

grid.arrange(p1, p2, p3, p4, p5, ncol=3)