نظرية التقدير

نظرية التقدير هي فرع من فروع الإحصاء يُعنى بتقدير قيم المعلمات بناءً على بيانات تجريبية مقاسة تتضمن عنصرًا عشوائيًا. تصف المعلمات بيئة فيزيائية كامنة بطريقة تؤثر قيمتها على توزيع البيانات المقاسة. يحاول المُقدِّر تقريب المعلمات المجهولة باستخدام القياسات. في نظرية التقدير، يُعتمد نهجان رئيسيان: [ 1 ]

  • يفترض النهج الاحتمالي (الموصوف في هذه المقالة) أن البيانات المقاسة عشوائية ولها توزيع احتمالي يعتمد على المعلمات محل الاهتمام
  • يفترض نهج الانتماء إلى المجموعة أن متجه البيانات المقاسة ينتمي إلى مجموعة تعتمد على متجه المعلمات.

أمثلة

على سبيل المثال، يُراد تقدير نسبة الناخبين الذين سيصوتون لمرشح معين. هذه النسبة هي المعلمة المطلوبة؛ ويستند التقدير إلى عينة عشوائية صغيرة من الناخبين. أو بدلاً من ذلك، يُراد تقدير احتمالية تصويت الناخب لمرشح معين، بناءً على بعض الخصائص الديموغرافية، كالعمر.

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

كمثال آخر، في نظرية الاتصالات الكهربائية، غالبًا ما ترتبط القياسات التي تحتوي على معلومات تتعلق بالمعلمات ذات الأهمية بإشارة مشوشة .

الأساسيات

بالنسبة لنموذج معين، يلزم توفر عدة "مكونات" إحصائية حتى يمكن تطبيق المُقدِّر. أولها عينة إحصائية - وهي مجموعة من نقاط البيانات المأخوذة من متجه عشوائي (RV) بحجم N. توضع في متجه ، x=[x[0]x[1]x[شمال-1]].{\displaystyle \mathbf {x} ={\begin{bmatrix}x[0]\\x[1]\\\vdots \\x[N-1]\end{bmatrix}}.} ثانيًا، هناك M معلمات θ=[θ1θ2θم]،{\displaystyle {\boldsymbol {\theta }}={\begin{bmatrix}\theta _{1}\\\theta _{2}\\\vdots \\\theta _{M}\end{bmatrix}},} والتي سيتم تقدير قيمها. ثالثًا، يجب تحديد دالة كثافة الاحتمال المستمرة (pdf) أو نظيرتها المتقطعة، دالة كتلة الاحتمال (pmf)، للتوزيع الأساسي الذي ولّد البيانات، وذلك بشرط معرفة قيم المعلمات: ص(x|θ).{\displaystyle p(\mathbf {x} |{\boldsymbol {\theta }}).\,} من الممكن أيضًا أن يكون للمعلمات نفسها توزيع احتمالي (مثل الإحصاءات البايزية ). عندئذٍ، يصبح من الضروري تحديد الاحتمال البايزي.π(θ).{\displaystyle \pi ({\boldsymbol {\theta }}).\,} بعد تشكيل النموذج، يكون الهدف هو تقدير المعلمات، والتي يُشار إليها عادةً بالرمز التالي:θ^{\displaystyle {\hat {\boldsymbol {\theta }}}}، حيث تشير "القبعة" إلى التقدير.

أحد المقدرات الشائعة هو مقدر الحد الأدنى لمتوسط ​​مربع الخطأ (MMSE)، والذي يستخدم الخطأ بين المعلمات المقدرة والقيمة الفعلية للمعلمات. هـ=θ^-θ{\displaystyle \mathbf {e} ={\hat {\boldsymbol {\theta }}}-{\boldsymbol {\theta }}} كأساس للأمثلية. ثم يتم تربيع حد الخطأ هذا، ويتم تقليل القيمة المتوقعة لهذا المربع لمقدر MMSE.

المقدرون

تشمل طرق التقدير (أساليب التقدير) والمواضيع ذات الصلة ما يلي:

أمثلة

ثابت غير معروف في الضوضاء البيضاء الغاوسية المضافة

لنفترض وجود إشارة منفصلة مستلمة ،x[ن]{\displaystyle x[n]}، لشمال{\displaystyle N}عينات مستقلة تتكون من ثابت غير معروفأ{\displaystyle A}مع ضوضاء غاوسية بيضاء مضافة (AWGN)w[ن]{\displaystyle w[n]}بمتوسط ​​صفر وتباين معروفσ2{\displaystyle \sigma ^{2}}( أي ،شمال(0،σ2){\displaystyle {\mathcal {N}}(0,\sigma ^{2})}بما أن التباين معروف، فإن المعلمة المجهولة الوحيدة هيأ{\displaystyle A}.

ثم يكون نموذج الإشارة x[ن]=أ+w[ن]ن=0،1،...،شمال-1{\displaystyle x[n]=A+w[n]\quad n=0,1,\dots ,N-1}

اثنان من المقدرات المحتملة (من بين العديد) للمعلمةأ{\displaystyle A}نكون:

  • أ^1=x[0]{\displaystyle {\hat {A}}_{1}=x[0]}
  • أ^2=1شمالن=0شمال-1x[ن]{\displaystyle {\hat {A}}_{2}={\frac {1}{N}}\sum _{n=0}^{N-1}x[n]}وهو متوسط ​​العينة

يبلغ متوسط ​​كلا هذين المقدرين 0.أ{\displaystyle A}ويمكن إثبات ذلك من خلال حساب القيمة المتوقعة لكل مقدر هـ[أ^1]=هـ[x[0]]=أ{\displaystyle \mathrm {E} \left[{\hat {A}}_{1}\right]=\mathrm {E} \left[x[0]\right]=A} و هـ[أ^2]=هـ[1شمالن=0شمال-1x[ن]]=1شمال[ن=0شمال-1هـ[x[ن]]]=1شمال[شمالأ]=أ{\displaystyle \mathrm {E} \left[{\hat {A}}_{2}\right]=\mathrm {E} \left[{\frac {1}{N}}\sum _{n=0}^{N-1}x[n]\right]={\frac {1}{N}}\left[\sum _{n=0}^{N-1}\mathrm {E} \left[x[n]\right]\right]={\frac {1}{N}}\left[NA\right]=A}

في هذه المرحلة، يبدو أن هذين المُقدِّرين يؤديان نفس الوظيفة. ومع ذلك، يصبح الفرق بينهما واضحًا عند مقارنة التباينات. vأر(أ^1)=vأر(x[0])=σ2{\displaystyle \mathrm {var} \left({\hat {A}}_{1}\right)=\mathrm {var} \left(x[0]\right)=\sigma ^{2}} و vأر(أ^2)=vأر(1شمالن=0شمال-1x[ن])=استقلال1شمال2[ن=0شمال-1vأر(x[ن])]=1شمال2[شمالσ2]=σ2شمال{\displaystyle \mathrm {var} \left({\hat {A}}_{2}\right)=\mathrm {var} \left({\frac {1}{N}}\sum _{n=0}^{N-1}x[n]\right){\overset {\text{independence}}{=}}{\frac {1}{N^{2}}}\left[\sum _{n=0}^{N-1}\mathrm {var} (x[n])\right]={\frac {1}{N^{2}}}\left[N\sigma ^{2}\right]={\frac {\sigma ^{2}}{N}}}

يبدو أن متوسط ​​العينة هو مقدر أفضل لأن تباينه أقل لكل N > 1.   

أقصى احتمال

استكمالاً للمثال باستخدام مُقدِّر الاحتمال الأقصى ، فإن دالة كثافة الاحتمال (pdf) للضوضاء لعينة واحدةw[ن]{\displaystyle w[n]}يكون ص(w[ن])=1σ2πخبرة(-12σ2w[ن]2){\displaystyle p(w[n])={\frac {1}{\sigma {\sqrt {2\pi }}}}\exp \left(-{\frac {1}{2\sigma ^{2}}}w[n]^{2}\right)} واحتماليةx[ن]{\displaystyle x[n]}يصبح (x[ن]{\displaystyle x[n]}يمكن اعتبارهاشمال(أ،σ2){\displaystyle {\mathcal {N}}(A,\sigma ^{2})}) ص(x[ن]؛أ)=1σ2πخبرة(-12σ2(x[ن]-أ)2){\displaystyle p(x[n];A)={\frac {1}{\sigma {\sqrt {2\pi }}}}\exp \left(-{\frac {1}{2\sigma ^{2}}}(x[n]-A)^{2}\right)} بالاستقلال ، احتمالx{\displaystyle \mathbf {x} }يصبح ص(x؛أ)=ن=0شمال-1ص(x[ن]؛أ)=1(σ2π)شمالخبرة(-12σ2ن=0شمال-1(x[ن]-أ)2){\displaystyle p(\mathbf {x} ;A)=\prod _{n=0}^{N-1}p(x[n];A)={\frac {1}{\left(\sigma {\sqrt {2\pi }}\right)^{N}}}\exp \left(-{\frac {1}{2\sigma ^{2}}}\sum _{n=0}^{N-1}(x[n]-A)^{2}\right)} أخذ اللوغاريتم الطبيعي لدالة كثافة الاحتمال lnص(x؛أ)=-شمالln(σ2π)-12σ2ن=0شمال-1(x[ن]-أ)2{\displaystyle \ln p(\mathbf {x} ;A)=-N\ln \left(\sigma {\sqrt {2\pi }}\right)-{\frac {1}{2\sigma ^{2}}}\sum _{n=0}^{N-1}(x[n]-A)^{2}} ومُقدِّر الاحتمال الأقصى هو أ^=argالأعلىlnص(x؛أ){\displaystyle {\hat {A}}=\arg \max \ln p(\mathbf {x} ;A)}

بأخذ المشتقة الأولى لدالة الاحتمال اللوغاريتمي أlnص(x؛أ)=1σ2[ن=0شمال-1(x[ن]-أ)]=1σ2[ن=0شمال-1x[ن]-شمالأ]{\displaystyle {\frac {\partial }{\partial A}}\ln p(\mathbf {x} ;A)={\frac {1}{\sigma ^{2}}}\left[\sum _{n=0}^{N-1}(x[n]-A)\right]={\frac {1}{\sigma ^{2}}}\left[\sum _{n=0}^{N-1}x[n]-NA\right]} وتعيينه على الصفر 0=1σ2[ن=0شمال-1x[ن]-شمالأ]=ن=0شمال-1x[ن]-شمالأ{\displaystyle 0={\frac {1}{\sigma ^{2}}}\left[\sum _{n=0}^{N-1}x[n]-NA\right]=\sum _{n=0}^{N-1}x[n]-NA}

ينتج عن ذلك مقدر الاحتمال الأقصى أ^=1شمالن=0شمال-1x[ن]{\displaystyle {\hat {A}}={\frac {1}{N}}\sum _{n=0}^{N-1}x[n]} وهو ببساطة متوسط ​​العينة. ومن هذا المثال، تبين أن متوسط ​​العينة هو مقدر الاحتمال الأقصى لـشمال{\displaystyle N}عينات من معلمة ثابتة وغير معروفة مشوهة بواسطة ضوضاء بيضاء غاوسية إضافية (AWGN).

الحد الأدنى لكريمر-راو

لإيجاد الحد الأدنى لكريمر-راو (CRLB) لمُقدِّر متوسط ​​العينة، من الضروري أولاً إيجاد رقم معلومات فيشر.أنا(أ)=هـ([أlnص(x؛أ)]2)=-هـ[2أ2lnص(x؛أ)]{\displaystyle {\mathcal {I}}(A)=\mathrm {E} \left(\left[{\frac {\partial }{\partial A}}\ln p(\mathbf {x} ;A)\right]^{2}\right)=-\mathrm {E} \left[{\frac {\partial ^{2}}{\partial A^{2}}}\ln p(\mathbf {x} ;A)\right]} والنسخ من الأعلى أlnص(x؛أ)=1σ2[ن=0شمال-1x[ن]-شمالأ]{\displaystyle {\frac {\partial }{\partial A}}\ln p(\mathbf {x} ;A)={\frac {1}{\sigma ^{2}}}\left[\sum _{n=0}^{N-1}x[n]-NA\right]}

بأخذ المشتقة الثانية 2أ2lnص(x؛أ)=1σ2(-شمال)=-شمالσ2{\displaystyle {\frac {\partial ^{2}}{\partial A^{2}}}\ln p(\mathbf {x} ;A)={\frac {1}{\sigma ^{2}}}(-N)={\frac {-N}{\sigma ^{2}}}} وإيجاد القيمة المتوقعة السالبة أمرٌ بديهي لأنها الآن ثابت حتمي -هـ[2أ2lnص(x؛أ)]=شمالσ2{\displaystyle -\mathrm {E} \left[{\frac {\partial ^{2}}{\partial A^{2}}}\ln p(\mathbf {x} ;A)\right]={\frac {N}{\sigma ^{2}}}}

وأخيرًا، وضع معلومات فيشر في vأر(أ^)1أنا{\displaystyle \mathrm {var} \left({\hat {A}}\right)\geq {\frac {1}{\mathcal {I}}}} النتائج في vأر(أ^)σ2شمال{\displaystyle \mathrm {var} \left({\hat {A}}\right)\geq {\frac {\sigma ^{2}}{N}}}

بمقارنة هذا بتباين متوسط ​​العينة (المحدد سابقًا)، يتضح أن متوسط ​​العينة يساوي الحد الأدنى لكريمر-راو لجميع قيمشمال{\displaystyle N}وأ{\displaystyle A}وبعبارة أخرى، فإن متوسط ​​العينة هو المقدر الفعال (الفريد بالضرورة) ، وبالتالي فهو أيضًا المقدر غير المتحيز ذو التباين الأدنى (MVUE)، بالإضافة إلى كونه مقدر الاحتمال الأقصى .

الحد الأقصى للتوزيع المنتظم

يُعدّ تقدير القيمة القصوى للتوزيع المنتظم أحد أبسط الأمثلة غير التافهة في مجال التقدير. ويُستخدم هذا المثال كتمرين عملي في الصف الدراسي، ولتوضيح المبادئ الأساسية لنظرية التقدير. علاوة على ذلك، في حالة التقدير القائم على عينة واحدة، يُبيّن هذا المثال قضايا فلسفية وسوء فهم محتمل في استخدام مُقدِّرات الاحتمال الأقصى ودوال الاحتمال .

بافتراض توزيع منتظم منفصل1،2،...،شمال{\displaystyle 1,2,\dots ,N}في حالة عدم معرفة القيمة القصوى، يُعطى مُقدِّر UMVU للقيمة القصوى بالصيغة التالية: ك+1كم-1=م+مك-1{\displaystyle {\frac {k+1}{k}}m-1=m+{\frac {m}{k}}-1} حيث m هي القيمة القصوى للعينة و k هو حجم العينة ، مع أخذ العينات بدون إرجاع. [ 2 ] [ 3 ] تُعرف هذه المسألة باسم مسألة الدبابات الألمانية ، وذلك بسبب تطبيق تقدير القيمة القصوى على تقديرات إنتاج الدبابات الألمانية خلال الحرب العالمية الثانية .

يمكن فهم الصيغة بشكل بديهي على النحو التالي:

"الحد الأقصى للعينة بالإضافة إلى متوسط ​​الفجوة بين الملاحظات في العينة"،

تُضاف الفجوة للتعويض عن الانحياز السلبي لأقصى قيمة في العينة كمُقدِّر لأقصى قيمة في المجتمع. [ ملاحظة 1 ]

وهذا له تباين قدره [ 2 ]1ك(شمال-ك)(شمال+1)(ك+2)شمال2ك2 للعينات الصغيرة كشمال{\displaystyle {\frac {1}{k}}{\frac {(N-k)(N+1)}{(k+2)}}\approx {\frac {N^{2}}{k^{2}}}{\text{ for small samples }}k\ll N} لذا فإن الانحراف المعياري يبلغ تقريبًاشمال/ك{\displaystyle N/k}، متوسط ​​حجم الفجوة بين العينات (في المجتمع)؛ قارنمك{\displaystyle {\frac {m}{k}}}أعلاه. يمكن اعتبار هذا حالة بسيطة للغاية لتقدير أقصى مسافة بين العناصر .

إن الحد الأقصى للعينة هو مقدر الاحتمال الأقصى للحد الأقصى للمجتمع، ولكن كما نوقش أعلاه، فهو متحيز.

التطبيقات

تتطلب العديد من المجالات استخدام نظرية التقدير. ومن هذه المجالات:

من المرجح أن تكون البيانات المقاسة عرضة للضوضاء أو عدم اليقين، ومن خلال الاحتمالية الإحصائية يتم البحث عن الحلول المثلى لاستخراج أكبر قدر ممكن من المعلومات من البيانات.

انظر أيضاً

ملحوظات

  1. لا تتجاوز القيمة القصوى للعينة القيمة القصوى للمجتمع، ولكنها قد تكون أقل، وبالتالي فهي مقدر متحيز : ستميل إلى التقليل من القيمة القصوى للمجتمع.

مراجع

الاقتباسات

  1. والتر، إي.؛ برونزاتو، ل. (1997). تحديد النماذج البارامترية من البيانات التجريبية . لندن، إنجلترا: سبرينغر-فيرلاغ.
  2. 1 2 جونسون، روجر (1994)، "تقدير حجم السكان"، تدريس الإحصاء ، 16 (2 (صيف)): 50-52 ، doi : 10.1111/j.1467-9639.1994.tb00688.x
  3. جونسون، روجر (2006)، "تقدير حجم السكان" ، الحصول على أفضل النتائج من تدريس الإحصاء ، مؤرشف من الأصل في 20 نوفمبر 2008

مصادر

  • إي إل ليمان وجي كاسيلا (1998). نظرية التقدير النقطي . سبرينغر. ISBN 0387985026.
  • ديل شيرمون (2009). هندسة تكلفة الأنظمة . دار نشر جوور. رقم ISBN 978-0-566-08861-2.
  • جون رايس (1995). الإحصاء الرياضي وتحليل البيانات . دار نشر دوكسبوري. رقم ISBN 0-534-209343.
  • ستيفن إم. كاي (1993). أساسيات معالجة الإشارات الإحصائية: نظرية التقدير . دار برنتيس هول للنشر. رقم ISBN 0-13-345711-7.
  • إتش. فنسنت بور (1998). مقدمة في كشف الإشارات وتقديرها . سبرينغر. ISBN 0-387-94173-8.
  • هاري ل. فان تريز (2001). الكشف والتقدير ونظرية التعديل، الجزء 1. وايلي. ISBN 0-471-09517-6تمت أرشفة النسخة الأصلية بتاريخ 28-04-2005.
  • دان سيمون (2006). تقدير الحالة الأمثل: مناهج كالمان، وH-infinity، والمناهج غير الخطية . وايلي. مؤرشف من الأصل بتاريخ 30-12-2010.
  • علي ح. سيد (2003). أساسيات الترشيح التكيفي . نيوجيرسي: وايلي. ISBN 0-471-46126-1.
  • يعقوب بار شالوم ؛ إكس. رونغ لي؛ ثياغالينغام كيروباراجان (2004). التقدير مع تطبيقات في التتبع والملاحة: النظرية والخوارزميات والبرمجيات . وايلي.
  • شعار ويكيميديا ​​كومنزالوسائط المتعلقة بنظرية التقدير على ويكيميديا ​​كومنز