boxp <- function(x,ylab){
five <- fivenum(x)
quart <- five[2:4]
fence.1 <- quart[1] - 3*(quart[2]-quart[1]) #1
fence.2 <- quart[3] + 3*(quart[3]-quart[2]) #2
fence.1.adj <- min(x[x > fence.1])
fence.2.adj <- max(x[x < fence.2])
y.range <- max(x)-min(x)
y.lim <- c(max(x)+0.1*y.range,min(x)-0.1*y.range)
outlier <- x[x < fence.1 | x > fence.2]
x.1 <- c(0,0)
par(xaxt="n")
plot(y.lim ~ x.1, type="n", xlim=c(-2,2), xlab="", ylab=ylab)
lines(c(-1,1),c(quart[3],quart[3]))
lines(c(-1,1),c(quart[2],quart[2]),lwd=3)
lines(c(-1,1),c(quart[1],quart[1]))
lines(c(-1,-1),c(quart[1],quart[3]))
lines(c(1,1),c(quart[1],quart[3]))
lines(c(0,0),c(quart[3],fence.2.adj))
lines(c(0,0),c(quart[1],fence.1.adj))
for (i in 1:length(outlier)) points(0,outlier[i],cex=1)
}
새로운 boxplot을 짜는 코드입니다. #1과 #2부분을 바꿔줌으로서 내가 원하는 boxplot을 만들 수 있습니다.
아래는 기존 플롯과 새로운 플롯의 비교 그림입니다.
이런 식으로 얼마든 바꿀 수 있습니다. 참 쉽죠?
2. 파이차트
원형차트를 말하는 건데요 누가 봐도 보기 쉽게 만들 수 있습니다.
기준이 되는 데이터는 제 하루 일과표로 해보겠습니다.
4시간 잠, 1시간 지하철, 학교에서 8시간, 1시간 지하철, 2시간 놀기, 4시간 공부, 2시간 오유 감상, 2시간 애니감상 으로 만들어보죠
\
위의 코드를 분석해봅시다. c(2,2,4,1,8,1,2,4)는 순서대로 넣을 숫자를(데이터를) 입력하였습니다.
col은 색깔입니다. 내가 원하는 색깔을 순서대로 저렇게 나오게 하겠다. 라는 뜻입니다.
label은 각 칸마다 이름을 붙이겠다는 뜻이구요
init.angle은 시작하는 각도입니다. 12시에서 시작했으니 90도로 잡았습니다.
clockwise는 시계방향이구요
main은 그래프의 이름입니다.
짠 쉽죠?
##파이 차트의 변형 - 반원차트입니다. 코드는 다음과 같습니다
half.pie <- function(x, title, color, names) {
I <- matrix(c(-1,0,0,1,1,0),3,2,byrow=T)
plot(I, xlab="", ylab="", axes=F, type="n", asp=1)
x.circle <- cos(seq(0,100,1)*pi/100)
y.circle <- sin(seq(0,100,1)*pi/100)
half.circle <- cbind(x.circle, y.circle)
k <- length(x)
cum.x <- cumsum(x)+1; cum.x.1 <- c(1,cum.x[1:(k-1)])
for (j in 1:k){
half.circle2<-rbind(half.circle[cum.x[j]:cum.x.1[j],],0.6*half.circle[cum.x.1[j]:cum.x[j],])
polygon(half.circle2,col=color[j],border=NA)
center <- round((cum.x[j]+cum.x.1[j])/2)
text(0.8*half.circle[center,1], 0.8*half.circle[center,2], names[j])
}
title(title)
}
이를 이용해 예제를 풀어보면
prop <- c(47,23,18,11,1)
names <- c("None","Buddhism","Protestant","Catholic","Others")
예제는 다음과 같은 데이터를 가져왔습니다. 사람들이 종교를 선호하는 선호도입니다.
## 또다른 변형입니다.
dots <- function(x, title, color, names){
theta <- seq(0, 2*pi, length=360)
circle <- cbind(cos(theta), sin(theta))
plot(circle, xlab="", ylab="", axes=F, type="n", asp=1)
k <- length(x)
cum.x <- cumsum(x); cum.x.1 <- c(0,cum.x[1:(k-1)])
for (j in 1:5){
points(cos(0.5*pi+((cum.x[j]-1):cum.x.1[j])/100*2*pi),sin(0.5*pi+((cum.x[j]-1):cum.x.1[j])/100*2*pi),col=color[j],pch=20)
center <- round((cum.x[j]+cum.x.1[j])/2)
text(0.8*cos(0.5*pi+center/100*2*pi), 0.8*sin(0.5*pi+center/100*2*pi),
names[j], col=color[j])
}
title(title, line=1)
}
이런 식으로 얼마든지 변형해서 사용하실 수 있습니다.
오늘은 여기까지구요 내일은 회귀식에 대해서 들고오겠습니다.
많은 분들이 과분한 관심을 주셨는데 제가 아는게 얼마 없어서 죄송합니다 ㅜㅜ
일단 계획한바로는 내일 회귀식, 내일 모레 비장의 무언가! 로 끝낼 예정입니다...
제가 더 공부하고 나서 다시 올릴게요 많이 기대하셨다면 죄송합니다 ㅜㅠㅠㅠ
