기본 콘텐츠로 건너뛰기

[data analysis]로그-노말 분포(Log-normal distribution)

[data analysis] 피벗테이블(Pivot table)

피벗테이블(Pivot table)

기준되는 항목을 축(행 또는 열)으로 지정하여 표로 작성하는 것을 피벗테이블(pivot table)이라하며 함수 pandas.pivot_table(data, values=None, index=None, columns=None, aggfunc='mean', …)를 적용합니다.

패키지 FinanceDataReader의 함수에 의해 호출한 주가 자료는 열이름은 Open, High, Low, Close, Volume, Change를 포함합니다. 각각 시가, 고가, 저가, 종가, 거래량을 의미하며 Change는 전날 종가에 비해 현재 시가의 비율입니다. 호출한 자료에서 Change를 목록화한후 그 목록에 따른 다른 값들의 평균을 나타내 봅니다. 이 표를 위해서 Change가 행 또는 열에 지정되어야 하고 나머지 중 타겟이 되는 부분을 values로 지정하여야 합니다.

import numpy as np 
import pandas as pd 
import FinanceDataReader as fdr
st=pd.Timestamp(2023,8, 26)
et=pd.Timestamp(2024, 8,8)
nme=['Open', 'High','Low','Close', 'Change', 'Volume']
da=fdr.DataReader('KS11', st, et)[nme]
da.head(2).round(3)
Open High Low Close Change Volume
Date
2023-08-28 2534.42 2543.41 2525.64 2543.41 0.010 483587568
2023-08-29 2550.76 2556.98 2545.27 2552.16 0.003 441365852

위 자료에서 그룹화하기 위한 "Change"를 목록변수로 변환합니다. pd.qcut()함수를 적용합니다.

chg, border=pd.qcut(da["Change"],[0, 0.25, 0.5, 0.75, 1], labels=[1,2,3,4], retbins=True)
chg.head(3)
Date
2023-08-28    4
2023-08-29    3
2023-08-30    3
Name: Change, dtype: category
Categories (4, int64): [1 < 2 < 3 < 4]
print(border)#각 그룹의 경계값
[-0.0877 -0.0073  0.0004  0.0077  0.0566]

위 객체들로 부터 다음 객체를 생성

da1=da.drop(labels=["Change"], axis=1)
da1=da1.pct_change()
da1['chg']=chg[0]
da1=da1.dropna()
da1.head(3).round(3)
Open High Low Close Volume chg
Date
2023-08-29 0.006 0.005 0.008 0.003 -0.087 3
2023-08-30 0.010 0.008 0.006 0.004 -0.175 3
2023-08-31 -0.004 -0.004 -0.005 -0.002 0.300 2

da1에서 chg를 기준 축으로 나머지 변수를 분류하고 각 그룹의 평균을 계산합니다.

allChg=pd.pivot_table(da1, columns="chg", aggfunc="mean")
allChg.round(3)
chg 1 2 3 4
Close -0.014 -0.003 0.003 0.014
High -0.008 -0.002 0.001 0.010
Low -0.010 -0.002 0.003 0.010
Open -0.004 -0.001 0.000 0.005
Volume 0.027 0.020 -0.006 0.070

위와 유사하지만 대상 values를 특정합니다.

closeChg=pd.pivot_table(da1, columns="chg", values=["Close", "Volume"], aggfunc="mean")
closeChg.round(3)
chg 1 2 3 4
Close -0.014 -0.003 0.003 0.014
Volume 0.027 0.020 -0.006 0.070

위 데이터 da1은 그룹화할 수 있는 변수가 1개입니다. 목록변수가 두개일 경우 각 변수에 따른 분류를 작성합니다. 즉 교차표를 작성합니다. 위 자료 중 Open을 목록화 합니다.

op, opBorder=pd.qcut(da1["Open"], 4, labels=range(1, 5), retbins=True)
op.head(3)
Date
2023-08-29    4
2023-08-30    4
2023-08-31    2
Name: Open, dtype: category
Categories (4, int64): [1 < 2 < 3 < 4]
print(opBorder.round(3))
[-0.04  -0.005  0.     0.006  0.032]
da2=da1
da2["Open"]=op
da2.head(3).round(3)
Open High Low Close Volume chg
Date
2023-08-29 4 0.005 0.008 0.003 -0.087 3
2023-08-30 4 0.008 0.006 0.004 -0.175 3
2023-08-31 2 -0.004 -0.005 -0.002 0.300 2

자료 da2에서 Open과 chg의 그룹에 대응하는 Close, Volume의 평균값에 대한 표를 작성합니다.

pd.pivot_table(da2, index="Open", values=["Close", "Volume"], columns="chg" )
Close Volume
chg 1 2 3 4 1 2 3 4
Open
1 -0.017533 -0.003300 0.002405 0.016712 0.045754 -0.005855 -0.083844 -0.131276
2 -0.010962 -0.002017 0.003407 0.014261 0.016789 0.004625 0.073962 0.043727
3 -0.011745 -0.002706 0.003402 0.012165 0.030558 0.114356 -0.052739 0.103737
4 -0.011536 -0.002113 0.003700 0.015344 -0.022143 -0.012375 0.005476 0.096930

자료 da2에서 Open과 chg의 각 그룹에 대응하는 Close의 빈도수에 대한 표를 작성합니다.

pd.pivot_table(da2, index="Open", values="Close", columns="chg" , aggfunc="count", margins=True)
chg 1 2 3 4 All
Open
1 26 17 9 6 58
2 16 19 18 5 58
3 10 12 18 18 58
4 7 12 11 28 58
All 59 60 56 57 232

댓글

이 블로그의 인기 게시물

유사변환과 대각화

내용 유사변환 유사행렬의 특성 대각화(Diagonalization) 유사변환(Similarity transformation) 유사변환 n×n 차원의 정방 행렬 A, B 그리고 가역 행렬 P 사이에 식 1의 관계가 성립하면 행렬 A와 B는 유사하다고 하며 이 변환을 유사 변환 (similarity transformation)이라고 합니다. $$\begin{equation}\tag{1} A = PBP^{-1} \Leftrightarrow P^{-1}AP = B \end{equation}$$ 식 1의 유사 변환은 다음과 같이 고유값을 적용하여 특성 방정식 형태로 정리할 수 있습니다. $$\begin{align} B - \lambda I &= P^{-1}AP – \lambda P^{-1}P\\ &= P^{-1}(AP – \lambda P)\\ &= P^{-1}(A - \lambda I)P \end{align}$$ 위 식의 행렬식은 다음과 같이 정리됩니다. $$\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)\\ &= \t

[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