中国少妇裸体淫交_女free性iy9v7交体内谢_女女闺蜜调教羞辱小说_韩国电影《少妇的激情》

學(xué)術(shù)刊物 生活雜志 SCI期刊 投稿指導(dǎo) 期刊服務(wù) 文秘服務(wù) 出版社 登錄/注冊 購物車(0)

首頁 > 精品范文 > 統(tǒng)計學(xué)變量的分類

統(tǒng)計學(xué)變量的分類精品(七篇)

時間:2023-08-10 16:50:29

序論:寫作是一種深度的自我表達(dá)。它要求我們深入探索自己的思想和情感,挖掘那些隱藏在內(nèi)心深處的真相,好投稿為您帶來了七篇統(tǒng)計學(xué)變量的分類范文,愿它們成為您寫作過程中的靈感催化劑,助力您的創(chuàng)作。

統(tǒng)計學(xué)變量的分類

篇(1)

【摘要】 目的 探討累積比數(shù)logit模型在有序分類資料中的正確應(yīng)用。方法 利用在陜西某地開展的一項(xiàng)隊列研究數(shù)據(jù),分別采用累積比數(shù)logit模型和偏比例優(yōu)勢模型進(jìn)行分析,對二者的結(jié)果進(jìn)行比較。結(jié)果 在資料不符合比例優(yōu)勢假定的情況下,直接采用累積比數(shù)logit模型不一定合適。結(jié)論 累積比數(shù)logit模型的應(yīng)用有一定的條件限制,當(dāng)不符合其應(yīng)用條件時,最好采用其他相應(yīng)方法。

【關(guān)鍵詞】 累積比數(shù)logit模型; 比例優(yōu)勢假定; 偏比例優(yōu)勢模型

Abstract: Objective To investigate the correct application of cumulative odds logit model on the ordinal data. Methods The cumulative odds logit model and partial proportional odds model were used for the data from a cohort study in Shaanxi Province, and the results were compared. Results The cumulative odds logit model might not be appropriate when the proportional odds assumption was violated. Conclusion The application of cumulative odds logit model is dependent on certain conditions. It is advisable to adopt other methods when the condition of cumulative odds logit model were violated.

Key words: cumulative odds logit model; proportional odds assumption; partial proportional odds model

累積比數(shù)logit模型是二分類logit模型的擴(kuò)展[1-2],主要用于處理反應(yīng)變量為有序分類變量的資料。該模型對資料要求不嚴(yán),解釋變量既可以是連續(xù)型變量,也可以是無序分類變量或有序分類變量。Amstrong等[3]經(jīng)模擬實(shí)驗(yàn)發(fā)現(xiàn),對于多分類有序反應(yīng)資料,如果采用一般的二分類logit模型而不是累積比數(shù)logit模型,則只能獲得50%~70%的檢驗(yàn)效能。說明在處理有序分類數(shù)據(jù)方面,累積比數(shù)logit模型要優(yōu)于二分類logit模型。但目前對該模型的使用條件、擬合優(yōu)度等內(nèi)容極少有探討。最近幾年已有不少文章對二分類logit模型的正確應(yīng)用進(jìn)行了探討[4-8],而對累積比數(shù)logit模型的研究卻不多見。本文通過實(shí)例對累積比數(shù)logit模型的使用條件和擬合優(yōu)度評價進(jìn)行了探討,為其正確應(yīng)用提供借鑒。

1 材料和方法

1. 1 資料來源 資料來自于陜西某地開展的胃黏膜病變與幽門螺桿菌關(guān)系的隊列研究。為探討幽門螺桿菌與胃黏膜病變進(jìn)展的關(guān)系,對400例淺表性胃炎患者檢測幽門螺桿菌感染狀態(tài),將其分為感染陽性和陰性2組,同時調(diào)查了2組患者的吸煙、飲酒情況,并對所有患者進(jìn)行隨訪。

1.2 胃黏膜病理診斷 用胃內(nèi)鏡觀察胃黏膜病變,并在胃體大小彎、胃角、胃竇大小彎、前后壁各取胃黏膜活檢組織1塊。病理切片按照全國胃、十二指腸活檢、病理診斷標(biāo)準(zhǔn)進(jìn)行診斷。每個受檢對象以最嚴(yán)重病變?yōu)榈谝辉\斷。病變嚴(yán)重程度由輕到重依次為淺表性胃炎(superficial gastritis, SG)、慢性萎縮性胃炎(chronic atrophic gastritis, CAG)、腸上皮化生(intestinal metaplasia, IM)和異型增生(dysplasia, DYS)。

1.3 統(tǒng)計學(xué)處理 本研究中的反應(yīng)變量為病變的不同嚴(yán)重程度,為有序分類資料。分析方法分別采用累積比數(shù)logit模型和偏比例優(yōu)勢模型進(jìn)行擬合。數(shù)據(jù)分析采用SAS 8.1統(tǒng)計軟件。P

2 結(jié) 果

2.1 累積比數(shù)logit模型分析 病例隨訪5年后,共失訪26例,剩余的374例中,135例仍為正?;騍G,139例進(jìn)展成為CAG,59例進(jìn)展為IM,41例進(jìn)展為DYS。多因素分析結(jié)果顯示,校正性別、年齡、吸煙、飲酒等因素后,幽門螺桿菌對胃黏膜進(jìn)展的影響無統(tǒng)計學(xué)意義(χ2=2.7439,P=0.0976)。見表1。

對該資料的比例優(yōu)勢假定條件進(jìn)行檢驗(yàn),結(jié)果發(fā)現(xiàn)該資料不滿足比例優(yōu)勢假定條件(χ2=38.85,P=0.0004),說明該資料采用累積比數(shù)logit模型分析可能并不合適。對模型的擬合優(yōu)度進(jìn)行分析,結(jié)果顯示,Deviance值和Pearson χ2值對應(yīng)的P值均小于0.05,說明模型擬合欠佳。表1 累積比數(shù)logit模型的參數(shù)估計

2.2 偏比例優(yōu)勢模型分析 偏比例優(yōu)勢模型也是基于累積logit的一種方法,但是它放寬了比例優(yōu)勢假定條件[9-11]。采用偏比例優(yōu)勢模型分析結(jié)果顯示(表2),幽門螺桿菌對胃黏膜進(jìn)展的影響有統(tǒng)計學(xué)意義(P=0.0113),這一結(jié)果與目前臨床和流行病學(xué)研究結(jié)果是一致的。而且從參數(shù)估計的標(biāo)準(zhǔn)誤來看,表2中的標(biāo)準(zhǔn)誤均小于表1中的標(biāo)準(zhǔn)誤。表2 偏比例優(yōu)勢模型的參數(shù)估計

3 討 論

本次研究分別采用累積比數(shù)logit模型和偏比例優(yōu)勢模型對同一組數(shù)據(jù)進(jìn)行分析,結(jié)果發(fā)現(xiàn),在數(shù)據(jù)違反比例優(yōu)勢假定條件的情況下,2種模型擬合結(jié)果差別較大。累積比數(shù)logit模型結(jié)果顯示幽門螺桿菌對胃黏膜病變的影響無統(tǒng)計學(xué)意義,而偏比例優(yōu)勢模型則表明幽門螺桿菌的影響有統(tǒng)計學(xué)意義。從統(tǒng)計學(xué)角度和實(shí)際角度綜合考慮,偏比例優(yōu)勢模型的誤差更小,且結(jié)果更為符合實(shí)際,提示本研究數(shù)據(jù)采用累積比數(shù)logit模型不合理。

實(shí)際中應(yīng)用累積比數(shù)logit模型時,主要注意的使用條件就是比例優(yōu)勢假定條件,即自變量的回歸系數(shù)應(yīng)與分割點(diǎn)k無關(guān)[12]。換句話說,無論從哪一點(diǎn)分類,對所有的累積logit,變量χk都有一個相同的βk估計。對于一個自變量χk而言,不同累積比數(shù)發(fā)生比的回歸線相互平行,只是截距參數(shù)有所差別。以前有研究認(rèn)為,累積比數(shù)logit模型對這一條件并不敏感,但本次研究發(fā)現(xiàn),當(dāng)這一基本假設(shè)條件不滿足時,結(jié)果的檢驗(yàn)效率會降低,容易產(chǎn)生假陰性錯誤。

對于累積比數(shù)logit模型的擬合效果,可通過Pearson χ2和Deviance值來判斷。Pearson χ2通過比較模型預(yù)測的和觀測的事件發(fā)生和不發(fā)生的頻數(shù)檢驗(yàn)?zāi)P统闪⒌募僭O(shè),其自由度為自變量不同水平的組合數(shù)目與參數(shù)數(shù)目之差,χ2統(tǒng)計量小就意味著預(yù)測值與觀測值之間沒有顯著性差別,可認(rèn)為模型較好地擬合了數(shù)據(jù)。Deviance通常稱為偏差,它是通過似然函數(shù)來測量所設(shè)模型與飽和模型之間的差異程度。當(dāng)Deviance值較小時,可認(rèn)為所設(shè)模型與飽和模型之間的差異小,即所設(shè)模型擬合較好。這2個指標(biāo)值在多數(shù)情況下是一致的,如果存在不一致,一般認(rèn)為,當(dāng)用最大似然值擬合logit模型時,Deviance值比Pearson χ2值更適用于測量擬合優(yōu)度[13]。需注意的是,當(dāng)模型中含有連續(xù)變量時,Pearson χ2和Deviance指標(biāo)將不再近似于χ2分布,從而也不再適用于擬合優(yōu)度評價。Bender等[14]建議,如果累積比數(shù)logit模型含有連續(xù)變量時,最好采用單獨(dú)的二分類logit模型利用Homser-Lemeshow指標(biāo)進(jìn)行評價。Pulksteins等[15]則提出了一種修正Pearson χ2和Deviance值,用于含連續(xù)變量的累積比數(shù)logit模型的擬合優(yōu)度評價等。

總之,累積比數(shù)logit模型的應(yīng)用不僅僅是參數(shù)求解,而是包含了一系列過程,如變量正確賦值、變量的初步描述、圖示檢查、變量選擇、參數(shù)求解、擬合優(yōu)度評價、結(jié)果的合理解釋等。它需要對專業(yè)和統(tǒng)計學(xué)知識的全面了解,也需要一個實(shí)事求是的態(tài)度。只有準(zhǔn)確掌握應(yīng)用條件,恰當(dāng)?shù)貙δP瓦M(jìn)行評價,并結(jié)合專業(yè)知識進(jìn)行合理的解釋,才能得出可靠的結(jié)論。

參考文獻(xiàn)

[1] Bender R,Gruven U.Ordinal logistic regression in medical research[J].J R Coll Physic Lond,1997,31(5):546-551.

[2] McCullagh P.Regression models for ordinal data [J].J R Stat Soc(B),1980,42(2):109-142.

[3] Amstrong BG,Sloan M.Ordinal Regression models for epidemiologic data[J].Am J Epidemiol,1989,129(1):191-204.

[4] 馮國雙,陳景武,周春蓮.logistic回歸應(yīng)用中容易忽視的幾個問題[J].中華流行病學(xué)雜志,2004:25(6):544-545

[5] 趙宇東,劉嶸,劉延齡,等.多元logistic回歸的共線性分析[J].中國衛(wèi)生統(tǒng)計,2001,17(5):259-261.

[6] 羅登發(fā),余松林.條件logistic回歸模型的殘差分析和影響診斷[J].中國衛(wèi)生統(tǒng)計,1997:14(1):13-15.

[7] 魏朝暉.logistic回歸診斷[J].中國衛(wèi)生統(tǒng)計,2001, 18(2):112-113.

[8] 趙清波,徐勇勇,夏結(jié)來.logistic回歸中高杠桿點(diǎn)的檢測[J].中國衛(wèi)生統(tǒng)計,1997,14(2):17-20.

