기본 콘텐츠로 건너뛰기

[ML] 결정트리(Decision Tree) 모델

한국과 외국의 시계열 자료의 합병

한국 금융자료와 외국 자료의 합병

시계열 자료는 날짜등의 시간에 따라 생산되는 자료입니다. 대표적으로 일일 주가자료는 거래일 기준으로 자료가 생산됩니다. 별도의 자료를 합치는 경우 시간이 이질적이라면 합쳐진 자료에 결측치가 포함됩니다. 이 경우 결측치를 조절하는 방법을 알아봅니다.

한국의 코스피와 환율, 외국의 지수 자료를 합쳐봅니다.

import numpy as np
import pandas as pd
import FinanceDataReader as fdr
지수명코스피환율다우나스탁필라델피아반도체변동성(미국)
코드명KS11USD/KRWDJIIXICSOXXVIX

코스피 지수의 모든 변수(Open, High, Low, Close, Volume, Change)와 다른 지수들의 종가(Close)를 호출하였습니다. 이 호출은 FinanceDataReader 패키지의 DataReader()함수를 적용하였습니다.

st=pd.Timestamp(2010,3, 1)
et=pd.Timestamp(2022, 6, 4)
nme={'exchg':'USD/KRW','dj':"DJI",'nasd':'IXIC','soxx':"SOXX", 'vix':'VIX'}
kos=fdr.DataReader('KS11', st, et)
stock={}
for i, j in zip(nme.keys(), nme.values()):
    stock[i]=fdr.DataReader(j, st, et)['Close']
stock=pd.DataFrame(stock.values(), index=stock.keys()).T
kos.tail(3)
Close Open High Low Volume Change
Date
2022-05-31 2685.90 2666.88 2685.90 2654.32 671680000.0 0.0061
2022-06-02 2658.99 2670.74 2674.00 2653.59 525260000.0 -0.0100
2022-06-03 2670.65 2679.57 2681.51 2663.00 562470000.0 0.0044
stock.tail(3)
exchg dj nasd soxx vix
Date
2022-06-01 1246.89 32813.10 11994.46 419.38 25.69
2022-06-02 1240.32 33248.02 12316.90 433.99 24.72
2022-06-03 1250.58 32898.91 12012.73 421.01 24.79

DataReader()에 의해 호출된 자료는 결측치가 0으로 작성됩니다. numpy.where()함수에 의해 감지할 수 있습니다.

np.where(kos==0), np.where(stock==0)
((array([ 110,  574,  597,  665,  695,  767,  813,  955, 1028, 1184, 1242,
         1397, 2129, 2205, 2249, 2274, 2500]),
  array([5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5])),
 (array([], dtype=int64), array([], dtype=int64)))

결측치는 DataFrame.mask()와 DataFrame.ffill()함수를 사용하여 처리할 수 있습니다.

kos1=kos.mask(kos==0, np.nan).ffill()
stock1=stock.mask(stock==0, np.nan).ffill()
np.where(kos1==0), np.where(stock1==0)
((array([], dtype=int64), array([], dtype=int64)),
 (array([], dtype=int64), array([], dtype=int64)))
kos1.shape, stock1.shape
((3027, 6), (3201, 5))

pandas.concat()함수를 사용하여 위의 두 객체를 결합합니다. 이 함수 중 매개변수 join은 결합방식을 지정하는 것으로 'outer' 두 객체의 모든 부분을 통합하고 'inner'은 공통부분만을 결합합니다.

Stock=pd.concat([kos1, stock1], join='outer', axis=1).ffill().dropna()
Stock.head(3)
Close Open High Low Volume Change exchg dj nasd soxx vix
Date
2010-03-02 1615.12 1612.43 1620.15 1608.97 415100000.0 0.0129 1145.00 10405.98 2280.79 47.51 19.06
2010-03-03 1622.44 1615.25 1622.63 1611.35 403470000.0 0.0045 1142.25 10396.76 2280.68 47.28 18.83
2010-03-04 1618.20 1625.35 1631.10 1612.47 383040000.0 -0.0026 1144.00 10444.14 2292.31 47.24 18.72
np.where(Stock==0), np.where(Stock==np.nan)
((array([], dtype=int64), array([], dtype=int64)),
 (array([], dtype=int64), array([], dtype=int64)))

댓글

이 블로그의 인기 게시물

[Linear Algebra] 유사변환(Similarity transformation)

