国产秋霞理论久久久电影-婷婷色九月综合激情丁香-欧美在线观看乱妇视频-精品国avA久久久久久久-国产乱码精品一区二区三区亚洲人-欧美熟妇一区二区三区蜜桃视频

最常用的的17種假設(shè)檢驗(附Python實例)

共 8520字,需瀏覽 18分鐘

 ·

2020-11-28 00:51

常用的假設(shè)檢驗

在統(tǒng)計學(xué)中,存在著數(shù)百種假設(shè)檢驗。而在數(shù)據(jù)分析、機器學(xué)習(xí)項目中,只有一小部分較常使用。本文介紹了17種常用的假設(shè)檢驗,包括適用場景及使用Python API的例子。讓我們一起來學(xué)習(xí)吧!

目錄

  • 1. 正態(tài)分布檢驗

    • 1.1 Shapiro-Wilk Test(W 檢驗)

    • 1.2 D’Agostino’s Test

    • 1.3 Anderson-Darling Test

  • 2. 相關(guān)性檢驗

    • 2.1 Pearson相關(guān)系數(shù)

    • 2.2 Spearman等級相關(guān)

    • 2.3 Kendall等級相關(guān)

    • 2.4 Chi-Squared Test卡方檢驗

  • 3. 參數(shù)檢驗

    • 3.1 學(xué)生t檢驗

    • 3.2 配對樣本的t檢驗

    • 3.3 方差分析ANOVA

    • 3.4 重復(fù)觀測的方差分析

  • 4. 非參數(shù)檢驗

    • 4.1 Mann-Whitney U Test

    • 4.2 Wilcoxon Signed-Rank Test

    • 4.3 Kruskal-Wallis H Test

    • 4.4 Friedman Test

  • 5. 時間序列的平穩(wěn)性檢驗

    • 5.1 Augmented Dickey-Fuller Unit Root Test

    • 5.2 Kwiatkowski-Phillips-Schmidt-Shin

1. 正態(tài)分布檢驗

1.1 Shapiro-Wilk Test(W 檢驗)

用于檢驗樣本數(shù)據(jù)是否來自服從正態(tài)分布的總體。

在實際應(yīng)用中,W 檢驗被認(rèn)為是一個可靠的正態(tài)性檢驗,但是也有人認(rèn)為該檢驗更適用于較小的數(shù)據(jù)樣本(數(shù)千個觀測值以內(nèi))。

使用前提:

  • 各樣本觀察值為獨立同分布的

原假設(shè):樣本數(shù)據(jù)服從正態(tài)分布

Python實現(xiàn)代碼:

#?Shapiro-Wilk?Normality?Test
from?scipy.stats?import?shapiro
data?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
stat,?p?=?shapiro(data)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布')
else:
????print('不服從正態(tài)分布')

1.2 D’Agostino’s Test

用于檢驗樣本數(shù)據(jù)是否來自服從正態(tài)分布的總體。

D’Agostino’s Test是通過計算樣本數(shù)據(jù)的峰度和偏度,來判斷其分布是否偏離正態(tài)分布。

  • 偏度是對數(shù)據(jù)分布對稱性的測度,衡量數(shù)據(jù)分布是否左偏或右偏。
  • 峰度是對數(shù)據(jù)分布平峰或尖峰程度的測度,它是一種簡單而常用的正態(tài)性統(tǒng)計檢驗量。

使用前提:

  • 各樣本觀察值為獨立同分布的

原假設(shè):樣本數(shù)據(jù)服從正態(tài)分布

Python實現(xiàn)代碼:

#?D'Agostino's?K^2?Normality?Test
from?scipy.stats?import?normaltest
data?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
stat,?p?=?normaltest(data)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布')
else:
????print('不服從正態(tài)分布')

1.3 Anderson-Darling Test

用于檢驗樣本數(shù)據(jù)是否服從某一已知分布。

該檢驗修改自一種更復(fù)雜的非參數(shù)的擬合良好的檢驗統(tǒng)計(Kolmogorov-Smirnov Test)。SciPy中的anderson()函數(shù)實現(xiàn)了Anderson-Darling檢驗,函數(shù)參數(shù)為樣本數(shù)據(jù)及要檢驗的分布名稱,默認(rèn)情況下,為'norm'正態(tài)分布,還支持對'expon'指數(shù)分布,'logistic'分布,以及'gumbel'耿貝爾分布的檢驗,它會返回一個包含不同顯著性水平下的p值的列表,而不是一個單一的p值,因此這可以更全面地解釋結(jié)果。

使用前提:

  • 各樣本觀察值為獨立同分布的

原假設(shè):樣本數(shù)據(jù)服從某一已知分布

Python實現(xiàn)代碼:

#?Anderson-Darling?Normality?Test
from?scipy.stats?import?anderson
data?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
result?=?anderson(data)
print('stat=%.3f'?%?(result.statistic))
for?i?in?range(len(result.critical_values)):
????sl,?cv?=?result.significance_level[i],?result.critical_values[i]
????if?result.statistic?????????print('顯著性水平為%.2f時,P值為%.1f,不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布'?%?(sl/100,?cv))
????else:
????????print('顯著性水平為%.2f時,P值為%.1f,?拒絕原假設(shè),樣本數(shù)據(jù)不服從正態(tài)分布'?%?(sl/100,?cv))

#?stat=0.424
#?顯著性水平為0.15時,P值為0.5,不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布
#?顯著性水平為0.10時,P值為0.6,不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布
#?顯著性水平為0.05時,P值為0.7,不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布
#?顯著性水平為0.03時,P值為0.8,不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布
#?顯著性水平為0.01時,P值為0.9,不能拒絕原假設(shè),樣本數(shù)據(jù)服從正態(tài)分布

2. 相關(guān)性檢驗

2.1 Pearson相關(guān)系數(shù)

用于檢驗兩樣本數(shù)據(jù)之間線性關(guān)系的強度。

