潛在變化分數模型與進階成長模型

20.1 把變化當成一個潛在變項

潛在變化分數模型由 McArdle 與同僚以一個長期的研究計畫發展出來 (McArdle, 2009; McArdle & Hamagami, 2001)，它建在一套讓初學者感到陌生的鷹架 (scaffolding) 上，因為其中幾乎每一條路徑都被固定為一。圖 20.1 把它畫了出來。每一個觀察分數的背後坐著一個潛在真分數 (latent true score)，已經去掉了測量誤差 (measurement error)；各個真分數之間以自我迴歸 (autoregression) 跨時間相連，而係數被固定為恰好一，因此每一個真分數原封不動地把前一個真分數整個帶過來。流進每一個真分數的是一個潛在變化分數 (latent change score)，路徑同樣固定為一，於是每一個時點的真分數字面上就是「前一個真分數加上變化」，這是一條恆等式 (identity)，不是一條迴歸式。模型化的工作發生在變化分數上。一個常數變化 (constant change) 因素以固定為一的負荷量 (loading) 負荷在每一個變化分數上，在每一步貢獻相同的增量 (increment)，而它的平均數就是加法變化參數。一條比例變化 (proportional change) 路徑則從每一個真分數走到下一個變化分數，係數在每個區間都相同，使每一次變化中有一部分是當下水準的一個比例。那些被固定的一並不是任意的：它們強制執行「變化就是相鄰兩個真分數之差」這條恆等關係，好讓那些自由參數，也就是常數變化的平均數與變異數以及比例係數，成為真正要推論 (inference) 的對象。

註：潛在真分數 \(\ell_t\) 以固定為一的負荷量撐起觀察分數 \(y_t\)；各真分數之間的自我迴歸固定為一；每一個潛在變化分數 \(\Delta_t\) 以固定為一的路徑進入它的真分數，強制執行 \(\ell_t = \ell_{t-1} + \Delta_t\) 這條恆等關係。自由參數都在變化上：常數變化因素 \(g\)（負荷量為一，平均數即加法變化）與從前一個水準過來的比例係數 \(\beta\)（紅）。正是那些被固定的一，才使變化成為一個被模型化的潛在變項。

把變化規則寫成 \(\Delta_t = g + \beta\,\ell_{t-1}\) 再把它展開，就看得出這套裝置為何值得費這番工夫：少少幾個參數生得出一整族的軌跡形狀。當只有常數變化時，軌跡是線性的，每一步都是同樣的增量。當只有比例變化時，軌跡是幾何式 (geometric) 的趨近，係數為負則朝零衰減 (decay)，為正則爆炸性成長。當兩者都在時，也就是雙變化 (dual change) 模型，軌跡以逐漸減速的方式朝 \(-g/\beta\) 這個漸近值 (asymptote) 上升或下降，正是第 14 章當作非線性混合模型 (nonlinear mixed model) 直接配適的那條負指數 (negative exponential) 形狀。圖 20.2 畫出這套分類，一條方程式依比例係數的正負與大小，生出常數、衰減、趨近漸近值與加速 (acceleration) 四種軌跡；表 20.1 則把變化分數的參數化對應到第 14 章的非線性曲線參數化。兩者描述的是同一批曲線；變化分數模型只是透過「生成它們的規則」抵達，而不是透過「它們所描出的函數」抵達。

註：由變化規則 \(\Delta_t = g + \beta\,\ell_{t-1}\) 生成的軌跡。比例係數為零給出線性變化；係數為負而沒有常數項給出衰減；正的常數項配上負的係數給出朝漸近值的趨近；係數為正則給出加速。生出這一族形狀的是那條生成規則，不是一條被配適出來的函數。

20.2 單變量變化分數模型的實作

