统计学入门-小白和诸葛先生的故事

David Sun原创,转载请注明出处。

古有诸葛先生夜观天象,可以知未来,有一叶知秋,见微知著的贤人。

国外有费米先生从纸片的飞落速度估算原子弹的爆炸能量。

今天我们有大数据的大师专家,从看似无奇枯燥的数据中,找到蛛丝马迹,预测趋势,判断分析根源和差异,决策商业未来。

小白一定很好奇,诸葛先生如何从天象猜测出未来若干几天的气候和风向,这个迷我们也先保密。

今天我们先从量化分析的角度,举一些案例,配上一些工具和算法,讲讲这个大千世界的数据和统计入门:

案例1: 小白在工厂管理5条生产线,某一天发现生产线1和5的故障率好像有点高,但是也不能有充分的依据证明,小白(产线质量主管)这时候来请教诸葛先生(统计学顾问)了。

小白:我不能瞎猜对吧,厂长会说我无凭无据,工人会说我找茬。

诸葛:对哦

小白:我们每天产线有30000个产品生产组装好,我每个都检查下,不得累死我。。。工人们会笑话死我。。。

诸葛:来来,我来给你把把脉,按照我的SOP走一遍先:

①你看看这些产品是否是同一批次,产品是否都规格一致

②找出关键的质量评判标准,例如产品的重量是否都严格控制在50±3g之间

③确定参数的分布类型,假设重量属于正态分布,那么选择抽样标准-随机抽样,水平抽样等;

④选择随机抽样方法,从5条生产线上分别抽选30个样品,各自代表了五条线五万个总体的一部分。

⑤采用t检验方法,检验150个样品各组之间组间的平均值是否有统计学显著性差异。

一共5组,每组之间都是相对独立的样本,随机抽样,没有配对。我们先找对比D,即2-5组间对比做案例。

两两对比


⑥建立假设:2组和5组的产线上,生产出来的产品平均来看都是一样的,差不多的。

即H0:μ2=μ5,H1:μ2≠μ5,α=0.05

⑦拿起我的钉耙,哦不,方天画戟,哦不,SPSS神器,或者Excel小刀一试。

⑨输入数据,用t检验算出sig(双侧),即显著性,看到方差相等,0.32大于0.05,方差不相等时候,也大于0.05。

⑩根据我们上面提到的α限值,我们得出统计学结论,即2组和5组生产的产品平均值相等,在统计学上有意义。我们认为2和5两天产线上的产品重量没有显著性的差异,基本上认为是一样的。小白:诸葛先生,我糊涂了,为什么有2个值呢,方差相等不相等怎么检查呢?

诸葛:两个样本方差相等与否,可以通过F检验来判定是否方差相等,我们从上面表格看到F值是6.617,SIG是0.013,可以认为两组之间方差不相等,有显著差异。SPSS里面的Levene检验是方差分析中检验方差相等与否的方法之一。

小白:这么多组和组之间,一个一个对比,好累人哦

诸葛:要不怎么说SPSS是神器呢,看下图


把这些组的数据全部加入检验变量

然后可以得出一系列组和组之间的t和F检验结果,下图仅为示例,和上图数据无关联性。

小白:好像懂一些了。但是如果我检查的重量到底是不是正态分布呢?

诸葛:神器就是神器

方法1:

看Q-Q图(不是腾讯的QQ)

如果数据点大部分都在直线上,可以近似认为是正态分布。

方法2:用两位苏联数学家Kolmogorov和Smirnov的名字命名的K-S检验

上面的显著性0.200,大于0.05,即认为和正态分布无显著性差异,属于正态分布。

方法3。。。

小白:先生住口,我先懂一些就好。但是如果我要看的不是产品重量,而是产品在市场上客户的评价呢?比如我在市场上找了100个用户,给产品的功能和性能打分,或者使用3个月以后的前后对比,这些数据的分布不一定是正态分布吧,那怎么分析呢?

诸葛:这时候就要用上非参数检验了,就是在总体方差未知或知道甚少的情况下,利用样本数据对总体分布形态等进行推断。

小白:比如我的产品是治疗高血压的药物,做临床实验,在上海北京几个城市,随机找了200个高血压的患者,分成2组,每组100人,分别用市场上已有的药物A和我们的产品药物B,跟踪记录患者使用药物前的血压,1个月以后,2个月以后和3个月以后的数据,然后能否对比出来,两种药物的效果差异呢?

诸葛:首先,我们不知道这些药物使用后的效果,血压和血压下降值是什么分布类型,所以要从非参数检验来分析。

其次,这200个样本,并没有配对研究,是随机的样本,所以是独立样本。把所有的数据输入到SPSS之后,就可以进行组件分布分析。

如果是2个样本组,那么用曼惠特尼检验,如果是3种或更多种药物,分成n组进行对比,那就要用Kruskal-Wallis检验。

如果要进一步分析,1-2-3月是否有血压递减的趋势,可以使用Jonckheere-Terpstra检验,看J-T值和显著性差异。

小白:好像有点晕了。。好像很有趣的样子。但是那么多数据,从哪来啊。

诸葛:数据统计博大精深,应用面极广,不局限于经济学,管理,金融股票投资期货分析,医药研究,生物检测,产品质量管控,气象学预测,地质学预测,社会调查和研究,人口和政策研究等等等等。大数据目前很热门啊,下次给你讲讲如何用python爬虫去互联网上找到一些神奇的数据吧。

小白:先生,我先去洗把脸。。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,539评论 6 497
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,911评论 3 391
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 161,337评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,723评论 1 290
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,795评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,762评论 1 294
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,742评论 3 416
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,508评论 0 271
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,954评论 1 308
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,247评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,404评论 1 345
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,104评论 5 340
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,736评论 3 324
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,352评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,557评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,371评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,292评论 2 352

推荐阅读更多精彩内容