第十章 藥學(xué)統(tǒng)計(jì)學(xué)

第十章藥學(xué)統(tǒng)計(jì)學(xué)第1頁(yè)，共46頁(yè)，2023年，2月20日，星期三統(tǒng)計(jì)學(xué)（statistics）：研究數(shù)據(jù)的收集、描述、分析、綜合和解釋，以獲得新信息、作出新推斷的學(xué)科。包括：（1）制定調(diào)查方案和設(shè)計(jì)科學(xué)實(shí)驗(yàn)（2）描述實(shí)驗(yàn)數(shù)據(jù)（3）進(jìn)行假設(shè)檢驗(yàn)（4）獲知變量間的定量關(guān)系概述第2頁(yè)，共46頁(yè)，2023年，2月20日，星期三第一節(jié)幾個(gè)基本統(tǒng)計(jì)學(xué)概念頻率：概率：一、必然事件與隨機(jī)事件必然事件：某條件實(shí)現(xiàn)后一定發(fā)生或一定不發(fā)生的事件隨機(jī)事件：一定條件下，不一定發(fā)生的事件二、頻率與概率第3頁(yè)，共46頁(yè)，2023年，2月20日，星期三1.總體與樣本的關(guān)系樣本參數(shù)用英文字母表示總體參數(shù)用希臘字母表示總體樣本

取樣觀測(cè)數(shù)據(jù)推斷三、總體與樣本第4頁(yè)，共46頁(yè)，2023年，2月20日，星期三

2.取樣的隨機(jī)性隨機(jī)性包括：①總體中個(gè)體的抽取必須是相互獨(dú)立的；②總體中所有個(gè)體被抽取的機(jī)會(huì)相等。滿足以上兩個(gè)要求的取樣，稱為簡(jiǎn)單隨機(jī)取樣（SRS，simplerandomsampling）這樣抽取的樣本稱為簡(jiǎn)單隨機(jī)樣本。第5頁(yè)，共46頁(yè)，2023年，2月20日，星期三四、觀測(cè)值的特征——集中位置與離散程度量度集中位置的統(tǒng)計(jì)量（1）均值樣本均值總體均值μ＝E（X）＝Σ取值·概率＝ΣXiPi均值的重要性質(zhì)：①觀測(cè)值與均值之差（偏差）之和②偏差的平方和最小第6頁(yè)，共46頁(yè)，2023年，2月20日，星期三（2）眾數(shù)頻數(shù)最大的觀測(cè)值叫眾數(shù)，常用于表示離散型隨機(jī)變量的集中位置。（3）中位數(shù)把變量的觀測(cè)值按大小順序排列，排在當(dāng)中的一個(gè)觀測(cè)值叫中位數(shù)。第7頁(yè)，共46頁(yè)，2023年，2月20日，星期三（1）極差（最大值和最小值之差）

R＝Xmax－Xmin（2）平均（絕對(duì)）偏差2.量度數(shù)據(jù)離散程度的統(tǒng)計(jì)量（3）方差式中n－1在統(tǒng)計(jì)學(xué)中叫自由度（degreeoffreedom），常用希臘字母ν表示第8頁(yè)，共46頁(yè)，2023年，2月20日，星期三五、兩個(gè)典型的概率分布1.二項(xiàng)分布是典型的離散型概率分布（1）特點(diǎn)：二項(xiàng)分布是一種每次試驗(yàn)只有兩種可能結(jié)果而不受以前實(shí)驗(yàn)影響的分布。在樣本容量相對(duì)于總體很小時(shí)，取樣試驗(yàn)后，返回不返回，對(duì)下一次試驗(yàn)的結(jié)果無(wú)影響。第9頁(yè)，共46頁(yè)，2023年，2月20日，星期三（2）二項(xiàng)分布的均值、方差、標(biāo)準(zhǔn)差均值＝np方差＝npq

