#导入数据
import pandas as pd
data = pd.read_excel(“第参数估计.xlsx”)
data[‘Age’].mean() #总体均值:25.46902269662129
#抽取100个样本
data_sample = data[‘Age’].sample(n=100)
样本均值的抽样分布为正态分布
import numpy as np
from scipy.stats import norm
import math
sample_mean = np.mean(data_sample)
sample_std = np.std(data_sample,ddof=1)
n = len(data_sample)
z = norm.interval(alpha=0.95, loc=0, scale=1)
print(“总体均值在0.95置信水平下的置信区间为:(%.4f,%.4f)” %(conf_intveral_lower,conf_intveral_upper))