الانحدار الذاتي المتكامل المتوسط المتحرك - أريما تعريف المتوسط المتحرك المتكامل الانتصاري - أريما نموذج تحليل إحصائي يستخدم بيانات السلاسل الزمنية للتنبؤ بالاتجاهات المستقبلية. وهو شكل من أشكال تحليل الانحدار الذي يسعى للتنبؤ بالتحركات المستقبلية على طول المشي العشوائي الذي يبدو من قبل الأسهم والسوق المالية من خلال دراسة الاختلافات بين القيم في سلسلة بدلا من استخدام قيم البيانات الفعلية. ويشار إلى التأخر في سلسلة الاختلاف باسم الانحدار الذاتي، ويشار إلى التأخر في البيانات المتوقعة كمتوسط متحرك. بريكينغ دون الانحدار التلقائي المتوسط المتحرك المتكامل - أريما يشار إلى هذا النوع من النماذج عموما باسم أريما (p، d، q)، مع الأعداد الصحيحة التي تشير إلى الانحدار الذاتي. متكاملة ومتحركة أجزاء من مجموعة البيانات، على التوالي. أريما النمذجة يمكن أن تأخذ في الاعتبار الاتجاهات والموسمية. (p، q) نماذج لتحليل السلاسل الزمنية - الجزء 3 هذا هو المنصب الثالث والأخير في سلسلة مصغرة على المتوسط المتحرك الانحدار الذاتي (p، q) نماذج لتحليل سلسلة زمنية - الجزء 3 أرما) لتحليل السلاسل الزمنية. قدمنا نماذج الانحدار الذاتي ونماذج المتوسط المتحرك في المقالات السابقة. الآن حان الوقت للجمع بينهما لإنتاج نموذج أكثر تطورا. في نهاية المطاف هذا سوف يقودنا إلى نماذج أريما و غارتش التي من شأنها أن تسمح لنا للتنبؤ عائدات الأصول وتوقع التقلبات. وستشكل هذه النماذج أساس إشارات التداول وتقنيات إدارة المخاطر. إذا كنت قد قرأت الجزء 1 والجزء 2 كنت قد رأيت أننا نميل إلى اتباع نمط لتحليلنا من نموذج سلسلة زمنية. سوء تكرار ذلك باختصار هنا: المبررات - لماذا نحن مهتمون في هذا النموذج معين تعريف - تعريف رياضي للحد من الغموض. كوريلوغرام - رسم عينة الرسم البياني لتصور سلوك النماذج. المحاكاة والمناسب - تركيب نموذج للمحاكاة، من أجل ضمان فهمنا النموذج بشكل صحيح. البيانات المالية الحقيقية - تطبيق نموذج لأسعار الأصول التاريخية الحقيقية. التنبؤ - توقعات القيم اللاحقة لبناء إشارات التداول أو الفلاتر. من أجل متابعة هذه المقالة فإنه من المستحسن أن نلقي نظرة على المواد السابقة على تحليل السلاسل الزمنية. ويمكن العثور عليها جميعا هنا. معيار معلومات بايزي في الجزء 1 من هذه المقالة سلسلة نظرنا في معيار المعلومات أكايك (إيك) كوسيلة لمساعدتنا على الاختيار بين أفضل نماذج أفضل سلسلة زمنية. وهناك أداة وثيقة الصلة هي معيار معلومات بايزي (بيك). أساسا لها سلوك مماثل ل إيك في أنه يعاقب نماذج وجود الكثير من المعلمات. وهذا قد يؤدي إلى الإفراط في الإمداد. والفرق بين بيك و إيك هو أن بيك أكثر صرامة مع فرض عقوبات إضافية على المعلمات. معيار معلومات بايزي إذا أخذنا وظيفة الاحتمال لنموذج إحصائي، الذي يحتوي على معلمات k، و L يزيد من احتمال. ثم يعطى معيار معلومات بايزي من قبل: حيث n هو عدد نقاط البيانات في السلاسل الزمنية. سنستخدم إيك و بيك أدناه عند اختيار نماذج أرما المناسبة (p، q). لتجونغ بوكس بوكس في الجزء 1 من هذه المقالة سلسلة راجان المذكورة في تعليقات ديسكوس أن اختبار لجونغ بوكس كان أكثر ملاءمة من استخدام معيار المعلومات أكايك لمعيار المعلومات بايزي في تقرير ما إذا كان نموذج أرما كان مناسبا لوقت سلسلة. اختبار يجونغ بوكس هو اختبار الفرضية الكلاسيكية التي تم تصميمها لاختبار ما إذا كانت مجموعة من أوتوكوريلاتيونس من نموذج سلسلة زمنية مجهزة تختلف اختلافا كبيرا عن الصفر. الاختبار لا يختبر كل تأخر الفردية عن العشوائية، وإنما اختبار العشوائية على مجموعة من التأخر. يجونغ-بوكس تيست نحدد الفرضية الفارغة على النحو التالي: إن بيانات السلاسل الزمنية عند كل تأخر هي i. i.d .. أي أن الارتباطات بين قيم السلسلة السكانية هي صفر. نحدد الفرضية البديلة على النحو التالي: إن بيانات السلاسل الزمنية ليست i. i.d. وتمتلك ارتباطا مسلسليا. نحسب إحصائية الاختبار التالية. س: حيث n هو طول عينة السلاسل الزمنية، فإن القبعة k هي الترابط الذاتي للعينة عند التأخر k و h هو عدد التأخيرات تحت الاختبار. وقاعدة القرار فيما يتعلق برفض الفرضية الصفرية هي التحقق مما إذا كانت Q غ تشي ch2، لتوزيع مربعات تشي مع h درجة من الحرية عند 100 (1 ألفا) من النسبة المئوية. في حين أن تفاصيل الاختبار قد تبدو معقدة قليلا، يمكننا في الواقع استخدام R لحساب الاختبار بالنسبة لنا، وتبسيط الإجراء إلى حد ما. المتوسط المتحرك المتحرك التلقائي (أرما) نماذج النظام p، q الآن بعد أن ناقشنا اختبار بيك واختبار بوكس، كنا مستعدين لمناقشة نموذجنا المختلط الأول، وهو المتوسط المتحرك للانحدار الذاتي للترتيب p أو q أو أرما (p، ف). وقد نظرنا حتى الآن في عمليات الانحدار الذاتي والمتوسط المتحرك. ويعتبر النموذج السابق سلوكه السابق كمدخلات للنموذج، وبهذه المحاولات للقبض على آثار المشاركين في السوق، مثل الزخم ومتوسط الانتعاش في تداول الأسهم. يستخدم هذا النموذج الأخير لتوصيف معلومات الصدمة لسلسلة، مثل إعلان مفاجئ للأرباح أو حدث غير متوقع (مثل انسكاب النفط بب ديبواتر هوريزون). وبالتالي، يحاول نموذج أرما التقاط كل من هذه الجوانب عند نمذجة السلاسل الزمنية المالية. لاحظ أن نموذج أرما لا يأخذ في الاعتبار تجميع التقلبات، وهو ظواهر تجريبية رئيسية للعديد من السلاسل الزمنية المالية. وهي ليست نموذجا غير متجانسة مشروطا. لذلك سنحتاج إلى الانتظار لنماذج أرش و غارتش. تعريف نموذج أرما (p، q) هو مزيج خطي من نموذجين خطيين، وبالتالي فهو في حد ذاته لا يزال خطي: الانحدار الذاتي المتوسط المتحرك للنموذج p، q نموذج السلاسل الزمنية، هو نموذج الانحدار الذاتي الانحداري للنظام p، q . أرما (p، q)، إف: ستارت alpha1 x alpha2 x لدوتس وت beta1 w beta2 w لدوتس بيتاق w إند حيث الضوضاء البيضاء مع E (وت) 0 والتباين sigma2. إذا نظرنا إلى مشغل التحول المتخلف. (انظر مقال سابق) ثم يمكننا إعادة كتابة أعلاه كدالة ثيتا و فاي من: يمكننا أن نرى بشكل مباشر أنه من خلال وضع p نيق 0 و q0 نحن استعادة أر (p) نموذج. وبالمثل إذا وضعنا p 0 و q نيق 0 نحن استرداد ما (q) نموذج. واحدة من السمات الرئيسية للنموذج أرما هو أنه شاذ ومزدوج في معلماته. وهذا يعني أن نموذج أرما غالبا ما يتطلب معلمات أقل من نموذج أر (p) أو ما (q) وحده. بالإضافة إلى ذلك إذا أعدنا كتابة المعادلة من حيث بسو، فإن ثيتا و فيي متعددة الحدود يمكن أن تشترك في بعض الأحيان عامل مشترك، مما يؤدي إلى نموذج أبسط. المحاكاة و كوريلوغرامز كما هو الحال مع نماذج الانحدار الذاتي والمتوسط المتحرك سنقوم الآن بمحاكاة مختلف سلسلة أرما ثم محاولة لتناسب نماذج أرما لهذه الإنجازات. نقوم بتنفيذ ذلك لأننا نريد أن نضمن أن نفهم الإجراء المناسب، بما في ذلك كيفية حساب فترات الثقة للنماذج، وكذلك التأكد من أن الإجراء فعلا استعادة تقديرات معقولة للمعلمات أرما الأصلية. في الجزء 1 والجزء 2 قمنا ببناء سلسلة أر و ما يدويا من خلال رسم N عينات من التوزيع الطبيعي ومن ثم صياغة نموذج سلسلة زمنية محددة باستخدام فترات تأخر هذه العينات. ومع ذلك، هناك طريقة أكثر مباشرة لمحاكاة أر، ما، أرما وحتى البيانات أريما، وذلك ببساطة عن طريق استخدام طريقة arima. sim في R. دعونا تبدأ مع أبسط نموذج أرما غير تافهة ممكن، وهي أرما (1،1 ) نموذج. وهذا هو، نموذج الانحدار الذاتي للنظام واحد جنبا إلى جنب مع نموذج المتوسط المتحرك للنظام واحد. مثل هذا النموذج له معاملين فقط، ألفا وبيتا، والتي تمثل الفواصل الأولى من السلسلة الزمنية نفسها وشروط الضوضاء البيضاء الصدمة. ويعطى هذا النموذج من قبل: نحن بحاجة إلى تحديد المعاملات قبل المحاكاة. يتيح أخذ ألفا 0.5 وبيتا -0.5: الإخراج هو كما يلي: يتيح أيضا رسم الرسم البياني: يمكننا أن نرى أنه لا يوجد ارتباط ذاتي كبير، والذي هو متوقع من نموذج أرما (1،1). وأخيرا، يتيح محاولة تحديد المعاملات والأخطاء القياسية باستخدام الدالة أريما: يمكننا حساب فترات الثقة لكل معلمة باستخدام الأخطاء القياسية: فترات الثقة لا تحتوي على قيم المعلمة الحقيقية لكلا الحالتين، ولكن يجب أن نلاحظ أن 95 فواصل الثقة واسعة جدا (نتيجة للأخطاء المعيارية الكبيرة المعقولة). يتيح الآن محاولة أرما (2،2) نموذج. وهذا هو، أر (2) نموذج جنبا إلى جنب مع ما (2) نموذج. نحن بحاجة إلى تحديد أربع معلمات لهذا النموذج: alpha1، ألفا 2، beta1 و beta2. دعونا تأخذ alpha1 0.5، alpha2-0.25 beta10.5 و beta2-0.3: إخراج أرما لدينا (2،2) نموذج على النحو التالي: و أوتوكوريلاتيون المقابلة: يمكننا الآن محاولة تركيب أرما (2،2) نموذج إلى البيانات: يمكننا أيضا حساب فترات الثقة لكل معلمة: لاحظ أن فترات الثقة لمعاملات العنصر المتوسط المتحرك (beta1 و beta2) لا تحتوي في الواقع على قيمة المعلمة الأصلية. ويوضح ذلك خطورة محاولة وضع النماذج على البيانات، حتى عندما نعرف قيم المعلمة الحقيقية ومع ذلك، فإننا نحتاج فقط لأغراض تجارية إلى أن تكون لها قدرة تنبؤية تتجاوز فرصة الإنتاج وتنتج ربحا كافيا فوق تكاليف المعاملات، لكي تكون مربحة في على المدى الطويل. الآن بعد أن رأينا بعض الأمثلة على نماذج أرما محاكاة نحن بحاجة إلى آلية لاختيار قيم p و q عند المناسب للنماذج إلى البيانات المالية الحقيقية. اختيار أفضل نموذج أرما (p، q) من أجل تحديد الترتيب p، q من نموذج أرما مناسب لسلسلة، نحتاج إلى استخدام إيك (أو بيك) عبر مجموعة فرعية من القيم p و q و ثم تطبيق اختبار لجونغ بوكس لتحديد ما إذا كان قد تم تحقيق تناسب جيد، لقيم معينة من p، س. لإظهار هذه الطريقة سنقوم أولا بمحاكاة عملية أرما (p، q) معينة. سنقوم ثم حلقة على جميع القيم الزوجية p في و q في وحساب إيك. وسوف نختار النموذج مع أدنى إيك ثم قم بتشغيل اختبار لجونغ بوكس على البقايا لتحديد ما إذا كنا قد حقق مناسبا. دعونا نبدأ من خلال محاكاة سلسلة أرما (3،2): سنقوم الآن بإنشاء كائن النهائي لتخزين أفضل نموذج تناسب وأدنى قيمة إيك. نحن حلقة على مختلف p، مجموعات q واستخدام الكائن الحالي لتخزين تناسب نموذج أرما (ط، ي)، لمتغيرات حلقة ط و j. إذا كان إيك الحالي أقل من أي إيك المحسوبة سابقا قمنا بتعيين إيك النهائي لهذه القيمة الحالية وحدد هذا الطلب. عند إنهاء حلقة لدينا ترتيب نموذج أرما المخزنة في final. order و أريما (p، د، ف) تناسب نفسها (مع مجموعة مكون المتكاملة ل 0) المخزنة كما نهائي.: لا يتيح إخراج إيك ، والنظام ومعاملات أريما: يمكننا أن نرى أن النظام الأصلي من نموذج أرما محاكاة تم استردادها، وهي P3 و Q2. يمكننا رسم مخطط المخلفات من نموذج لمعرفة ما إذا كانت تبدو وكأنها تحقيق الضوضاء البيضاء منفصلة (دون): و كوريلوغرام تبدو فعلا مثل تحقيق دون. وأخيرا، نحن إجراء اختبار يجونغ بوكس لمدة 20 تأخر لتأكيد هذا: لاحظ أن قيمة P أكبر من 0.05، التي تنص على أن المخلفات مستقلة على مستوى 95 وبالتالي أرما (3،2) نموذج يوفر نموذج جيد صالح. ومن الواضح أنه يجب أن يكون هذا هو الحال منذ أن تم محاكاة البيانات أنفسنا ومع ذلك، هذا هو بالضبط الإجراء الذي سوف نستخدم عندما نأتي لتناسب أرما (ص، ف) نماذج إلى مؤشر SampP500 في القسم التالي. البيانات المالية الآن بعد أن حددنا الإجراء لاختيار نموذج السلسلة الزمنية المثلى لسلسلة محاكاة، فمن السهل إلى حد ما لتطبيقه على البيانات المالية. لهذا المثال سوف نختار مرة أخرى مؤشر الأسهم الأمريكية SampP500. يتيح تحميل أسعار الإغلاق اليومية باستخدام كوانتمود ثم إنشاء سجل عوائد تيار: يتيح تنفيذ الإجراء المناسب نفسه كما في محاكاة أرما (3،2) سلسلة أعلاه على سجل يعود سلسلة من SampP500 باستخدام إيك: أفضل نموذج المناسب لديه أمر أرما (3،3): يتيح مؤامرة بقايا النموذج المجهزة ل SampP500 سجل تيار العوائد اليومية: لاحظ أن هناك بعض قمم كبيرة، وخاصة في فترات تأخر أعلى. وهذا يدل على سوء صالح. دعونا إجراء اختبار لجونغ بوكس لمعرفة ما إذا كان لدينا أدلة إحصائية لهذا: كما نشتبه، قيمة P أقل من 0.05 وعلى هذا النحو لا يمكننا أن نقول أن بقايا هي تحقيق الضوضاء البيضاء منفصلة. وبالتالي هناك علاقة ذاتية إضافية في المخلفات التي لم يتم تفسيرها من قبل أرما المجهزة نموذج (3،3). الخطوات التالية كما ناقشنا على طول في هذه المقالة سلسلة شهدنا أدلة على التغايرية المشروط (تجميد التقلب) في سلسلة SampP500، وخاصة في الفترات 2007-2007. عندما نستخدم نموذج غارتش في وقت لاحق في سلسلة المقال سوف نرى كيفية القضاء على هذه أوتوكوريلاتيونس. في الممارسة العملية، نماذج أرما هي عادة لا تناسب بشكل جيد لعائدات الأسهم سجل. نحن بحاجة إلى أن نأخذ بعين الاعتبار عدم التفاوت المشروط واستخدام مزيج من أريما و غارتش. ستنظر المقالة التالية أريما وتبين كيف يختلف المكون المتكامل عن نموذج أرما الذي كنا ننظر فيه في هذه المقالة. انقر أدناه لمعرفة المزيد حول. المعلومات الواردة في هذا الموقع هو رأي المؤلفين الفرديين استنادا إلى ملاحظاتهم الشخصية، وبحوثهم، وسنوات الخبرة. الناشر ومؤلفيه ليست مسجلة مستشارين الاستثمار، والمحامين، كباس أو غيرها من المهنيين الخدمات المالية ولا تقدم القانونية والضريبية والمحاسبية، وتقديم المشورة الاستثمارية أو غيرها من الخدمات المهنية. المعلومات التي يقدمها هذا الموقع هو التعليم العام فقط. ولأن كل حالة من الحالات الواقعية تختلف عن ذلك، ينبغي للقارئ أن يلتمس مستشاره الشخصي. لا يتحمل المؤلف أو الناشر أي مسؤولية أو مسؤولية عن أي أخطاء أو سهو، ولا يتحمل أي مسؤولية أو مسؤولية تجاه أي شخص أو كيان فيما يتعلق بالأضرار التي يتسبب فيها أو يزعم أنها ناجمة بشكل مباشر أو غير مباشر عن المعلومات الواردة في هذا الموقع. استخدام على مسؤوليتك الخاصة. بالإضافة إلى ذلك، قد يتلقى هذا الموقع تعويضا ماليا من الشركات المذكورة من خلال الإعلانات، والبرامج التابعة لها أو غير ذلك. تتغير الأسعار والعروض المقدمة من المعلنين الذين يظهرون على هذا الموقع بشكل متكرر، وأحيانا دون إشعار. في حين أننا نسعى جاهدين للحفاظ على المعلومات في الوقت المناسب ودقيقة، قد تكون تفاصيل العرض قديمة. ولذلك ينبغي للزائرين التحقق من شروط أي من هذه العروض قبل المشاركة فيها. يتحمل المؤلف وناشره المسؤولية عن تحديث المعلومات وإخلاء المسؤولية عن محتوى الطرف الثالث ومنتجاته وخدماته بما في ذلك عند الوصول إليه من خلال الارتباطات التشعبية والإعلانات على هذا الموقع. نموذج المتوسط المتحرك للخروج في الإحصاءات. (أرما). تسمى أحيانا نماذج بوكس-جينكينز بعد جورج بوكس و G. M. جينكينز. يتم تطبيقها عادة على بيانات السلاسل الزمنية. نظرا لسلسلة زمنية من البيانات X t. فإن نموذج أرما هو أداة لفهم القيم المستقبلية في هذه السلسلة وربما التنبؤ بها. ويتكون النموذج من جزأين، جزء الانحدار الذاتي (أر) ومتوسط متحرك (ما). وعادة ما يشار إلى النموذج باسم نموذج أرما (p، q) حيث p هو ترتيب جزء الانحدار الذاتي و q هو ترتيب جزء المتوسط المتحرك (كما هو موضح أدناه). نموذج الانحدار الذاتي يشير الرمز أر (p) إلى نموذج الانحدار الذاتي للترتيب p. ويكتب نموذج أر (p) نموذج الانحدار الذاتي هو في الأساس مرشح استجابة النبضات اللانهائية مع بعض التفسيرات الإضافية الموضوعة عليه. بعض القيود ضرورية على قيم معلمات هذا النموذج لكي يبقى النموذج ثابتا. على سبيل المثال، العمليات في نموذج أر (1) مع 1 غ 1 ليست ثابتة. مثال: أر (1) - عملية تحرير تعطى أر (1) - بروسيس التي ينتج عنها ملف لورنتزيان للكثافة الطيفية: حساب معلمات أر تحرير المعادلة أر (p) تعطى بواسطة المعادلة لأن آخر جزء من المعادلة غير صفري إلا إذا كانت m 0، عادة ما تحل المعادلة من خلال تمثيلها كمصفوفة ل غ 0، وبالتالي الحصول على المعادلة اشتقاق تحرير المعادلة التي تحدد عملية أر هي ضرب كلا الجانبين بواسطة X تم واتخاذ المتوقع القيمة التي تنتج معادلات يول-ووكر: نموذج المتوسط المتحرك إديت يشير الرمز (q) إلى نموذج المتوسط المتحرك للنظام q. حيث 1. q هي معلمات النموذج و t. t-1. هي مرة أخرى، وشروط الخطأ. نموذج المتوسط المتحرك هو في الأساس مرشح استجابة النبض المحدود مع بعض التفسيرات الإضافية الموضوعة عليه. نموذج المتوسط المتحرك للإنحدار الذاتي. يشير الرمز أرما (p. q) إلى النموذج مع عبارات الانحدار الذاتي p و q متوسط المصطلحات المتحركة. يحتوي هذا النموذج على نماذج أر (p) و ما (q)، ملاحظة حول مصطلحات الخطأ تحرير N (0، 2) حيث 2 هو التباين. قد تضعف هذه الافتراضات ولكن القيام بذلك سيغير خصائص النموذج. على وجه الخصوص، تغيير في i. i.d. فإن الافتراض سيحدث فارقا جوهريا نوعا ما. مواصفة من حيث عامل التأخر في بعض النصوص، تحدد النماذج من حيث عامل التأخر L. وفي هذه المصطلحات، يعطى النموذج أر (p) حيث يمثل تعدد الحدود (ما) يمثل نموذج ما (q) حيث يمثل الحدود متعدد الحدود وأخيرا، يعطى نموذج أرما (p. q) مجتمعا بواسطة أو أكثر بإيجاز، يمكن أن النماذج أرما بشكل عام، بعد اختيار p و q، يتم تركيبها من قبل أقل المربعات الانحدار للعثور على قيم المعلمات التي تقلل من خطأ المدى. ويعتبر عموما من الممارسات الجيدة العثور على أصغر قيم p و q التي توفر ملاءمة مقبولة للبيانات. لنموذج أر نقية ثم المعادلات يول ووكر يمكن استخدامها لتوفير مناسبا. التعميمات يعد اعتماد X على القيم السابقة وشروط الخطأ t خطيا ما لم يحدد خلاف ذلك. وإذا كان الاعتماد غير خطي، فإن النموذج يسمى على وجه التحديد المتوسط المتحرك غير الخطري (نما)، أو الانحدار الذاتي غير الخطية (نار)، أو نموذج المتوسط المتحرك غير الخطي للانحدار الذاتي (نارما). ويمكن تعميم نماذج المتوسط المتحرك للانحدار الذاتي بطرق أخرى. أنظر أيضا نماذج الانحدار الذاتي المشروط (أرش) ونماذج الانحدار الذاتي المتكامل (أريما). وفي حالة تركيب سلاسل زمنية متعددة، يمكن تركيب نموذج أريما (أو فريما) فيكتوريد. إذا كانت السلسلة الزمنية المعنية تظهر ذاكرة طويلة ثم كسور أريما (فريما، وتسمى أحيانا أرفيما) النمذجة هو المناسب. إذا كان يعتقد أن البيانات تحتوي على تأثيرات موسمية، يمكن أن يكون نموذجها نموذج ساريما (الموسمية أريما). التعميم آخر هو نموذج الانحدار الذاتي متعدد (مار). يتم فهرسة نموذج مار بواسطة العقد من شجرة، في حين يتم فهرسة نموذج الانحدار الذاتي القياسي (الوقت المنفصل) بواسطة الأعداد الصحيحة. انظر نموذج الانحدار الذاتي متعدد اللغات للحصول على قائمة المراجع. انظر أيضا تحرير المراجع إديت جورج بوكس أند F. M. جنكينز. تحليل سلسلة الوقت: التنبؤ والتحكم. الطبعة الثانية. أوكلاند، كاليفورنيا: هولدن-داي، 1976. ميلز، تيرينس C. تقنيات سلسلة الوقت للاقتصاديين. مطبعة جامعة كامبريدج، 1990. برسيفال، دونالد B. أندرو T. والدن. التحليل الطيفي للتطبيقات الفيزيائية. كامبريدج ونيفرزيتي بريس، 1993.
Comments
Post a Comment