当前位置:网站首页>置信区间的画法
置信区间的画法
2022-07-04 21:14:00 【小蜗笔记】
err_inter = result.conf_int()
#置信区间的中心
err_mean = (err_inter[1]+err_inter[0])/2
# 置信区间的长度
err_long = (err_inter[1]-err_inter[0])/2
#置信区间画图
import matplotlib.pyplot as plt
from scipy import stats
import numpy as np
plt.style.use('seaborn-whitegrid')
fig, ax = plt.subplots(figsize=(50, 20),dpi=60)
ax.spines['bottom'].set_linewidth('2.0')#设置边框线宽为2.0
ax.spines['bottom'].set_color('black')
ax.spines['top'].set_linewidth('2.0')#设置边框线宽为2.0
ax.spines['top'].set_color('black')
ax.spines['right'].set_linewidth('2.0')#设置边框线宽为2.0
ax.spines['right'].set_color('black')
ax.spines['left'].set_linewidth('2.0')#设置边框线宽为2.0
ax.spines['left'].set_color('black')
plt.xticks(fontsize = 50)
plt.yticks(fontsize = 50)
for i in range(1,len(err_inter)):
ax.errorbar(i-10,err_mean[i] , yerr=err_long[i],
fmt='.', markersize=10, capsize=15, color='black')
xticks = [-9,-8, -7, -6, -5, -4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11]
xticklabes =['-9','-8','-7','-6','-5','-4','-3','-2','-1','0','1','2','3','4','5','6','7','8','9','10','11']
plt.xticks(xticks,xticklabes,size=50)
yticks = [-2,-1,0,1,2,3,4,5]
yticklabes =['-2','-1','0','1','2','3','4','5']
ax.plot([-9.5,11.5],[0,0], color='black',linestyle = '--',linewidth=2, markersize=12)
plt.xticks(xticks,xticklabes,size=50)
ax.set_ylim(-2.5, 5.5)
ax.set_xlim(-9.5, 11.5)
#ax.fill_betweenx([1, 0], facecolor='lightgrey', alpha=0.3)
#ax.legend(loc='best', fontsize=11, framealpha=1, frameon=True)
#ax.set_xticklabels()
#ax.set_ticklabels()
#ax.set_yticklabels(labels=labels)
ax.set_xlabel('time', fontsize=50)
ax.set_ylabel('coefficient', fontsize=50)
#ax.yaxis.set_major_formatter(plt.NullFormatter())
#fig.tight_layout();
#plt.show()
#plt.savefig(r'C:\Users\sc35131\Desktop\wangyao\Cell_consistence_alogrithm\picture9\\'+'1234',bbox_inches = 'tight',pad_inches = 0,dpi =350)
plt.savefig('平行趋势检验.png')
边栏推荐
- MP3是如何诞生的?
- gtest从一无所知到熟练使用(3)什么是test suite和test case
- Delphi SOAP WebService 服务器端多个 SoapDataModule 实现相同的接口方法,接口继承
- Jerry's ad series MIDI function description [chapter]
- Redis transaction
- redis发布订阅的使用
- 每日一题-LeetCode1200-最小绝对差-数组-排序
- redis RDB AOF
- Le module minidom écrit et analyse XML
- new IntersectionObserver 使用笔记
猜你喜欢
At the right time, the Guangzhou station of the city chain science and Technology Strategy Summit was successfully held
Three or two things about the actual combat of OMS system
做BI开发,为什么一定要熟悉行业和企业业务?
开源之夏专访|Apache IoTDB社区 新晋Committer谢其骏
超详细教程,一文入门Istio架构原理及实战应用
gtest从一无所知到熟练使用(3)什么是test suite和test case
如何借助自动化工具落地DevOps
Case sharing | integrated construction of data operation and maintenance in the financial industry
Daily question -leetcode1200- minimum absolute difference - array - sort
How is the entered query SQL statement executed?
随机推荐
Maidong Internet won the bid of Beijing life insurance
超详细教程,一文入门Istio架构原理及实战应用
创客思维在高等教育中的启迪作用
Solve the problem of data disorder caused by slow asynchronous interface
Shutter WebView example
Jerry's ad series MIDI function description [chapter]
OMS系统实战的三两事
QT—绘制其他问题
[weekly translation go] how to code in go series articles are online!!
更强的 JsonPath 兼容性及性能测试之2022版(Snack3,Fastjson2,jayway.jsonpath)
Bookmark
new IntersectionObserver 使用笔记
At the right time, the Guangzhou station of the city chain science and Technology Strategy Summit was successfully held
Jerry's ad series MIDI function description [chapter]
Redis 排查大 key 的3种方法,优化必备
案例分享|金融业数据运营运维一体化建设
Cadeus has never stopped innovating. Decentralized edge rendering technology makes the metauniverse no longer far away
Jerry added the process of turning off the touch module before turning it off [chapter]
Arcgis 10.2.2 | arcgis license server无法启动的解决办法
淘宝商品评价api接口(item_review-获得淘宝商品评论API接口),天猫商品评论API接口