예제 4.8 변환법 - 스튜던트 t 분포

표준정규분포와 카이제곱 분포를 활용하여 student t 분포를 따르는 난수를 생성할 수 있다. 4.7과 함께 변환법의 다른 예시이다.

# Z~ N(0,1) 과 V~Chisq(m)이 독립이면
# T=Z/(V/m)**(1/2) 는 자유도 m의 student t 분포를 따름.

n <- 1000
z <- rnorm(n,0,1) #분포 앞에 r 붙이면 난수 추출한다는 뜻
v <- rchisq(n,m)
x <- z/(v/m)**(1/2) # R은 1/2제곱 연산 되네...?

q <- qt(ppoints(n),m)
qqplot(q,x,cex=0.25,xlab="t(2)",ylab="Sample")
qqline(q,col=2)
#qqline: 정규분포 qq plot에서 1Q와 3Q 지나는 선 그리는 함수
#col : 색깔. 1: 검정, 2:빨강 ...

++ qqplot함수 더 알아봐야겠다.

x <- c(1,2,3) #벡터 만드려면 c()함수를 써야 합니다.
y <- c(2,1,3) #x에세 순서만 바꿨다.
qqplot(x,y) #0, 0.5, 1 점에 각각 점 1개씩 찍힌다.

x <- c(1,2,4)
y <- c(1,4,2)
qqplot(x,y) #자동으로 오름차순 정렬된다.