# Packages
import pandas as pd
import numpy as np
from statsmodels.stats.weightstats import DescrStatsW
# Load Data
year_cur=2022
FilePath = "Data/scfp" + np.str_(year_cur) + "s/rscfp" + np.str_(year_cur) + ".dta"
df = pd.read_stata(FilePath)
# Median
wq = DescrStatsW(data=df.loc[:,'networth'], weights=df.loc[:,'wgt'])
NW_Median = wq.quantile(probs=0.50, return_pandas=False)
print("Median = " + np.str_(NW_Median[0]))
Median = 192700.0
# Mean
NW_Mean = np.average(a=df.loc[:,'networth'], weights=df.loc[:,'wgt'])
print("Mean = " + np.str_(NW_Mean.round(0)))
Mean = 1059470.0