[9] Bender R,Benner A.Calculating ordinal regression models in SAS and S-Plus[J]. Biom J,2000,42(6):677-699.

[10] Peterson B, Harrell F.Partial proportional ordinal model for ordinal response variables[J].Appl Statist,1990,39(2):205-217.

[11] Lall R,Campbell MJ,Walters SJ,et al.A review of ordinal regression models applied to health-related quality of life assessments[J].Stat Methods Med Res,2002,11(1):49-67.

[12] Fahrmeir L, Tutz G.Multivariate statistical modeling based on generalized linear models [M].2nd ed. Berlin:Springer-Verlag,1998:75-79.

[13] 王濟(jì)川, 郭志剛.Logistic回歸模型——方法與應(yīng)用[M].北京:高等教育出版社,2001:58-65.

篇(2)

【關(guān)鍵詞】高齡老年人;生命質(zhì)量評價;總體健康狀況自評

【中圖分類號】R195 【文獻(xiàn)標(biāo)識碼】A 【文章編號】1004-7484(2012)13-0449-02

1947 年世界衛(wèi)生組織給健康下的定義為: “健康不僅僅是沒有疾病和病痛, 而且還包括身體、心理和社會方面的完好狀態(tài)[1] ”。1990 年WHO提出了“健康老齡化”,1993 年第15 屆國際老年學(xué)會提出了“科學(xué)為健康老齡化服務(wù)”[2]的人口老齡化應(yīng)對目標(biāo)。國內(nèi)外對老年人生命質(zhì)量進(jìn)行了大量研究,取得了較多的研究成果。本研究以SF-36量表為測量工具,在山東省選擇70歲及以上老年人進(jìn)行生命質(zhì)量狀況的調(diào)查,以了解其生命質(zhì)量的總體狀況及影響因素,為提高老年人生活質(zhì)量,促進(jìn)健康老齡化的制度安排提供參考。

1 對象與方法

1.1 研究對象 本次調(diào)查的對象是山東省居家養(yǎng)老和機(jī)構(gòu)養(yǎng)老的70歲及以上老年人。采用分層方法把山東省的老年人劃分為沿海、中部、西部3種地區(qū)類型,每個類型隨機(jī)抽取3個地級市,在被抽取的地級市中以市(縣)為單位進(jìn)行隨機(jī)抽樣。

1.2 調(diào)查方法 調(diào)查問卷包括一般情況表(養(yǎng)老方式、性別、年齡、受教育程度、婚姻狀況、過去從事的職業(yè)、目前的收入來源、生活狀況等)和健康狀況調(diào)查問卷(SF-36),分為居家養(yǎng)老卷與機(jī)構(gòu)養(yǎng)老卷,調(diào)查者是通過專門培訓(xùn)的研究生及本科生采用統(tǒng)一導(dǎo)語入戶對老年人進(jìn)行調(diào)查,能獨(dú)立填寫問卷的由其獨(dú)立填寫,對于沒有閱讀能力或不能填寫問卷的老年人,由調(diào)查人員朗讀問卷內(nèi)容,請其作答,由調(diào)查人員代為填寫問卷。共發(fā)放問卷590份,回收有效問卷558份,回收率94.6%。

1.3 統(tǒng)計學(xué)方法 應(yīng)用Epidata3.1數(shù)據(jù)庫雙機(jī)錄入,SPSS16.0軟件包進(jìn)行數(shù)據(jù)處理,對一般人口學(xué)特征等進(jìn)行統(tǒng)計描述;單因素分析應(yīng)用?2檢驗(yàn),多因素分析應(yīng)用有序分類變量的Logistic回歸分析。

2 結(jié)果與分析

2.1 人口社會學(xué)特征的描述分析