該檢驗將兩個變量之間的協(xié)方差進(jìn)行歸一化處理以給出可解釋的分?jǐn)?shù),為一個介于-1到1之間的值,-1表示完全負(fù)相關(guān),1表示完全正相關(guān),0表示沒有相關(guān)性。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 樣本數(shù)據(jù)服從正態(tài)分布
  • 每個樣本觀測值的方差相等

原假設(shè):兩變量相互獨立

Python實現(xiàn)代碼:

#?Pearson's?Correlation?test
from?scipy.stats?import?pearsonr
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[0.353,?3.517,?0.125,?-7.545,?-0.555,?-1.536,?3.350,?-1.578,?-3.537,?-1.579]
stat,?p?=?pearsonr(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('兩變量相互獨立')
else:
????print('兩變量可能存在線性相關(guān)關(guān)系')

2.2 Spearman等級相關(guān)

用于檢驗兩變量是否具有單調(diào)關(guān)系。

當(dāng)兩變量因非線性關(guān)系相關(guān),或者不服從正態(tài)分布時,Spearman相關(guān)系數(shù)可以用來反映變量間的相關(guān)性強度。如果存在線性關(guān)系,也可以使用這種方法來檢驗,但是可能導(dǎo)致計算出的相關(guān)系數(shù)較低。

Spearman不是使用樣本數(shù)據(jù)本身的協(xié)方差和標(biāo)準(zhǔn)差來計算相關(guān)系數(shù)的,而是根據(jù)樣本值的相對秩次來計算統(tǒng)計量,這是非參數(shù)統(tǒng)計中常用的方法。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 各樣本數(shù)據(jù)可定序

原假設(shè):兩變量相互獨立

Python實現(xiàn)代碼:

#?Spearman's?Rank?Correlation?Test
from?scipy.stats?import?spearmanr
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[0.353,?3.517,?0.125,?-7.545,?-0.555,?-1.536,?3.350,?-1.578,?-3.537,?-1.579]
stat,?p?=?spearmanr(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('兩變量相互獨立')
else:
????print('兩變量可能存在相關(guān)關(guān)系')

2.3 Kendall等級相關(guān)

用于檢驗兩變量是否具有單調(diào)關(guān)系。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 各樣本數(shù)據(jù)可定序

原假設(shè):兩變量相互獨立

Python實現(xiàn)代碼:

#?Kendall's?Rank?Correlation?Test
from?scipy.stats?import?kendalltau
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[0.353,?3.517,?0.125,?-7.545,?-0.555,?-1.536,?3.350,?-1.578,?-3.537,?-1.579]
stat,?p?=?kendalltau(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('兩變量相互獨立')
else:
????print('兩變量可能存在相關(guān)關(guān)系')

2.4 Chi-Squared Test卡方檢驗

用于檢驗兩分類變量是否相關(guān)。

卡方檢驗的零假設(shè)是一個分類變量的觀測頻數(shù)與該分類變量的期望頻數(shù)相吻合。檢驗統(tǒng)計量服從卡方分布。

使用前提:

  • 用于計算兩變量的列聯(lián)表的觀測值是獨立的
  • 列聯(lián)表的每個單元格的期望計數(shù)不小于5

原假設(shè):兩變量相互獨立

Python實現(xiàn)代碼:

#?Chi-squared?test?with?similar?proportions
from?scipy.stats?import?chi2_contingency
from?scipy.stats?import?chi2

#?列聯(lián)表
table?=?[[10,?20,?30],
????????[6,??9,??17]]
print('列聯(lián)表')
print(table)

stat,?p,?dof,?expected?=?chi2_contingency(table)
print('自由度dof=%d'?%?dof)
print('期望分布')
print(expected)

#?[[10.43478261?18.91304348?30.65217391]
#??[?5.56521739?10.08695652?16.34782609]]

#?采用統(tǒng)計量推斷
prob?=?0.95
critical?=?chi2.ppf(prob,?dof)
print('probability=%.3f,?critical=%.3f,?stat=%.3f'?%?(prob,?critical,?stat))
if?abs(stat)?>=?critical:
????print('拒絕原假設(shè),兩變量存在相關(guān)關(guān)系')
else:
????print('不能拒絕原假設(shè),兩變量相互獨立')

#?采用p值推斷
alpha?=?1.0?-?prob
print('significance=%.3f,?p=%.3f'?%?(alpha,?p))
if?p?<=?alpha:
????print('拒絕原假設(shè),兩變量存在相關(guān)關(guān)系')
else:
????print('不能拒絕原假設(shè),兩變量相互獨立')

3. 參數(shù)檢驗

3.1 學(xué)生t檢驗

用于檢驗兩獨立樣本數(shù)據(jù)的均值是否存在顯著差異。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 樣本數(shù)據(jù)服從正態(tài)分布
  • 每個樣本觀測值的方差相等

原假設(shè):兩樣本集的均值相同

Python實現(xiàn)代碼:

#?Student's?t-test
from?scipy.stats?import?ttest_ind
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
stat,?p?=?ttest_ind(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),兩樣本集分布相同')
else:
????print('拒絕原假設(shè),樣本集分布可能不同')

3.2 配對樣本的t檢驗

用于檢驗兩配對樣本的均值是否存在顯著差異。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 樣本數(shù)據(jù)服從正態(tài)分布
  • 每個樣本觀測值的方差相等
  • 樣本的數(shù)據(jù)是一一對應(yīng)的關(guān)系

原假設(shè):各樣本集的均值相同

Python實現(xiàn)代碼:

#?Paired?Student's?t-test
from?scipy.stats?import?ttest_rel
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
stat,?p?=?ttest_rel(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),兩樣本集分布相同')
else:
????print('拒絕原假設(shè),兩樣本集分布可能不同')

3.3 方差分析ANOVA

用于檢驗兩個及以上獨立樣本數(shù)據(jù)的均值是否存在顯著差異。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 樣本數(shù)據(jù)服從正態(tài)分布
  • 每個樣本觀測值的方差相等

原假設(shè):各樣本集均值相同

Python實現(xiàn)代碼:

#?Analysis?of?Variance?Test
from?scipy.stats?import?f_oneway
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
data3?=?[-0.208,?0.696,?0.928,?-1.148,?-0.213,?0.229,?0.137,?0.269,?-0.870,?-1.204]
stat,?p?=?f_oneway(data1,?data2,?data3)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),樣本集分布相同')
else:
????print('拒絕原假設(shè),樣本集分布可能不同')

3.4 重復(fù)觀測的方差分析

用于檢驗兩或多配對樣本的均值是否存在顯著差異。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 樣本數(shù)據(jù)服從正態(tài)分布
  • 每個樣本觀測值的方差相等
  • 樣本的數(shù)據(jù)是一一對應(yīng)的關(guān)系

Python中Scipy暫不支持該檢驗

4. 非參數(shù)檢驗

4.1 Mann-Whitney U Test

用于檢驗兩獨立樣本分布是否相同。該檢驗的原理是,將兩樣本集合并起來,一起排序,然后計算兩樣本的值是隨機混合排序的,還是聚集的。秩序隨機則代表兩樣本集分布無差異,而出現(xiàn)樣本值聚集情況則表明有差異。為了使檢驗有效,要求每個樣本數(shù)據(jù)至少有20個觀測值。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 各樣本集中數(shù)據(jù)可定序

原假設(shè):兩樣本集分布相同

Python實現(xiàn)代碼:

#?Mann-Whitney?U?Test
from?scipy.stats?import?mannwhitneyu
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
stat,?p?=?mannwhitneyu(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),兩樣本集分布相同')
else:
????print('拒絕原假設(shè),樣本集分布可能不同')

4.2 Wilcoxon Signed-Rank Test

用于檢驗兩配對樣本分布是否相同。

樣本不獨立,因此不能使用Mann-Whitney U檢驗,而是使用Wilcoxon符號秩檢驗。這相當(dāng)于配對樣本的學(xué)生t檢驗,但是適用于定序數(shù)據(jù),而不是服從正態(tài)分布的實測數(shù)據(jù)。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 各樣本集中數(shù)據(jù)可定序
  • 樣本的數(shù)據(jù)是一一對應(yīng)的關(guān)系

原假設(shè):兩樣本集分布相同

Python實現(xiàn)代碼:

#?Wilcoxon?Signed-Rank?Test
from?scipy.stats?import?wilcoxon
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
stat,?p?=?wilcoxon(data1,?data2)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),兩樣本集分布相同')
else:
????print('拒絕原假設(shè),樣本集分布可能不同')

4.3 Kruskal-Wallis H Test

用于檢驗兩個及以上獨立樣本的分布是否相同。

Kruskal-Wallis 檢驗是單因素方差分析檢驗的非參數(shù)版本,可以用來確定兩個以上的獨立樣本是否有不同的分布,可以被認(rèn)為是Mann-Whitney U檢驗的一般化。該檢驗的零假設(shè)是“所有樣本數(shù)據(jù)從同一個分布中抽取的”,若拒絕原假設(shè),則表明有足夠證據(jù)證明至少有一個樣本與其他樣本不同,但是該假設(shè)并不能確定具體是哪幾個不同。因此,如果需要確定,可再使用Mann-Whitney檢驗進(jìn)行兩兩對比。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 各樣本集中數(shù)據(jù)可定序

原假設(shè):各樣本集的分布相同 備擇假設(shè):至少有一個樣本集與其他樣本不同

Python實現(xiàn)代碼:

#?Kruskal-Wallis?H?Test
from?scipy.stats?import?kruskal
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
data3?=?[-0.208,?0.696,?0.928,?-1.148,?-0.213,?0.229,?0.137,?0.269,?-0.870,?-1.204]
stat,?p?=?kruskal(data1,?data2,?data3)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),樣本集分布相同')
else:
????print('拒絕原假設(shè),樣本集分布可能不同')

