민감도 분석

Python으로 배우는 Monte Carlo 시뮬레이션

Izzy Weber

Curriculum Manager, DataCamp

민감도 분석

입력 범위의 영향 파악에 도움
표나 플롯 요약 시 패턴·추세를 보여줌

Monte Carlo 시뮬레이션으로 bmi와 hdl 값을 늘리거나 줄이면, 예측 y 값(질병 진행)은 어떻게 변할까요?

매개변수 정의

cov_dia = dia[["age", "bmi", "bp", "tc", "ldl", "hdl", "tch", "ltg", "glu"]].cov()
mean_dia = dia[["age", "bmi", "bp", "tc", "ldl", "hdl", "tch", "ltg", "glu"]].mean()

시뮬레이션 함수 정의

def simulate_bmi_hdl(cov_dia, mean_list):

    list_ys = []
    for i in range(50):
        simulation_results = st.multivariate_normal.rvs(mean=mean_list, 
                                                        size=500, cov=cov_dia)
        df_results = pd.DataFrame(simulation_results, 
                                  columns=["age","bmi","bp","tc","ldl","hdl","tch","ltg","glu"])
        predicted_y = regr_model.predict(df_results)
        df_y = pd.DataFrame(predicted_y, columns=["predicted_y"])
        df_summary = pd.concat([df_results, df_y], axis=1)
        y = np.mean(df_summary["predicted_y"]) 
        list_ys.append(y)

    return(np.mean(list_ys))

입력 매개변수 범위로 시뮬레이션 수행

hdl = []
bmi = []
simu_y = []
for mean_hdl_inc in np.arange(-20, 50, 30): 
    for mean_bmi_inc in np.arange(-7, 11, 3):

        mean_list = mean_dia + np.array([0, mean_bmi_inc, 0, 0, 0, mean_hdl_inc, 0, 0, 0])
        hdl.append(mean_hdl_inc)
        bmi.append(mean_bmi_inc)

        mean_y = simulate_bmi_hdl(cov_dia, mean_list)

        simu_y.append(mean_y)

df_sa = pd.concat([pd.Series(hdl), pd.Series(bmi), pd.Series(simu_y)], axis=1)

df_sa.columns = ["hdl_inc", "bmi_inc", "y"]

민감도 분석 결과의 스타일된 DataFrame

df_sa.sort_values(by=['hdl_inc', 'bmi_inc']).pivot(index='hdl_inc',
                                             columns='bmi_inc',
                                             values='y').style.background_gradient(
                                             cmap=sns.light_palette("red", as_cmap=True))

정렬, 피벗, 스타일 적용 후의 df_sa

민감도 분석 결과용 헥스빈 플롯

df_sa.plot.hexbin(x='hdl_inc',y='bmi_inc', C='y',
                    reduce_C_function=np.mean,
                    gridsize=10, cmap="viridis",
                    sharex=False)

df_sa의 헥스빈

조밀한 매개변수 공간의 헥스빈 플롯

연습해 봅시다!

Python으로 배우는 Monte Carlo 시뮬레이션