نموذج التأثيرات الثابتة

في الإحصاء ، يُعرف نموذج التأثيرات الثابتة بأنه نموذج إحصائي تكون فيه معلمات النموذج ثابتة أو غير عشوائية. وهذا يختلف عن نماذج التأثيرات العشوائية والنماذج المختلطة التي تكون فيها جميع معلمات النموذج أو بعضها متغيرات عشوائية. في العديد من التطبيقات، بما في ذلك الاقتصاد القياسي [ 1 ] والإحصاء الحيوي [ 2 ] [ 3 ] [ 4 ] [ 5 ] [ 6 ] ، يشير نموذج التأثيرات الثابتة إلى نموذج انحدار تكون فيه متوسطات المجموعات ثابتة (غير عشوائية)، على عكس نموذج التأثيرات العشوائية الذي تكون فيه متوسطات المجموعات عينة عشوائية من المجتمع الإحصائي. [ 7 ] [ 6 ] عمومًا، يمكن تجميع البيانات وفقًا لعدة عوامل مُلاحَظة. ويمكن نمذجة متوسطات المجموعات كتأثيرات ثابتة أو عشوائية لكل مجموعة. في نموذج التأثيرات الثابتة، يكون كل متوسط ​​مجموعة كمية ثابتة خاصة بتلك المجموعة.

في بيانات اللوحات التي تتضمن مشاهدات طولية لنفس الشخص، تمثل التأثيرات الثابتة المتوسطات الخاصة بكل شخص. في تحليل بيانات اللوحات، يُستخدم مصطلح مُقدِّر التأثيرات الثابتة (المعروف أيضًا باسم مُقدِّر التداخل ) للإشارة إلى مُقدِّر معاملات نموذج الانحدار الذي يتضمن تلك التأثيرات الثابتة (مُعامل تقاطع ثابت مع الزمن لكل شخص).

الوصف النوعي

تساعد هذه النماذج في التحكم في تحيز المتغيرات المحذوفة الناتج عن عدم التجانس غير الملحوظ عندما يكون هذا التباين ثابتًا بمرور الوقت. ويمكن إزالة هذا التباين من البيانات عن طريق التفاضل، على سبيل المثال عن طريق طرح متوسط ​​مستوى المجموعة بمرور الوقت، أو عن طريق أخذ الفرق الأول الذي من شأنه إزالة أي مكونات ثابتة بمرور الوقت من النموذج.

هناك افتراضان شائعان حول التأثير الفردي المحدد: افتراض التأثيرات العشوائية وافتراض التأثيرات الثابتة. يفترض افتراض التأثيرات العشوائية أن التأثيرات الفردية المحددة غير مرتبطة بالمتغيرات المستقلة. أما افتراض التأثيرات الثابتة فيفترض أن التأثيرات الفردية المحددة مرتبطة بالمتغيرات المستقلة. إذا تحقق افتراض التأثيرات العشوائية، يكون مقدِّر التأثيرات العشوائية أكثر كفاءة من مقدِّر التأثيرات الثابتة. أما إذا لم يتحقق هذا الافتراض، فإن مقدِّر التأثيرات العشوائية يكون غير متسق . يُستخدم اختبار دوربين-وو-هاوسمان غالبًا للتمييز بين نموذجي التأثيرات الثابتة والعشوائية. [ 8 ] [ 9 ]

النموذج الرسمي والافتراضات

ضع في اعتبارك نموذج التأثيرات غير الملحوظة الخطية لـشمال{\displaystyle N}الملاحظات وتي{\displaystyle T}الفترات الزمنية:

yأنات=Xأناتβ+αأنا+uأنات{\displaystyle y_{it}=X_{it}\mathbf {\beta } +\alpha _{i}+u_{it}}لت=1،...،تي{\displaystyle t=1,\dots ,T}وأنا=1،...،شمال{\displaystyle i=1,\dots ,N}

أين:

  • yأنات{\displaystyle y_{it}}المتغير التابع الذي تمت ملاحظته للفردأنا{\displaystyle i}في ذلك الوقتت{\displaystyle t}.
  • Xأنات{\displaystyle X_{it}}هو المتغير الزمني1×ك{\displaystyle 1\times k}(عدد المتغيرات المستقلة) متجه الانحدار.
  • β{\displaystyle \beta }هوك×1{\displaystyle k\times 1}مصفوفة المعاملات.
  • αأنا{\displaystyle \alpha _{i}}هو التأثير الفردي غير الملحوظ والثابت بمرور الوقت. على سبيل المثال، القدرة الفطرية للأفراد أو العوامل التاريخية والمؤسسية للدول.
  • uأنات{\displaystyle u_{it}}هو حد الخطأ .