유사변환(Similarity transformation) n×n 차원의 정방 행렬 A, B 그리고 가역 행렬 P 사이에 식 1의 관계가 성립하면 행렬 A와 B는 유사행렬(similarity matrix)이 되며 행렬 A를 가역행렬 P와 B로 분해하는 것을 유사 변환(similarity transformation) 이라고 합니다. $$\tag{1} A = PBP^{-1} \Leftrightarrow P^{-1}AP = B $$ 식 2는 식 1의 양변에 B의 고유값을 고려한 것입니다. \begin{align}\tag{식 2} B - \lambda I &= P^{-1}AP – \lambda P^{-1}P\\ &= P^{-1}(AP – \lambda P)\\ &= P^{-1}(A - \lambda I)P \end{align} 식 2의 행렬식은 식 3과 같이 정리됩니다. \begin{align} &\begin{aligned}\textsf{det}(B - \lambda I ) & = \textsf{det}(P^{-1}(AP – \lambda P))\\ &= \textsf{det}(P^{-1}) \textsf{det}((A – \lambda I)) \textsf{det}(P)\\ &= \textsf{det}(P^{-1}) \textsf{det}(P) \textsf{det}((A – \lambda I))\\ &= \textsf{det}(A – \lambda I)\end{aligned}\\ &\begin{aligned}\because \; \textsf{det}(P^{-1}) \textsf{det}(P) &= \textsf{det}(P^{-1}P)\\ &= \textsf{det}(I)\end{aligned}\end{align} 유사행렬의 특성 유사행렬인 두 정방행렬 A와 B는 'A ~ B' 와 같

[matplotlib] 히스토그램(Histogram)

히스토그램(Histogram) 히스토그램은 확률분포의 그래픽적인 표현이며 막대그래프의 종류입니다. 이 그래프가 확률분포와 관계가 있으므로 통계적 요소를 나타내기 위해 많이 사용됩니다. plt.hist(X, bins=10)함수를 사용합니다. x=np.random.randn(1000) plt.hist(x, 10) plt.show() 위 그래프의 y축은 각 구간에 해당하는 갯수이다. 빈도수 대신 확률밀도를 나타내기 위해서는 위 함수의 매개변수 normed=True로 조정하여 나타낼 수 있다. 또한 매개변수 bins의 인수를 숫자로 전달할 수 있지만 리스트 객체로 지정할 수 있다. 막대그래프의 경우와 마찬가지로 각 막대의 폭은 매개변수 width에 의해 조정된다. y=np.linspace(min(x)-1, max(x)+1, 10) y array([-4.48810153, -3.54351935, -2.59893717, -1.65435499, -0.70977282, 0.23480936, 1.17939154, 2.12397372, 3.0685559 , 4.01313807]) plt.hist(x, y, normed=True) plt.show()

R 미분과 적분

내용 expression 미분 2차 미분 mosaic를 사용한 미분 적분 미분과 적분 R에서의 미분과 적분 함수는 expression()함수에 의해 생성된 표현식을 대상으로 합니다. expression expression(문자, 또는 식) 이 표현식의 평가는 eval() 함수에 의해 실행됩니다. > ex1<-expression(1+0:9) > ex1 expression(1 + 0:9) > eval(ex1) [1] 1 2 3 4 5 6 7 8 9 10 > ex2<-expression(u, 2, u+0:9) > ex2 expression(u, 2, u + 0:9) > ex2[1] expression(u) > ex2[2] expression(2) > ex2[3] expression(u + 0:9) > u<-0.9 > eval(ex2[3]) [1] 0.9 1.9 2.9 3.9 4.9 5.9 6.9 7.9 8.9 9.9 미분 D(표현식, 미분 변수) 함수로 미분을 실행합니다. 이 함수의 표현식은 expression() 함수로 생성된 객체이며 미분 변수는 다음 식의 분모의 변수를 의미합니다. $$\frac{d}{d \text{변수}}\text{표현식}$$ 이 함수는 어떤 함수의 미분의 결과를 표현식으로 반환합니다. > D(expression(2*x^3), "x") 2 * (3 * x^2) > eq<-expression(log(x)) > eq expression(log(x)) > D(eq, "x") 1/x > eq2<-expression(a/(1+b*exp(-d*x))); eq2 expression(a/(1 + b * exp(-d * x))) > D(eq2, "x") a * (b * (exp(-d * x) * d))/(1 + b