from numpy import mean,median,ptp,var,stdimport numpy as npimport pandas as pd#算术平均值data = [1,2,3,4,5,5,5,6,7,8,8,9]data_mean = mean(data)print("算术平均数 = ")print(data_mean)#中位数data_median = median(data)print("中位数 = ")print (data_median)#众数data_mode = mode(data)print("众数 = ")print (data_mode)#极差data_range = max(data)-min(data)print("极差 = ")print (data_range)#方差和标准差data_var = var (data)print ("方差 = ")print (data_var)#标准差data_sd = std(data)print ("标准差 = " )print (data_sd)#变异系数data_cv = mean(data) / std(data)print ("变异系数 = " )print (data_cv)#四分卫极差dt = pd.Series(np.array([1,2,3,4,5,5,5,6,7,8,8,9]))print('Q1:', dt.quantile(.25))print('Q3:', dt.quantile(.75))print('Q3 - Q1 :', dt.quantile(.75) - dt.quantile(.25) )

#整理成描述性统计的函数my_describe <- function(x){  options(digits = 3)  N = length(x);  Min = min(x, na.rm = TRUE);  Q1 = quantile(x, probs = 0.25, na.rm = TRUE);  Median = median(x, na.rm = TRUE);  Q3 = quantile(x, probs = 0.75, na.rm = TRUE);  Max = max(x, na.rm = TRUE);  Mean = mean(x, na.rm = TRUE);  Var = var(x, na.rm = TRUE);  Sd = sd(x, na.rm = TRUE);  Range = abs(diff(range(x)));  #返回结果  return(data.frame(N = N, Min = Min, Q1 = Q1, Median = Median, Q3 = Q3, Max = Max, Mean = Mean, Var = Var, Sd = Sd, Range = Range))}#可对data数据集批量使用sapply(data, my_describe)