على عكسXأنات{\displaystyle X_{it}}،αأنا{\displaystyle \alpha _{i}}لا يمكن ملاحظتها مباشرة.

بخلاف نموذج التأثيرات العشوائية حيث غير الملاحظαأنا{\displaystyle \alpha _{i}}مستقل عنXأنات{\displaystyle X_{it}}للجميعت=1،...،تي{\displaystyle t=1,...,T}يسمح نموذج التأثيرات الثابتة (FE)αأنا{\displaystyle \alpha _{i}}أن يكون مرتبطًا بمصفوفة الانحدارXأنات{\displaystyle X_{it}}. الاستقلالية التامة فيما يتعلق بحد الخطأ الخاصuأنات{\displaystyle u_{it}}لا يزال مطلوبًا.

التقدير الإحصائي

مقدر التأثيرات الثابتة

منذαأنا{\displaystyle \alpha _{i}}لا يمكن ملاحظته، ولا يمكن التحكم فيه بشكل مباشر. يقوم نموذج العناصر المحدودة بإلغاءαأنا{\displaystyle \alpha _{i}}عن طريق إزالة معنى المتغيرات باستخدام التحويل الداخلي :

yأنات-y¯أنا=(Xأنات-X¯أنا)β+(αأنا-α¯أنا)+(uأنات-u¯أنا)y¨أنات=X¨أناتβ+u¨أنات{\displaystyle y_{it}-{\overline {y}}_{i}=\left(X_{it}-{\overline {X}}_{i}\right)\beta +\left(\alpha _{i}-{\overline {\alpha }}_{i}\right)+\left(u_{it}-{\overline {u}}_{i}\right)\implies {\ddot {y}}_{it}={\ddot {X}}_{it}\beta +{\ddot {u}}_{it}}

أينy¯أنا=1تيت=1تيyأنات{\displaystyle {\overline {y}}_{i}={\frac {1}{T}}\sum \limits _{t=1}^{T}y_{it}}،X¯أنا=1تيت=1تيXأنات{\displaystyle {\overline {X}}_{i}={\frac {1}{T}}\sum \limits _{t=1}^{T}X_{it}}، وu¯أنا=1تيت=1تيuأنات{\displaystyle {\overline {u}}_{i}={\frac {1}{T}}\sum \limits _{t=1}^{T}u_{it}}.

منذαأنا{\displaystyle \alpha _{i}}ثابت،αأنا¯=αأنا{\displaystyle {\overline {\alpha _{i}}}=\alpha _{i}}وبالتالي يتم التخلص من التأثير. مقدر التأثيرات الثابتةβ^Fهـ{\displaystyle {\hat {\beta }}_{FE}}ثم يتم الحصول عليها من خلال انحدار المربعات الصغرى العادية لـy¨{\displaystyle {\ddot {y}}}علىX¨{\displaystyle {\ddot {X}}}.

توجد ثلاثة بدائل على الأقل للتحويل الداخلي مع اختلافات:

  • أحدها هو إضافة متغير وهمي لكل فردأنا>1{\displaystyle i>1}(مع استبعاد الفرد الأول بسبب الارتباط الخطي المتعدد ). هذا مكافئ عدديًا، وليس حسابيًا، لنموذج التأثيرات الثابتة، ولا ينجح إلا إذا كان مجموع عدد السلاسل وعدد المعلمات العامة أقل من عدد المشاهدات. [ 10 ] يُعدّ أسلوب المتغيرات الوهمية مُستهلكًا كبيرًا لذاكرة الحاسوب، ولا يُنصح به للمسائل الأكبر من سعة ذاكرة الوصول العشوائي المتاحة، وقدرة برنامج التجميع المُستخدم.
  • Second alternative is to use consecutive reiterations approach to local and global estimations.[11] This approach is very suitable for low memory systems on which it is much more computationally efficient than the dummy variable approach.
  • The third approach is a nested estimation whereby the local estimation for individual series is programmed in as a part of the model definition.[12] This approach is the most computationally and memory efficient, but it requires proficient programming skills and access to the model programming code; although, it can be programmed including in SAS.[13][14]

