数据探查

#… 省略部分代码
#计算统计量
minv=app.api.WorksheetFunction.Min(y)    #最小值
maxv=app.api.WorksheetFunction.Max(y)    #最大值
meanv=app.api.WorksheetFunction.Average(y)    #均值
median=app.api.WorksheetFunction.Median(y)    #中值
rng=maxv-minv    #极差
gmean=app.api.WorksheetFunction.GeoMean(y)    #几何均值
hmean=app.api.WorksheetFunction.HarMean(y)    #调和均值
tmean=app.api.WorksheetFunction.TrimMean(y,0.05)    #截尾均值
pt25=app.api.WorksheetFunction.Percentile(y,0.25)    #25%分位数
pt75=app.api.WorksheetFunction.Percentile(y,0.75)    #75%分位数
iqr=pt75-pt25    #内四分极值

#画横线和标注
draw_line(cht,0,minv,68,minv,'最小值')
draw_line(cht,0,maxv,68,maxv,'最大值')
draw_line(cht,0,meanv,60,meanv,'均值')
draw_line(cht,0,median,68,median,'中值')
draw_line(cht,0,gmean,45,gmean,'几何平均值')
draw_line(cht,0,hmean,45,hmean,'调和平均值')
draw_line(cht,0,tmean,45,tmean,'截尾平均值')
draw_line(cht,0,pt25,55,pt25,'25%分位数')
draw_line(cht,0,pt75,55,pt75,'75%分位数')
#... 省略部分代码

root=os.getcwd()    #获取当前工作路径
app=xw.App(visible=True,add_book=False)    #创建Excel应用
wb=app.books.open(root+r'/data.xlsx',read_only=False)    #打开数据文件返回工作簿对象
sht=wb.sheets('Sheet1')    #获取指定工作表对象
rng=[0 for _ in range(4)]
rng[0]=sht.range('A2:A201')    #各组数据对应的单元格区域对象
rng[1]=sht.range('B2:B201')
rng[2]=sht.range('C2:C201')
rng[3]=sht.range('D2:D201')
for i in range(4):
    cht=draw_hist(rng[i],200)    #绘制直方图
    set_style(cht)    #设置样式

root=os.getcwd()    #获取当前工作路径
app=xw.App(visible=True,add_book=False)    #创建Excel应用
wb=app.books.open(root+r'/data.xlsx',read_only=False)    #打开数据文件返回工作簿对象
sht=wb.sheets('Sheet1')    #获取指定工作表对象
data=[0 for _ in range(4)]
data[0]=sht.range('A2:A201').value    #获取数据
data[1]=sht.range('B2:B201').value
data[2]=sht.range('C2:C201').value
data[3]=sht.range('D2:D201').value
app.kill()    #退出应用
#从comtypes包中导入CreateObject函数
from comtypes.client import CreateObject
app2=CreateObject("Excel.Application")    #创建Excel应用
app2.Visible=True    #应用窗口可见
app2.ScreenUpdating=False
wb2=app2.Workbooks.Open(root+r'/data.xlsx')    #添加工作簿
sht2=wb2.Sheets('Sheet1')    #获取第1个工作表
minx=[0 for _ in range(4)]
maxx=[0 for _ in range(4)]
miny=[0 for _ in range(4)]
maxy=[0 for _ in range(4)]
minx[0]=-30    #4个图表横轴和纵轴的取值范围
maxx[0]=30
minx[1]=0
maxx[1]=1
minx[2]=-3
maxx[2]=4
minx[3]=-6
maxx[3]=4
maxy[0]=0.4
maxy[1]=1.1
maxy[2]=0.5
maxy[3]=0.5
for i in range(4):
    shp=sht2.Shapes.AddChart2()    #创建空白图表
    shp.Left=20+(i+1)*10
    cht=shp.Chart    #获取图表
    cht.ChartType=-4169    #散点图
    ax1=cht.Axes(1)    #获取横轴
    ax2=cht.Axes(2)    #获取纵轴
    ax1.MinimumScale=minx[i]    #设置横轴和纵轴的取值范围
    ax1.MaximumScale=maxx[i]
    ax2.MinimumScale= 0
    ax2.MaximumScale=maxy[i]
    cht.HasLegend=False    #不显示图例

    set_style(cht)    #设置样式

    ax1.CrossesAt=ax1.MinimumScale    #坐标轴相交位置
    ax2.CrossesAt=ax2.MinimumScale

    #绘制核密度估计曲线图
    draw_kde(cht,data[i],0,0,0,255,minx[i],maxx[i])