統(tǒng)計學(xué)-方差分析實踐

題目
數(shù)據(jù)背景:有A删掀、B翔冀、C、D四個地區(qū)披泪,不同地區(qū)的銷售量不一樣纤子,現(xiàn)抽取了不同時間段內(nèi)每個地區(qū)的銷售量,試解決:

1款票、每個地區(qū)間的銷售量是否相同控硼?

2、不同月份的銷售量是否相同艾少?

3卡乾、不同時間與地區(qū)的銷售量是否相同?

       date     A     B     C     D
0 2019-12-01  8087  5923  9299  6287
1 2019-12-02  7851  5617  8984  6078
2 2019-12-03  4679  5481  8461  6735
3 2019-12-04  4296  6473  8407  1749
4 2019-01-04  7875  1688  5371  6377
5 2019-01-05  7271  6078  5923  6161
6 2019-01-06  9299  6735  5617  6861
7 2019-01-07  1749  2588  8407  1782
8 2019-03-08  3196  2595  9299  6820
9 2019-03-09  1220  1132  9627  6591

1缚够、每個地區(qū)間的銷售量是否相同幔妨?
用單因素方差分析鹦赎,先假設(shè),再構(gòu)造統(tǒng)計量误堡,最后進(jìn)行決策
(1)假設(shè)H0钙姊,即假設(shè)樣本A,B,C,D的均值相同,不同地區(qū)無顯著性影響埂伦。H1:樣本A,B,C,D均值不全相等煞额,有顯著性影響。
?(2)構(gòu)造統(tǒng)計量沾谜,因素為地區(qū)膊毁,水平為4,因變量為銷售量基跑。

import numpy as np
import scipy.stats as stats
import pandas as pd
#讀取數(shù)據(jù)
data=pd.read_excel("/Downloads/1銷售數(shù)據(jù).xlsx")

print(data.describe())  #統(tǒng)計數(shù)量婚温、均值、標(biāo)準(zhǔn)差媳否、上下四分位
print(data.mean().sum()/4)  #輸出對應(yīng)的總體均值
                A            B            C            D
count    11.000000    11.000000    11.000000    11.000000
mean   5150.454545  4329.909091  7457.909091  5588.727273
std    3049.346368  2079.064090  2234.142809  1912.142311
min    1132.000000  1132.000000  2642.000000  1749.000000
25%    2472.500000  2591.500000  5770.000000  6056.500000
50%    4679.000000  5481.000000  8407.000000  6287.000000
75%    7863.000000  6000.500000  9141.500000  6663.000000
max    9299.000000  6735.000000  9627.000000  6861.000000
5631.75

總平方和SST=95533231.1875+61867762.1875+86597368.1875+36583242.6875
組間平方和SSA=57894573.5196352
組內(nèi)平方和SSE=222687030.727273
MSA=SSA/(4-1)=19298191.17 , MSE =SSE/(N-K)=31812432.96,
檢驗統(tǒng)計量F=MSA/MSE=0.606624184 , 在給定顯著性水平a=0.05,在F分布表中查找分子自由度df1=4-1=3,分母自由度df2=n-k=11-4=7,相應(yīng)的臨界值Fa(3,7)=4.347,F<Fa 則不拒絕原假設(shè)H0,即認(rèn)為各地區(qū)對銷售量影響不顯著栅螟。
問題2. 不同月份銷售量不同?
1)假設(shè)H0:各月的均值相等篱竭,即各個月份對銷售量無顯著影響力图;H1:各個月份均值不完全相同,即不同月份對銷售量有顯著影響掺逼。

(2)構(gòu)造和計算統(tǒng)計量
(所有數(shù)據(jù)按照月份分組,2019年3月數(shù)據(jù)少4行吃媒,因此保留前12行分析數(shù)據(jù))
1)計算總體均值:5898.5,水平為3
2)總平方和SST=205887565
3)組間平方和SSA=39687746
4)組內(nèi)平方和SSE=SST?SSA=205887565?39687746=166199819
5)MSA=SSA/(3-1)=19843873
MSE=SSE/(n-k)=SSE/(12-3)=18466646.56
6)F=MSA\MSE=1.074579131
(3)決策分析吕喘,F(xiàn)<Fa=Fa(2,9)=4.256赘那,不拒絕原假設(shè),即各個月份影響不顯著氯质。

