4주

From Biocourse

Jump to: navigation, search
간단한 그림 그리기..

>  x= -5:5                                                                                      ## x에 -5부터 5까지 한개씩 증가하면서 할당.
>  y= x^2                                                                                      ## y에 x를 제곱하여 할당.
>  par(mfrow=c(3, 2))                                                                  ## 그림을 파티션함. 3*2로
>  plot(x, y)                                                                                  ## plot을 이용하여 그림을 그린다. 기본값이면 dot형태의 그림
                                                                                                      ## 아래 그림에서 1,1위치
>  plot(x, y, type="l")                                                                  ## plot에서 type="l", line으로 그린다. 아래 그림에서 1,2 위치
>  plot(x, y, type="b")                                                                 ## plot에서 type="b", line과 dot를 같이 그린다. 2,1위치
>  plot(x, y, type="h")                                                                 ## plot에서 type="h", 수직선으로 그린다. 2,2 위치
>  plot(x, y, type="o")                                                                 ## plot에서 type="o", line과 dot를 같이 그린다. 3,1위치
>  plot(x, y, type="s")                                                                 ## plot에서 type="s", line으로 계단형식으로 그린다. 3,2위치

>  mtext("Different options for the plot parameter type", side=3, outer=T, line=0.5)
## 그림안에 text를 삽입시 사용된다. side는 그림에서의 위치, outer는 sub 그림여부, line은 전체 그림에서의 위치를 나타낸다.

image:r-4.jpeg

위의 코드를 이용하여 다음과 같은 그림을 그릴수 있다.
사용자가 적당한 option을 바꿔 가면서 사용을 하면서, 익숙해 지길 바란다.



>  x= seq( -5, 5, 1)                                                                      ## -5부터 5까지 1개씩 증가하여 x에 할당한다.
>  y= x^2                                                                                    ## x를 제곱하여 y에 할당한다.
>  plot(x, y, pch="X", main="Main Title", sub="Subtitle",  xlab="X Axis Label", ylab="Y Axis Label", xlim=c(-8, 8)
+  , type="o", lty=2)
 

## x와 y를 이용하여 그림을 그린다. main, sub, xlab, ylab, xlim,type,lty와 같은 옵셥을 사용한다.


image:r-4-1.jpeg


pch : 그림에서 point을 나타낸다.
main : Main Title 그림 상단에 글씨를 보인다.
sub :  Subtitle     그림 하단에 글씨를 보인다.
xlab : X Axis Label 그림의 x축에 글씨를 보인다.
ylab : Y Axis Label 그림의 y축에 글씨를 보인다.
xlim : x 축의 정도를 정한다. (-8부터 8까지)
type : 그림의 형태를 정한다.
lty : 선의 type을 정한다.



>  x= seq(0, 2*pi, length=21)                                                                                      ## x를 다음과 같이 할당한다.
>  y= sin(x)                                                                                                                   ## y를 sin 함수를 이용하여 할당한다.
>  plot(x, y, axes=F, type="b",pch="x",xlab="",ylab="")                                      ## x,y를 그린다. x,y축은 그리지 않는다.
>  axis(1, at=c(0, 1, 2, pi, 4, 5, 2*pi),   labels=c(0, 1, 2, "Pi", 4, 5, "2 Pi"), pos=0)  ## x축을 그린다.
>  axis(2, at=c(-1, -0.5, 0, 0.25, 0.5, 0.75, 1))                                                              ## y축을 그린다.
>  abline(h=c(-1, -0.5, 0.5, 1), lty=3)                                                                           ## 추가로 선을 그린다.
>  text(pi, 0.1, "sin(pi)=0", adj=0)                                                                              ## 글씨를 추가한다.
>  title("The sine function\nfrom 0 to 2 Pi")                                                             ## 그림에 title을 추가한다.

image:r-4-2.jpeg

## 사용자가 다양한 형태의 그림을 그릴수 있다.
## 자기가 원하는 x,y축을 표현할 수 있고, 추가로 선을 넣을수 있으며, 글씨를 추가할 수 있다.




## 같은 자료로 다양한 그림 그리기


> par(mfrow=c(2,2))                                                                                ## 2 * 2 로 part를 나눈다.
> x<-1:100/100:1                                                                                     ## x를 할당한다.
> plot(x)                                                                                                   ## plot을 그린다.
> plot(x,type="l")                                                                                   ## plot에서 option을 이용해서 그린다.
> hist(x)                                                                                                   ## 히스토 그램을 그린다.
> boxplot(x)                                                                                            ## Boxplot을 그린다.

image:r-4-3.jpeg




  R-code 3




실습


1. Read => homework-2week.txt                        ==> 3week : solution1


1.1  Boxplot ( Sample 별 )                                      ==> 4week : solution1

 
1.2 Gene 별로 Missing이 없는 자료 선택. 

     => Boxplot ( by samples ) , Boxplot( by genes )  ==> 4week : solution2


     => 평균 Histogram                                              ==> 4week : solution3


     => 분산 Histogram                                              ==> 4week : solution4


     => 평균 vs. 분산 산점도                                      ==> 4week : solution5


     => max 평균 gene vs. min 평균 gene 산점도       ==> 4week : solution6


1.3 Sample별missing=0 삽입.                                 ==> 4week : solution7


     => 평균 Histogram                                              ==> 4week : solution8


     => 분산 Histogram                                              ==> 4week : solution9