虛擬變量

虛擬變量（dummy variable），又稱啞變量，指在統計學和計量經濟學，尤其是迴歸分析中僅取0或1的值，以此表示某些可能改變結果的定性效應之有無的變量。可以認為，虛擬變量是迴歸模型中用數字來替代並表示定性事實，從而將數據分為互斥的類別（如吸菸者及非吸菸者）。
虛擬自變量（又稱虛擬解釋變量）取值為0時，該變量的係數對因變量沒有影響，而當虛擬變量取值為1時，其係數會改變截距。例如，假定某人是否為某組的成員是與迴歸相關的定性變量之一，如果屬於該組成員以1代表，則其他所有人獲得0值，此時截距對非成員是常數項，對成員則是常數項加上“成員”虛擬變量的係數。
虛擬變量常用於時間序列分析，包括狀態轉換、季節性分析和定性數據應用。

**引入虛擬自變量**
1=wage = α0 + δ0female + α1education + U的圖像，其中δ0 0可以理解為教育條件一樣時，女性和男性之間的工資差異。因此，δ0有助於確定男女之間的工資是否存在差異。例如，如果δ0 > 0（正係數），則女性的工資高於男性（其他因素不變）。附加在虛擬變量前的係數稱為差別截距係數（differential intercept coefficient）。在模型的圖像上，它可視為女性及男性之間截距的差異。圖中顯示的為δ0 iD2i = 1, D3i = 0) = α1 + α2
南部地區公立學校教師平均工資：
 1=E(YiD2i = 0, D3i = 1) = α1 + α3
西部地區公立學校教師平均工資：
 1=E(YiD2i = 0, D3i = 0) = α1
（期望值中不含誤差項，因為我們假設它滿足通常的普通最小二乘法條件，即1=E(ui) = 0。）
期望值可以解釋如下：西部公立學校教師的平均工資等於多元迴歸方程中的截距項α1，加上差別截距係數α2和α3，後者解釋為南部、北部地區的教師工資均值與西部地區相差多少。因此，北部和南部教師的平均工資與西部教師的平均工資進行比較。西部地區因此成為基準組，即與之進行比較的組。省略的組別（即沒有分配虛擬變量的組別）作為基準組別。
使用這些數據進行迴歸的結果為：
 1=Ŷi = 26,158.62 − 1734.473D2i − 3264.615D3i
 1=se =           (1128.523)   (1435.953)               (1499.615)
 1=t  =            (23.1759)    (−1.2078)                (−2.1776)
 1=p  =             (0.0000)     (0.2330)                 (0.0349)
 1=R2 = 0.0901
其中，se = 標準誤差，t = t-statistic，p = p值。
迴歸結果可以解釋為：西部教師（基準組）的平均工資約為$26,158；與之相比，北部教師的工資低約$1734（$26,158.62 - $1734.473 = $24,424.14，即平均工資北部教師的工資），南部教師的工資低約$3265（$26,158.62 - $3264.615 = $22,894，即南部教師的平均工資）。
要確定南、北部教師的平均工資與西部教師的平均工資（比較組別）在統計學上是否有差異，我們必須找出迴歸結果的斜率係數是否具有統計學意義。為此，我們需要考慮p值。北部的估計斜率係數在統計上並不顯著，因為它的p值為23%；然而，南部的p值僅在3.5%左右，在5%的水平上具有統計顯著性。因此，總體結果是：西部和北部教師的平均工資在統計上沒有顯著差異，但南部教師的平均工資在統計上比西部低$3265左右。模型如圖2所示。該模型是一個方差分析模型，其中，一個定性變量具有3個類別。

**有兩個定性變量的方差分析模型**
考慮具有兩個定性變量的方差分析模型，每個變量都有兩個類別：因變量“時薪”使用定性變量“婚姻狀況”（已婚/未婚）和“地理區域”（北部/非北部）來解釋。其中，婚姻狀況和地理區域是兩個解釋性虛擬變量。
假設基於某些給定數據的迴歸輸出如下所示：
 1=Ŷi = 8.8148 + 1.0997D2 − 1.6729D3
其中，
 Y為時薪（單位：$）
 D2為婚姻狀況，1 = 已婚，0 = 未婚
 D3為地理區域，1 = 北部，0 = 非北部
該模型為每個定性變量分配了一個虛擬變量，比每個變量中包含的類別數少一個。
這裡，基準組是省略的組別：未婚、非北方地區的人。所有比較都與該基準組或省略的組別相關。基準組的平均時薪約為$8.81（截距項）。相比之下，已婚人士的平均時薪高出約$1.10，約為$9.91（$8.81 + $1.10）。相比之下，居住在北方的人的平均時薪低約$1.67，約為$7.14（$8.81 - $1.67）。
一般來說，如果迴歸中包含一個以上的定性變量，則應選擇一個省略的組別作為基準組。所有比較都與該組相關。截距項將顯示基準組的期望值，斜率係數顯示其他組別與基準（省略）組別的差異程度。

