getwd()
## [1] "/Users/coop2711/Dropbox/Works/Research_Data/Doorae/R"
library(pander)
library(knitr)
library(RCurl)
## Loading required package: bitops
library(data.table)
## data.table 1.9.6 For help type ?data.table or https://github.com/Rdatatable/data.table/wiki
## The fastest way to learn (by data.table authors): https://www.datacamp.com/courses/data-analysis-the-data-table-way
# panderOptions()
panderOptions("table.split.table", Inf)
options(width=180)
githubURL <- "https://github.com/kwlee58/Doorae/R/doorae_new1015.rda"
# load(url(githubURL))
# tmpfile <- tempfile()
# download.file(githubURL, destfile=tmpfile, method="curl")
# url.data <- tmpfile
# load(url.data)
# getURLContent(githubURL, binary=TRUE)
load("~/Dropbox/Works/Research_Data/Doorae/R/doorae_new.rda")
(tapply.ys <- tapply(doorae$"업력", doorae$"공동체상태", summary))
## $경영개선
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.40 2.65 4.60 5.00 6.20 12.50
##
## $창업준비
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.200 1.000 1.500 1.694 1.600 11.800
tapply.ys[[1]]
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.40 2.65 4.60 5.00 6.20 12.50
tapply.ys[1]
## $경영개선
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 1.40 2.65 4.60 5.00 6.20 12.50
str(tapply.ys[1])
## List of 1
## $ 경영개선:Classes 'summaryDefault', 'table' Named num [1:6] 1.4 2.65 4.6 5 6.2 12.5
## .. ..- attr(*, "names")= chr [1:6] "Min." "1st Qu." "Median" "Mean" ...
str(tapply.ys[[1]])
## Classes 'summaryDefault', 'table' Named num [1:6] 1.4 2.65 4.6 5 6.2 12.5
## ..- attr(*, "names")= chr [1:6] "Min." "1st Qu." "Median" "Mean" ...
names(tapply.ys[[1]])
## [1] "Min." "1st Qu." "Median" "Mean" "3rd Qu." "Max."
names(tapply.ys[[1]]) <- c("최소값", "제1사분위수", "중위수", "산술평균", "제3사분위수", "최대값")
names(tapply.ys[[2]]) <- c("최소값", "제1사분위수", "중위수", "산술평균", "제3사분위수", "최대값")
tapply.ys
## $경영개선
## 최소값 제1사분위수 중위수 산술평균 제3사분위수 최대값
## 1.40 2.65 4.60 5.00 6.20 12.50
##
## $창업준비
## 최소값 제1사분위수 중위수 산술평균 제3사분위수 최대값
## 0.200 1.000 1.500 1.694 1.600 11.800
pander(tapply.ys)
doorae.ag <- aggregate(업력~공동체상태, data=doorae, summary)
names(aggregate(업력~공동체상태, data=doorae, summary))
## [1] "공동체상태" "업력"
xtabs(업력~공동체상태, doorae.ag)
##
## 공동체상태 Min. 1st Qu. Median Mean 3rd Qu. Max.
## 경영개선 1.400 2.650 4.600 5.000 6.200 12.500
## 창업준비 0.200 1.000 1.500 1.694 1.600 11.800
doorae.xtab <- xtabs(업력~공동체상태, doorae.ag)
str(doorae.xtab)
## xtabs [1:2, 1:6] 1.4 0.2 2.65 1 4.6 ...
## - attr(*, "dimnames")=List of 2
## ..$ 공동체상태: chr [1:2] "경영개선" "창업준비"
## ..$ : chr [1:6] "Min." "1st Qu." "Median" "Mean" ...
## - attr(*, "class")= chr [1:2] "xtabs" "table"
## - attr(*, "call")= language xtabs(formula = 업력 ~ 공동체상태, data = doorae.ag)
rownames(doorae.xtab)
## [1] "경영개선" "창업준비"
colnames(doorae.xtab)
## [1] "Min." "1st Qu." "Median" "Mean" "3rd Qu." "Max."
dimnames(doorae.xtab)<-list("공동체상태"=rownames(doorae.xtab), "업력"=names(tapply.ys[[1]]))
str(doorae.xtab)
## xtabs [1:2, 1:6] 1.4 0.2 2.65 1 4.6 ...
## - attr(*, "dimnames")=List of 2
## ..$ 공동체상태: chr [1:2] "경영개선" "창업준비"
## ..$ 업력 : chr [1:6] "최소값" "제1사분위수" "중위수" "산술평균" ...
## - attr(*, "class")= chr [1:2] "xtabs" "table"
## - attr(*, "call")= language xtabs(formula = 업력 ~ 공동체상태, data = doorae.ag)
kable(xtabs(업력~공동체상태, doorae.ag), caption="<공동체상태와 업력>")
<공동체상태와 업력>
| 경영개선 |
1.4 |
2.65 |
4.6 |
5.000 |
6.2 |
12.5 |
| 창업준비 |
0.2 |
1.00 |
1.5 |
1.694 |
1.6 |
11.8 |
kable(doorae.xtab, caption="<공동체상태와 업력>")
<공동체상태와 업력>
| 경영개선 |
1.4 |
2.65 |
4.6 |
5.000 |
6.2 |
12.5 |
| 창업준비 |
0.2 |
1.00 |
1.5 |
1.694 |
1.6 |
11.8 |