# Ceteris Lab downloadable Python script
# Course: Fundamentals of Python for Financial Econometrics

import numpy as np

rng = np.random.default_rng(1401)
draws = rng.normal(0, 1, 100_000)
probability = np.mean(draws > 1.96)
print(round(probability, 4))

import numpy as np
from scipy import stats

sample = np.array([2.1, 2.4, 2.0, 2.7, 2.3, 2.5])
mean = sample.mean()
se = stats.sem(sample)
interval = stats.t.interval(0.95, df=len(sample)-1, loc=mean, scale=se)
print(round(mean, 3), tuple(round(v, 3) for v in interval))

import numpy as np
rng = np.random.default_rng(1401)
draws = rng.normal(0, 1, 100_000)
probability = np.mean(draws > 1.96)
print(round(probability, 4))
import numpy as np