居家養(yǎng)老老人324人,機(jī)構(gòu)養(yǎng)老老人234人,共558人,其中男288人(51.6%),女270人(48.4%)。年齡70~79歲393人(70.4%),80歲及以上的高齡老人165人(29.6%)?;橐鰻顩r:有配偶234人(41.9%),喪偶306人(54.8%),離異4人(0.8%),未婚14人(2.5%)。文化程度:未受過教育269人(48.2%),小學(xué)174人(31.2),初中48人(8.6%),高中、中專及技校36人(6.5%),大專及以上31人(5.5%)。居家養(yǎng)老老人和機(jī)構(gòu)養(yǎng)老老人年齡、性別,受教育程度分布差異無統(tǒng)計學(xué)意義,過去從事的職業(yè)、婚姻狀況、生活狀況差異有統(tǒng)計學(xué)意義(P

2.2總體健康狀況自評及相關(guān)因素差異性分析

本研究將老年人的總體健康狀況分為5個等級,分別是“1=差、2=一般、3=好、4=很好、5=非常好”??傮w健康狀況自評的居家養(yǎng)老老人與機(jī)構(gòu)養(yǎng)老老人差異無統(tǒng)計學(xué)意義(P>0.05)?!靶詣e”在老年人總體健康狀況之間的差異性無統(tǒng)計學(xué)意義(p>0.05),可以認(rèn)為老年人總體健康狀況在性別分布上無差異?!澳挲g”、“受教育程度”、“婚姻狀況”與老年人總體健康狀況之間的差異性有統(tǒng)計學(xué)意義(P

2.3總體健康狀況自評影響因素的Logistic回歸結(jié)果及分析

單因素分析僅對總體健康狀況自評的差異性進(jìn)行檢驗(yàn),而沒有對影響因素進(jìn)行歸因分析,因此需要進(jìn)一步做多因素分析。以總體健康狀況自評為因變量,運(yùn)用SPSS16.0統(tǒng)計軟件中的有序分類變量回歸方法進(jìn)行分析,納入模型中的自變量有年齡、性別、受教育程度、目前婚姻狀況、職業(yè)、養(yǎng)老方式以及生活狀況7個自變量。模型檢驗(yàn)及回歸結(jié)果如下。

2.3.1 模型檢驗(yàn)

表2為對模型中是否所有自變量偏回歸系數(shù)全為0進(jìn)行似然比檢驗(yàn),結(jié)果P

2.3.2 方程中的有效變量及參數(shù)檢驗(yàn)

在納入模型的7個變量中有年齡、養(yǎng)老方式及生活狀況3個自變量對總體健康狀況自評的影響具有統(tǒng)計學(xué)意義(P

3 討論

3.1 年齡與總體健康狀況評價的程度呈負(fù)相關(guān)關(guān)系?!澳挲g”因素在老年人對生命質(zhì)量總體滿意度評價的回歸分析中具有統(tǒng)計學(xué)意義。可以認(rèn)為老人的年齡越大,對生命質(zhì)量的總體滿意度越低,本結(jié)論與景睿、劉曉東等[3]的研究結(jié)果一致。原因在于,年齡越高,身體機(jī)能越低,日?;顒幽芰ο陆?。同時對死亡的恐懼加大,故自我健康的評價較低??梢?,年齡是健康狀況的風(fēng)險因素,不僅表現(xiàn)在日?;顒幽芰Γ彝渡涑鲂睦頎顩r的變化。

3.2 經(jīng)濟(jì)條件對總體健康狀況評價具有積極作用?!吧顮顩r與當(dāng)?shù)匾话慵彝サ谋容^”對70歲及以上老人總體滿意度影響因素的回歸分析中,生活狀況“很富?!薄氨容^富裕”“一般”的老年人,對生命質(zhì)量的總體滿意度更高。原因可能出于以下幾個方面:一是生活狀況水平高,老年人就可能有更大的選擇空間和余地,不為經(jīng)濟(jì)所累,過自己想要的生活,生活比較自由[4]。二是生活狀況較好的老年人有更好的條件享受醫(yī)療保健,更多地參加社交活動,在心理上能夠獲得更多的滿足感,因此對生命質(zhì)量的總體滿意度會較高。

3.3 居家養(yǎng)老老人的健康自評好于機(jī)構(gòu)養(yǎng)老老人。 “養(yǎng)老方式” 因素在老人對生命質(zhì)量總體滿意度評價的的回歸分析中,居家養(yǎng)老模式下的老人對自我生命質(zhì)量評價的總體滿意度高于機(jī)構(gòu)養(yǎng)老老人。首先,老年人觀念相對保守,接受新事物的能力較弱。居家養(yǎng)老是我國傳統(tǒng)的養(yǎng)老方式,在這種環(huán)境中老年人往往具有更高的歸宿感。其次,老年人年齡大,易產(chǎn)生孤獨(dú)感與陌生感。居家養(yǎng)老模式能使老年人更容易獲得子女以及鄰里生活照料、精神慰藉的家庭和社會支持,從而減少孤獨(dú)感與陌生感[5],因此對生命質(zhì)量的總體滿意度高。

3.4 受教育程度與總體健康狀況自評之間的差異性有待于進(jìn)一步研究。在描述性分析的差異性檢驗(yàn)中,受教育程度在總體健康狀況自評的差異有統(tǒng)計學(xué)意義。受教育程度較高的老年人對生命質(zhì)量的總體滿意度高于教育程度較低的老年人。然而在“受教育程度”對總體滿意度自評的影響因素的回歸結(jié)果顯示無統(tǒng)計學(xué)意義(見表3),可能的原因在于混雜因素的影響,此問題有待于進(jìn)一步研究。

參考文獻(xiàn)

[1] Velarde-Jurado E et al.Salud Publica Mex,2002,44:349-361.

[2] 李秀燕, 郭繼志. 老年人生命質(zhì)量評價的現(xiàn)狀及展望[J]. 國外醫(yī)學(xué)?社會醫(yī)學(xué)分冊, 2003,20(4):154-158.

[3] 景睿,劉曉冬等.山東省農(nóng)村老年人生命質(zhì)量評價及影響因素分析[J].中國農(nóng)村衛(wèi)生事業(yè)管理,2008-8 ,8(28)

[4] 林江,楊繼峰,劉強(qiáng). 健康狀態(tài)認(rèn)知理論的概述[R]. 廣西中醫(yī)學(xué)院院報 2010,13(1):74-75

[5] 王生鋒,齊玉梅.中等城市社區(qū)人群生命質(zhì)量評價及影響因素調(diào)查[J].山西醫(yī)藥雜志,2008-10,10(37).

作者簡介:

(1973.03-)碩士研究生,濰坊醫(yī)學(xué)院講師,研究方向:老年社會學(xué)

篇(3)

關(guān)鍵詞 地統(tǒng)計學(xué);生態(tài)學(xué);尺度;時空特征;應(yīng)用

中圖分類號 S153 文獻(xiàn)標(biāo)識碼 A 文章編號 1007-5739(2014)13-0245-01

Application of Geostatistics in Ecology

LI Xiu-mei 1 ZHOU Shi-xue 2 LUO Sheng-jun 2 LI Chang-zhou 2 LIU Li-ping 2 *

(1 Forestry Ecology Department,Hubei Ecology Vocational College,Wuhan Hubei 430200; 2 Macheng Forestry Bureau)

Abstract Natural phenomena has spatial variability and autocorrelation characteristics.As the theory and methods in researching space distribution,goestatistics became an effective method for quantitative analysis of spatial characteristics,and gradually introduced to ecology.Based on the basic theory and method,the paper introduced the application of geostatistics in ecology.

Key words geostatistics;ecology;scale;temporal and spatial characteristics;application

地統(tǒng)計學(xué)起源于20世紀(jì)60年代,是以區(qū)域化變量為理論基礎(chǔ),以半變異函數(shù)和插值分析為主要工具的一種地質(zhì)統(tǒng)計學(xué)方法;是通過分析空間數(shù)據(jù)探索空間過程的信息分析技術(shù)[1]。與傳統(tǒng)生態(tài)學(xué)方法相比,地統(tǒng)計學(xué)變異函數(shù)中的變程(a)包含了距離和方向2種含義,即量化的空間尺度信息,比依靠經(jīng)驗(yàn)直接確定空間取樣尺度更合理;比通過尺度推繹方式間接獲得目標(biāo)現(xiàn)象的特征更精確;修正了傳統(tǒng)景觀格局研究方法中小尺度生態(tài)學(xué)現(xiàn)象在區(qū)域內(nèi)均質(zhì)性的不合理假設(shè)。

由于該理論充分考慮了樣點(diǎn)的位置、方向和彼此間距離等空間結(jié)構(gòu)信息,為實(shí)現(xiàn)參數(shù)的離散化與空間化提供了一種有效工具,可以定量化區(qū)域變量的空間特征,進(jìn)而對未知樣點(diǎn)進(jìn)行無偏最優(yōu)估值,以直接反映自然現(xiàn)象的隨機(jī)性和結(jié)構(gòu)性,廣泛用于地理學(xué)、環(huán)境科學(xué)、土壤學(xué)等諸多研究領(lǐng)域[2-3]。揭示了經(jīng)典統(tǒng)計方法難以發(fā)現(xiàn)的規(guī)律,有利于融合格局、尺度、過程關(guān)系,完善生態(tài)學(xué)理論與方法。

1 地統(tǒng)計學(xué)的基本理論與方法

1.1 區(qū)域化變量理論

地統(tǒng)計學(xué)處理的對象為區(qū)域化變量,區(qū)域化變量的兩大特點(diǎn)是隨機(jī)性和結(jié)構(gòu)性。基于此,地統(tǒng)計學(xué)引入隨機(jī)函數(shù)及其概率分布模型為理論基礎(chǔ),對區(qū)域化變量加以研究[4]。

1.2 變異函數(shù)

變異函數(shù)是地統(tǒng)計學(xué)方法的基礎(chǔ),根據(jù)已知樣本點(diǎn)來確定變量在空間上的變化規(guī)律,推算未知點(diǎn)的屬性。其優(yōu)點(diǎn)在于根據(jù)已知樣本點(diǎn)計算某未知點(diǎn)的屬性值時,考慮了不同距離、不同方向空間點(diǎn)位間的相關(guān)性,使估計值更精確。

1.3 空間插值

與傳統(tǒng)的插值方法不同,地統(tǒng)計學(xué)考慮樣點(diǎn)的方向、位置和彼此間的距離,可以研究既有一定隨機(jī)性又有一定結(jié)構(gòu)性的各種變量的空間分布規(guī)律[5]??死锔瘢↘riging)插值是一種最優(yōu)、無偏的估值方法,在生態(tài)學(xué)中的應(yīng)用最廣泛,可以給出每一估計樣點(diǎn)的不確定性(即產(chǎn)生誤差的幾率和大小),并利用多種附屬信息填補(bǔ)采樣不足的缺陷。這種在誤差允許范圍內(nèi)的空間差值既節(jié)省工作量,又彌補(bǔ)因資料不足帶來的困難[6]。

1.4 空間模擬

空間模擬方法是當(dāng)今地統(tǒng)計領(lǐng)域中由已知推斷未知的最活躍的一種方法。由于Kriging方法具有平滑作用,不適宜用于獲取變量極值的空間分布,例如通過重金屬含量的極值來尋找湖泊水體污染源;而空間模擬方法即能模擬變量的空間變化趨勢又能保留變量的極值;此外,空間模擬通過多種實(shí)現(xiàn)(realization)系統(tǒng)的表現(xiàn),進(jìn)行各種情況下的模擬,具有較好的統(tǒng)計效用[1]。

2 地統(tǒng)計學(xué)在生態(tài)學(xué)中的應(yīng)用

地統(tǒng)計學(xué)檢驗(yàn)、模擬和估計空間特征的作用,對認(rèn)識不同尺度生態(tài)學(xué)功能與過程具有重要意義[7]。20世紀(jì)80年代初,引起生態(tài)學(xué)者關(guān)注,廣泛應(yīng)用于描述生態(tài)因子的空間自相關(guān)性、繪制生態(tài)因子分布圖以及設(shè)計抽樣方案,分析自然因子普遍存在的空間相關(guān)程度、距離和方向等。利用其對可信程度和誤差的評價,解決了定量地測度空間尺度和更精確地繪制自然因子的空間分布[8]。生態(tài)現(xiàn)象所涉及的任何屬性,如植被類型、生物量、土壤化學(xué)元素含量、污染物濃度等,均可以作為地統(tǒng)計模型的變量。

地統(tǒng)計學(xué)方法可以用于研究離散現(xiàn)象的空間特征(群落、種群的格局分布),在處理具有空間連續(xù)性特征的變量(土壤性質(zhì))更占優(yōu)勢。20世紀(jì)80年代初,區(qū)域化變量理論和地統(tǒng)計方法成為量化土壤物理化參數(shù)空間變異的有效方法[9]。尤其是地統(tǒng)計學(xué)方法中半方差圖和Kriging插值法適合于土壤特性的空間預(yù)測[8]。目前廣泛用于土壤養(yǎng)分[10]、水分[11]等的分布、污染物擴(kuò)散、土壤肥力質(zhì)量評價[12]、土壤分類制圖、試驗(yàn)設(shè)計和采樣方法探討等[13]。

地統(tǒng)計學(xué)在生態(tài)學(xué)的應(yīng)用主要集中在生態(tài)系統(tǒng)尺度以下,彌補(bǔ)了景觀生態(tài)學(xué)在處理物種、種群和群落空間信息的不足。區(qū)域化隨機(jī)變量耦合了地形(海拔、坡度、坡向)、土壤、水分等環(huán)境因子,是分析種群、群落的空間異質(zhì)性與環(huán)境因子關(guān)系的有力工具[14];既是對生態(tài)系統(tǒng)尺度空間異質(zhì)性和格局問題的研究,也從生態(tài)系統(tǒng)的組成成分入手,同樣是對生態(tài)系統(tǒng)功能過程的細(xì)化[7]。

3 結(jié)語

地統(tǒng)計學(xué)無論是空間異質(zhì)性分析法還是空間模擬差值法,通常都是對因子的靜態(tài)研究,對處理時間變化問題存在缺陷。有學(xué)者嘗試建立基于地統(tǒng)計學(xué)原理的時間動態(tài)模型[15]。鑒于地統(tǒng)計學(xué)方法側(cè)重于空間分析,而傳統(tǒng)統(tǒng)計學(xué)方法注重于時序分析,二者相結(jié)合,將時間因子融入空間變異中是今后需要解決的問題[6]。非線性多維時間序列分析預(yù)測模型(GS-SVR)就是半變異函數(shù)基礎(chǔ)上,結(jié)合SVR和新的定階方法發(fā)展的一種新的多維時間序列最優(yōu)階數(shù)判斷法[16]。

地統(tǒng)計學(xué)與其他格局分析方法相結(jié)合,將有助于解決現(xiàn)象的時空發(fā)展過程與機(jī)理問題。近期“3S”技術(shù)與地統(tǒng)計學(xué)相結(jié)合,通過地理數(shù)據(jù)確定樣點(diǎn)之間的距離,通過屬性數(shù)據(jù)計算出變量之間的差異,二者結(jié)合得到地統(tǒng)計學(xué)所需要的步長和半方差函數(shù)關(guān)系,強(qiáng)化了大尺度空間信息的分析[10,17]。但目前地統(tǒng)計學(xué)應(yīng)用于景觀尺度的研究較少,并且尚未形成較成熟的定量化研究方法。

4 參考文獻(xiàn)

[1] 王勁峰,武繼磊,孫英君,等.空間信息分析技術(shù)[J].地理研究,2005,24(3):464-472.

[2] 唐濤,蔡慶華,潘文斌.地統(tǒng)計學(xué)在淡水生態(tài)學(xué)中的應(yīng)用[J].湖泊科學(xué),2000,12(3):280-288.

[3] 張坤,洪偉,吳承禎,等.基于地統(tǒng)計學(xué)和GIS的福建省降雨侵蝕力空間格局[J].山地學(xué)報,2009,27(5):538-540.

[4] 孫英君,王勁峰,柏延臣.地統(tǒng)計學(xué)方法進(jìn)展研究[J].地球科學(xué)進(jìn)展,2004,19(2):268-274.

[5] 王景雷,孫景生,張寄陽,等.基于GIS和地統(tǒng)計學(xué)的作物需水量等值線圖[J].農(nóng)業(yè)工程學(xué)報,2004,20(5):51-54.

[6] 王正軍,李典謨,商晗武,等.地質(zhì)統(tǒng)計學(xué)理論與方法及其在昆蟲生態(tài)學(xué)中的應(yīng)用[J].昆蟲知識,2002,39(6):405-411.

[7] 王政,王慶成,李哈濱.紅松老齡林主要樹種的空間異質(zhì)性特征與比較的定量研究[J].植物生態(tài)學(xué)報,2000,24(6):718-723.

[8] 葛劍平,郭海燕,仲莉娜.地統(tǒng)計學(xué)在生態(tài)學(xué)中的應(yīng)用(Ⅰ)―基本理論和方法[J].東北林業(yè)大學(xué)學(xué)報,1995,23(2):88-94.

[9] 廖桂堂,李廷軒,王永東,等.基于GIS和地統(tǒng)計學(xué)的低山茶園土壤肥力質(zhì)量評價[J].生態(tài)學(xué)報,2007,27(5):1978-1986.

[10] 祝錦霞,許紅衛(wèi),王珂,等.基于GIS和地統(tǒng)計學(xué)的低丘紅壤地區(qū)三種土壤性質(zhì)空間變異性研究[J].土壤,2008,40(6):960-965.

[11] 張繼光,陳洪松,蘇以榮.喀斯特山區(qū)洼地表層土壤水分的時空變異[J].生態(tài)學(xué)報,2008,28(12):6334-6343.

[12] 王曉春,韓士杰,鄒春靜,等.長白山岳樺種群格局的地統(tǒng)計學(xué)分析[J].應(yīng)用生態(tài)學(xué)報,2002(7):781-784.

[13] 李艷,史舟,徐建明,等.地統(tǒng)計學(xué)在土壤科學(xué)中的應(yīng)用及展望[J].水土保持學(xué)報,2003,17(1):178-183.

[14] 尚占環(huán),姚愛興,龍瑞軍.干旱地區(qū)山地荒漠草原陰坡植物群落空間異質(zhì)性[J].生態(tài)學(xué)報,2005,25(2):312-318.

[15] 夏敏,趙小敏,張佳寶,等.基于GIS和地統(tǒng)計學(xué)的土壤養(yǎng)分時空變異分析[J].江西農(nóng)業(yè)大學(xué)學(xué)報,2007,29(2):312-317.

篇(4)

關(guān)鍵詞:冠心病心絞痛 蒙醫(yī)證型 診斷標(biāo)準(zhǔn) 臨床研究

心血管疾病已成為影響人類健康的頭號殺手,其中絕大多數(shù)是由冠狀動脈粥樣硬化所致。隨著人們生活水平的迅速提高,本病近年來呈增長趨勢。冠心病心絞痛屬于蒙醫(yī)學(xué)“心刺痛”范疇,目前對冠心病心絞痛蒙醫(yī)證型的診斷尚未完全一致。本研究收集410例包含蒙醫(yī)信息的冠心病心絞痛患者的臨床資料進(jìn)行統(tǒng)計學(xué)分析,并結(jié)合專家意見,初步建立了冠心病心絞痛蒙醫(yī)證型的診斷標(biāo)準(zhǔn)。

一、臨床資料

1、病例入選標(biāo)準(zhǔn)

冠心病心絞痛入選標(biāo)準(zhǔn)依據(jù)這幾年中西醫(yī)結(jié)合冠心病心絞痛會議制定的診斷標(biāo)準(zhǔn)。

2、病例來源

本研究對象為2009年10月-2011年6月期間在本所蒙醫(yī)心病科、中醫(yī)心病科及CCU住院診治的部分患者。

二、研究方法

1、制定病例觀察表

通過回顧文獻(xiàn),確定本研究的病例入選標(biāo)準(zhǔn)和觀察指標(biāo)(包括西醫(yī)和蒙醫(yī)內(nèi)容),并據(jù)此制定統(tǒng)一的病例觀察表。病例觀察表包括與冠心病心絞痛患者蒙醫(yī)診斷相關(guān)信息的指標(biāo)。

2、總體方法

按照西醫(yī)冠心病心絞痛診斷標(biāo)準(zhǔn)確定入選病例并進(jìn)行臨床觀察,填寫病例觀察表,請蒙西醫(yī)結(jié)合專家對每份觀察表信息進(jìn)行辨證并做出相應(yīng)的證型診斷;根據(jù)證型分組,進(jìn)行統(tǒng)計學(xué)分析,篩選對證型診斷有意義的指標(biāo);再請蒙西醫(yī)結(jié)合專家對指標(biāo)進(jìn)行討論,最后制定冠心病心絞痛蒙醫(yī)證型的診斷標(biāo)準(zhǔn)。

3、觀察項(xiàng)目及方法

按照入選標(biāo)準(zhǔn)入選冠心病心絞痛病例,應(yīng)用統(tǒng)一的病例觀察表前瞻性地收集臨床資料。主要觀察指標(biāo)包括:性別、年齡、住院轉(zhuǎn)歸、冠心病心絞痛發(fā)病誘因、既往健康狀況、發(fā)生冠心病心絞痛的持續(xù)時間及發(fā)作次數(shù)、住院天數(shù),入組72 h內(nèi)每12 h及心絞痛發(fā)作時的心電圖、心肌酶[肌酸激酶同工酶(CK-MB)、肌紅蛋白(Mb)、肌鈣蛋白(cTnT)]、血常規(guī)、血生化及血脂[總膽固醇(TC)、三酰甘油(TG)、高密度脂蛋白膽固醇(HDL-C)、低密度脂蛋白膽固醇(LDL-C)]。記錄72 h內(nèi)各指標(biāo)的值。蒙醫(yī)指標(biāo)包括癥狀、舌象、脈象、體態(tài)、根源及因素、疼痛、累及部位等。研究終點(diǎn)為患者在住院期間死亡或放棄治療。

4、統(tǒng)計學(xué)分析

所有觀察資料應(yīng)用SPSS13.0統(tǒng)計軟件進(jìn)行統(tǒng)計分析。對計量資料進(jìn)行兩獨(dú)立樣本t檢驗(yàn),結(jié)果以x±s表示;對分類變量資料采用卡方檢驗(yàn)。

5、質(zhì)控措施

采取參研人員培訓(xùn);研究過程中抽查一定數(shù)量(觀察病例總例數(shù)的10%)病例觀察表與原始住院病歷進(jìn)行核對,復(fù)印一定數(shù)量(觀察病例總例數(shù)的5%)的原始病歷予以存檔等方法控制研究質(zhì)量。

6、數(shù)據(jù)分析方法

進(jìn)行統(tǒng)計分析后,找出各證型中對證型診斷有意義的指標(biāo)。請蒙西醫(yī)結(jié)合專家進(jìn)行討論,確定診斷指標(biāo)。

三、結(jié)果

冠心病心絞痛蒙醫(yī)證型分布經(jīng)專家對部分包含蒙醫(yī)信息的冠心病心絞痛患者資料進(jìn)行分析。

1、粘邪型的指標(biāo)篩選

分析比較粘邪型和非粘邪型的計量資料及兩證型間的分類變量資料,結(jié)果見表1、表2。表1 粘邪型與非粘邪型血脂比較(略) 表2 粘邪型與非粘邪型分類變量資料統(tǒng)計結(jié)果(略)注:癥狀:1=無,2=有;持續(xù)時間:1≥5 min,2≥10 min,3≥15 min,4≥30 min;郝衣希拉體態(tài):1=無,2=有;疼痛程度:1=輕度,2=中度,3=重度; CK-MB、Mb、cTnT:1=正常,2=升高。

2、楚斯型的指標(biāo)篩選

分析比較楚斯型和非楚斯型的計量資料及兩證型間的分類變量指標(biāo),結(jié)果見表3、表4。表3 楚斯型與非楚斯型血脂比較(略)表4 楚斯型與非楚斯型分類變量資料統(tǒng)計結(jié)果(略)注:癥狀:1=無,2=有;持續(xù)時間:1≥5 min,2≥10 min,3≥15 min,4≥30 min;CK-MB、Mb、cTnT:1=正常,2=升高;單純T波改變:1=正常, 2=單純T波改變。

3、胃痧型的指標(biāo)篩選

分析比較胃痧型組和非胃痧型組的計量資料及兩證型間的分類變量資料,結(jié)果見表5、表6。表5 胃痧型與非胃痧型分類變量資料統(tǒng)計結(jié)果(略)注:癥狀:1=無,2=有;CK-MB、Mb、cTnT:1=正常,2=升高;心電圖下壁ST-T改變:1=無,2=有;體態(tài):1=巴達(dá)干赫依,2=非巴達(dá)干赫依表6 胃痧型與非胃痧型血脂比較(略)。

篇(5)

關(guān)鍵詞:獨(dú)立性檢驗(yàn);假設(shè)檢驗(yàn);臨界值

中圖分類號:G632 文獻(xiàn)標(biāo)識碼:A 文章編號:1002-7661(2012)12-178-02

本文將就以下幾個問題展開闡述:

一、獨(dú)立性檢驗(yàn)的形成

獨(dú)立性檢驗(yàn)的基本過程是根據(jù)客觀實(shí)踐情況和經(jīng)驗(yàn),提出原假設(shè),選好統(tǒng)計量,進(jìn)行抽樣、試驗(yàn)、計算、檢驗(yàn),進(jìn)行判斷.也就是說,整個過程貫穿著通過實(shí)踐提出假設(shè)理論,再通過實(shí)踐進(jìn)行檢驗(yàn).假設(shè)的過程其實(shí)就是類似與數(shù)學(xué)證明中的反證法,其基本步驟如下:

假設(shè): :兩分類變量沒有關(guān)系,用A、B表示兩個分類變量,若 成立 事件A與事件B獨(dú)立 .這單純是從概率的角度衡量兩個分類變量的是否有關(guān).我們需要更進(jìn)一步對相關(guān)程度進(jìn)行檢驗(yàn),就是在假設(shè) 下,如果出現(xiàn)一個與 相矛盾的小概率事件,就可以推斷 不成立,且該推斷犯錯誤的概率不超過這個小概率.

二、獨(dú)立性檢驗(yàn)的基本思想

在新課程標(biāo)準(zhǔn)數(shù)學(xué)2-3第三章第二節(jié)對獨(dú)立性檢驗(yàn)進(jìn)行了明確的闡述,課本首先通過對分類變量進(jìn)行定義,分類變量也稱屬性變量或定性變量,它們的不同取值僅表示個體所屬的類別,其取值是離散的.如性別變量,只能取男、女兩個值,商品的等級變量只取一級、二級……,是否吸煙,,國籍等等都是分類變量.分類變量的均值和方差沒有實(shí)際意義,所以不做研究.接著定義列聯(lián)表:一般為兩個或兩個以上分類變量的匯總統(tǒng)計表.在我們的教材中僅限于兩個分類變量的列聯(lián)表,并且每個分類變量只取兩個值,這樣的列聯(lián)表稱為2×2列聯(lián)表,如下:

總計

c

總計

在假設(shè) 成立的條件下,A表示 ,B表示 ,可以通過求分類變量 占總數(shù)與分類變量 占的總數(shù)的概率(用頻率估計概率), 恰好為事件AB發(fā)生的頻數(shù); 和 恰好分別為事件A和B發(fā)生的頻數(shù).由于我們可以利用頻率估計概率,所以在 成立的條件下應(yīng)該有:

, ,可得:

.

即: .

因此, 越小,說明常上網(wǎng)與不及格之間的關(guān)系越弱,否則,關(guān)系越強(qiáng).從這個角度這能說明兩個變量間關(guān)系的強(qiáng)弱,而不能判斷它們具體有多大程度上有關(guān),在此基礎(chǔ)上為了使不同樣本的數(shù)據(jù)有一個統(tǒng)一而又合理的評判標(biāo)準(zhǔn),統(tǒng)計學(xué)家們經(jīng)過研究后構(gòu)造了一個隨機(jī)變量(卡方) = ,并且統(tǒng)計學(xué)家們通過實(shí)踐還得到了如下的卡方臨界值表:

P(K2≥k0) 0.50 0.40 0.25 0.15 0.10

k0 0.455 0.708 1.323 2.072 2.706

P(K2≥k0) 0.05 0.025 0.010 0.005 0.001

k0 3.841 5.024 6.635 7.879 10.828

在表格中的數(shù)字與上述式子能夠說明一個什么問題呢?上面的表格中的第一行是作為檢驗(yàn)的犯錯的上界(上界也是我們要找的

那個小概率),下面的 是取值的臨界值,接下來我們就從一個具體實(shí)例中做以分析:

例 1為了考察高中生的性別與是否喜歡數(shù)學(xué)課程之間的關(guān)系,在我校學(xué)生中隨機(jī)抽取300名學(xué)生,得到如下列聯(lián)表:

喜歡數(shù)學(xué)課程 不喜歡數(shù)學(xué)課程 總計

男 37 85 122

女 35 143 178

總計 72 228 300

由表中數(shù)據(jù)計算 的觀測值.能夠以95%的把握認(rèn)為高中生的性別與是否喜歡數(shù)學(xué)課程之間有關(guān)系嗎?

解:可以有95%以上的把握認(rèn)為“性別與喜歡數(shù)學(xué)課程之間有關(guān)系”.

= 4.514

因此應(yīng)該斷定“性別與喜歡數(shù)學(xué)課程之間有關(guān)系”成立,并且這種判斷結(jié)果出錯的可能性約為5%,所以,約有95%的把握認(rèn)為“性別與喜歡數(shù)學(xué)課程之間有關(guān)系”。

這種利用隨機(jī)變量 來確定是否能以一定的把握認(rèn)為“兩個分類變量之間有關(guān)系”的方法,稱為兩個分類變量的獨(dú)立性檢驗(yàn)。

在上述題目做判斷兩分類變量是否有關(guān)時出現(xiàn)了“斷定”一詞,“斷定”一詞在獨(dú)立性檢驗(yàn)中的含義是指檢驗(yàn)判斷,“斷定為A,B有關(guān)系”就是檢驗(yàn)判斷為A,B有關(guān)系,也就是拒絕A,B無關(guān)系,即拒絕原假設(shè) (接受假設(shè) 的對立面)。“約有95%的把握”中“把握”一詞在獨(dú)立性檢驗(yàn)中的含義是指不犯錯誤的可信度,“有95%的把握”就是有95%可信度(可能性)。換而言之,應(yīng)該是在原假設(shè) 成立的條件下,檢驗(yàn)判斷接受原假設(shè) 犯錯誤的概率不超過5%,而不犯錯誤的概率超過95%。換句話說,就是在原假設(shè) 成立的條件下,不犯錯誤接受對立假設(shè) 錯誤的概率超過95%。對與求出的 的觀測值 越大說明可信度越高,犯錯誤的概率就越小

三、獨(dú)立性檢驗(yàn)的做題步驟

通過以上的分析我們可以知道對于獨(dú)立性檢驗(yàn)問題如何去分析,接下來我們就要從實(shí)際操作中研究怎么去處理這部分問題。首先我們知道從2×2列聯(lián)表的角度來說,我們對列聯(lián)表的中概率的分析可以在直觀上看出它們的概率關(guān)系,而這種直觀判斷不足之處在于不能給出推斷“兩個分類變量有關(guān)系”犯錯概率,但是獨(dú)立性檢驗(yàn)就可以彌補(bǔ)這個不足.即首先直觀上判斷兩分類變量是否有關(guān)系,然后獨(dú)立性檢驗(yàn)主要從是否有關(guān)和有多大的把握認(rèn)為它們有關(guān)這兩個方面來考查,這樣以來就可以比較清晰的看出變量關(guān)系以及相關(guān)程度。那么這時候就需要借助隨機(jī)變量 來求值,進(jìn)而判斷,即要推斷“X與Y有關(guān)系”,可以通過頻率估計概率進(jìn)行直觀判斷,再按下面的步驟進(jìn)行:

1.根據(jù)實(shí)際需要確定容許推斷“兩個分類變量有關(guān)系”犯錯概率的上界 ,然后查表確定臨界值 ;

2.根據(jù)2×2列聯(lián)表與公式計算 的觀測值K;

3.如果 ,就可以推斷“兩個分類變量有關(guān)系”,這種推斷犯錯誤的概率不超過 ;否則就犯錯誤的概率不超過 的前提下不能推斷“兩個分類變量有關(guān)系”,最后做出判斷。

例2在500人身上試驗(yàn)?zāi)撤N血清預(yù)防感冒的作用,把他們一年中的感冒記錄與另外500名未用血清的人的感冒記錄作比較,結(jié)果如表所示。問:該種血清能否起到預(yù)防感冒的作用?