標(biāo)準(zhǔn)差＝（3）二項(xiàng)分布的實(shí)際意義一些只有兩個(gè)結(jié)果的互斥事件都與二項(xiàng)分布有關(guān)，在藥學(xué)方面，可用于分析從制劑批抽出n個(gè)個(gè)體中不合格個(gè)體數(shù)的概率。第10頁(yè)，共46頁(yè)，2023年，2月20日，星期三2.正態(tài)分布是一種連續(xù)型概率分布，各種分布都以它為中心在一定條件下相互轉(zhuǎn)化。

設(shè)連續(xù)性隨機(jī)變量X的概率密度為

則稱X服從參數(shù)為、（-∞<<+∞，＞0）的正態(tài)分布，記為X～N（，

）第11頁(yè)，共46頁(yè)，2023年，2月20日，星期三（1）特點(diǎn)：①極大值在μ；②以X＝μ為對(duì)稱軸；③X軸是漸進(jìn)線；④拐點(diǎn)在X＝μ±σ。第12頁(yè)，共46頁(yè)，2023年，2月20日，星期三（2）標(biāo)準(zhǔn)正態(tài)分布①一般的正態(tài)分布取決于均值和標(biāo)準(zhǔn)差②計(jì)算概率時(shí)，每一個(gè)正態(tài)分布都需要有自己的正態(tài)概率分布表，這種表格是無(wú)窮多的若能將一般的正態(tài)分布轉(zhuǎn)化為標(biāo)準(zhǔn)正態(tài)分布，計(jì)算概率時(shí)只需要查一張表第13頁(yè)，共46頁(yè)，2023年，2月20日，星期三任何一個(gè)一般的正態(tài)分布，可通過(guò)下面的線性變換轉(zhuǎn)化為標(biāo)準(zhǔn)正態(tài)分布xms一般正態(tài)分布=1Z標(biāo)準(zhǔn)正態(tài)分布第14頁(yè)，共46頁(yè)，2023年，2月20日，星期三（3）標(biāo)準(zhǔn)正態(tài)分布表的使用①將一個(gè)一般的轉(zhuǎn)換為標(biāo)準(zhǔn)正態(tài)分布計(jì)算概率時(shí)，查標(biāo)準(zhǔn)正態(tài)概率分布表②對(duì)于負(fù)的x，可由(-x)x得到③對(duì)于標(biāo)準(zhǔn)正態(tài)分布，即X~N(0,1)，有P(aXb)baP(|X|a)2a1④對(duì)于一般正態(tài)分布，即X~N(,)，有第15頁(yè)，共46頁(yè)，2023年，2月20日，星期三六、均值的分布和中心極限定理均值的分布

如一個(gè)隨機(jī)變量X呈總體均值為μ、總體方差為σ2的正態(tài)分布，則其容量為n的樣本均值也呈正態(tài)分布，其總體均值仍為μ，但總體方差為σ2/n，即X第16頁(yè)，共46頁(yè)，2023年，2月20日，星期三2.中心極限定理設(shè)從均值為，方差為2的一個(gè)任意總體中抽取容量為n的樣本，當(dāng)n充分大時(shí)，樣本均值的抽樣分布近似服從均值為μ、方差為σ2/n的正態(tài)分布一個(gè)任意分布的總體X當(dāng)樣本容量足夠大時(shí)(n

30)，樣本均值的抽樣分布逐漸趨于正態(tài)分布第17頁(yè)，共46頁(yè)，2023年，2月20日，星期三第二節(jié)取樣

按與調(diào)查目的有關(guān)的某個(gè)主要性標(biāo)志將總體單位劃分為若干層（也稱類、組或子總體），然后從多層中按隨機(jī)原則分別抽取一定數(shù)目的單位構(gòu)成樣本。一、隨機(jī)取樣與隨機(jī)數(shù)表簡(jiǎn)單隨機(jī)抽樣是最基本的抽樣組織方式。是對(duì)總體單位不進(jìn)行任何劃分或排隊(duì)，完全隨機(jī)地直接從總體中抽取樣本單位，使每個(gè)總體單位都有完全均等的機(jī)會(huì)被抽中。二、分層取樣第18頁(yè)，共46頁(yè)，2023年，2月20日，星期三

