Expanding window functions ด้วย pandas

การจัดการข้อมูล Time Series ใน Python

Stefan Jansen

Founder & Lead Data Scientist at Applied Artificial Intelligence

Expanding windows ใน pandas

จาก rolling window สู่ expanding window
คำนวณค่าสถิติตั้งแต่ต้นจนถึงวันปัจจุบัน
ชุดข้อมูลเวลาใหม่สะท้อนค่าประวัติศาสตร์ทั้งหมด
เหมาะสำหรับผลตอบแทนสะสม และ min/max สะสม
สองวิธีใน pandas:
- .expanding() — ใช้เหมือน .rolling()
- .cumsum(), .cumprod(), cummin()/max()

แนวคิดพื้นฐาน

df = pd.DataFrame({'data': range(5)})

df['expanding sum'] = df.data.expanding().sum()

df['cumulative sum'] = df.data.cumsum()

df

   data  expanding sum  cumulative sum
0     0            0.0               0
1     1            1.0               1
2     2            3.0               3
3     3            6.0               6
4     4           10.0              10

โหลดข้อมูล S&P 500

data = pd.read_csv('sp500.csv', parse_dates=['date'], index_col='date')

DatetimeIndex: 2519 entries, 2007-05-24 to 2017-05-24
Data columns (total 1 columns):
SP500    2519 non-null float64

ch3_2_v2 - Expanding Window Functions with Pandas.013.png

วิธีคำนวณผลตอบแทนสะสม

ผลตอบแทนรายงวด $r_t$: ราคาปัจจุบันหารด้วยราคางวดก่อนลบ 1:

$$r_t = \frac{P_t}{P_{t-1}} - 1$$
- ผลตอบแทนหลายงวด: ผลคูณของ $(1 + r_t)$ ทุกงวด ลบ 1:
$$R_T = (1 + r_1)(1 + r_2)...(1 + r_T) - 1$$
- ผลตอบแทนรายงวด: .pct_change()
- คณิตศาสตร์พื้นฐาน .add(), .sub(), .mul(), .div()
- ผลคูณสะสม: .cumprod()

ผลตอบแทนสะสมในทางปฏิบัติ

pr = data.SP500.pct_change() # period return

pr_plus_one = pr.add(1)

cumulative_return = pr_plus_one.cumprod().sub(1)

cumulative_return.mul(100).plot()

ch3_2_v2 - Expanding Window Functions with Pandas.021.png

การหา min & max สะสม

data['running_min'] = data.SP500.expanding().min()

data['running_max'] = data.SP500.expanding().max()

data.plot()

ch3_2_v2 - Expanding Window Functions with Pandas.023.png

ผลตอบแทนรายปีแบบ rolling

def multi_period_return(period_returns):
    return np.prod(period_returns + 1) - 1

pr = data.SP500.pct_change() # period return

r = pr.rolling('360D').apply(multi_period_return)

data['Rolling 1yr Return'] = r.mul(100)

data.plot(subplots=True)

ผลตอบแทนรายปีแบบ rolling

data['Rolling 1yr Return'] = r.mul(100)

data.plot(subplots=True)

ch3_2_v2 - Expanding Window Functions with Pandas.027.png

มาฝึกกันเถอะ!

การจัดการข้อมูล Time Series ใน Python