未感冒 感冒 總計

使用血清 258 242 500

未使用血清 216 284 500

總計 474 526 1000

分析:在使用該種血清的人中,有 的人患過感冒;在沒有使用該種血清的人中,有 的人患過感冒,使用過血清的人與沒有使用過血清的人的患病率相差較大。從直觀上來看,使用過血清的人與沒有使用過血清的人的患感冒的可能性存在差異。

解:提出假設(shè) :感冒與是否使用該種血清沒有關(guān)系。由列聯(lián)表中的數(shù)據(jù),求得:

當(dāng) 成立時, 的概率約為0.01,

我們有99%的把握認(rèn)為:該種血清能起到預(yù)防感冒的作用。

篇(6)

關(guān)鍵詞:計量統(tǒng)計模型;向量自回歸模型;經(jīng)濟(jì)增長;脈沖響應(yīng)函數(shù);方差分析

中圖分類號:F224 文獻(xiàn)識別碼:A 文章編號:1001-828X(2016)033-0000-01

隨著計量統(tǒng)計學(xué)與經(jīng)濟(jì)學(xué)知識的相互融合,各種計量統(tǒng)計模型在經(jīng)濟(jì)增長問題中的應(yīng)用越來越廣泛,這里主要通過向量自回歸模型來探討計量統(tǒng)計模型在其中的應(yīng)用。