問題3. 不同時間與地區(qū)的銷售量是否相同募舟?
這里按照無交互作用的雙因素方差分析方法
(1)假設(shè)(2)構(gòu)造統(tǒng)計量和計算(3)決策

(1)行因素:H0:假設(shè)日期對銷售量無顯著影響;H1:日期對銷售量有顯著影響;
列因素:H0:假設(shè)地區(qū)對銷售量無顯著影響H1:地區(qū)對銷售量有顯著影響列因素闻察。
(2)構(gòu)造統(tǒng)計量并計算:




(3)從上表中可以看出Fr<Fa拱礁,F(xiàn)c<Fa,分別從地區(qū)蜓陌、日期都對銷售量影響不顯著觅彰。

參考:《統(tǒng)計學(xué)(第七版)》賈俊平
https://blog.csdn.net/long636/article/details/104437179

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末吩蔑,一起剝皮案震驚了整個濱河市钮热,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌烛芬,老刑警劉巖隧期,帶你破解...
    沈念sama閱讀 212,884評論 6 492
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件飒责,死亡現(xiàn)場離奇詭異,居然都是意外死亡仆潮,警方通過查閱死者的電腦和手機(jī)宏蛉,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,755評論 3 385
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來性置,“玉大人拾并,你說我怎么就攤上這事∨羟常” “怎么了嗅义?”我有些...
    開封第一講書人閱讀 158,369評論 0 348
  • 文/不壞的土叔 我叫張陵,是天一觀的道長隐砸。 經(jīng)常有香客問我之碗,道長,這世上最難降的妖魔是什么季希? 我笑而不...
    開封第一講書人閱讀 56,799評論 1 285
  • 正文 為了忘掉前任褪那,我火速辦了婚禮,結(jié)果婚禮上式塌,老公的妹妹穿的比我還像新娘博敬。我一直安慰自己,他們只是感情好峰尝,可當(dāng)我...
    茶點故事閱讀 65,910評論 6 386
  • 文/花漫 我一把揭開白布冶忱。 她就那樣靜靜地躺著,像睡著了一般境析。 火紅的嫁衣襯著肌膚如雪囚枪。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 50,096評論 1 291
  • 那天劳淆,我揣著相機(jī)與錄音链沼,去河邊找鬼。 笑死沛鸵,一個胖子當(dāng)著我的面吹牛括勺,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播曲掰,決...
    沈念sama閱讀 39,159評論 3 411
  • 文/蒼蘭香墨 我猛地睜開眼疾捍,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了栏妖?” 一聲冷哼從身側(cè)響起乱豆,我...
    開封第一講書人閱讀 37,917評論 0 268
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎吊趾,沒想到半個月后宛裕,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體瑟啃,經(jīng)...
    沈念sama閱讀 44,360評論 1 303
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 36,673評論 2 327
  • 正文 我和宋清朗相戀三年揩尸,在試婚紗的時候發(fā)現(xiàn)自己被綠了蛹屿。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 38,814評論 1 341
  • 序言:一個原本活蹦亂跳的男人離奇死亡岩榆,死狀恐怖错负,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情勇边,我是刑警寧澤湿颅,帶...
    沈念sama閱讀 34,509評論 4 334
  • 正文 年R本政府宣布,位于F島的核電站粥诫,受9級特大地震影響油航,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜怀浆,卻給世界環(huán)境...
    茶點故事閱讀 40,156評論 3 317
  • 文/蒙蒙 一谊囚、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧执赡,春花似錦镰踏、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,882評論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至首懈,卻和暖如春绊率,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背究履。 一陣腳步聲響...
    開封第一講書人閱讀 32,123評論 1 267
  • 我被黑心中介騙來泰國打工滤否, 沒想到剛下飛機(jī)就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人最仑。 一個月前我還...
    沈念sama閱讀 46,641評論 2 362
  • 正文 我出身青樓藐俺,卻偏偏與公主長得像,于是被迫代替她去往敵國和親泥彤。 傳聞我的和親對象是個殘疾皇子欲芹,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 43,728評論 2 351

推薦閱讀更多精彩內(nèi)容