Finally, each of the above alternatives can be improved if the series-specific estimation is linear (within a nonlinear model), in which case the direct linear solution for individual series can be programmed in as part of the nonlinear model definition.[15]

First difference estimator

An alternative to the within transformation is the first difference transformation, which produces a different estimator. For t=2,,T{\displaystyle t=2,\dots ,T}:

yityi,t1=(XitXi,t1)β+(αiαi)+(uitui,t1)Δyit=ΔXitβ+Δuit.{\displaystyle y_{it}-y_{i,t-1}=\left(X_{it}-X_{i,t-1}\right)\beta +\left(\alpha _{i}-\alpha _{i}\right)+\left(u_{it}-u_{i,t-1}\right)\implies \Delta y_{it}=\Delta X_{it}\beta +\Delta u_{it}.}

The FD estimator β^FD{\displaystyle {\hat {\beta }}_{FD}} is then obtained by an OLS regression of Δyit{\displaystyle \Delta y_{it}} on ΔXit{\displaystyle \Delta X_{it}}.

When T=2{\displaystyle T=2}, the first difference and fixed effects estimators are numerically equivalent. For T>2{\displaystyle T>2}, they are not. If the error terms uit{\displaystyle u_{it}} are homoskedastic with no serial correlation, the fixed effects estimator is more efficient than the first difference estimator. If uit{\displaystyle u_{it}} follows a random walk, however, the first difference estimator is more efficient.[16]

Equality of fixed effects and first difference estimators when T=2

For the special two period case (T=2{\displaystyle T=2}), the fixed effects (FE) estimator and the first difference (FD) estimator are numerically equivalent. This is because the FE estimator effectively "doubles the data set" used in the FD estimator. To see this, establish that the fixed effects estimator is: FET=2=[(xi1x¯i)(xi1x¯i)+(xi2x¯i)(xi2x¯i)]1[(xi1x¯i)(yi1y¯i)+(xi2x¯i)(yi2y¯i)]{\displaystyle {FE}_{T=2}=\left[(x_{i1}-{\bar {x}}_{i})(x_{i1}-{\bar {x}}_{i})'+(x_{i2}-{\bar {x}}_{i})(x_{i2}-{\bar {x}}_{i})'\right]^{-1}\left[(x_{i1}-{\bar {x}}_{i})(y_{i1}-{\bar {y}}_{i})+(x_{i2}-{\bar {x}}_{i})(y_{i2}-{\bar {y}}_{i})\right]}

Since each (xi1x¯i){\displaystyle (x_{i1}-{\bar {x}}_{i})} can be re-written as (xi1xi1+xi22)=xi1xi22{\displaystyle (x_{i1}-{\dfrac {x_{i1}+x_{i2}}{2}})={\dfrac {x_{i1}-x_{i2}}{2}}}, we'll re-write the line as:

FET=2=[i=1Nxi1xi22xi1xi22+xi2xi12xi2xi12]1[i=1Nxi1xi22yi1yi22+xi2xi12yi2yi12]{\displaystyle {FE}_{T=2}=\left[\sum _{i=1}^{N}{\dfrac {x_{i1}-x_{i2}}{2}}{\dfrac {x_{i1}-x_{i2}}{2}}'+{\dfrac {x_{i2}-x_{i1}}{2}}{\dfrac {x_{i2}-x_{i1}}{2}}'\right]^{-1}\left[\sum _{i=1}^{N}{\dfrac {x_{i1}-x_{i2}}{2}}{\dfrac {y_{i1}-y_{i2}}{2}}+{\dfrac {x_{i2}-x_{i1}}{2}}{\dfrac {y_{i2}-y_{i1}}{2}}\right]}