把整座設定階梯 (specification ladder) 配適到一份已知為「減速 (decelerating) 上升」的資料上，選擇就變得具體。在本章的示範資料上，那是一個由雙變化規則生成的模擬歷程，只有常數變化的模型配適不良，因為它硬把一條直線穿過一條彎曲的軌跡；只有比例變化的模型配適得更差，因為少了固定的增量，它捕捉不住那段持續的上升；而雙變化模型配適良好，並且還原了生成參數，加法變化接近它的真值，比例係數也接近它為負的真值。圖 20.3 把三條隱含軌跡 (implied trajectory) 與觀察平均數放在一起：只有雙變化模型跟得上那段減速。比較的依據是訊息準則 (information criteria)，而在模型互相巢套 (nested) 之處還加上卡方差異 (chi-square difference)，雙變化模型以很大的差距勝出。表 20.2 是參數的名詞表，表 20.3 則把參數的組合對應到它們所隱含的動態。有兩種估計上的病症會一再出現，它們值得的是冷靜的回應而不是驚慌：當資料接近線性時，比例係數可能漂到邊界 (boundary) 上；而變化分數的變異數可能變成略為負值，那又是第 13 章那個邊界問題 (boundary problem)，回應的方式是簡化模型，或如第 17 章所示改用貝氏估計 (Bayesian estimation)。

註：一個減速歷程的觀察平均數（灰）對三種潛在變化分數設定的隱含軌跡。只有常數變化與只有比例變化的模型都抓不到那個形狀；雙變化模型把固定的增量與朝漸近值拉的比例項結合起來，跟得上它。訊息準則決定性地偏好雙變化模型。

比例係數帶著一項在此首度登場、而後主宰整個連續時間章節的告誡。它是以觀測區間為單位 (per observation interval) 定義的：係數為 \(-0.2\) 的意思是，在相鄰兩個被測量到的時點之間，與漸近值之間的差距縮小了五分之一；而如果這兩個時點的間隔拉成兩倍，係數就會不同，因為在比較長的區間裡縮小的幅度比較大。圖 20.4 把同一個歷程分別以「每一波」與「每隔一波」重新配適，展示了這件事：比例係數變動得相當可觀，由大約 \(-0.2\) 變成大約 \(-0.4\)，不是因為歷程變了，而是因為區間變了。這種區間依賴 (interval dependence) 意味著，來自不同測量時程的研究，它們的比例係數彼此不可比較；而這也是支持第 27 章那個連續時間重新參數化最強的一項論據，因為在那裡，速率是歷程本身的性質，而不是取樣方式的性質。

註：同一個歷程分別以每一個時點與每隔一個時點取樣時，估計出來的比例係數。係數的絕對值幾乎翻了一倍，不是因為動態變了，而是因為區間變了。離散的比例係數是「每個區間」的量，不能拿來跨測量時程比較；只有連續時間的速率可以。

20.3 雙變量變化分數模型與耦合

引起最多興趣、也惹出最多麻煩的那項擴充，是把兩套變化分數系統耦合起來，讓其中一個歷程的水準驅動另一個歷程的變化。圖 20.5 勾勒出這個架構：兩套雙變化系統並排，中間有耦合路徑 (coupling paths) 把每一個歷程在某一個時點的真分數帶進另一個歷程在下一個時點的變化。從第一個歷程到第二個歷程的變化有一條正的耦合，意思是在第一個歷程上分數較高，可以預測第二個歷程隨後的上升，這就是「第一個歷程領先第二個歷程」這項主張的結構性表達。本章的示範資料是以一組不對稱的耦合生成的，其中一個歷程驅動另一個而反向不成立，而雙變量模型還原的正是這個結構：在生成的方向上有一條可觀的耦合，在另一個方向上則微不足道，正確地辨識出領先落後的結構。此時的誘惑是把這樣的結果讀成「因果問題已經解決」，而這份誘惑必須被抵抗。