在大規(guī)模生產(chǎn)中，表明供應(yīng)方產(chǎn)品質(zhì)量特征的隨機(jī)變量應(yīng)不超出某規(guī)定范圍。三、系統(tǒng)取樣

常用于質(zhì)量控制，每第n個(gè)個(gè)體選取一個(gè)，但最初選取的個(gè)體是隨機(jī)的。四、驗(yàn)收取樣第19頁(yè)，共46頁(yè)，2023年，2月20日，星期三第三節(jié)數(shù)據(jù)制圖

與統(tǒng)計(jì)方法有關(guān)的圖解：①以傳遞信息為目的的圖解：如直方圖、條形圖②在坐標(biāo)中描述變量之間關(guān)系的圖解一、引言第20頁(yè)，共46頁(yè)，2023年，2月20日，星期三在許多研究數(shù)據(jù)中，響應(yīng)Y的對(duì)數(shù)而不是Y本身與獨(dú)立變量X呈線性關(guān)系，這樣就要進(jìn)行半對(duì)數(shù)標(biāo)繪，即變量X的坐標(biāo)是自然數(shù)而響應(yīng)（函數(shù)）Y的是它的對(duì)數(shù)logY。二、描述頻數(shù)的圖解直方圖、條形圖、扇形圖三、描述變量關(guān)系的標(biāo)繪圖第21頁(yè)，共46頁(yè)，2023年，2月20日，星期三第四節(jié)統(tǒng)計(jì)估計(jì)和假設(shè)檢驗(yàn)描述統(tǒng)計(jì)推斷統(tǒng)計(jì)統(tǒng)計(jì)估計(jì)假設(shè)檢驗(yàn)統(tǒng)計(jì)方法第22頁(yè)，共46頁(yè)，2023年，2月20日，星期三樣本統(tǒng)計(jì)量均值和方差S2分別是其總體均值μ和總體方差σ2的最佳估計(jì)。置信區(qū)間（confidenceinterval）：是我們相信統(tǒng)計(jì)量如總體均數(shù)所在的區(qū)間，由總體的性質(zhì)、參數(shù)的樣本估計(jì)值和想達(dá)到的置信度決定。以樣本均值為對(duì)稱中心的雙側(cè)置信區(qū)間：

P％置信區(qū)間＝X一、統(tǒng)計(jì)估計(jì)第23頁(yè)，共46頁(yè)，2023年，2月20日，星期三1.用t分布確定置信區(qū)間P％置信區(qū)間＝2.幾個(gè)不同置信區(qū)間的構(gòu)造（1）σ未知，需要由樣本估計(jì)（2）σ已知，直接代入公式（3）比率的置信區(qū)間一、統(tǒng)計(jì)估計(jì)第24頁(yè)，共46頁(yè)，2023年，2月20日，星期三二、假設(shè)檢驗(yàn)1.引言（1）概念①事先對(duì)總體參數(shù)或分布形式作出某種假設(shè)②然后利用樣本信息來(lái)判斷原假設(shè)是否成立（2）類型①參數(shù)假設(shè)檢驗(yàn)②非參數(shù)假設(shè)檢驗(yàn)（3）特點(diǎn)采用邏輯上的反證法依據(jù)統(tǒng)計(jì)上的小概率原理第25頁(yè)，共46頁(yè)，2023年，2月20日，星期三零假設(shè)（nullhypothesis）：關(guān)于兩個(gè)總體在某一參數(shù)如均值上一致（實(shí)際上差異為零）的假設(shè)H0。如H0：μ1＝μ2備擇假設(shè)（alternativehypothesis）：如果舍棄H0則接受Ha，Ha也可用H1表示。如Ha：μ1≠μ2第26頁(yè)，共46頁(yè)，2023年，2月20日，星期三例：表觀測(cè)單一總體均值的實(shí)驗(yàn)示例同類藥物的治愈率新藥的治愈率值得臨床研究重視得血壓降低標(biāo)準(zhǔn)值臨床前研究中n只大鼠血壓降低均值藥典片劑通則關(guān)于崩解時(shí)限的規(guī)定藥片片劑平均崩解時(shí)限標(biāo)示量藥片片劑的平均含量假設(shè)值或標(biāo)準(zhǔn)值樣本均值對(duì)此例中每一類問(wèn)題的假設(shè)檢驗(yàn)，可提出三個(gè)假設(shè)檢驗(yàn)：①H0:μ=μ0,Ha:μ

