본문 바로가기

공부 : 금융공학, 금융, 통계, 공학, 경제 등

(39)
예제 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"..
예제 4.9 합과 혼합: 컨볼루션 convolution, chisq의 사례로. 정규분포의 합성곱(convolution)을 통해 카이제곱분포를 만들어내는 문제다. Zv~N(0,1) 의 제곱의 합 sum(Zv^2 ... Zv^2)은 Chisq(v)을 따른다. convolution 개념이 여기서도 나오는구나...? 그럼 코드 보자. 예제에서는 v=2, n(표본 개수)=10000을 다룬다. 1000은 적으니 10000으로 합시다. # 내가 직접 짠 방법 n
예제 4.8 변환법 - 스튜던트 t 분포 표준정규분포와 카이제곱 분포를 활용하여 student t 분포를 따르는 난수를 생성할 수 있다. 4.7과 함께 변환법의 다른 예시이다. # Z~ N(0,1) 과 V~Chisq(m)이 독립이면 # T=Z/(V/m)**(1/2) 는 자유도 m의 student t 분포를 따름. n
예제 4.7 [qqplot() 이해가 잘 안 감.] 베타 분포를 생성. 감마분포 2개를 만들어서 베타분포를 만들어냄. 정의에 따라 베타분포가 두 개의 감마분포에 의해 형성됨. cf. $$ U \sim Gamma(r,\lambda), V\sim Gamma(s,\lambda) \\ \rightarrow X= U/(U+V) \sim Beta(r,s) ; $$ 라고 합니다. (단, U, V는 독립) 이걸 바탕으로 코드를 짜 보면 #Beta(3,2)를 따르는 난수 추출 n
예제 4.5 [베타 분포 공부 필요] 베타 분포에서의 난수 발생을 채택-기각법을 사용하여 수행. 채택-기각법으로 Beta(2,2) 분포를 따르는 1,000개의 변량을 생성하려면? 코드 확인해보자. n
예제 3.8 장난감 5개 중 한 명의 아이가 가지고 노는 장난감의 수 X의 분포는 아래와 같음. 0 0.03 1 0.16 2 0.30 3 0.23 4 0.17 5 0.11 이때 평균은 2.68이고 표준편차는 1.3106이다. 이를 모의실험을 통해 추정해보자. 모의실험은 총 4단계로 나뉜다. 1단계: 가설 설정 2단계: 각 균일난수에 대해 숫자 할당. 3단계: 임의의 균일난수 추출. 이를 할당된 범위에 대입함. 4단계: 3단계를 충분히 반복함. 이를 통해 모집한 표본에 대해 평균과 표준편차를 구한다. 이를 R 코드로 구현하면 아래와 같다. n
예제 3.7 [코드만 있음. 글자 살짝 깨짐.] 1500명 성인, 비율 p=0.15. p를 추정하기 위해 표본비율 $\hat{p}$을 이용하고자 함. $\hat{p}$ 은 근사적으로 평균 $\mu = 0.15$이고 $\sigma=0.0092$인 정규분포를 따르는 확률변수이다. Solution. m
R을 이용한 통계계산 문제풀이 여기다가 링크 남긴다. 시험 끝나고 푼 만큼 코드 적을게요. 안뇽.
Hierarchical Risk Parity (HRP) 8. Linkage Matrix Linkage Matrix는 Hierarchical Clustering에서 사용되는 특별한 행렬의 일종이다. Matrix의 각 행이 하나의 link를 의미하며, 각 열은 link에 대한 정보를 담고 있다. array([[ 0. , 1. , 1. , 2. ], [ 3. , 4. , 1. , 2. ], [ 6. , 7. , 1. , 2. ], [ 9. , 10. , 1. , 2. ], [ 2. , 12. , 1.29099445, 3. ], [ 5. , 13. , 1.29099445, 3. ], [ 8. , 14. , 1.29099445, 3. ], [11. , 15. , 1.29099445, 3. ], [16. , 17. , 5.77350269, 6. ], [18. , 19. , 5.77350269, 6. ..
Hierarchical Risk Parity (HRP) 7. Euclidean Distance Matrix 오랜만에 수식을 봐서 머리가 지끈할 수도 있으니, n차원 좌표계에서의 직선 거리(=Euclidean Distance)를 피타고라스 정리에서부터 유도해보자. 우리가 아는 피타고라스 정리는 다음과 같다. 직각삼각형의 빗변 c와 나머지 두 변 a, b에 대해, c**2 = a**2 + b**2 … (1), 또한 c = (a**2 + b**2)**(1/2) …(2) 도 자연스레 유도할 수 있을 것이다. ( ‘**’은 제곱 기호이다.) 또한 이를 응용한 좌표 상의 거리 공식인 두 점 A(x1,y1), B(x2,y2)에 대해 Distance(A,B) = ((x1-x2)^2 + (y1-y2)^2)^(1/2)…(3) 또한 유도할 수 있을 것이다. (x축, y축 그리고, 두 점 찍고, 두 점 사이에 직선을 긋고, 나머지..