R-code 4
From Biocourse
## T분포에서 기각역 그림. ( 자료 : barley.txt )
> barley = read.table("c://barley.txt",sep="\t",header=T)
## 자료를 C 디렉토리에서 읽어서 barley로 저정한다.
> attach(barley) ## barley 자료를 활성화 시킨다.
> x<-yield[site=="Crookston"] ## site 가 Crookston인 자료의 yield값을 x에 저장
> y<-yield[site=="Morris"] ## site 가 Morris인 자료의 yield값을 y에 저장
> df<-length(x)-1 ## x의 길이 -1 을 df에 저장한다.
> alpha<-0.05 ## 0.05를 alpha에 저장한다.
> bound.left<- -3 ## -3을 bound.left에 저장한다.
> bound.right<- 3 ## 3을 bound.right에 저장한다.
> xaxis<-seq(bound.left,bound.left,length=1000)
## bound.left에서 bound.left까지를 1000등분한다.
> yaxis<-dt(xaxis,df) ## xaxis의 값을 t 분포 자유도가 df인 density를 구한다.
> plot(xaxis,yaxis,type="l",xlab="",ylab="t density") ## xaxis와 yaxis로 그림을 그린다.
> critical.left<-qt(alpha/2,df)
## alpha/2를 t분포에서 자유도가 df인 qantile 값을 구함
> critical.right<-qt(1-alpha/2,df)
## 1 - alpha/2를 t분포에서 자유도가 df인 qantile 값을 구한다.
> xaxis<-seq(critical.left,critical.right,length=100)
## critical.left에서 critical.right까지를 100등분한다.
> yaxis<-c(dt(xaxis,df),0,0) ## yaxis의 값을 정의한다.
> xaxis<-c(xaxis,critical.left,bound.left) ## xaxis의 값을 정의한다.
> polygon(xaxis,yaxis,density=25)
## polygon을 이용하여 아래 그림 처럼 구간을 색으로 채운다.
> xaxis<-seq(critical.right,bound.right,length=100)
## critical.right에서 critical.right까지를 100등분한다.
> yaxis<-c(dt(xaxis,df),0,0) ## yaxis의 값을 정의한다.
> xaxis<-c(xaxis,bound.right,critical.right) ## xaxis의 값을 정의한다.
> polygon(xaxis,yaxis,density=25)
##polygon을 이용하여 아래 그림 처럼 구간을 색으로 채운다.
> test.stat<-t.test(x,y,paired=T,conf.level=1-alpha)$statistic
##.test 쌍 자료 분석 결과에서 신뢰구간 1-alpha값을 저장한다.
> points(test.stat,0.01,cex=1.5,pch="*") ## test.stat를 점으로 표현한다.
> text(test.stat,0.04,"paired test statistics",adj=1)
##그림에 paired test statistics를 삽입한다
> test.stat<-t.test(x,y,paired=F,conf.level=1-alpha)$statistic
##.test 두 자료 분석 결과에서 신뢰구간 1-alpha값을 저장한다.
> points(test.stat,0.01,cex=1.5,pch="*") ## test.stat를 점으로 표현한다.
> text(test.stat,0.04,"unpaired",adj=0) ## 그림에 unpaired를 삽입한다.
> text(bound.left, 0.08,"rejection \n region",adj=0) ## 그림에 rejection region를 삽입한다.
> text(bound.right,0.08,"rejection \n region",adj=1) ## 그림에 rejection region를 삽입한다.
> text((bound.left+bound.right)/2,0.16,"acceptance region")
##그림에 acceptance region를 삽입한다.
> xaxis<-c(rep(critical.left,2),rep(critical.right,2))
> yaxis<-c(0.12,0.14,0.14,0.12)
> lines(xaxis,yaxis)