**協方差分析模型**
包含定量和定性變量的混合迴歸模型稱為Analysis of covariance（ANCOVA）模型。 協方差分析模型是方差分析模型的擴展。它們統計控制定量解釋變量（也稱為協變量或控制變量）的影響。
為了說明如何包含定性和定量回歸變量來形成協方差分析模型，假設我們考慮在方差分析模型中使用的相同示例，其中有一個定性變量：A國三個地理區域的公立學校教師的平均年薪。如果我們包括一個定量的變量：州政府對每名公立學校學生的財政支出，則這一回歸中，可得到以下模型：
圖3：協方差分析模型示例的迴歸結果：公立學校教師的年薪（Y）與公立學校每名學生的州財政支出之間的關係。
 1=Yi = α1 + α2D2i + α3D3i + α4Xi + Ui
其中，
 Yi = i州公立學校教師的平均年薪
 Xi = 公立學校每名學生的州財政支出
 D2i = 1，若i州在北部地區
: 否則D2i = 0
 D3i = 1，若i州在南部地區
: 否則D3i = 0
假設這個模型的迴歸輸出是
 1=Ŷi = 13,269.11 − 1673.514D2i − 1144.157D3i + 3.2889Xi
結果表明，公立學校每名學生的州財政支出每增加$1，公立學校教師的平均工資就會增加約$3.29。此外，對於北部地區的一個州，教師的平均工資比西部地區低約$1673，對於南部地區的一個州，教師的平均工資比西部地區低約$1144。圖3描繪了該模型。假設支出係數不因州而異，則平均工資線彼此平行。圖中分別顯示了每個類別中，兩個定量變量之間的關係，即公立學校教師的工資（Y）與公立學校每名學生的州財政支出（X）之間的關係。

**虛擬變量間的交互作用**
迴歸模型中的定量回歸變量通常存在交互作用。同樣，定性的迴歸變量，即虛擬變量之間也可能有交互作用，這些交互作用可在迴歸模型中描述。例如，在涉及工資確定的迴歸中，如果考慮兩個定性變量：性別和婚姻狀況，則婚姻狀況和性別之間可能存在交互作用。這些交互作用可以在迴歸方程中顯示，如下例所示。
由於兩定性變量是性別和婚姻狀況，定量的解釋變量是受教育年限，解釋變量純線性的迴歸將是
 1=Yi = β1 + β2D2,i + β3D3,i + αXi + Ui
其中
 i表示某個特定個人
 Y = 時薪（單位：$）
 X = 受教育年限
 D2 = 1若為女性，否則為0
 D3 = 1若已婚，否則為0
這一模型無法表達兩個定性變量D2和D3之間發生交互作用的可能性。例如，已婚女性的工資與未婚男性的工資差異，其數額與單獨為女性和單獨為已婚的差值之和不同。為了考慮這種可能性，工資的確定可以指定為：
 1=Yi = β1 + β2D2,i + β3D3,i + β4(D2,iD3,i) + αXi + Ui
其中，
 β2 = 作為女性的差別效應
 β3 = 已婚的差別效應
 β4 = 女性且已婚的進一步差別效應
根據該等式，在沒有非零誤差的情況下，未婚男性的工資為β1 + αXi，未婚女性的工資為β1 + β2 + αXi，已婚男性的工資為β1 + β3 + αXi，而已婚女性的比例是β1 + β2 + β3 + β4 + αXi（其中任何虛擬變量係數的估計值都可能是為正、零，或負）。
由此，交互虛擬變量（兩個虛擬變量的乘積）會改變因變量，使之與單獨考慮兩個虛擬變量時得到的值有所不同。
或者，可以通過使用不同的數據分類方案來避免使用虛擬變量的乘積來表達交互作用，即使用根據特徵組合指定組別的方案。如果令
 D4 = 1若為未婚女性，否則為0
 D5 = 1若為已婚男性，否則為0
 D6 = 1若為已婚女性，否則為0
那麼僅需指定如下回歸：
 1=Yi = δ1 + δ4D4,i + δ5D5,i + δ6D6,i + αXi + Ui
那麼在零誤差項下，因變量的值是對於基準組未婚男性為δ1 + αXi，對於未婚女性為δ1 + δ4 + αXi，對於已婚男性為δ1 + δ5 + αXi，對於已婚女性為δ1 + δ6 + αXi。該方法在等號右側使用的變量數與先前具有交互項的方法相同，並且使用本方法還是指定交互項目的方法，對於Xi下任意定性特徵組合下的因變量預測值，給出的迴歸結果是相同的。

