예제 4.10: 합성곱

이제 슬슬 R이 손에 익는다. 모르는 문제 위주로만 풀어도 괜찮을 것 같다.

# 예제 4.10
# X1 ~ N(0,1) 이고
# X2 ~ N(3,1) 일 때,
# S = X1 + X2의 분포를 알아보고자 한다.
# 이때 X1과 X2는 독립이다.

# 수리적으로 S ~ N(3,2)를 따른다고 하는데,
# X1 + X2를 통해 구해보자.

n=10000 #10000번으로 고정하자 그냥
X.1 = rnorm(n,0,1)
X.2 = rnorm(n,3,1)
S = X.1 + X.2
N.3.2 = rnorm(n,3,2^0.5) #분산 말고 표준편차 써야 하는 것 유의.
qqplot(S,N.3.2,cex=0.25)
plot(density(S),xlim=c(-5,10),ylim=c(0.,0.4),lty=1,xlab="x",main="") #lty가 선 모양
lines(density(X.1),lty=2)
lines(density(X.2),lty=3)
lines(density(N.3.2), lty=4)
legend(4.5, 0.4, c("X1+X2","X1","X2","N(3,2)"), lty=1:4)
#legend: '범례'라는 뜻. 여태까지 쓰면서도 몰랐네..
#lty:선 종류 (line type?) #lwd:선 굵기
#lty 1:4 lty 1번부터 4번까지 쓴다는 뜻.

추가로 정규혼합 푸는 문제 있긴 한데, 얘는 나중에 풀자.

시간이 많이 없어서...