一、VAR模型在經(jīng)濟(jì)增長問題中的應(yīng)用

首先抽取統(tǒng)計數(shù)據(jù),確定相關(guān)變量,然后建立向量自回歸VAR[1]模型,并用Johansen協(xié)整檢驗(yàn)、脈沖響應(yīng)函數(shù)與方差分解[2]等方法對多個自變量的相關(guān)性進(jìn)行了實(shí)證研究。VAR模型是進(jìn)行多個自變量的分析時最易操作的計量統(tǒng)計模型,經(jīng)常用于預(yù)測相關(guān)時間序列,并分析隨機(jī)擾動項(xiàng)對變量的影響,進(jìn)而說明各種因素對變量造成的影響。此外運(yùn)用VAR 模型的優(yōu)勢在于可合理地描述變量間的互動關(guān)系,而無需事先對變量的外生性和內(nèi)生性進(jìn)行辨別。本文以高等教育的投入和教育規(guī)模與經(jīng)濟(jì)增長之間的關(guān)系為例進(jìn)行說明。

1.自變量選擇、數(shù)據(jù)處理

通常以地區(qū)生產(chǎn)總值(GDP)來代表經(jīng)濟(jì)增長,而收集到的數(shù)據(jù)是絕對數(shù)據(jù),因此會受通貨膨脹等因素的影響。將數(shù)據(jù)換算成某年的不變價格,然后以某年為基期來計算,用以剔除通貨膨脹等因素的影響。為了消除異方差,從而更好的分析某省教育投入(EDU)、教育規(guī)模(SCA)、地區(qū)生產(chǎn)總值(GDP)之間的關(guān)系,對EDU、SCA、GDP取自然對數(shù),分別用InEDU、LnSCA和InGDP[3]表示。

2.平穩(wěn)性檢驗(yàn)

在估計VAR 回歸模型之前,因?yàn)楝F(xiàn)實(shí)中的經(jīng)濟(jì)時間序列一般均為非平穩(wěn)的,采用傳統(tǒng)的計量統(tǒng)計學(xué)知識去建立模型,較容易產(chǎn)生虛假回歸[3]的現(xiàn)象。因此我們需要對各個變量進(jìn)行ADF檢驗(yàn),若各個變量存在單整階數(shù),而且階數(shù)相等時,這樣才可以有效地確定多個變量間的相互關(guān)系。

3.VAR模型建立

運(yùn)用通過平穩(wěn)性檢驗(yàn)的變量建立VAR模型。傳統(tǒng)的計量統(tǒng)計方法不足以說明變量之間的動態(tài)關(guān)系,同時內(nèi)生變量在方程兩端均有出現(xiàn),將會使判斷和估計變得特別復(fù)雜。為了解決這些問題,Christopher Sims[1]于1980年提出了VAR模型,并運(yùn)用到經(jīng)濟(jì)學(xué)中,使經(jīng)濟(jì)系統(tǒng)動態(tài)性分析加速發(fā)展,進(jìn)而解釋各種自變量對經(jīng)濟(jì)增長形成的影響。

VAR(p)模型的數(shù)學(xué)表達(dá)式[1]是:

Yt=a+p1Yt-1+p2Yt-2+…+ppYt-p+Ut

其中,Yt是包含多個內(nèi)生變量的列向量,這里內(nèi)生變量就是指教育投入、教育規(guī)模和經(jīng)濟(jì)增長,p是滯后階數(shù),t是樣本個數(shù)。pj是n行n列的系數(shù)矩陣。Ut是隨機(jī)擾動項(xiàng),j=1,2,…,n。

4.VAR模型的穩(wěn)定性檢驗(yàn)

對上述建立的VAR模型進(jìn)行穩(wěn)定性檢驗(yàn),如果模型中AR特征多項(xiàng)式的全部根的倒數(shù)值都在單位圓內(nèi)[4],說明VAR模型是穩(wěn)定的,則可以說這個VAR模型是有效可靠的,能夠進(jìn)行接下來的脈沖響應(yīng)和方差分解分析。

5.Johansen協(xié)整檢驗(yàn)

Juselius與Johansen [1]共同提出了以VAR模型為基礎(chǔ)的檢驗(yàn)回歸系數(shù)的方法。進(jìn)行多個變量協(xié)整檢驗(yàn)采用這種方法更為合適。此方法是為了說明各變量之間是否存在長期穩(wěn)定的均衡關(guān)系。運(yùn)用Johansen檢驗(yàn)法對各變量進(jìn)行協(xié)整檢驗(yàn),建立標(biāo)準(zhǔn)化的協(xié)整方程,通過系數(shù)看出各變量間的均衡關(guān)系。

6.Granger因果關(guān)系檢驗(yàn)

由Johansen協(xié)整檢驗(yàn)可得,某自變量與經(jīng)濟(jì)增長存在長期的協(xié)整關(guān)系,但這種關(guān)系是否是一種因果關(guān)系,需要進(jìn)行Granger因果關(guān)系檢驗(yàn)。Granger因果關(guān)系檢驗(yàn)[1]是檢驗(yàn)一個變量的滯后變量是否能夠引入到其他變量中。若某個變量受到其他變量的滯后影響,則稱它們具有Granger因果關(guān)系[5]。

7.脈沖響應(yīng)函數(shù)與方差分解

為了更加明確多個自變量與經(jīng)濟(jì)增長之間的動態(tài)關(guān)系,在分析VAR模型時,實(shí)際上是考慮一個脈沖項(xiàng)發(fā)生變化進(jìn)而對整個系統(tǒng)的影響程度,這就是脈沖響應(yīng)函數(shù)分析法。方差分解幫助我們了解每一個結(jié)構(gòu)沖擊對內(nèi)生變量變化的貢獻(xiàn)度[6],用以分析不同結(jié)構(gòu)沖擊的重要程度。

二、小結(jié)

本文以計量統(tǒng)計學(xué)中典型的VAR模型為例,說明了計量統(tǒng)計模型應(yīng)用在經(jīng)濟(jì)增長問題中的具體步驟。在計量統(tǒng)計學(xué)中,還有很多的統(tǒng)計方法適用于經(jīng)濟(jì)問題,比如主成分分析法、層次分析法等等,大家都可以此為例進(jìn)行展開學(xué)習(xí)與應(yīng)用。運(yùn)用計量統(tǒng)計學(xué)分析各因素如何影響經(jīng)濟(jì)增長,為制定經(jīng)濟(jì)政策的參與者提供好的建議,提升計量統(tǒng)計學(xué)的現(xiàn)實(shí)意義。

參考文獻(xiàn):

[1]王少平.計量經(jīng)濟(jì)學(xué)[M]高等教育出版社,2011:1-992.

[2]張曉峒.Eviews使用指南與案例[M].北京:機(jī)械工業(yè)出版社.2007:84-87,121-122

[3]鄧媛,李瑞光.VAR模型實(shí)證分析云南省教育投入與經(jīng)濟(jì)增長的關(guān)系[J].2009(2):1-6.

[4]馬鵬媛,米紅.高等教育規(guī)模與經(jīng)濟(jì)增長關(guān)系演變的實(shí)證研究[J].2012(2):1-5.

[5]蔡文伯,程璐.新疆教育增L與高等教育規(guī)模VAR模型[J].2012(12):1-3.

篇(7)

中圖分類號:B844.2 文獻(xiàn)標(biāo)識碼:A 文章編號:1000-6729(2008)008-0628-02

焦慮是青少年最常見的心理問題之一,其發(fā)病率甚至高于成人[1],但對青少年群體不同亞組的焦慮分布研究少見。本文以高中學(xué)生為樣本,調(diào)查其焦慮水平與一些人口學(xué)資料之間的關(guān)系。

1 對象與方法

1.1 對象

采用分層隨機(jī)整群抽樣方法,調(diào)查長沙和岳陽各一所中學(xué)(城市)、瀏陽一所中學(xué)(農(nóng)村)的高一至高三年級學(xué)生。發(fā)放問卷1612份,收回有效問卷1561份,其中長沙中學(xué)558人,岳陽中學(xué)41人,瀏陽中學(xué)592人;男生755人,女 生779 人,缺失27人;獨(dú)生子女 726 人,非獨(dú)生子女 791 人,缺失44人;父親受教育程度:文盲60人,小學(xué)359人,初中371人,高中516人,大學(xué)或以上139人,缺失116人;母親受教育程度:文盲89人,小學(xué)410人,初中424人,高中458人,大學(xué)或以上66人,缺失114人;家庭月收入:< 1000 元249 人, 1000-1500元386 人,1500-2000元294 人,2000-3000元245 人,3000-4000元119 人,4000-5000元46 人,> 5000元42 人,缺失180人。

1.2 工具

兒童多維焦慮量表中文版(The Multidimensional Anxiety Scale for Children-Chinese Version,MASC-C)[2]為自評量表。共包括39個條目,分為4個因子:軀體癥狀、躲避傷害、社交焦慮、分離焦慮/恐怖。要求受試者說明最近一周內(nèi)焦慮癥狀或感覺出現(xiàn)的頻度,采用0-3級評分。

1.3統(tǒng)計方法

進(jìn)行獨(dú)立樣本t檢驗(yàn)、單因素協(xié)方差分析、Logistic 回歸分析。

2 結(jié) 果

2.1 不同人口學(xué)特征學(xué)生兒童多維焦慮量表中文版的評分比較