≠

μ0②H0:μ=μ0,Ha:μ<μ0③H0:μ=μ0,Ha:μ>μ0第27頁(yè)，共46頁(yè)，2023年，2月20日，星期三2.檢驗(yàn)假設(shè)的步驟（1）提出檢驗(yàn)統(tǒng)計(jì)量T并進(jìn)行實(shí)驗(yàn)（2）提出零假設(shè)H0和備擇假設(shè)Ha（3）規(guī)定顯著性水平α顯著性水平：即根據(jù)統(tǒng)計(jì)檢驗(yàn)的結(jié)果舍棄零假設(shè)H0（存在顯著差異）而實(shí)際上不存在的概率。（4）計(jì)算檢驗(yàn)統(tǒng)計(jì)量T并進(jìn)行顯著性檢驗(yàn)第28頁(yè)，共46頁(yè)，2023年，2月20日，星期三3.幾種具體的假設(shè)檢驗(yàn)（1）σ已知的單樣本雙側(cè)均值檢驗(yàn)（2）σ已知的單樣本單側(cè)均值檢驗(yàn)（3）σ未知的單樣本雙側(cè)均值檢驗(yàn)先求t值，再計(jì)算置信區(qū)間，進(jìn)行t檢驗(yàn)（4）獨(dú)立樣本方差S12和S22的比較：F檢驗(yàn)第29頁(yè)，共46頁(yè)，2023年，2月20日，星期三1.什么是回歸分析：①?gòu)囊唤M樣本數(shù)據(jù)出發(fā)，確定變量之間的數(shù)學(xué)關(guān)系式②對(duì)這些關(guān)系式的可信程度進(jìn)行各種統(tǒng)計(jì)檢驗(yàn)，并從影響某一特定變量的諸多變量中找出哪些變量的影響顯著，哪些不顯著③利用所求的關(guān)系式，根據(jù)一個(gè)或幾個(gè)變量的取值來(lái)預(yù)測(cè)或控制另一個(gè)特定變量的取值，并給出這種預(yù)測(cè)或控制的精確程度第五節(jié)回歸分析第30頁(yè)，共46頁(yè)，2023年，2月20日，星期三2.回歸模型的類型一個(gè)自變量?jī)蓚€(gè)及兩個(gè)以上自變量回歸模型多元回歸一元回歸線性回歸非線性回歸線性回歸非線性回歸第31頁(yè)，共46頁(yè)，2023年，2月20日，星期三3.回歸模型（1）回答“變量之間是什么樣的關(guān)系？”（2）方程中運(yùn)用①1個(gè)數(shù)字的因變量(響應(yīng)變量)——被預(yù)測(cè)的變量②1個(gè)或多個(gè)數(shù)字的或分類的自變量(解釋變量)——用于預(yù)測(cè)的變量（3）主要用于預(yù)測(cè)和估計(jì)第32頁(yè)，共46頁(yè)，2023年，2月20日，星期三概念要點(diǎn)（1）當(dāng)只涉及一個(gè)自變量時(shí)稱為一元回歸，若因變量y與自變量x之間為線性關(guān)系時(shí)稱為一元線性回歸（2）對(duì)于具有線性關(guān)系的兩個(gè)變量，可以用一條線性方程來(lái)表示它們之間的關(guān)系（3）描述因變量y如何依賴于自變量x的方程稱為回歸模型第33頁(yè)，共46頁(yè)，2023年，2月20日，星期三5.簡(jiǎn)單線性回歸方程的形式如下

Y=0+1x