4.4 Friedman Test

用于檢驗兩個及以上配對樣本的分布是否相同。Friedman Test是重復(fù)測量方差分析的非參數(shù)版本,需要每組至少有10個樣本,當(dāng)?shù)贸鲲@著性結(jié)果,表明至少有一個樣本與其他樣本存在差異。

使用前提:

  • 各樣本觀察值為獨立同分布的
  • 各樣本集中數(shù)據(jù)可定序
  • 樣本的數(shù)據(jù)是一一對應(yīng)的關(guān)系

原假設(shè):各樣本集的分布相同

Python實現(xiàn)代碼:

#?Friedman?Test
from?scipy.stats?import?friedmanchisquare
data1?=?[0.873,?2.817,?0.121,?-0.945,?-0.055,?-1.436,?0.360,?-1.478,?-1.637,?-1.869]
data2?=?[1.142,?-0.432,?-0.938,?-0.729,?-0.846,?-0.157,?0.500,?1.183,?-1.075,?-0.169]
data3?=?[-0.208,?0.696,?0.928,?-1.148,?-0.213,?0.229,?0.137,?0.269,?-0.870,?-1.204]
stat,?p?=?friedmanchisquare(data1,?data2,?data3)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('不能拒絕原假設(shè),樣本集分布相同')
else:
????print('拒絕原假設(shè),樣本集分布可能不同')

5. 時間序列的平穩(wěn)性檢驗

下面是針對時間序列數(shù)據(jù)的平穩(wěn)性檢驗。直、當(dāng)數(shù)據(jù)沒有呈現(xiàn)明顯的模式特征的話(趨勢性、季節(jié)性),可認(rèn)為它是平穩(wěn)的。平穩(wěn)序列中各觀察值在某個固定水平上波動,雖然在不同時間段波動程度不同,但是并不存在某種規(guī)律,可認(rèn)為是隨機的。

5.1 Augmented Dickey-Fuller Unit Root Test