農(nóng)村組的焦慮癥狀總分、軀體癥狀及躲避傷害因子分均高于城市組(47.1±15.9/44.7±16.5、11.0±6.5/9.5±6.7、15.3±3.9/14.8±4.3,t=2.819、4.091、2.243,P=0.005、0.000、0.025)。女性組的焦慮癥狀總分、軀體癥狀、躲避傷害、社交焦慮及分離焦慮因子分均高于男性組(49.9±16.1/41.5±15.4、11.3±6.9/8.8±6.2、15.5±4.1/14.4±4.2、15.0±5.2/13.4±5.5、8.1±4.4/4.9±3.8,t=10.470、7.315、5.257、5.727、15.272,均P

父親受教育水平較低的高中生的焦慮癥狀總分及4個因子分均較高,母親受教育水平較低的高中生的焦慮癥狀總分及軀體癥狀、社交焦慮、分離焦慮因子分均較高(見表1);按家庭月收入分組的各組被試在焦慮癥狀總分及各因子分上的差異均未達(dá)到統(tǒng)計學(xué)顯著水平。單因素協(xié)方差分析結(jié)果顯示:控制其他人口學(xué)變量后,不同父親教育程度組別之間的焦慮癥狀總分及各因子分差異均未達(dá)到統(tǒng)計學(xué)顯著水平,不同母親教育程度組別之間在社交焦慮因子分上仍存在顯著差異(F=2.942,P=0.019)。

2.2 焦慮影響因素的多分類有序變量Logistic回歸分析

將被試的MAS-C焦慮癥狀總分為因變量,按百分等級分為3組:57分為3組(百分等級75以上)。將居住地、性別、同胞狀況、父母受教育水平引入多分類有序反應(yīng)變量Logistic 逐步回歸模型,最終僅性別進(jìn)入Logistic 回歸方程(OR 值為0.436,P

3 討 論

本研究結(jié)果顯示女性是青少年產(chǎn)生焦慮癥狀的危險因素,與Twenge JM、王芳芳等的研究一致[3-4]。居住地為農(nóng)村、非獨(dú)生子女以及父母受教育程度低的高中生,其焦慮水平較高,但進(jìn)行協(xié)方差分析后上述差異就消失了。提示居住地、同胞狀況及社會經(jīng)濟(jì)狀況等變量并不是獨(dú)立對焦慮水平起作用的,可能相互間存在作用的重疊。

參考文獻(xiàn)

[1] 劉賢臣,孫良民,唐茂芹,等. 2462名青少年焦慮自評量表測查結(jié)果分析. 中國心理衛(wèi)生雜志,1997,11(2): 75-77.

[2] 鄒濤, 姚樹橋, 朱熊兆,等. MASC 中文版信效度的分析及與西方樣本的比較. 中國臨床心理學(xué)雜志,2007, 15(5): 452-455.

[3] Twenge JM, Nolen-Hoeksema S.Age, gender, race, socioecono-mic status, and birth cohort diffe-rences on the Children's Depression Inventory: A meta-analysis. J Abnorm Psychol,2002, 111 (4): 578-588.

    欧美日韩一区二区三区自拍| 久久综合亚洲鲁鲁五月天69堂| 欧美巨大丰满少妇xxxx人| 伊人久久艹| 性少妇videosexfr老师| brazzers精品女教师| 波多野结衣不打码视频| 91精品国产乱码久久久久久蜜臀 | 精品亚洲欧美无人区乱码| 色偷偷av一区二区| 野外震车做爰视频播放| 免费看日本zzzwww色| 日本r级未删减在线观看| 试衣间做爰过程| 正在播放淫亚洲| 国产成a人片在线观看视频下载| 性片免费看| 肉体粗喘娇吟国产91| 娇妻出轨放荡呻吟| 久久久久69| 国产片侵犯亲女| 成年人黄色片| free性欧美| 秘密基地电影免费观看| av一区久久| 狠狠色视频| 黄色软件成人在线| 亚洲a∨精品一区二区三区| 稚嫩的玉腿在稚嫩h| 日韩三级麻豆| 他的硕大进出粉嫩h校园| 国产乱国产乱老熟300部网站| 亚洲无极品| 综合精品欧美日韩国产在线| xxxx性开放按摩hd| 久久久久欧美精品| 国产自产c区| 亚洲精品乱码久久久久久v| 亚洲精品久久久久久久久久吃药| 性猛交乱淫av| 男人桶进女人屁股| 韩国免费黄色| 啪啪床吻戏| 噜噜噜av| 苍井空电影在线观看| 日本女人毛耸耸| 白丝校花让我c| 夜夜躁很很躁日日躁麻豆| 国产真实伦在线观看| 亚洲一线二线三线写真| av集中营| 我的错误电影在线观看| fre性欧美69巨大| 国产精品免费看久久久8精臀av| 娇妻被交换粗又大又硬产品视频| 久久久久久久久福利| 日本特黄特色aaa大片免费| china女人ass裸体欣赏| 欧美高大丰满少妇xxxx| 19禁大尺度做爰| 中国av在线播放| 女主人享受舌奴用嘴伺候| 男男重囗味sm群虐| 99精品国产丝袜在线拍国语| 欧美精产国品一二三产品特点| 老师好紧开裆蕾丝内| 欧美人狂配大交3d怪物一区| 美女和帅哥鸡鸡对鸡鸡| 韩国免费黄色片| 久久香蕉精品| 日本最新伦理片| 校花大胸奶汁乳流奶水出来| 2区3区在线涩网涩| 免费av网址在线观看| 精跪趴灌满h室友4p动漫| 丝袜ol办公室秘书系列番号| 韩国吻胸激情视频| 泰国做爰视频| 久久午夜神器| rylskyart极品少妇人体印度| 亚洲人成网站色www| 亚洲乱码国产乱码精品精软件| 明星后宫艳荡史| 公交车猛烈进出婷婷2的视频| 中国少妇和黑人老外做爰疼痛感受| a级片日本| 老翁h狠狠躁死你h| 妹妹我要操| 成年人免费观看黄色| 亚洲欧美日韩综合一区二区| 欧美色图片区| 男人和女人日比视频| 精品久久久久久无码中文字幕| 亚洲另类春色综合婷婷| 国产精品久久久久久福利| 久久久久激情电影| 双性荡夫冰山美人h| 婷婷精品国产一区二区三区日韩| 教练你好大轻点练车网站| 久久综合伊人中文字幕| 精品精品欲天堂导航| 啊轻点灬大巴太粗太长了小熊| 打扑克一边摸一边叫床一边爽| 男生用鸡鸡插女生| 亚洲欧美另类色图| 丰满老女人高潮呻吟| 99精品无人区乱码在线观看| 双性h被强迫各种姿势侵犯| ass白嫩少妇ass巨大戾| 和护士长疯狂作爱| 免费黄色亚洲| 樱桃视频免费观看视频大全| 国产精品国产亚洲精品看不卡| www.日韩毛片| 小sao货撅起屁股扒开c打视频| 趴下我要c你校园h| 黄色裸视频| 国产精品igao视频网网址| 欧美影院二区| 日韩av专区| 亚洲成aⅴ人在线观看| 九色porny丨国产精品| 双性男仆胸大又放荡的h| 天天躁夜夜躁天干天干2022| chinamodel私拍大尺度| 伊人精品久久| 亚洲一区精品二人人爽久久| beegxxxxxxxhd日本老师| 欧美激情偷拍| 三级午夜理伦三级交换的一天| 出轨闺蜜交换高h| 亚洲欧美乱综合图片区| 精品偷拍一区二区三区在线看| 国产一在线精品一区在线观看| 久久99亚洲精品久久99| 狠狠干毛片| 女帝在龙椅上受辱沦陷h| 人乳喂奶hd无中字| 日本人性生活视频| 欧美极品少妇感bbbbbb| 日本aaaa级毛片| 欧美做爰猛烈床戏大尺度电影视频| 色爽一av| 91九色在线视频| 尤物视频在线观看| 美丽姑娘高清版在线电影| 最新国产精品自在线观看| 哦┅┅快┅┅用力啊┅美艳后宫| 濑亚美莉种子| 天堂躁躁人人躁婷婷视频ⅴ| 亚洲天堂激情| 国产夫妇交换群4p大战| 亚洲精品国偷自产在线91正片| 绿帽淫传h| 偷拍大众浴池精品视频在线播放| 亚洲精品国产一区二区| 奴性少妇下跪被踩头女虐女 | 老女人裸体视频| 夫妻日b视频| 国产农村妇女精品一二区| 艳妇乳肉潘金莲1—5在线看 | 日本夜爽爽一区二区三区| 久久99精品久久久久国产越南| 又色又爽又黄又免费的照片| 扒开腿挺进湿润的花苞电影| 香港黄色一级片| 女生下面是什么味道| 隣の若妻さん波多野结衣 | 天天躁日日躁狠狠的躁天龙影院| 是用玉势不是用玉柱| av漫画在线观看| 欧美国产亚洲日韩在线一区 | 福利cosplayh裸体の视频| 中国黄色三级| 男女生搞逼| 国语精品对白露脸少妇网站| 女性爱爱视频| 久色视频网站| 免费一二二区视频| 免费无遮挡男女交性视频| av片中文| 少妇亲子伦av| 亚洲偷怕自拍| 青青草97国产精品免费观看| 国产精品免费一区二区三区四区| 秋霞一区二区三区| 一级色网站| 男女裸体无遮挡激情床吻视频 | 亚洲国产精品久久亚洲精品| 亚洲精品永久视频| 啊轻点灬太粗嗯太深了用力的视频| 国产a三级| 娇妻被恩人征服h版| 国产精品久久久久久爽爽爽| 边吃奶边添下面好爽| 远古野人粗壮h灌满3p| 男女日bb视频| 一区二区偷拍| 欧美激情福利| 日韩视频黄色| 硅胶人做爰全过程| 看很黄很黄的细节小说| 裸体精品bbbbbbbbb美国| 揉捏她饱满的双乳呻吟视频| 欧美另类高潮| 国产精品一区二区久久精品| 成年片人尤蜜免费视频| 国产亚洲成aⅴ人片在线观看麻豆| 娇喘呻吟柳腰雪臀紧湿| 壮攻被下药成骚受| 广东少妇aaaaaa毛片| 男操女人逼| 亚洲国产午夜精品理论片在线播放| 把极品白丝班长啪到腿软| 嗯~啊~快点死我男男视频| 亚洲欧美日韩人成在线播放| 麻豆91av| 日本做爰猛烈床戏大尺度电影网盘| 久久久久国产a免费观看rela| 特黄一级视频| 91视频黄色软件| 女明星三级做爰无极| 久久综合久| av电影在线观看| 婷婷综合久久中文字幕蜜桃三电影| 国模欢欢大尺度| 日本黄色电影网站| 欧美亚洲另类激情另类| 操你妹影视| 亚洲天堂色2017| 国产男女免费视频| 久久99精品国产麻豆婷婷洗澡| 久久综合久久美利坚合众国| 欧美最顶级丰满的aⅴ艳星| 91视频下载污| 伊人久久av| 男女无遮挡动态图| 白丝班长扒开腿让我爽了一夜| 一极黄色片| 麻豆一区二区三区精品视频| 日产精品久久久一区二区| 日韩伦理毛片| 亚洲伊人天堂| 操操网av| 黄页网站视频| 3d小舞被黄动漫免费视频| 一个色综合av| 免费做爰姿势视频| av在线操| 人操人人爽| 日本免费人成在线观看网站| 欧美偷拍色图| 完美世界在线观看免费播放| 好姑娘4影视剧在线观看| 萍萍的性荡生活第二部| 丰满的护士hd高清| 99av国产精品欲麻豆| 麻麻张开腿让我桶| 色一区二区三区四区| 操孕妇的逼| 国产真人无遮挡作爱免费视频| 黑人干日本少妇| 丰满放荡岳乱妇91www| 黄色在线小视频| 亚洲小说欧美激情另类| a级情欲片在线观看| 国产一区二区精品久久岳| 刘亦菲高潮水多又紧又爽| 亚洲国产精品成人久久久| 男人添女人囗交做爰曰本| 天天操 夜夜操| 老校长胯下的呻吟h| 国产人久久人人人人爽按摩 | 制服.丝袜.亚洲.另类.中文| 免费看男人躁女人过性生活视 | 国产在线精品一区二区不卡麻豆 | 情趣蕾丝内衣少妇啪啪av| 中文字幕精品一区二区2021年| 国产做爰又粗又大又深人物| 国产精品美女久久久网av| 亚洲国产精品第一区二区| 扒开美女狂揉视频在酒店| 国产乱人激情h在线观看| 好姑娘4影视剧在线观看| 亚洲一区天堂| 少妇被黑人4p嗷嗷叫| 亚洲劲爆av| 亚洲精品久久久中文| 葫芦里面不卖药千片万片你需要| 欧美精品老牛影视| 日日夜夜精品视频免费观看| 日本电影少妇| 黄色录像片子| 91老师片黄在线观看| av片在线观看免费| 欧美亚洲日韩国产区三| 巨ru麻麻被调教成肉奴小说| 黄色片网站大全| 日本护士做爰取精视频| 男女插插插| 别揉我胸啊嗯~喷水了日本| 久久99精品久久久久久水蜜桃| 亚洲欧美久久久| 日本一级裸体洗澡| 解开老师裙子猛烈进入| 亚洲精品99久久久久中文字幕| 91香蕉视频污污| 边走边做h挺进闺蜜娇妻季红豆| 8拔插拔插x8免费视频x_拔插拔插更快乐| av网站免费线看| 摸着她得双乳好兴奋| ass中国艳妇裸体欣赏| 亚洲色图网友自拍| 久久久综合网| 亚洲人成777| 我爱av好色| 污黄视频在线免费观看| 性猛交ⅹ×××乱大交| 农村少妇bbwbbw精品hd| 欧美一级淫片丝袜脚交| 呻吟耸动撞击女同学| 青青影院av| 91av麻豆| 国产乱码av| 亚洲精品无播放器在线播放| 久久久日韩精品一区二区| 天天色天天操天天射| 免费三级伦理片| 久久久精品欧美一区二区免费| 黄色免费影片| 激情小说亚洲图片伦| 国产性色强伦免费视频| 日韩av中文字幕第一页| 和女房东啪啪呻吟不断| 国内精品久久久久久影院| 久久日韩乱码一二三四区别| 丰满少妇毛片| 每天裸体跪着被主人调教| 久久精品久久久久观看99水蜜桃| 少妇自拍一区| 是偶然吗电视剧全集免费播放 | 国产999精品久久久久久| 国产欧美日韩| 天天天天躁天天爱天天碰2018| 亚洲人成亚洲精品| 第1章双乳喷奶水| 日韩伦理免费看| 日韩夜夜爽| 日韩av免费片| 甘雨到爽高潮痉挛动漫| 纯爱无遮挡h肉动漫在线播放| 国产福利视频一区二区| 性感美女视频一二三| 日韩性色视频| 麻豆伦理电影| 免费国产黄线在线播放| 性福利视频| 国产乱淫av片免费看| 乱lun乱lun合集| 国内精品国产三级国产a久久| free性护士vidos猛烈主演| 跨坐湿润晃动h小三| 五月婷婷一区二区| 看电影来5566一区二区黑夜| 久久精品国产久精国产一老狼| 一区二区三区免费看视频| 国产激情婷婷| 特级淫片裸体免费看视频| 做受高潮aaaa动漫| 挺进美妇玉蚌深处名器田野的春天| 久久久黄色免费网站| 黄色二级网站| 久久久久久亚洲精品中文字幕| 慈溪职高吃乳门图无遮挡| 天堂av网手机版| 黄色av网址大全| 美女黄色网页| 波多野结衣福利视频| 娇妻被朋友在客厅呻吟| 很详细的肉肉床文过程片段视频| 给处女破处的过程| 国产伦一区二区三区色一情| 国产麻豆精品一区二区三区v视界| 日韩精品色哟哟| 黄色av入口| 主人~别揉了~尿了~啊哈网站| 欧美三级香港三级日本三级| 伊人av超碰久久久麻豆| 少妇被粗大猛进进出出s视频| 亚洲日韩在线中文字幕综合| 亚洲国产精品成人久久| 日本色黄视频| 欧美乱色视频| 美女被吸乳得到大胸app| 国产精品线在线精品| 放荡的女邻居4伦理| 韩国三级毛片| 国产精品日日摸夜夜摸av| 优衣库试衣间未删减视频| 国产精品日日摸夜夜添夜夜av| 日本少妇交换做爰| 少妇一线天| 国产粉嫩一区二区三区| 亚洲综合激情五月久久| 折磨小男生性器羞耻的故事| 天堂а在线地址8最新版| 国产精品福利视频萌白酱app| 亚洲欧美乱日韩乱国产| 日本护士囗交片| 免费黄av| 正在播放白嫩极品美乳视频动漫| 在线观看伦理片| jzzijzzij日本成熟少妇| 欧美狠狠操| 日本少妇与黑人| 一级做a爰片性色毛片16| 百合h啪肉np文〈肉| 黑人一级淫片40厘米| 伊人av影院| 精品一区二区三区在线视频| 厨房娇妻被朋友胯下挺进视频| 少妇被黑人操| 激情综合五月丁香亚洲| 国产男女裸体做爰爽爽| 曰韩黄色片| 黄色日漫视频| 小少妇av| 同桌把中指伸进了我的下面h| 亚洲精品久久激情国产片| 久久久二区| 美艳老师的性福生活| 爽好大快水蜜桃| 揉我胸啊嗯~出水了国产精品| 挺进高潮颤抖h| 美少年高潮h跪趴扩张调教喷水| 欧洲精品免费一区二区三区| 青草精品国产福利在线视频 | 丰满少妇被猛烈进入高清视频| 波多野结衣福利视频| 小的p好紧好紧小说| 麻豆影院免费夜夜爽日日澡| 快穿被各种男人啪h男男视频| 影音先锋亚洲一区| 日日摸日日操| 999久久久国产精品| 欧美日韩精品| 做爰丰满少妇内衣1313图片| 最新日韩av| 美女搞黄在线观看| 诱人的乳峰奶水2| ass少妇猛性bbwbbw| 浑圆乳房大力地揉搓着h| 任你躁久久久久久妇女av| 潘金莲裸体高潮| av黄色图片| av免费资源在线观看| 日韩色黄大片| 国产亚洲精品久久久久久久软件| 欧美性调教| 懂色av蜜臀av粉嫩avdnlt| 美女黄网站| 天堂国产一区二区三区| 女人张开腿让男人桶爽30分钟| 日韩两性视频| 日产精品久久久一区二区 | 男人亲女人下面的视频| 激乳放荡的护士2| 欧美三级韩国三级日本三斤| 夹住小核h颤抖| 校花上面好大下面好湿| 调教少妇小说| 亚洲精品日日夜夜| 巨熟乳波霸av田中瞳| 国产精品成久久久久三级6二k| 日本亲子乱子伦xxxx| 一区二区三区四区av| 日本少妇丰满做爰图片| 天干夜天干天天天爽视频| 国产69久久精品成人看| 淫辱的秘书调教15| 亚洲超丰满肉感bbw| 亚洲无人区码一码二码三码| 欧美一区二区网站| 日韩污视频网站| 综影视肉文| 亚洲国产欧美日韩精品一区二区三区 | 国产喷白浆精品一区二区豆腐 | bwbwbwbwbw大高潮| 男人呻吟双腿大开男男h互攻视频| 免费成人黄色网| 国产在线码观看清码视频| 国语精品对白露脸少妇网站| 亚洲欧洲日韩| 91黄瓜视频| 和女邻居做爰伦理| 91欧美日韩麻豆精品| 成人区精品一区二区婷婷| 日韩av资源| 停不了的爱在线观看| 国产女人高潮合集特写| 国产白嫩美女在线观看| 亚洲另类制服| ass人间极品裸体ass| 免费黄色视屏| 韩国三级中文字幕hd爱的色放 | 日韩一区av在线| jizz精品| 91黄瓜视频| 日韩av卡一卡二| 乱肉妇精品av| 校霸双腿打开揉弄高潮h视频| 叶玉卿三级露全乳视频| 糈精国产xxxx在线观看| 很详细的肉肉床文过程片段视频| 青青影院av| 国产电影一区二区三区爱妃记| 全篇肉高h秘书被c办公室有| 免费看裸体网站黑桃| 亚洲女人久久久| 男生jj插女生jj| 美女被男生猛戳视频| 年轻的私人女教师5| 久久综合久久久久| 关晓彤下面好紧进不去h| 美女屁股无遮挡图片| 久久久久久av| free性欧美69巨大| 穿乳环蒂环露出高h| 黄色91大片| 抱着娇妻让粗大玩3p| 香港黄a三级三级三级看三级 | 稚嫩校花高耸双乳h| 美女喷白浆| 久久午夜电影网| 人妖一区二区三区| 武则天黄色片| 少妇性色av| 亚洲综合精品香蕉久久网| 少妇无内裤下蹲露大唇视频| 多人一起做人爱视频| 啊轻点灬大巴太粗太长了视频| 全黄h全肉禁乱公交车| 亚洲精品女人| 尻逼免费视频| 亚洲黄色av电影| 国产又粗又长又黄又猛小说推荐 | 两个男生做酿酿酱酱的视频| 黄色精品久久| 欧美黑人一级视频| 中文字幕久一区二区| 国产乱人偷精品视频免费| 亚洲无线av| 无码中文字幕日韩专区| 适合夫妻晚上看的电视剧大全| 亚洲精品久久久久久一区二区| 男女互操视频| 国产女烈大尺度sm受刑视频| 国产a三级| 欧美二区乱c黑人| 小荡货腿张开让我cao视频| 精品久久久久久国产| 小受被用各种姿势进入np| ysl口红水蜜桃色号1258| 久久伊人精品| 自拍偷拍视频二区| 嫩模被强到高潮不断视频| 浴室嗯啊啪呻吟h| 亚洲视频福利| av网站国产| 最新国产精品自在线观看| 黑人壮汉把娇妻曰到高潮| 免费av在线网站| 特级毛片在线| 国产绿奴调教91| 日韩欧美a级片| 办公室双腿打开揉弄高潮淑芬| 放荡娇妻肉交换h文章| 久久爱一区二区三区电影| 久久久黄色大片| 成在线人免费视频| 久久免费福利视频| 把睡裙往上撩嗯啊h| 国产精品爱久久久久久久电影蜜臀 | 《头等舱》电影大尺度| 女s调教vk| 少妇装睡让我滑了进去| 水蜜桃涩涩| 男生和女生一起差差差app| 亚洲欧美日韩精品色xxx| 91pony九色丨交换| 亚洲成在人线在线播放| 亚洲一区二区久久| 亚洲欧洲日产国码久在线| 岳把我用嘴含进满足我视频| julia乱码中文一二三区| 无遮无挡爽爽免费视频软件| 免播放器看av| 少妇做爰免费直播| 三级毛片av| 肌肉猛男的巨大粗爽gⅴ猛男自| 日韩一区二区在线观看视频| sm男男调教〈h丿海棠| 我和岳交换夫妇爽4p晓娟小说 | 性爱电影免费看| 欧美另类高清zo欧美| 在线看不卡的av| 日本va欧美va精品发布| 亚洲厕所偷拍| 扒开双腿猛烈进入床戏gif| 邻居少妇张开双腿让我爽一夜| 亚洲成年人电影| 亚洲第二区| 国产精品igao为爱做激情| 老师露出强行让男生揉小说| 精品爱爱视频| 日韩精品久久久肉伦网站| 国产精品视频免费一区二区| 国模大尺度视频| bbbbb女女女女女bbbb| sm人妖前列腺高潮视频| 操日本女人逼视频| 妻欲迷城h未删减| 午夜激情爱爱| 成人欧美日韩一区二区三区| 巴西少妇xxb大毛又多| 丰满放荡岳乱妇91ww| 亚洲一区二区丝袜| 精品偷拍被偷拍在线观看| 国产99久久精品一区二区| 黄色录象片| 《制服肉奴隷》在线| 拍戏床戏被肉高h文| 法国艳星性猛交| 狠狠操夜夜爱| 被cao哭高h奶水h| 少妇口述高潮| 在公车被弄得高潮呻吟| 宿舍里的贱奴女虐女被荡秋千玩| 午夜久久久久| 浓精灌孕h校园h乱小视频| 男男稚嫩小受哭叫粗大撑破h动漫 日本三级日本三级日本三级极 | 新素女艳谭| 绿帽h啪肉np辣文| 男插女青青影院| 小雪被老外黑人撑破了视频| 秀婷程仪公欲息肉婷在线观看| 欧美人狂配大交3d怪物一区| 国产精品第一国产精品| 国产精品自在线拍国产| 成人免费一区二区三区视频| 日本免费视频| 日本女人日b| 男女视频h| 亚洲国产精品一区二区久久hs| 调教小奴高潮惩罚play道具| 老汉压着稚嫩的小身体h| 小莹与公翁办公室激情| 嗯...嗯啊~嗯...轻一点岳| 亚洲综合图色40p| 337p嫩模大胆色肉噜噜噜| 和女房东啪啪呻吟不断| 少妇口述最爽的3p交换经历| 叶修被狂c躁到高潮失禁| 日本一级性生活| 欧美日韩高清一区二区国产亚洲免费看| 一区二区精品视频在线| 香港黄a三级三级三级看三级| av色影院| 亚洲国产中文在线二区三区免| 亚洲日韩一中文字暮| 少妇一级片免费看| 黑人咬住花蒂猛吸高潮电影| 性侏儒videosfree另类| h漫在线看| 男人干女人逼| chinese高潮videos叫床| 亚洲精品少妇久久久久久| 久久国产精品嫩草影院介绍| 猛烈顶弄h禁欲老师h春潮视频| 天天想夜夜爽| 日本大乳子极品少妇| 久久久久女教师免费一区| 扒开美女猛进爽爽的视频网址大全| 日本大尺度做爰呻吟舌吻| 久久一级黄色片| 黑人干白人| 婷婷久久香蕉五月综合加勒比| gogogo高清免费完整版有一首歌 | 男~女~乱~交| 又黄又免费视频| 黄色搞鸡视频| 亚洲精品乱码久久久久久蜜桃| 大狼拘与少妇牲小说| 和邻居寡妇做爰免费观看| av不卡一区| 白浆av在线| 日本韩国亚洲欧美在线| 久久精品亚洲日本佐佐木明希| 无毛一区二区| 亚洲美女激情视频| 国产精品自在拍在线拍| 欧美人与动xxxxz0oz| 极品美女啪啪| 国产精品另类激情久久久免费| 女教师~淫辱の教室动漫网站| 天天鲁一鲁摸一摸爽一爽视频| av片网站在线观看| 情趣内衣a∨片在线观看| 国产suv精品一区二区| 老师的圣水重口慎入h| 亚洲gay猛男无套网站| 欧美国产激情二区三区| 男人操女人的鸡鸡| 台湾性裸体古装电影| 波多野结衣福利视频| 国产激情艳情在线看视频| 国产精品视频免费一区二区| 四季av一区二区凹凸精品| 日本午夜寂寞影院| 日本少妇一级| 日本的黄色一级片| 久久超碰99| 把雪白的岳弄怀孕了av| 中文字幕乱码亚洲精品一区 | 久久国产精品99精品国产| 秋霞一级黄色片| 免费黄av| av片网站在线观看| 回民做爰xxxx性视频| 亚洲综合图色40p| 亚洲欧美日韩综合一区| 亚洲一码二码三码四码| 亚洲国产精品美女久久久久| 被女同桌高h了女百合扶她 | 男人和女人在曰批40分钟| 久久久综合久久| 女人裸体做爰免费视频| 久久久久久五月天| 欧美三级韩国三级日本一级| 男人扒开我下面舌头伸进去了| 《情欲电车2》日本| a级黄色片网站| 日韩三级国产| 日韩黄色一级| 邻居少妇张开双腿让我爽一夜| 亚洲欧美激情四射| av在线dvd| 日日碰狠狠躁久久躁婷婷| 美女灌肠浣肠无遮挡| 美女被下药脱了个精光| cijilu在线视频| 精东影视传媒mv国产剧能看不| 国产精品成人免费一区久久羞羞| 日韩久久久久久久久| 巨胸喷奶水wwww贱多视频| 天干天干天啪啪夜爽爽av小说| 日日偷换短篇高h合集溪夕汐 | 名人明星裸体hdv三级合集| 日本少妇激三级做爰在线| 亚洲国产日韩精品一区二区三区| 粉嫩av虎白女粉嫩| 宝贝你好湿h| 久久久久久久久久性| 强扒开腿挺进湿润的花苞漫画| 中文字幕亚洲一区二区va在线| 狂c躁到翻白眼高潮失禁| 骚乱np浪荡肉h| 班长哭着说再c就坏掉了| 粗大的内捧猛烈进出| 中文字幕亚洲综合久久| 黑人巨大精品欧美黑寡妇| 电影《曼娜艳荡性史》免费观看| 《艳妇荡乳》在线观看| 国产吃瓜黑料一区二区| 97色精品视频| 夜夜爽日日柔柔日日人人| 欧美性少妇videosex| 啊嗯轻一点进去了第一次网站| 欧美性xxxx极品少妇| 欧美高h视频| 狠狠色狠狠色综合日日五| 少妇口述疯狂的性经历| 小东西你好湿| 日韩精品久久久久久免费| 日本激情按摩| 一边摸一边爽一边叫床免费| 泰国一级片| 欧美三级午夜理伦三级| 对白刺激国产对白精品城中村| vvvv99日韩精品亚洲| 全文都是肉到高潮的小说推荐| 一本一本久久aa综合精品| 脱岳裙子从后面挺进去电影| 欧美日韩久久中文字幕| 日本大尺度r级情欲片我| 人妖黄色片| 日韩一区二区三区射精| 国产丰满麻豆| a一区二区三区乱码在线 | 欧洲| 她双腿大开被绑在床调教ktv| 《无极》床戏视频| 亚洲欧美激情精品一区二区| 三上悠亚与上司在被窝在线观看| 催眠侵犯の奶水授乳羞羞漫画| 欧美黄色片在线播放| 欧美寡妇偷汉性猛交| free性满足hd国产精品久| 又大又紧又爽水又多18p| 日本伦理片情欲电车| 自拍偷拍视频二区| 久久久精品波多野结衣| 捆绑凌虐一区二区三区| 美女mm131爽爽爽| 国产一性一交一伦| 日韩二级片| 亚洲老妈激情一区二区三区| 日韩av三级在线| 97久久精品| 澳门黄色一级片| 国产精品高清一区二区三区不卡| av生活片| 国产精品久久久久久久久久黑人| youjizz国产精品| 办公室勾搭新婚少妇| 欧美老肥妇做.爰bbww| 国内自拍视频一区二区三区| 91麻豆产精品久久久| 女人脱精光免费直播app| 女人与公狼做交十配视频v1| 淫辱的秘书调教| 墨雨云间免费观看完整版| 娇妻在厨房被朋友挺进| 扒开双腿猛进入校花高h的作文| 欧美黄色动态图| 日本体内she精视频| 少妇毛茸茸bbw高清| 给处女破处的过程| 日本三级欧美三级人妇视频黑白配| 日日骚一区二区| 男攻把男受弄到高潮视频| 午夜黄色一级片| av岬奈奈美一区二区三区| 中日韩av电影| 男女互被舔69式免费视频| 揉捏花蒂出绵软呻吟小核嗯啊| 12—13女人做爰视频| 日日摸天天爽天天爽视频| 小说h调教| 亚洲精品国产在线观看| 日本美女的性生活| 99精品国产在热久久婷婷| 久久久综合电影| 黄色无毒视频| 中国女人av| 欧美精品看片| 调教小太正黄腐文h| 男女aa视频| 人人干av| 美女草逼网站| 少妇的性事| 国内精品久久久久影院优| 欧美亚洲二区| 精东影业毛片| 《貂蝉艳史》三级| juliaann欧美二区三区| 揉到高潮嗯啊~喷水求饶免费视频| 特黄特色特刺激视频免费播放| 国产黑色丝袜在线播放| 亚洲精品粉嫩美女一区| 爽到失禁(h)黄文| 337p欧洲亚洲日本大胆色噜噜| 456亚洲人成在线播放网站| 国产精品免费看久久久| 黄色污污视频| 毛片一级黄色| 女女同性百合av大片在线播放| 小伙强行挺进护士体内| 久久婷婷久久| 黄色a在线| 交h粗暴调教小公主h| 久久伊人精品| 婷婷综合av| 国产高潮国语对白精品视频网站| 亚洲偷怕自拍| 黄色片网络| 网红一级片| 黄视频在线观看免费视频| 一级生活黄色片| 国产精品人妖在线| 亚洲久久久| 国产欧美综合一区二区三区| 国产精品福利视频萌白酱app| 亚洲日韩中文字幕一区| 性生生活实战片| 日本黄色片一级片| 欧美黑人一区二区| 91黄瓜视频| 懂色av粉嫩av蜜臀av| 91操人视频| 友田真希三级理伦电影| 被黑人扒开双腿猛进的小说| 多强被c到爽调教sm视频| 娇妻被两黑人按摩韩国电影| 亚洲区日韩| 英语老师穿肉色丝袜让我玩| 丰满少妇被粗大猛烈进| 精品国产一区二区三区观看不卡| 国产亚洲精品久久久久久鸭绿欲| 国产亚洲精品综合一区| 少妇中出视频| youjizzxxxx国语对白| 三上悠亚被老人绝伦| 国产黄色片在线观看| 少妇口述18厘米受不了| 看很黄很黄的细节小说| 99久久99久久久精品齐齐| 《年轻的丝袜老师2》| 口技姿势男人爽到叫| 丰满少妇毛片| 亚洲国产操| 日韩av色| 久久福利电影| 欧美人与禽zozzo性伦交| 综合久久久久久| 国产1区2区| 三上悠亚被邻居疯狂侵犯小说| 大尺度私拍视频| 两个女人互添下身视频在线观看| 国产精品免费视频二三区| 久久精品女人天堂| 日韩做受高潮电影o| 亚洲精品欧美二区三区中文字幕| 亚洲精品国产精品乱码不99按摩| 我揉着护士白嫩的大乳漫画| 亚洲午夜久久久影院| 国产激情麻豆| 精品视频一区在线| 亚洲第一狼人区| 国产精品爽爽久久久久久| 久久久久久久久久久身材好| 免费在线观看av的网站| chinese性旺盛的老妇| 欧美伦理电影一区二区三区| 爆操小护士av| 法国内衣秀| 中文字幕亚洲日韩无线码| 皇上灌到你怀孕为止h| 精品不卡av| 免费的av| 粗大的内捧猛烈进出小视频| 日本老师扒裙子打屁股网站| 我和亲女在客厅作爱h在| 宿舍里的贱奴女虐女被荡秋千玩| 亚洲精品久久久久久中文| 欧美性受xxxx黑人猛交| 大尺度做爰床戏裸体视频| 色狠久久av北条麻妃081| av综合一区| 小泽玛利亚强忍高潮| 2019午夜75福利不卡片在线| 欧美人与禽猛交乱配| 波多野结衣作品番号| 91黄91| 天天操夜夜爽| 国产真实露脸乱子伦| 黄视频在线免费| 裸着身体给老头吃奶|