註：兩套雙變化系統，分別屬於歷程 \(x\)（上）與歷程 \(y\)（下），各自有自己的比例式自我回饋。耦合路徑把一個歷程在前一個時點的水準帶進另一個歷程的變化：\(\gamma_{yx}\)（紅）是 \(x\) 的水準對 \(y\) 的變化的效果，\(\gamma_{xy}\)（紫）則是反方向。兩條耦合之間的不對稱，就是一項領先落後主張的結構性表達。

解讀一個耦合參數需要一套逐條列舉的紀律，因為耦合文獻一再過度宣稱。第一，耦合是以整個系統為條件的，因此加進或拿掉另一條路徑時，它的數值甚至正負都可能改變，這使得「逐一拿掉一條路徑」的敏感度分析 (sensitivity analysis) 成為必要。第二，它有區間依賴，繼承了比例係數的問題。第三，除非模型把真分數與誤差分開，否則耦合會把測量誤差吸收 (absorb) 進去，這是採用二階版本的一項論據。第四，一個耦合的意義取決於與它並排的自我回饋 (self-feedback) 係數，因為同一個水準效果，放在自我調節 (self-regulating) 的歷程旁與放在自我放大 (self-amplifying) 的歷程旁，讀起來並不一樣。第五，某些雙變量變化分數模型在統計上等價於加了限制的交叉延宕追蹤模型 (cross-lagged panel model)，因此一項耦合結果可能有一個等價模型 (equivalent model) 的影子，一個不同但配適得一樣好的設定會用另一套說法來描述它，這正是第 21 章要處理的統整工作。本書對此的標準是一整套敏感度檢驗 (sensitivity battery)：在一系列站得住腳的設定之下重新配適那個耦合，報告估計值的穩定程度，而不是報告單一個點估計 (point estimate)。圖 20.6 對本章的耦合跑了這套檢驗，估計值令人放心地穩定；表 20.4 則是那份檢核表。通過這套檢驗的耦合，才配得上一個謹慎的領先落後讀法；沒有通過的，就只是某個特定設定造出來的假象。

註：在幾種替代設定之下估計出來的耦合，對照它的真值。估計值是穩定的，而這正是一個謹慎的領先落後讀法所倚靠的東西。一個在這些設定之間大幅移動的耦合，是所選模型造出來的假象而不是一項發現，而只報告單一個點估計會把這個風險藏起來。

耦合系統是一組差分方程，而這樣讀它就預告了第六部。那一對變化規則在兩個水準的每一種組合上，都定義出各自預期的一步變化量，把它們畫成箭頭就得到一個向量場 (vector field)，也就是一個動力系統 (dynamical system) 的流 (flow) 在離散時間中的影像。圖 20.7 為配適出來的耦合系統畫出這個場：箭頭顯示這個歷程不論從哪裡出發都朝一個平衡點 (equilibrium) 流動，而流的形狀，是螺旋、直接收斂還是發散，就是那些參數所編碼的定性動態 (qualitative dynamics)。這正是第六部的向量自我迴歸模型 (vector autoregressive model) 與連續時間模型要形式化的那幅圖，而在此處把它畫出來，也就坐實了變化分數模型的橋樑地位。

註：每一支箭頭都是配適出來的耦合模型在該組水準之下所預期的一步變化量。這個場顯示系統朝它的平衡點流動。把變化分數模型做這樣的動力系統式解讀，就是通往第六部那些向量自我迴歸模型與連續時間模型的橋樑。

20.4 二階成長模型

