[摘 要]傳統(tǒng)的抽樣調查課程教學過于注重理論推導，缺乏實際應用，以至于學生在學習過程中未能真正掌握其中方法，不懂得如何運用所學的方法解決實際問題。課題組對實際抽樣時整群抽樣的抽樣與估計步驟進行了系統(tǒng)的講解，并結合R軟件完成了整群抽樣的樣本抽取與總體參數估計。通過案例教學，幫助學生掌握整群抽樣的理論知識，并提高了學生運用R軟件解決實際問題的能力。

[關鍵詞]案例教學；抽樣調查；整群抽樣；R軟件

[中圖分類號] G642 [文獻標識碼] A [文章編號] 2095-3437（2018）07-0065-03

○、引言

抽樣調查作為一門應用性、實用性很強的統(tǒng)計學科，是統(tǒng)計學類專業(yè)的專業(yè)必修課程。然而，傳統(tǒng)的抽樣調查課程教學過于注重理論推導，缺乏實際應用，以至于學生在學習過程中未能真正掌握其中方法，不懂得如何運用所學的方法解決實際問題。為此，候震梅（2016）[1]探討了抽樣調查課程實驗教學改革的思路與途徑；張學新（2015）[2]提出一種抽樣調查課程貫穿統(tǒng)計軟件使用的教學改革方法；盧玉桂（2015）[3]提出應運用案例教學法進行抽樣調查課程教學。

整群抽樣具有實施調查便利、節(jié)省費用的優(yōu)點，因而被廣泛應用于調查中。但筆者在教授該知識點時，發(fā)現大部分學生無法真正理解與掌握整群抽樣方法。因此，筆者擬在前人研究基礎[4][5]上，結合多年教授抽樣調查課程的經驗，探討如何利用R軟件幫助學生更好的理解與掌握“整群抽樣”抽樣與估計過程。

一、整群抽樣簡介

整群抽樣（cluster sampling）[5]是將總體劃分為若干個群，然后以群為抽樣單元，從總體N（總體群數）抽取n個樣本群，并對樣本群中的所有單元進行調查的一種抽樣方法。整群抽樣分為群規(guī)模相等和群規(guī)模不等兩種情況。所謂群規(guī)模是指組成群的單元數量，群規(guī)?？纱罂尚　Ｈ阂?guī)模大，則估計精度差但費用省；群規(guī)模小，則可提高估計精度但費用也高。因此，一般而言，群規(guī)模不宜過大。本文主要介紹群規(guī)模相等的整群抽樣（也稱為等概率整群抽樣），等概率整群抽樣是指在總體N個群（初級抽樣單元）中，每個群所包含的單元（二級單元）數皆等于M。

等概率整群抽樣一般用簡單隨機抽樣方法抽取群，此時抽樣比為f=n / N，則總體均值[y]的無偏估計為：

二、群抽樣的R軟件實現

筆者在進行整群抽樣教學時，發(fā)現多數學生學起來吃力，無法理解整群抽樣與簡單[?]隨機抽樣的關系，以至于產生厭學、不學等一些負面情緒。為了幫助學生更好的掌握整群抽樣方法，筆者將通過文獻[5]的例子，講解整群抽樣過程。

例 某郵局該轄區(qū)共有5000戶，并劃分為500個群，則每個群有10戶居民。為了解郵局管轄區(qū)內每個家庭的月平均定報份數及其95%的置信區(qū)間，運用整群抽樣方法從500個群中抽取5個群進行抽樣調查。

本例中，每個群的規(guī)模均為10，故為等概率整群抽樣。因此，抽樣方法為運用簡單隨機抽樣從500個群中抽取5個群。因為是按照簡單隨機抽樣抽取群的，所以每個群的入樣概率均為5/500=1/100，又因為總體中的某個群一旦被抽中，則群內的所有單元全部入樣。因此，每個總體單元的入樣概率都為1/100，這也驗證了群規(guī)模相等的整群抽樣為何也可稱為等概率整群抽樣。

實施抽樣調查時，抽樣與估計是分開進行的，且先抽樣后估計。整群抽樣的抽樣步驟可分為以下兩步：

第一步：編制抽樣總體數據框，即編制抽樣框。

假設事先已獲取管轄區(qū)內5000戶住戶的名錄，并已劃分為500個群。注意，本例中為簡單起見，住戶名錄用數字編號代替，實際抽樣時，應為真實名錄，這樣完成抽樣后，才能根據抽中的名錄清單，開展調查。

其R代碼如下：

> a=rep（1：500，each=10）

> b=1：5000

> data=data.frame（“群號”=a，“住戶ID”=b）

第二步：調用整群抽樣函數cluster，完成抽樣。