**虛擬因變量**

**因變量是虛擬變量的情景**
具有虛擬因變量（也稱為定性因變量）的模型中，因變量受解釋變量的影響是定性的。例如，一些關於執行“多少”行為的決定涉及事先決定是否執行該行為；“事前決策”的迴歸在迴歸模型中有一個因虛擬變量。
例如，潛在工人成為勞動力一部分的決策是一個虛擬因變量。該決策是二分的，即有兩種可能的結果：是和否。因此，虛擬因變量Participation在“參與”時取值為1，“不參與”則取值為0。下面給出其他一些二元的虛擬因變量的例子：- 決策：職業選擇。虛擬因變量：若選擇成為高管，則 高管=1，否則為0。- 決策：隸屬於某政黨。虛擬因變量：若隸屬於某政黨，則 隸屬=1，不加入則為0。- 決策：退休。虛擬因變量：若已退休，則 退休=1，0表示未退休。
當虛擬因變量具有兩個以上的值時（例如隸屬於多個政黨），它就變成了多響應（或多項式、多叉分枝）模型。

**虛擬因變量模型**
虛擬因變量模型的分析有不同的方法。其中一種是通常的普通最小二乘法，在該情景下稱為linear probability model。另一種方法是假設存在一個不可觀察的連續潛變量Y*，如果Y* > 0，則觀察到的二分變量Y = 1，否則為0。這是邏輯斯諦和Probit model的基本理念。這些模型將在下面簡要討論。

**線性概率模型**
普通最小二乘法模型中，若因變量Y是二分的虛擬變量，取值0或1，則稱該模型為linear probability model（LPM）。假設考慮以下回歸：
 
其中。
  = 家庭收入
 若住房為家庭自有，0若住房非家庭自有
該模型稱為線性概率模型，是因為迴歸是線性的。給定Xi條件下的Yi的條件均值記作，可解釋為在相應Xi值下，事件發生的條件概率，即1=Pr(Yi = 1  Xi) 。在這個例子中，表示的是當家庭收入為Xi時，一個家庭自有住房的概率。
此時，使用普通最小二乘法的假設，可得到
 
LPM模型中存在一些固有問題：
# 迴歸線無法做到擬合良好，因此諸如R2之類的顯著性度量將不可靠。
# 使用LPM方法分析的模型有異方差干擾。
# 誤差項有非正態分佈。
# LPM可能會給出大於1或小於0的因變量的預測值。這將很難作出解釋，因為預測值是概率，必須介於0和1之間。
# LPM模型的變量之間可能存在非線性關係，在這種情況下，線性迴歸將無法準確擬合數據。

**線性概率模型的替代品**
圖 4：累積分佈函數
為了避免LPM的侷限性，我們需要一個模型，其中隨著解釋變量Xi的增加，1=Pi = E(Yi = 1  Xi)應恆在0到1之間。因此，自變量和因變量之間的關係必然是非線性的。
為此，可以使用累積分佈函數（CDF） 來估計虛擬因變量回歸。圖4顯示了一條S形曲線，它類似於隨機變量的累積分佈函數。該模型中，概率在0和1之間，並且已經捕獲了非線性關係。此時的問題就是如何選用累積分佈函數。
可以使用兩種替代性的CDF：邏輯斯諦CDF和正態CDF。邏輯斯諦CDF產生邏輯斯諦模型，正態CDF產生Probit model。

**邏輯斯諦模型**
LPM的缺點促使人們開發了一種改良的模型，稱為logit模型。該模型中，迴歸方程中誤差項的累積分佈是Logistic distribution。因為它是非線性的，其迴歸更接近現實。
使用最大似然法估計logit模型。該模型中， ，即在給定自變量的情況下，因變量取值為1的概率為：
 
然後模型以odds ratio的形式表示：邏輯斯諦迴歸（對數幾率迴歸）中建模的是幾率的自然對數，幾率定義為 。取幾率的自然對數，logit（Li）表示為
 
這種關係表明，Li與Xi呈線性關係，但概率與Xi呈非線性關係。

**概率單位模型**
概率單位（probit）模型是為彌補LPM的缺點而提出的另一個模型。Probit模型使用與logit模型相同的非線性方法。但是，它使用正態CDF而不是邏輯CDF。

**參見**
 Binary regression- 鄒檢驗- 假說檢定- 指示函數- 線性判別分析- 多重共線性

**參考文獻**

**延伸閱讀**
 pages=209–230- isbn=90-237-2919-6

**外部鏈接**
 archive-date=December 16, 2011- dead-url=no - archive-date=March 1, 2006