第 18 章之所以要確立測量恆等性 (measurement invariance)，就是為了讓那張支票在這裡兌現。二階成長模型 (second-order growth model) 不是把成長結構放在觀察分數上，也不是放在複合分數 (composite score) 上，而是放在由多個指標所測量的潛在構念上：一個帶恆等負荷量與恆等截距的一階測量模型 (first-order measurement model)，在每個時點定義出一個因素，而一個二階的成長模型再把截距 (intercept) 與斜率 (slope) 配適到那些因素上。這樣的建構需要整條恆等性鏈條，其中截距恆等性 (scalar invariance) 是因素平均數可比較的必要條件，因此也是「因素平均數的成長有意義」的必要條件；表 20.5 說明每一級各自許可了什麼。回報有兩重：成長被配適在一個無誤差的構念上，因此測量誤差不再把它衰減 (attenuate)；而且恆等性是可以被檢定的，不是被假定的。圖 20.8 在第 18 章那份投入程度資料上呈現這份回報，該資料的第二個指標 (indicator) 有漂移 (drift)：二階模型把成長配適在部分恆等的潛在構念上，還原了真實的斜率，而複合分數默默假定了嚴格恆等性 (strict invariance)，於是高估了它，這正是第 18 章那個後果模擬在真實資料上的回聲。把成長配適到複合分數上，就是把它配適到一份意義可能已經位移的測量上；把成長配適到一個恆等的因素上，才是把它配適到那個構念上。二階模型有兩種多變量的形式，回答的是不同的問題：因素曲線模型 (curve-of-factors model) 把一個成長歷程配適到「每一波都被測量到的一個因素」上，而曲線因素模型 (factor-of-curves model) 則把一個高階因素 (higher-order factor) 配適到「好幾個各自獨立的成長歷程」上；兩者之間的選擇取決於那些構念在實質上的結構。

註：投入程度構念的平均斜率，分別由一個建在恆等潛在因素上的二階成長模型與一個建在複合分數上的成長模型估計出來，並與真實斜率對照。二階模型還原了真值；複合分數默默假定了那個會漂移的題目所違反的嚴格恆等性，因而高估了成長。第 18 章那道測量檢查，對改變速率的估計有實實在在的後果。

20.5 結構化的非線性成長

在變化分數模型與二階模型之外，還有一族結構化潛在曲線 (structured latent curve) 模型，它們透過「由某個函數的形式推導出來的受限負荷量」，強加一個特定的非線性目標函數，可以是指數 (exponential)、羅吉斯 (logistic) 或樣條 (spline) (Blozis, 2004)。這些模型兼有「一個具名函數形式所帶來的解讀精確度」與「潛在變項架構的彈性」，位置就落在第 19 章那個不強加任何形狀的放開負荷量 (freed-loading) 潛在基底 (latent basis)、與第 14 章那個強加一個僵硬形狀的完全參數化非線性混合模型之間。本章把它們當成地圖而不是住所：當理論指明了一條具有可解讀參數的特定曲線時，用結構化潛在曲線或非線性混合模型去配適它；當形狀未知時，用潛在基底去估計它；而當個別不同的測量時點 (individually varying times)、或結構方程架構處理起來不順手的多層次巢套佔了主導時，第 14 章與第 30 章的混合模型與加法模型 (additive model) 才是更好的歸宿。重點在於，這些架構構成一條關於「改變的形狀」的假設連續體，從什麼都不假設一直到完全指定一個函數，而分析者要挑的是理論與資料所擔保得起的那個位置。

20.6 在 R 中執行變化分數模型

超過幾個波次之後，手寫潛在變化分數的語法 (syntax) 就容易出錯，因此隨書的分析提供了一個產生器，它由變項名稱與波次數目生出 lavaan 語法，是一件讀者可以重複使用並加以擴充的配套產物。

library(lavaan)
# lcs_syntax(v, W, dual, prop) 產生單變量模型（見 ch20_analysis_V01.R）
# dual = 雙變化；const = 只有常數
dual  <- sem(lcs_syntax("x", W = 5, dual = TRUE,  prop = TRUE),  data = d)
const <- sem(lcs_syntax("x", W = 5, dual = TRUE,  prop = FALSE), data = d)
anova(const, dual)                      # 比較兩種設定

雙變量的產生器把兩套單變量系統合起來再加上耦合路徑，而二階成長模型則把截距因素與斜率因素放到第 18 章那些恆等的波次因素上。