方程的圖示是一條直線，因此也稱為直線回歸方程0是回歸直線在y軸上的截距，是當(dāng)x=0時(shí)y的期望值1是直線的斜率，稱為回歸系數(shù)，表示當(dāng)x每變動(dòng)一個(gè)單位時(shí)，y的平均變動(dòng)值第34頁(yè)，共46頁(yè)，2023年，2月20日，星期三6.最小二乘法

回歸直線是對(duì)若干個(gè)數(shù)據(jù)對(duì)（X,Y）繪制的，但由于遵守正態(tài)分布的實(shí)驗(yàn)誤差，數(shù)據(jù)對(duì)中兩個(gè)以上的點(diǎn)，不能恰好在一條直線上，即恰好通過(guò)所有數(shù)據(jù)點(diǎn)的直線是找不到的，只能找到一條離所有點(diǎn)都近的直線，所用的客觀方法叫最小二乘法（methodofleastsquares）。第35頁(yè)，共46頁(yè)，2023年，2月20日，星期三7.最小二乘線（1）定義給定n個(gè)數(shù)據(jù)對(duì)（X,Y），找到一條規(guī)定X，Y關(guān)系的直線，使所有數(shù)據(jù)對(duì)（X，Y）與擬規(guī)定的直線在Y軸方向的距離平方和最小——用數(shù)學(xué)語(yǔ)言表示就是最小，這樣的一條直線就是描述X，Y間關(guān)系的最佳直線，稱為最小二乘線。第36頁(yè)，共46頁(yè)，2023年，2月20日，星期三（2）圖示xy(xn,yn)(x1,y1)(x2,y2)(xi,yi)｝ei=yi-yi＾第37頁(yè)，共46頁(yè)，2023年，2月20日，星期三（3）和的計(jì)算公式第38頁(yè)，共46頁(yè)，2023年，2月20日，星期三二、回歸分析在藥物研究中的應(yīng)用例：把6個(gè)濃度按兩倍遞增的純青霉素溶液（1~32U/ml）置于進(jìn)行生物測(cè)定的杯碟中。下表給出每一濃度以mm表示的抑菌圈直徑，在這個(gè)長(zhǎng)度的測(cè)量中發(fā)生較大的誤差。如果對(duì)青霉素溶液的濃度C取以2為底的對(duì)數(shù)（筆，binit），則抑菌圈的直徑與青霉素溶液的濃度成線性關(guān)系。試求Y在X上的回歸直線方程。24.7723.1321.3519.5217.7815.87Y＝抑菌圈直徑/mm543210X＝log2C32168421C/ml表青霉素溶液濃度的對(duì)數(shù)變換第39頁(yè)，共46頁(yè)，2023年，2月20日，星期三解進(jìn)行Y在X上的回歸計(jì)算

ΣX＝15n＝6ΣY＝122.42＝2.50＝20.40

ΣX2＝55ΣXY＝337.24ΣY2＝2553.3880于是Y在X上的回歸直線方程是＝15.94＋1.782X

第40頁(yè)，共46頁(yè)，2023年，2月20日，星期三第六節(jié)實(shí)驗(yàn)設(shè)計(jì)1.基本原理（1）對(duì)照：齊同可比（2）重復(fù)：①獲得總體標(biāo)準(zhǔn)差的估計(jì)值②用大容量樣本重復(fù)多次獲得的均值，能把因素在實(shí)驗(yàn)中的效應(yīng)估計(jì)得更精密。（3）隨機(jī)：有助于把可能存在得外部因素效應(yīng)平均化而減免。一、實(shí)驗(yàn)設(shè)計(jì)的基本原理第41頁(yè)，共46頁(yè)，2023年，2月20日，星期三2.步驟①確認(rèn)并陳述問(wèn)題②選擇因素和水平③選擇響應(yīng)變量④選擇實(shí)驗(yàn)設(shè)計(jì)⑤進(jìn)行實(shí)驗(yàn)⑥分析數(shù)據(jù)⑦做結(jié)論并提建議第42頁(yè)，共46頁(yè)，2023年，2月20日，星期三二、方差分析（analysisofvariance，ANOVA）方差分析能確定引起生產(chǎn)和實(shí)驗(yàn)