670339.py

import numpy as np
import pandas as pd

install = np.concatenate([np.ones(19), np.zeros(1)])
repair = np.concatenate([np.ones(19), [np.nan]])
outage = np.concatenate([np.ones(17), np.zeros(1), [np.nan, np.nan]])
df = pd.DataFrame({'install': install, 'repair': repair, 'outage': outage})
print(df)

proportions = []

for i in range(1_000_000):
    indexes = np.random.choice(a=len(df), size=len(df), replace=True)
    p = np.nanmean(df.iloc[indexes].to_numpy())
    proportions.append(p)

proportions = np.array(proportions)
mean = np.mean(proportions)
std = np.std(proportions)
print(f'{mean=!s}, {std=!s}')