# 雙變量耦合：一個歷程的水準預測另一個歷程的變化
biv <- sem(lcs_bivariate(W = 5, coupling = "both"), data = d)
# 二階成長：成長因素建在潛在構念上（測量恆等）
# I =~ 1*f1 + 1*f2 + 1*f3 + 1*f4 ; S =~ 0*f1 + 1*f2 + 2*f3 + 3*f4

完整的分析，包括設定階梯、區間依賴的示範、帶敏感度檢驗的雙變量耦合、隱含的向量場，以及二階成長與複合分數的比較，都寫在隨書腳本 ch20_analysis_V01.R 之中。本章的各張圖形由中文版的 ch20_figures_zh_V01.R 繪製，而所用的那份資料，則由 gen_lcs_V01.R 產生。

20.7 常見的誤解

關於變化分數模型，有幾種信念會把人帶偏。第一種是認為變化分數模型不過是多繞幾道手續的成長模型；它提出的是一項不同的科學主張，講的是生成變化的規則而不是軌跡的形狀，而且用「規則配得上、形狀配不上」的資料就能把兩者區分開。第二種是認為負的比例係數是迴歸趨中 (regression to the mean) 造成的假象；它可以是測量誤差造成的假象，而二階版本把那個誤差移除掉，但它同樣可以是貨真價實的自我調節，兩者要靠模型化來分辨，不是靠假定。第三種是認為雙變量變化分數模型解決了因果方向；敏感度檢驗之所以存在，正是因為它沒有解決，而一個耦合只有在各種設定之下都穩定時才可信。第四種是認為二階模型是心理計量純粹主義者的玩意；關於偏誤 (bias) 與檢定力 (power) 的證據說，凡是要估計一個潛在構念之變化的人，都該在意它。第五種是認為比例係數可以跨研究比較；它有區間依賴，因此不能，除非先做連續時間的重新參數化。

接下來要去哪裡

本章把「描述變化」轉了個彎，走向「模型化變化的動態」，而這個彎道通往本書其餘的部分。第 21 章證明潛在變化分數模型、交叉延宕追蹤模型與隨機截距 (random intercept) 的追蹤模型是同一個統整 (unified) 家族的成員，解決了耦合那一節提出的等價模型影子問題，對於在此處初次遇到耦合的讀者，那是最自然的下一步。第 27 章解決區間依賴的問題，估計一個屬於歷程本身而不屬於取樣時程的連續時間速率，並且把本章那個差分方程式的讀法落到實處。耦合系統的向量場圖，是第 24 章與第 25 章那些向量自我迴歸模型與動態結構方程模型 (dynamic structural equation model, DSEM) 的離散預告，那些模型在密集測量之中模型化的是個體之內的同一種流。二階成長模型完成了第 18 章開始的那條測量優先的鏈條，而第 34 章的發展連鎖 (developmental cascade) 應用，則把雙變量耦合這套做法投到「一個發展領域如何把另一個領域拉著走」這種實質問題上。

習題

1. 20.1 推導那條軌跡。對一個給定常數變化參數與比例參數的四波次雙變化模型，把遞迴關係式展開算出隱含的平均數，並拿去對照一個 lavaan 的配適結果加以驗證。
2. 20.2 把整座階梯配一遍。在所提供的資料上配適只有常數、只有比例與雙變化三種模型，比較它們，並寫出一段關於變化規則、站得住腳的結論。
3. 20.3 區間依賴實作。模擬一個雙變化歷程，分別在每一個時點與每隔一個時點上重新配適，示範比例係數會改變而軌跡不會。
4. 20.4 耦合的敏感度。在雙變量資料上估計那個耦合，跑完整套敏感度檢驗，然後寫出這套檢驗所支持的、帶保留條件的領先落後結論。
5. 20.5 二階成長。由一個恆等的測量模型建出二階成長模型，把它的斜率推論與建在複合分數上的成長模型相比，並解讀兩者的差異。
6. 20.6 使用產生器。讀懂隨書的語法產生器，用它配適一個五波次的模型，並為它加上一個新的選項。