進行抽樣前，首先需要下載并加載sampling包（抽樣程序包），然后才可以調用整群抽樣函數。整群抽樣函數cluster（ ）的第一個參數為總體數據框，第二個參數為分群變量，參數size設定樣本群數，參數method設定抽樣方法，包括不放回簡單隨機抽樣（srswor）、放回簡單隨機抽樣（srswr）、泊松抽樣（poisson）、系統(tǒng)抽樣（systematic）4種抽樣方法，description為邏輯型向量，取TRUE值，表示輸出抽樣信息。

其R代碼如下：

> install.packages（‘sampling） #下載sampling包

> library（sampling）#加載sampling包

> n=5#樣本群數

#調用整群抽樣函數

>data.c=cluster（data，“群號”，size=n，method=“srswor”，description=TRUE）

Number of selected clusters：5

Number of units in the population and number of selected units：5000 50

#從總體數據框中提取樣本數據

> data.c=getdata（data，data.c）

>data.c[1：20，]#顯示前20個樣本的信息

結果顯示：整群抽樣從總體群500中抽取了5個群，入樣的群編號依次為72、140、146、430、460，且總體單元數為5000，樣本單元數為50；整群抽樣函數cluster（ ）返回被抽中單元的住戶ID、群號、抽樣單元編號（ID_unit）和入樣概率（Prob）。

完成抽樣后，將對入樣的住戶展開調查，獲取相關樣本數據，完成樣本數據的收集。假設經過調查后，獲取的樣本數據為表2。

根據獲取的樣本數據，可開始進行總體參數的估計。整群抽樣的估計步驟也可以分為兩步：

第一步：估計的前期工作。

進行抽樣估計前，需要下載并加載survey包（估計程序包），另外還需要加載一個grid基礎包，才能開始進行抽樣估計。另外，抽樣估計前，需要調用svydesign函數定義抽樣設計，但在使用該函數前，還需要設定每個樣本的樣本權重pw和fpc。其中，pw為入樣概率的倒數N/n，即為500/5=100，fpc為總體群數N，即為500。其R代碼如下：

> install.packages（‘survey）#下載survey包

> library（survey） #加載survey包

> library（grid）#加載grid包

> N=500#設定總體群數

> n=5#設定樣本群數

> pw=rep（N/n，nrow（data1））#設定樣本權重

> fpc=rep（N，nrow（data1））#設定fpc變量

> data1.c=as.data.frame（cbind（data1，pw，fpc））#合并樣本數據框

> data1.c[1：5，]#顯示前5行樣本數據框

第二步：調用svydesign函數定義抽樣設計，并完成估計。抽樣設計函數svydesign（ ）參數id定義群變量，參數weights定義樣本權重，參數data定義樣本數據框，參數fpc定義fpc變量。

> d.c<-svydesign（id=～群號，weights=～pw，data=data1.c，fpc=～fpc）#抽樣設計

> summary（d.c）#查看抽樣設計

1 - level Cluster Sampling design

With （5） clusters.

svydesign（id = ～群號，weights = ～pw，data = data1.c，fpc = ～fpc）

Probabilities：

Min. 1st Qu. Median Mean 3rd Qu. Max.

0.01 0.01 0.01 0.01 0.01 0.01

Population size （PSUs）：500

Data variables：

[1] “住戶ID”“群號”“ID_unit”“Prob”“訂報數” “pw”“fpc”

> svymean（～訂報數，d.c）#訂報數的均值估計和標準誤差

mean SE

訂報數 2.02 0.1062

結果顯示：訂報數的均值估計值為2.02份，標準誤差為0.1062，則訂報數的均值估計值的95%置信區(qū)間為[1.8119，2.2282]。

三、結論

本文以一個簡單的抽樣案例，對整群抽樣的抽樣與估計過程進行了系統(tǒng)的講解，并結合R軟件完成了整群抽樣的抽樣與估計。通過運用案例教學方法進行整群抽樣的課堂教學，不僅可以幫助學生理解與掌握整群抽樣的理論知識，同時還可以提高學生靈活運用R軟件的能力，以及運用R軟件解決實際抽樣調查的能力。為此，在進行抽樣調查課程教學時，應注意運用實際案例和R軟件演示相結合進行教學，這將有利于學生對抽樣理論知識理解與掌握，同時提高學生運用R軟件解決實際問題的能力。

[ 參 考 文 獻 ]

[1] 侯震梅. 《抽樣調查》實驗教學創(chuàng)新性研究[J]. 現代商貿工業(yè)，2016（10）：168-169.

[2] 張學新. 《抽樣調查》課程的統(tǒng)計軟件教學方法實踐[J]. 寧夏師范學院學報，2015（6）：83-91.

[3] 盧玉桂. 案例教學法在《抽樣調查》課程教學中的應用[J]. 亞太教育，2015（4）：104+93.

[4] 王偉，陳志軍，徐辰武. 基于R語言的隨機抽樣方法及其應用[J]. 揚州大學學報（農業(yè)與生命科學版），2014（2）：77-81.

[5] 金勇進等.抽樣技術：第四版[M].北京：中國人民大學出版社，2015.

[責任編輯：林志恒]