=[i=1N2xi2xi12xi2xi12]1[i=1N2xi2xi12yi2yi12]{\displaystyle =\left[\sum _{i=1}^{N}2{\dfrac {x_{i2}-x_{i1}}{2}}{\dfrac {x_{i2}-x_{i1}}{2}}'\right]^{-1}\left[\sum _{i=1}^{N}2{\dfrac {x_{i2}-x_{i1}}{2}}{\dfrac {y_{i2}-y_{i1}}{2}}\right]}
=2[i=1N(xi2xi1)(xi2xi1)]1[i=1N12(xi2xi1)(yi2yi1)]{\displaystyle =2\left[\sum _{i=1}^{N}(x_{i2}-x_{i1})(x_{i2}-x_{i1})'\right]^{-1}\left[\sum _{i=1}^{N}{\frac {1}{2}}(x_{i2}-x_{i1})(y_{i2}-y_{i1})\right]}
=[i=1N(xi2xi1)(xi2xi1)]1i=1N(xi2xi1)(yi2yi1)=FDT=2{\displaystyle =\left[\sum _{i=1}^{N}(x_{i2}-x_{i1})(x_{i2}-x_{i1})'\right]^{-1}\sum _{i=1}^{N}(x_{i2}-x_{i1})(y_{i2}-y_{i1})={FD}_{T=2}}

Chamberlain method

Gary Chamberlain's method, a generalization of the within estimator, replaces αi{\displaystyle \alpha _{i}} with its linear projection onto the explanatory variables. Writing the linear projection as:

αi=λ0+Xi1λ1+Xi2λ2++XiTλT+ei{\displaystyle \alpha _{i}=\lambda _{0}+X_{i1}\lambda _{1}+X_{i2}\lambda _{2}+\dots +X_{iT}\lambda _{T}+e_{i}}

this results in the following equation:

yit=λ0+Xi1λ1+Xi2λ2++Xit(λt+β)++XiTλT+ei+uit{\displaystyle y_{it}=\lambda _{0}+X_{i1}\lambda _{1}+X_{i2}\lambda _{2}+\dots +X_{it}(\lambda _{t}+\mathbf {\beta } )+\dots +X_{iT}\lambda _{T}+e_{i}+u_{it}}

which can be estimated by minimum distance estimation.[17]

Hausman–Taylor method

Need to have more than one time-variant regressor (X{\displaystyle X}) and time-invariant regressor (Z{\displaystyle Z}) and at least one X{\displaystyle X} and one Z{\displaystyle Z} that are uncorrelated with αi{\displaystyle \alpha _{i}}.

Partition the X{\displaystyle X} and Z{\displaystyle Z} variables such that X=[X1itTN×K1X2itTN×K2]Z=[Z1itTN×G1Z2itTN×G2]{\displaystyle {\begin{array}{c}X=[{\underset {TN\times K1}{X_{1it}}}\vdots {\underset {TN\times K2}{X_{2it}}}]\\Z=[{\underset {TN\times G1}{Z_{1it}}}\vdots {\underset {TN\times G2}{Z_{2it}}}]\end{array}}} where X1{\displaystyle X_{1}} and Z1{\displaystyle Z_{1}} are uncorrelated with αi{\displaystyle \alpha _{i}}. Need K1>G2{\displaystyle K1>G2}.

Estimating γ{\displaystyle \gamma } via OLS on di^=Ziγ+φit{\displaystyle {\widehat {di}}=Z_{i}\gamma +\varphi _{it}} using X1{\displaystyle X_{1}} and Z1{\displaystyle Z_{1}}حيث أن الأدوات تعطي تقديرًا متسقًا.

التعميم مع عدم اليقين في المدخلات

عندما يكون هناك عدم يقين في المدخلات لـy{\displaystyle y}بيانات،دلتاy{\displaystyle \delta y}ثم الـχ2{\displaystyle \chi ^{2}}ينبغي تقليل قيمة القيمة، بدلاً من مجموع مربعات البواقي. [ 18 ] ويمكن تحقيق ذلك مباشرة من خلال قواعد الاستبدال:

yأناتدلتاyأنات=βXأناتدلتاyأنات+αأنا1دلتاyأنات+uأناتدلتاyأنات{\displaystyle {\frac {y_{it}}{\delta y_{it}}}=\mathbf {\beta } {\frac {X_{it}}{\delta y_{it}}}+\alpha _{i}{\frac {1}{\delta y_{it}}}+{\frac {u_{it}}{\delta y_{it}}}}،

ثم القيم والانحرافات المعيارية لـβ{\displaystyle \mathbf {\beta } }وαأنا{\displaystyle \alpha _{i}}يمكن تحديدها من خلال تحليل المربعات الصغرى العادية الكلاسيكية ومصفوفة التباين والتباين المشترك .

يُستخدم لاختبار الاتساق

قد تكون تقديرات التأثيرات العشوائية غير متسقة أحيانًا في حالة السلاسل الزمنية الطويلة، إذا كانت التأثيرات العشوائية غير محددة بشكل صحيح (أي أن النموذج المختار للتأثيرات العشوائية غير مناسب). مع ذلك، قد يظل نموذج التأثيرات الثابتة متسقًا في بعض الحالات. على سبيل المثال، إذا لم تكن السلسلة الزمنية قيد الدراسة مستقرة، فقد لا تكون نماذج التأثيرات العشوائية التي تفترض الاستقرار متسقة في حالة السلاسل الزمنية الطويلة. ومن الأمثلة على ذلك وجود اتجاه تصاعدي في السلسلة الزمنية. فمع ازدياد طول السلسلة، يُعيد النموذج تقييم متوسط ​​الفترات السابقة بالزيادة، مما يُعطي تنبؤات متحيزة بشكل متزايد للمعاملات. مع ذلك، لا يجمع النموذج ذو التأثيرات الزمنية الثابتة المعلومات عبر الزمن، وبالتالي لن تتأثر التقديرات السابقة.

في مثل هذه الحالات، حيث يُعرف أن نموذج التأثيرات الثابتة متسق، يمكن استخدام اختبار دوربين-وو-هاوسمان لاختبار ما إذا كان نموذج التأثيرات العشوائية المختار متسقًا.ح0{\displaystyle H_{0}}صحيح، كلاهماβ^Rهـ{\displaystyle {\widehat {\beta }}_{RE}}وβ^Fهـ{\displaystyle {\widehat {\beta }}_{FE}}متسقة، ولكن فقطβ^Rهـ{\displaystyle {\widehat {\beta }}_{RE}}فعال. إذاحأ{\displaystyle H_{a}}صحيح أن اتساقβ^Rهـ{\displaystyle {\widehat {\beta }}_{RE}}لا يمكن ضمان ذلك.

انظر أيضاً

ملحوظات

  1. غرين، دبليو إتش، 2011. التحليل الاقتصادي القياسي ، الطبعة السابعة، برنتيس هول
  2. ^ ديجل ، بيتر ج. هيجرتي، باتريك؛ ليانغ، كونغ يي؛ زيجر، سكوت ل. (2002). تحليل البيانات الطولية (  الطبعة الثانية). مطبعة جامعة أكسفورد. ص 169 – 171. ISBN  0-19-852484-6.
  3. فيتزموريس، غاريت م.؛ ليرد، نان م.؛ وير، جيمس هـ. (2004). التحليل الطولي التطبيقي . هوبوكين: جون وايلي وأولاده. ص 326-328 . ISBN  0-471-21487-6.
  4. ليرد، نان م.؛ وير، جيمس هـ. (1982). "نماذج التأثيرات العشوائية للبيانات الطولية". القياسات الحيوية . 38 (4): 963-974 . doi : 10.2307/2529876 . JSTOR 2529876 . 
  5. غاردينر، جوزيف سي؛ لو، زيهوي؛ رومان، لي آن (2009). "التأثيرات الثابتة، والتأثيرات العشوائية، ونموذج المعادلات التقديرية المعممة: ما هي الاختلافات؟". الإحصاء في الطب . 28 (2): 221-239 . doi : 10.1002/sim.3478 . PMID 19012297. S2CID 16277040 .  
  6. 1 2 غوميز، ديلان جي إي (20 يناير 2022). "هل أستخدم التأثيرات الثابتة أم التأثيرات العشوائية عندما يكون لدي أقل من خمسة مستويات لعامل التجميع في نموذج التأثيرات المختلطة؟" . PeerJ . 10 e12794 . doi : 10.7717/peerj.12794 . PMC 8784019. PMID 35116198 .  
  7. رامزي، ف.، وشافر، د.، 2002. المحقق الإحصائي: دورة في أساليب تحليل البيانات ، الطبعة الثانية. دار نشر دوكسبوري.
  8. كاميرون، أ. كولين؛ تريفيدي، برافين ك. (2005). الاقتصاد القياسي الجزئي: الأساليب والتطبيقات . مطبعة جامعة كامبريدج. ص 717-719 . ISBN  978-0-521-84805-3.
  9. نيرلوف، مارك (2005). مقالات في الاقتصاد القياسي لبيانات اللوحات . مطبعة جامعة كامبريدج. ص 36-39 . ISBN  978-0-521-02246-0.
  10. غارسيا، أوسكار. (1983). "نموذج معادلة تفاضلية عشوائية لنمو ارتفاع الأشجار في الغابات". القياسات الحيوية . 39 (4): 1059-1072 . doi : 10.2307/2531339 . JSTOR 2531339 . 
  11. تايت، ديفيد؛ سيزوسكي، كريس جيه؛ بيلا، إيمري إي. (1986). "ديناميكيات غابات الصنوبر اللودج بول". المجلة الكندية لأبحاث الغابات 18 ( 10): 1255-1260 . doi : 10.1139/x88-193 .
  12. ستروب، مايك؛ سيزوسكي، كريس جيه. (2006). "خصائص ثبات العمر الأساسي لتقنيتين لتقدير معلمات نماذج مؤشر الموقع". علوم الغابات . 52 (2): 182-186 . doi : 10.1093/forestscience/52.2.182 .
  13. ستروب، مايك؛ سيزوسكي، كريس جيه. (2003). بوركهارت، إتش إيه (محرر). ملاءمة معلمات مؤشر الموقع العالمي عندما يُعامل مؤشر موقع القطعة أو الشجرة كمعلمة إزعاج محلية . وقائع ندوة الإحصاء وتكنولوجيا المعلومات في مجال الغابات؛ 8-12 سبتمبر 2002؛ بلاكسبيرغ، فيرجينيا: معهد فيرجينيا للفنون التطبيقية وجامعة ولاية فيرجينيا. ص 97-107 . 
  14. Cieszewski, Chris J.; Harrison, Mike; Martin, Stacey W. (2000). "Practical methods for estimating non-biased parameters in self-referencing growth and yield models"(PDF). PMRC Technical Report. 2000 (7): 12. Archived from the original(PDF) on 2016-03-04. Retrieved 2015-12-24.
  15. Schnute, Jon; McKinnell, Skip (1984). "A biologically meaningful approach to response surface analysis". Can. J. Fish. Aquat. Sci. 41 (6): 936–953. doi:10.1139/f84-108.
  16. Wooldridge, Jeffrey M. (2001). Econometric Analysis of Cross Section and Panel Data. MIT Press. pp. 279–291. ISBN 978-0-262-23219-7.
  17. Chamberlain, Gary (1984). Chapter 22 Panel data. Handbook of Econometrics. Vol. 2. pp. 1247–1318. doi:10.1016/S1573-4412(84)02014-6. ISBN 978-0-444-86186-3. ISSN 1573-4412.
  18. Ren, Bin; Dong, Ruobing; Esposito, Thomas M.; Pueyo, Laurent; Debes, John H.; Poteet, Charles A.; Choquet, Élodie; Benisty, Myriam; Chiang, Eugene; Grady, Carol A.; Hines, Dean C.; Schneider, Glenn; Soummer, Rémi (2018). "A Decade of MWC 758 Disk Images: Where Are the Spiral-Arm-Driving Planets?". The Astrophysical Journal Letters. 857 (1): L9. arXiv:1803.06776. Bibcode:2018ApJ...857L...9R. doi:10.3847/2041-8213/aab7f5. S2CID 59427417.

References

  • Christensen, Ronald (2002). Plane Answers to Complex Questions: The Theory of Linear Models (Third ed.). New York: Springer. ISBN 0-387-95361-2.
  • Gujarati, Damodar N.; Porter, Dawn C. (2009). "Panel Data Regression Models". Basic Econometrics (Fifth international ed.). Boston: McGraw-Hill. pp. 591–616. ISBN 978-007-127625-2.
  • Hsiao, Cheng (2003). "Fixed-effects models". Analysis of Panel Data (2nd ed.). New York: Cambridge University Press. pp. 95–103. ISBN 0-521-52271-4.
  • وولدريدج، جيفري م. (2013). "تقدير التأثيرات الثابتة". الاقتصاد القياسي التمهيدي: منهج حديث (  الطبعة الدولية الخامسة). ماسون، أوهايو: ساوث ويسترن. ص 466-474 . ISBN  978-1-111-53439-4.