用于檢驗時間序列是否具有單位根,例如是否具有趨勢或者更普遍的自回歸。單位根是時間序列的一種特性,使時間序列具有非平穩(wěn)性。單位根的存在意味著時間序列是非平穩(wěn)的。

使用前提:

  • 觀察值在時間上是有序的

原假設(shè):存在一個單位根,序列是非平穩(wěn)的 備擇假設(shè):不存在單位根,序列平穩(wěn) (注意:該檢驗的原假設(shè)與下一個檢驗的假設(shè)是不同的;如果p值小于顯著性水平,認(rèn)為該時間序列平穩(wěn))

Python實現(xiàn)代碼:

#?Augmented?Dickey-Fuller?unit?root?test
from?statsmodels.tsa.stattools?import?adfuller
data?=?[0,?1,?2,?3,?4,?5,?6,?7,?8,?9]
stat,?p,?lags,?obs,?crit,?t?=?adfuller(data)
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('時間序列是不平穩(wěn)的')
else:
????print('時間序列平穩(wěn)')

5.2 Kwiatkowski-Phillips-Schmidt-Shin

用于檢驗時間序列在確定性趨勢上的平穩(wěn)性。

使用前提:

  • 觀察值在時間上是有序的

Python實現(xiàn)代碼:

注意,在默認(rèn)情況下,該函數(shù)參數(shù)regression='c',代表 原假設(shè):時間序列是平穩(wěn)的。(圍繞一個常數(shù)) 備擇假設(shè):時間序列是不平穩(wěn)的。

#?Kwiatkowski-Phillips-Schmidt-Shin?test
#?The?data?is?stationary?around?a?constant?(default).
from?statsmodels.tsa.stattools?import?kpss
data?=?[0,?1,?2,?3,?4,?5,?6,?7,?8,?9]
stat,?p,?lags,?crit?=?kpss(data,?regression='c')
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('時間序列是平穩(wěn)的')
else:
????print('時間序列是不平穩(wěn)的')

#?stat=0.410,?p=0.073
#?時間序列是平穩(wěn)的

而使參數(shù)regression='ct',可檢驗時間序列在確定趨勢上的平穩(wěn)性。

#?Kwiatkowski-Phillips-Schmidt-Shin?test
#?The?data?is?stationary?around?a?trend
from?statsmodels.tsa.stattools?import?kpss
data?=?[0,?1,?2,?3,?4,?5,?6,?7,?8,?9]
stat,?p,?lags,?crit?=?kpss(data,?regression='ct')
print('stat=%.3f,?p=%.3f'?%?(stat,?p))
if?p?>?0.05:
????print('時間序列是趨勢穩(wěn)定的')
else:
????print('時間序列是趨勢不穩(wěn)定的')

#?stat=0.433,?p=0.010
#?時間序列是趨勢不穩(wěn)定的

參考資料:

https://machinelearningmastery.com/statistical-hypothesis-tests-in-python-cheat-sheet/

·················END·················

推薦閱讀

  1. 說說心里話

  2. 寫給所有數(shù)據(jù)人。

  3. 從留存率業(yè)務(wù)案例談0-1的數(shù)據(jù)指標(biāo)體系

  4. NB,真PDF神處理工具!

  5. 超級菜鳥如何入門數(shù)據(jù)分析?


歡迎長按掃碼關(guān)注「數(shù)據(jù)管道」

瀏覽 310
點贊
評論
收藏
分享

手機掃一掃分享

分享
舉報
評論
圖片
表情
推薦
點贊
評論
收藏
分享

手機掃一掃分享

分享
舉報

感谢您访问我们的网站,您可能还对以下资源感兴趣:

国产秋霞理论久久久电影-婷婷色九月综合激情丁香-欧美在线观看乱妇视频-精品国avA久久久久久久-国产乱码精品一区二区三区亚洲人-欧美熟妇一区二区三区蜜桃视频 国产最新在线视频| 中文字幕精品视频在线观看| 高H网站| 国产精品色婷婷99久久精品| 91成人做爰A片| 人人妻人人澡人人爽人人DVD| 国内精品卡一卡二卡三| 久操视频免费在线观看| 日韩精品成人av| 欧美一区二区精品| 久久国产精品在线| 国产熟女一区二区久久| jizz亚洲| 影音先锋男人站| 啪啪视频免费观看| www四虎| 熟妇槡BBBB槡BBBB图| 成人无码网站在线观看| 嗯啊在线视频| 亚洲AV无码成人精品区东京热| 久久天堂| 樱桃av| 38t6电影网logo高清图片 | 影音先锋成人资源站| 一曲二曲三曲在线观看中文字| 亚洲午夜免费视频| 东北女人操逼视频| 一级黄色电影在线观看| 亚洲AV成人片色在线观看麻豆 | 免费一区二区三区四区| 91日本| 成人黄色av| 精品国产久久久久久| 色噜噜狠狠一区二区三区| 国产精品熟女| 欧美日韩视频| 中文字幕福利| 国产麻豆一区二区三区| 欧美日韩一区二区三区视频| 午夜黄色操逼视频| 欧美性爱网址| 天堂网在线视频| 亚洲AV成人精品一区二区三区| 亚洲天堂日本| 俺去俺来也在线www色情网| 国产精品v欧美精品v日韩精品| 大鸡吧在线观看| 久久骚| 精品无码专区| 亚洲特级毛片| 国产伦精一品二品三品app| 中文字幕手机在线视频| 国产精品99视频| 欧美老女人操逼视频| 国产h在线播放| 清清草在线视频| 先锋资源国产| 爱爱亚洲| 黄色av网| 一级a一级a爰片免费免免中国A片| 97色色婷婷五月天| 欧美在线日韩在线| 99久久综合国产精品二区| 欧美肏屄网| 大鸡巴视频在线| 日韩AV无码免费| 国产婷婷色一区二区| 风情万种AV| 成人大香蕉| 一区在线观看| 一区二区三区精品| 免费视频一区二区| 久久精品大屁股| 在线免费毛片| 亚洲免费毛片| 欧美日韩黄色片| 久久久久无码国产精品不卡| 欧美日韩免费视频| 亚洲三级网| 中文字字幕在线中文| 视频二区| 美女福利在线| 久操福利视频| 不卡视频一区二区| 99热这里有精品| 欧美丰满老熟妇XXXXX性| 黄色视频久久| 婷婷五月天黄色| 色色a| 操啊操| 久久夜色精品国产欧美乱极品| 一本色道久久综合熟妇| 欧美另类激情| 18XXX亚洲HD护士JD| 久久蜜桃| 91人人澡| 国产中文字幕亚洲综合欧美| 一级做a爰片毛片A片| 国产熟女一区二区三区五月婷| 日韩无码久| 午夜视频99| 亚洲欧洲免费看| 影音先锋乱伦电影| 鲁鲁鲁鲁鲁鲁鲁777777| 亚洲中文字幕成人| 亚洲免费小电影| 日本人妻在线观看| 久久婷婷婬片A片AAA| 午夜福利视频网站| 国产69精品久久久久久| 又黄又爽无遮挡| 一本色道久久综合狠狠躁的推荐| 日批免费网站| 欧美一级A| 人人澡人人澡| 无码成人AV在线看免费| 欧美激情色色| 亚洲小电影在线观看| 大鸡巴操B视频| 怡红院成人在线| 特级特黄AAAA免费看| 成人免费操| 刘玥一区二区| 久久久久久久麻豆| 爆乳尤物一区二区三区| 99久久婷婷国产综合精品hsex,亚 91小宝寻花一区二区三区三级 | 亚洲jiZZjiZZ日本少妇| 国产视频123区| 波多野结衣av在线观看| 99精品免费在线观看| 国产超碰| 91人人妻人人澡人人爽人人| 99人妻人人爽人人添人人精品| 嫩草av| 尻屄视频| 农村一级婬片A片| 中文字幕特黄A片| 一区免费在线观看| 小明成人免费视频| 岛国av免费看| 色婷婷Av| 日韩欧美中文字幕在线观看| 99视频精品| 日韩无码砖区| 婷婷日逼| 欧美在线视频a| 国产操b视频| www.激情| 久久99精品视频| 国产成人久久| 干欧美美女| 乱伦91视频| 亚洲无码中文视频| 日韩中文无码字幕| 中文字幕23页| 中文av字幕| 另类老妇奶性BBWBBw| 91丨九色丨蝌蚪丨肥女| 成人在线三级| 欧美日韩一道本| 91吴梦梦一区二区传媒| 成人av一区| 大香蕉福利视频| 欧美毛视频| 大鸡吧大香蕉| 中韩日美免费看的电影| 人人爱人人摸人人操| 国产无码在线看| 日屄视频在线观看| 免费岛国av大片| 亚洲在线无码播放| 国产一级a毛一级a毛观看视频网站www.jn | 操b视频在线播放| 免费无码一级A片大黄在线观看 | 老师机性爱视频在线播放| 日本特黄一级片| 久久久久久久9999| 日韩高清无码毛片| 操啊操| 日韩成人电影| 亚洲欧美婷婷五月色综合| 亚洲伊人影院| 麻豆日韩| 免费看黃色AAAAAA片| 亚洲无码一区二区三区| 秋霞A片| 国产激情欧洲在线观看一区二区三区 | 97香蕉久久夜色精品国产| av久操| 99热在线观看免费| 国产乱伦电影| 一级做a爰片毛片A片| 精品久草| 自拍偷拍精品视频| 日本不卡中文字幕| 色九九| 国产精品一区二区黑人巨大| 青青免费在线视频| 婷婷亚洲色| 黄色精品视频| 亚洲美女网站在线观看| 欧美天天撸| 无码日韩视频| 国产黄色自拍视频| 黄色录像一级带| 蜜臀久久精品久久久久| 伊人成年网| 亚洲第二页| 国产91在线中日| 久久久久亚洲精品| 黄片免费看网站| 黄片aaa| 久久中文字幕电影| 翔田千里无码免费播放| 午夜爽爽爽| 亚洲精品内射| 91丨牛牛丨国产人妻| 天天日日日干| 亚洲视频免费在线播放| 天天干夜夜骑| 伊人性视频| 久久久波多野结衣| 亚洲综合天堂| 精品乱子伦一区二区三区,亚洲国产成| 天天干天天射天天| AA丁香综合激情| 亚洲一级在线观看| 搡BBBB搡BBB搡五十粉嫩| 一边做一边说国语对白| 中文字幕不卡无码| 乱伦专区| 国产AV综合网| 中文在线资源| 一道本在线视频| 亚洲精品无码在线播放| a√免费看| 亚州精品人妻一二三区| 91爱爱| 国产精品一区二区免费| 岛国无码破解AV在线播放 | 久久影院三级片| 久久精品久久久久久久| 日韩一二三区| 熟女人妻一区二区三区| 黄色电影av| 成人免费在线| 怡春院中文字幕| jizz丝袜| 九哥草逼网| 雾水情缘电影港片| 国产欧美日韩综合在线视频| 午夜激情五月天| 99久久精品国产毛片| 亚洲欧美日韩成人| 欧美三级黄色| 波多野结衣视频在线播放| 91乱了伦国产乱子伦| 久草热视频| 精品内射| 毛片a级| 日韩中文字幕免费| 欧美一级黄色片| 蜜桃Av噜噜| 黄片日逼| 午夜精品18码视频国产17c| 亚洲日韩精品无码| 中文字幕在线有码| 婷婷丁香人妻天天爽| 国产毛片久久久久久久| 久草综合在线| 久热只有精品| 豆花视频免费观看| 亚洲欧美成人在线| 国产婬片lA片www777| 国产精品a久久久久| 日韩不卡| 久久看片| 国产TS丝袜人妖系列视频| 俺也去com| 久久一区二区三区四区五区| 亚洲骚妇| 免费看黄色毛片| 中文字幕黄色电影| 日韩不卡精品| 欧美久久大香蕉| 操逼99| 99精品亚洲| 欧美成人性爱网站| 强波多野结衣黑人| 青青操成人在线视频| 蜜桃av秘无码一区二区三区| 国产狼友| 欧美在线免费观看| 蕉蕉视| 久久久久久久久久国产精品| 爽好紧别夹喷水网站| 亚洲国产精品成人综合色在线婷婷 | 欧美色色色色色| 操逼日韩欧美| 日韩精品视频一区二区三区| 美女被操免费网站| 日皮视频| 午夜AV影院| AV天堂小说网| 中文字幕无码日韩| 亚洲女与黑人正在播放| 俺也来www俺也色com| 色一区二区三区| 欧美系列在线| 日韩AV免费在线| 亚洲国产精品成人网站| 欧美V∧| 蜜臀av在线| 97国产精品视频人人做人人爱| 国产91白丝在一线播放| 日韩AV网站在线观看| 欧美成人激情视频| 亚洲高清福利| 精品欧美乱码久久久久久| 性感成人在线| 四虎成人无码A片观看| 欧美成人黄色电影| 色婷婷国产精品视频| 青青草免费公开视频| 高潮91PORN蝌蚪九色| 成人性爱av| 黄色的视频网站| 久久不卡视频| 青操AV| 国产精品自拍偷拍| 无码不卡视频在线观看| 精品人妻午夜| 嫩BBB嗓BBBB榛BBBB| 一级性爱毛片| 一区二区久久| 欧美国产性爱| 丁香五月天婷婷| 国产精品美女久久久久AV爽| 亚洲欧美性爱| 88在线无码精品秘入口九色| 成人自拍偷拍视频| 成人黄网站免费视频| 国产中文在线| 男人V天堂| 日韩人妻精品中文字幕| 老太老熟女城中层露脸60| 一级a片免费看| 四虎影成人精品A片| 影音先锋男人站| 波多野结衣久久精品| 69AV在线观看| 一区二区经典| 一级片免费| 欧美BBWBBWBBWBBWBBwBBW| 亚洲中文字幕第一页| 日本爱爱免费播放视频| 婷婷五月丁香五月| 又黄又湿的视频| 精品第一页| 日本aa视频| 日本女优婷婷青青草原| 中文字幕一二三区| 少妇高潮av久久久久久| 91人妻人人澡人人爽人人精| 青娱乐国产精品| 黄片免费看视频| 先锋成人影音| 三级成人无码| 欧美一级免费视频| 夜夜操影院| 亚洲A∨无码无在线观看| 亚洲天堂久久| 日韩三级片在线视频| 成人在线免费观看国产| 欧美日韩激情视频| 69av天堂| 青娱乐久久| 午夜性爱网址| 天堂网在线视频| 天天色天| 青误乐在线播放| 欧美级毛片一进一出夜本色| 美女久草| 日逼电影网| 一本久久综合亚洲鲁鲁五月天 | 人人操人人摸人人爽| 天天爱天天操| 永久免费无码中文字幕| 成人在线伊人| 丁香六月综合| 国产精品人人人人| 操逼网站免费| 欧美性猛交XXXX乱大交HD | 免费观看A级毛片| 日韩天堂在线| 国产一级a一片成人AV| 成人aV免费观看| 亚洲人妻一区二区| 日本黄色大全| 亚洲精品免费视频| 国产成人午夜高潮毛片| 永井玛丽亚av无码中出流出| 国产一二区| 久久久久久无码精品亚洲日韩麻豆 | 久久久久亚洲AV无码专区| 人人妻人人摸| 毛片久久久| 四季AV一区二区凹凸懂色桃花| 亚洲免费视频播放| 成人精品一区二区无码| 奇米久久| 欧美在线视频一区| 日韩黄色视频在线观看| a天堂视频| 99精品视频国产| 国产欧美高清在线| 色婷婷一级A片AAA毛片| 波多野结衣精品无码| 国产无码一区| 99xav| 中文字幕成人网站中文字幕| 欧美黄色性爱视频| 国产三级片无码| 超碰在线进入| 91超碰在线免费观看| 亚洲黄色一级电影| 中文视频免费播放| 日韩久久视频| 免费看黄色的网站| 亚洲777| 久久99久久99久久99| av电影在线免费观看| 精产国品一区二区区别| 懂色AV一区二区三区国产中文在线| 夜色福利网| 综合成人| 另类Av| 国产A片免费视频| 北条麻妃人妻中文无码| 无码免费在线观看视频| 国产精品18进进出出17c| 九哥草逼网| 99久热在线精品| 日韩天堂av| 中文字幕国产精品| 精品成人AV| 浮力影院久久| 久久性爱网| 骚逼逼影院| 特级西西444WWW高清大视频| 青青草操逼视频| 天天日狠狠操| 亚洲码AV波多野| 国产精品视频播放| 亚洲无码高清一区| 97超碰人人摸| 美女三片| 97成人人妻一区二区三区| a片视频免费| 国产男女av| 欧美去干网| 大香蕉在线75| 三级毛片网站| 国产操逼片| 欧美丰满老熟妇XXXXX性| 亚洲av网站在线观看| 国产十欧洲十美国+亚洲一二三区在线午夜 | 伊人网综合| 福利三区| 国产久久性爱| 亚洲国产97| 久久久久成人片免费观看蜜芽| 日韩一级免费电影| 成人免费一级视频| 91在线视频播放| 岛国免费AV| 五月丁香婷婷啪啪| 性毛片| 欧美成人性爱网| 国产欧美日韩一区二区三区| 亚洲加勒比久久88色综合| 午夜福利高清在线观看| 无码少妇视频| 国产农村妇女精品一二区| 成人三级视频在线| 91香蕉| 国产精品无码AV| 日B视频在线观看| 久草视频资源| 水蜜桃网站| 日本久久久| 在线免费观看视频黄| 国产乱伦视屏| 中文字幕免费视频在线播放| 久久中文字幕视频| 国产一级免费在线观看| 免费V片| 午夜天堂精品久久久| 高清毛片AAAAAAAAA片| 8050午夜网| 男人天堂资源网| 在线观看av网站中文字幕| 99久久99九九99九九九| 国产aa片| 免费国产黄色视频网站| 加勒比日韩在线| 大香蕉伊在线观看| 国产一卡二卡在线| 天天天天天天天操| 久久天堂一区| 久久精品电影| 香蕉毛片| 日韩a在线| 豆花成人视频在线观看| 日韩欧美国产成人| 蜜桃91精品秘入口| 九色影院| 亚洲视频一区二区三区四区娇小视频在线观看视频 | 无码国产99精品久久久久网站| 日都一级A片| 丁香色婷婷| 日韩久久视频| 特级西西444WWW无码视频兔费看 | 国产精品久久久久久久久久久久久久久久 | 中文字幕日本人妻| 亚洲AV永久无码国产精品久久 | 国内无码| 日韩人妻无码一区| 午夜性爱AV| 久艹综合| 欧美V在线| 蜜桃av秘一区二区三区| 亚洲第一影院| 国产欧美综合一区二区三区| 免费无码婬片AAAA片老婦| 亚洲高清无码中字| 国产寡妇亲子伦一区二区三区四区| 国产美女裸体网站| 无码人妻一区二区三区免费n狂飙 性猛交AAAA片免费看蜜桃视频 | 亚洲AV观看| 大香蕉欧美在线| 麻豆AV片| 国产日韩欧美成人| 中文字幕免费看高清| 国产午夜福利视频| 国产1区2区3区| 97精品在线观看| 另类老妇性BBwBBw图片| 少妇人妻一级A毛片| 自拍无码视频| 亚洲色图15p| 一级免费黄片| 国产成人一级片| 激情内射网站| 亚洲一级黄色大片| 无码人妻丰满熟妇区17水蜜桃| 日韩精品一区二区三| 黄色福利视频| 国精产品一区一区三区| 夜夜操网站| 牛牛精品一区| 成人性爱免费网站| 台湾色综合| 伊人三级| 69久久久久久久久久| 亚州成熟少妇视频在线观看| 久久嫩草精品久久久久| 无码不卡在线| 91麻豆免费视频| 无码一区二区高清| 伊人五月丁香| 国产热视频| 色色激情视频| 黄色电影免费在线观看| 欧美日韩性| 狠狠干干| 中文电视剧字幕在线播放免费视频| 四川美人搡BBw搡BBw| 伊人影院在线看| 簧片网站免费| 亚洲AV色香蕉一区二区三区| 91麻豆大奶巨乳一区白虎| 亚洲资源在线| 亚洲精品乱码| 国产一区二区视频在线观看| 超碰人人爽| 亚洲午夜免费视频| 亚洲AV成人网| 色婷婷老师| 人妻爽爽| 在线观看亚洲无码视频| 黑人AV在线播放| 国产黄色AV片| 亚洲秘无码一区二区三区欧美| 麻豆AV96熟妇人妻| 91大鸡巴| 欧美亚洲天堂网| 亚洲成人AV在线观看| 在线免费观看亚洲| 欧美久久久久久| 精品国产女人| 九色PORNY丨自拍蝌蚪| 亚洲无码人妻在线| 91视频美女| 欧美激情xxx| 日逼网站国产| 一级黄色AV片| 超碰自拍私拍二区三区区| 久久三级| 精品国产AV| 亚洲狠狠干| 日韩中文字幕AV| 91人妻无码成人精品一区二区| 日本一区二区视频| 国产免费看片| 国产夫妻在线视频| aaa片| 熟女人妻ThePorn| 99久久精品国产一区二区三区| 亚洲黄色一级电影| 屁屁影院CCYYCOM国产| 亚洲日本无码50p| 美女福利导航| 一级a一级a爰片免费免免中国A片| 中文字幕无码高清| 91国产视频在线播放| 中文字幕在线观看网站| 亚洲国产高清无码| 熟女人妻在线| 国产免费AV片在线无码| AA精品| 日韩肏屄网| 国产色天使| 91熟女丰满原味| 国内精品一区二区| 北条麻妃久久| 欧美黄片一区| 熟女综合网| 超碰一区| 免费视频久久| 成人久久久久一级大黄毛片中国| 高潮流水视频| 国产成人久久精品麻豆二区 | 揉BBB搡BBB搡BBB| 无码欧美人XXXXX日本无码 | 亚洲免费观看高清视频| 国产中文人人国际| 成人免费毛片果冻日本| 淫荡97| 乱子伦一区二区三区视频在线观看 | 操一炮在线视频| 国产麻豆AⅤMDMD0071| MAD033_后宫秘密陶子.| 午夜色色福利| 国产午夜91人妻| 久久伊人春色| 乱子伦一区二区三区视频在线观看 | 无码AV高清| 九九九九九九国产| 99精品视频北条麻妃国产版 | 国产精品无码毛片| 欧美精产国品一| 亚洲女人天堂AV| 日本午夜影院| 日韩AV手机在线观看| 香蕉国产在线| 黄色AV免费观看| 超清无码在线| 中文字幕巨肉乱码中文乱码| 日韩人妻无码电影| 男女啪啪国产| 99黄片| 蜜臀久久99精品久久久兰草影视 | 操操干| 91人人在线| 五十路AV| 一本一道无码| 免费av大全| 国产18| 在线看片a| 成人在线观看AV| 亚洲国产精品成人va在线观看| 人妻丝袜蕾丝高跟双飞| 蜜桃av无码| 7777精品伊人久久7777| 天天干天天干天天操| 骚骚网| 波多野结衣一二三区| 大香蕉96| 亚人精品中文字幕在线观看| 午夜成人无码| 玩弄人妻少妇500系列视频| 久久精品视频18| 亚洲午夜AV久久乱码| 操骚逼视频| 亚洲综合一二三区| 亚洲wwwwww| 91国产免费视频| 11一12周岁女毛片| 三级成人视频| 老太色HD色老太HD-百度| 波多无码在线| 精品国产黄色| 刘玥91精一区二区三区| 无码人妻一区二区一牛影视| 一级欧美一级日韩| 日韩精品高清中文| 青草视频精品| 日韩无码一卡二卡| 欧美MV日韩MV国产网站| 精品黄色电影| 国产操逼免费视频| 久久免费视频播放| 一区二区三区高清| 免费在线观看A片| 日韩成人片| 成人女人18女人毛片| 国产乱子伦-区二区| 国产综合无码| 大香蕉五月丁香| 谁有毛片网址| 五月天开心网| 欧美男人天堂| 北条麻妃在线无码| 国产嫩草久久久一二三久久免费观看| 丁香花在线小说免费阅读| 日韩色图在线观看| 久久精品女同亚洲女同13| 中文在线无码| 亚洲小电影在线| 人妻FrXXeeXXee护士| 国产成人精品一区二区| AV天堂免费播放| 思思热这里只有精品| 黄色三级片视频| 中文字幕北条麻妃在线| 五月色丁香| 免费看一区二区三区| 午夜激情视频在线观看| 日韩无码成人| 特猛特黄AAAAAA片| 北条麻妃成人视频| 影音先锋男人资源网| 99热5| 天堂中文网| 精品国产99久久久久久www| 99艹艹| 欧美三级视频在线| 91麻豆一区| 淫荡人妻视频| 黑人无码在线| 久久AV电影| 日本色情网| 欧美激情视频一区| 久久永久免费视频| 黑人大荫蒂女同互磨| 亚洲va国产va天堂va久久| 欧美精品在线观看| 四川妇搡BBBB搡BBBB| av三级网站| 亚洲网站在线播放| 欧美午夜在线| 日日av| 久久日av| 日本黄色片在线播放| 天天爽夜夜| 亚欧免费视频| 久久精品视频免费| 精品中文在线视频| 国产剧情一区二区| 精品国产天线2024| 日韩黄色无码视频| 中文字幕伊人| 精产国品一区二区| 天天操夜夜操狠狠| 久久无码一区二区| 3344gc在线观看入口| 一级片免费观看| 欧美亚洲综合在线观看| 日批无码| 伊人蕉久| 欧美精品在线播放| 脓肿是什么原因引起的,该怎么治疗 | 99久久伊人| 俺去久久| 成年人毛片| 西西www444无码大胆| 操逼网址大全| 亚洲免费一级| 91伊人网| 精品一区二区三区毛片| 中文字幕一二三四| 免费无码婬片AAAA片老婦| 久久嫩草精品久久久久| 麻豆三级| 成人黄网站在线观看| 欧美黄色性爱| 色五月天婷婷| 免费看无码网站成人A片| 亚洲成人人妻| 免费欧美黄色| 黄色视频网站国产| 91免费在线| 四虎av在线播放| 免费国产在线视频| 久久黑人| 免费操逼网站| 成人性生交片无码免费看人| 免费看操逼逼| 四川少妇BBBB槡BBBB槡| 91国产精品| 无码国精品一区二区免费蜜桃| 坏男人内射老太太| 人妻无码专区| 91一区二区| 一区高清| 无码在线观看免费视频| 免费成人黄片| 国产免费观看av| 翔田千里中文字幕无码| 桃色一区| 国产无码一区| 男人的天堂婷婷| 欧一美一婬一伦一区二区三区自慰,| 91亚洲国产成人精品一区| 天天操天天操免费视频| 色情综合| 奇米99| 免费看操逼视频| 安徽妇搡BBBB搡BBB| 青青色视频| 日韩欧美爱爱| 免费一级婬片AA片观看| 日韩AV自拍| 亚洲AV无码久久精品色无码蜜桃 | 国产人妻人伦精品1国产丝袜| yw视频在线观看| 91色噜噜狠狠色婷婷| 小泽玛利亚一区二区免费| 亚洲综人网| 色天使亚洲| 亚洲狠狠| 99久久人妻无码中文字幕系列| 亚洲成人Av| 欧美四区| 日韩激情视频在线观看| 亚洲乱码中文字幕| www.水蜜桃| 三级片网站国产| 日韩在线成人视频| 午夜在线视频| 天天天天色| AV无码一区二区| 老司机免费视频| 人妻少妇中文字幕久久牛牛| 四川BBB嫩BBBB爽BBBB| av天堂亚洲| 亚洲无码激情在线| 亚洲欲色| 亚洲无码中文字幕在线播放| 蜜臀av一区二区三区| 久草久热| av免费在线播放| 亚洲精品鲁一鲁一区二区三区| 啪啪网站免费观看| 91站街农村熟女露脸| 女女女女女女BBBBBB手| 亲子伦视频一区二区三区| 婷婷五月花| 亚洲播播在线视频| 综合色网站| 欧美日韩三级| 久久免费视频3| 特级西西人体444.444人体聚色| 亚洲猛男操逼欧美国产视频| 国产亚洲精品成人a| 操少妇| 欧美第一网站| 香蕉伊人视频| 久久久性爱视频| 国产一级乱伦| 欧美性爱福利视频| 日韩黄色电影| 丰滿人妻一区二区三| 自拍av在线| 婷婷五月在线播放| 久操青青| 亚洲专区免费| 尹人大香蕉网| 丰满人妻一区二区三区精品高 | 亚洲一线视频| 亚洲无码黄片| 国产女人精品视频| 99er在线观看视频| 无码人妻一区二区三区精品不付款 | 日逼大香蕉| 亚洲无码精品久久| 久久久久久三级电影| 日本成人A片| 无码人妻一区二区三区|