مشكلة الأنواع غير المرئية

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

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

يمكن تحليل مشكلة الأنواع غير المرئية رياضياً على النحو التالي: إذان{\displaystyle n}يتم أخذ عينات مستقلة،XنX1،...،Xن{\displaystyle X^{n}\triangleq X_{1},\ldots ,X_{n}}ثم إذام{\displaystyle m}تم أخذ المزيد من العينات المستقلة، ويُعطى عدد الأنواع غير المرئية التي سيتم اكتشافها بواسطة العينات الإضافية بالمعادلة التالية: يويو(Xن،Xن+1م+ن)|{Xن+1م+ن}{Xن}|،{\displaystyle U\triangleq U(X^{n},X_{n+1}^{m+n})\triangleq \left|\{X_{n+1}^{m+n}\}\setminus \{X^{n}\}\right|,} معXن+1م+نXن+1،...،Xن+م{\displaystyle X_{n+1}^{m+n}\triangleq X_{n+1},\ldots ,X_{n+m}}كونها المجموعة الثانية منم{\displaystyle m}عينات.

تاريخ

في أوائل أربعينيات القرن العشرين، أمضى ألكسندر ستيفن كوربيت  عامين في مالايا البريطانية يصطاد الفراشات. [ 2 ] وقد دوّن عدد الأنواع التي رصدها، وعدد أفراد كل نوع التي تم اصطيادها. فعلى سبيل المثال، رصد 74  نوعًا مختلفًا، لم يصطد منها سوى  فراشتين فقط.

عندما عاد كوربيت إلى المملكة المتحدة، تواصل مع عالم الإحصاء الحيوي رونالد فيشر وسأله عن عدد أنواع الفراشات الجديدة التي يمكن أن يتوقع اصطيادها إذا ذهب إلى نصب الفخاخ لمدة عامين آخرين؛ [ 3 ] .

ردّ فيشر بتقدير بسيط: فمع استمرار  كوربيت في نصب الفخاخ لمدة عامين إضافيين، يمكنه أن يتوقع اصطياد 75  نوعًا جديدًا. وقد استخدم في ذلك عملية جمع بسيطة (البيانات مقدمة من أورليتسكي [ 3 ] في الجدول من المثال أدناه). يو=أنا=1ن(-1)أنا+1φأنا=118-74+44-24+-12+6=75.{\displaystyle U=\sum _{i=1}^{n}(-1)^{i+1}\varphi _{i}=118-74+44-24+\cdots -12+6=75.} هناφأنا{\displaystyle \varphi _{i}}يتوافق مع عدد الأنواع الفردية التي تمت ملاحظتهاأنا{\displaystyle i}مرات. تم تأكيد مجموع فيشر لاحقًا بواسطة جود-تولمين. [ 2 ]

المقدرون

لتقدير عدد الأنواع غير المرئية، لنفترضتم/ن{\displaystyle t\triangleq m/n}ليكن عدد العينات المستقبلية (م{\displaystyle m}) مقسومًا على عدد العينات السابقة (ن{\displaystyle n})، أوم=تن{\displaystyle m=tn}. يتركφأنا{\displaystyle \varphi _{i}}ليكن عدد الأنواع الفردية التي تمت ملاحظتهاأنا{\displaystyle i}مرات (على سبيل المثال، إذا كان هناك 74 نوعًا من الفراشات مع وجود فردين مُلاحظين في جميع العينات، فإنφ2=74{\displaystyle \varphi _{2}=74}).

مقدر جود-تولمين

طُوِّر مُقدِّر غود-تولمين (GT) بواسطة غود وتولمين عام 1953. [ 4 ] ويُعطى تقدير الأنواع غير المرئية بناءً على مُقدِّر غود-تولمين بالصيغة التالية: يوجي تييوجي تي(Xن،ت)-أنا=1(-ت)أناφأنا.{\displaystyle U^{\text{GT}}\triangleq U^{\text{GT}}(X^{n},t)\triangleq -\sum _{i=1}^{\infty }(-t)^{i}\varphi _{i}.} لقد ثبت أن مقدر جود-تولمين هو تقدير جيد لقيمت1.{\displaystyle t\leq 1.}كما أن مقدر جود-تولمين يحقق تقريبًا ما يلي: هـ(يوجي تي-يو)2نت2.{\displaystyle \operatorname {\mathbb {E} } (U^{\text{GT}}-U)^{2}\lesssim nt^{2}.} هذا يعني أنيوجي تي{\displaystyle U^{\text{GT}}}تقديراتيو{\displaystyle U}إلى الداخلنت،{\displaystyle {\sqrt {n}}\cdot t,}طالمات1.{\displaystyle t\leq 1.}

لكن بالنسبة لـت>1،{\displaystyle t>1,}يفشل مُقدِّر غود-تولمين في الحصول على نتائج دقيقة. وذلك لأنه إذات>1،{\displaystyle t>1,}يوجي تي{\displaystyle U^{\text{GT}}}زيادة بنسبة(-ت)أناφأنا{\displaystyle (-t)^{i}\varphi _{i}}لأنا{\displaystyle i}معφأنا>0،{\displaystyle \varphi _{i}>0,}بمعنى أنه إذاφأنا>0،{\displaystyle \varphi _{i}>0,}يوجي تي{\displaystyle U^{\text{GT}}}ينمو بشكل فائق الخطية فيت،{\displaystyle t,}لكنيو{\displaystyle U}يمكن أن ينمو بشكل خطي على الأكثر معت.{\displaystyle t.}لذلك، عندمات>1،{\displaystyle t>1,}يوجي تي{\displaystyle U^{\text{GT}}}ينمو أسرع منيو{\displaystyle U}ولا يقترب من القيمة الحقيقية. [ 3 ]

وللتعويض عن ذلك، أظهر إيفرون وثيستيد في عام 1976 [ 1 ] أن تحويل أويلر المقتطع يمكن أن يكون أيضًا تقديرًا قابلاً للاستخدام (تقدير "ET"): يوETأنا=1نححETφأنا،{\displaystyle U^{\text{ET}}\triangleq \sum _{i=1}^{n}h_{h}^{\text{ET}}\cdot \varphi _{i},} مع حأناET(-ت)أنا+1P(Xأنا)،{\displaystyle h_{i}^{\text{ET}}\triangleq (-t)^{i+1}\cdot \mathbb {P} (X\geq i),} أينXسلة المهملات(ك،11+ت)،{\displaystyle X\sim \operatorname {Bin} \left(k,{\frac {1}{1+t}}\right),}و P(Xأنا)={ج=أناك(كج)تك-ج(1+ت)ك ل أناك،0 ل أنا>ك،{\displaystyle \mathbb {P} (X\geq i)={\begin{cases}\displaystyle \sum _{j=i}^{k}{\binom {k}{j}}{\frac {t^{kj}}{(1+t)^{k}}}&{\text{ لـ }}i\leq k,\\0&{\text{ لـ }}i>k,\end{cases}}} أينك{\displaystyle k}هو الموقع الذي تم اختياره لقطع تحويل أويلر.

مقدر جود-تولمين المُعدَّل

على غرار منهج إيفرون وثيستيد، طوّر كلٌّ من ألون أورليتسكي ، وأناندا ثيرثا سوريش، وييهونغ وو، مُقدِّر غود-تولمين السلس. وقد أدركوا أن فشل مُقدِّر غود-تولمين يعود إلى النمو الأسي، وليس إلى تحيُّزه. [ 3 ] ولذلك، قدّروا عدد الأنواع غير المرئية عن طريق اقتطاع السلسلة. يول-أنا=1ل(-ت)أناφأنا.{\displaystyle U^{l}\triangleq -\sum _{i=1}^{l}(-t)^{i}\varphi _{i}.} وأشار كل من أورليتسكي وسوريش وو أيضًا إلى أن التوزيعات ذاتت>1{\displaystyle t>1}، الحد الرئيسي في تقدير المجموع هول-ذ{\displaystyle l-{\text{th}}}المصطلح، بغض النظر عن قيمةل{\displaystyle l}[ 2 ] لحل هذه المسألة، اختاروا عددًا صحيحًا عشوائيًا غير سالب .ل{\displaystyle L}، واختصر السلسلة عندل{\displaystyle L}ثم أخذ المتوسط ​​على توزيع حولل{\displaystyle L}[ 3 ] المُقدِّر الناتج هو يول=هـل[-أنا=1ل(-ت)أناφأنا].{\displaystyle U^{L}=\operatorname {E} _{L}\left[-\sum _{i=1}^{L}(-t)^{i}\varphi _{i}\right].} تم اختيار هذه الطريقة بسبب تحيزيول{\displaystyle U^{l}}تغير الإشارات بسبب(-ت)أنا{\displaystyle (-t)^{i}}معامل. حساب المتوسط ​​على توزيع منل{\displaystyle L}وبالتالي يقلل من التحيز. وهذا يعني أنه يمكن كتابة المُقدِّر على شكل توليفة خطية للانتشار: [ 2 ]يول=هـل[-أنا1(-ت)أناφأنا1أنال]=-أنا1(-ت)أنابرو(لأنا)φأنا.{\displaystyle U^{L}=\operatorname {E} _{L}\left[-\sum _{i\geq 1}(-t)^{i}\varphi _{i}\mathbf {1} _{i\leq L}\right]=-\sum _{i\geq 1}(-t)^{i}\Pr(L\geq i)\varphi _{i}.} اعتمادًا على توزيعل{\displaystyle L}ستختلف النتائج باختلاف الخيار المُختار. باستخدام هذه الطريقة، يمكن وضع تقديرات لـتlnن{\displaystyle t\propto \ln n}وهذا هو الأفضل على الإطلاق. [ 3 ]

منحنى اكتشاف الأنواع

يمكن أيضًا استخدام منحنى اكتشاف الأنواع . يربط هذا المنحنى عدد الأنواع الموجودة في منطقة ما كدالة للزمن. ويمكن إنشاء هذه المنحنيات باستخدام مُقدِّرات (مثل مُقدِّر جود-تولمين) ورسم عدد الأنواع غير المرئية عند كل قيمة.ت{\displaystyle t}[ 5 ]

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

مثال: فراشات كوربيت

على سبيل المثال، لننظر إلى البيانات التي قدمها كوربيت لفيشر في أربعينيات القرن العشرين. [ 3 ] باستخدام نموذج جود-تولمين، يتم إيجاد عدد الأنواع غير المرئية باستخدام يو=-أنا=1(-ت)أناφأنا.{\displaystyle U=-\sum _{i=1}^{\infty }(-t)^{i}\varphi _{i}.} ويمكن استخدام ذلك بعد ذلك لإنشاء علاقة بينت{\displaystyle t}ويو{\displaystyle U}.

البيانات المقدمة إلى فيشر من قبل كوربيت [ 3 ]
عدد الأعضاء الذين تمت ملاحظتهم،أنا{\displaystyle i}123456789101112131415
عدد الأنواع،φأنا{\displaystyle \varphi _{i}}11874442429222019201512146126

تظهر هذه العلاقة في الرسم البياني أدناه.

عدد الأنواع غير المرئية كدالة لـ t، ونسبة العينات الجديدة إلى العينات السابقة.

يتضح من الحبكة أنه فيت=1{\displaystyle t=1}، وهو ما يمثل قيمةت{\displaystyle t}ذلك الذي قدمه كوربيت إلى فيشر، التقدير الناتج عنيو{\displaystyle U}تبلغ القيمة 75، وهو ما يتطابق مع ما توصل إليه فيشر. كما يمثل هذا الرسم البياني منحنى اكتشاف الأنواع لهذا النظام البيئي، ويحدد عدد الأنواع الجديدة التي سيتم اكتشافها.ت{\displaystyle t}يزداد (ويتم أخذ المزيد من العينات).

استخدامات أخرى

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

مثال: كم عدد الكلمات التي كان يعرفها شكسبير؟

استنادًا إلى بحثٍ أجراه ثيستيد وإيفرون حول أعمال شكسبير المعروفة، يبلغ إجمالي عدد الكلمات 884,647 كلمة. [ 1 ] كما وجد البحث أن هناك إجماليًا قدرهشمال=864{\displaystyle N=864}الكلمات المختلفة التي تظهر أكثر من 100 مرة. وبناءً على ذلك، وُجد أن العدد الإجمالي للكلمات الفريدة هو 31,534 كلمة. [ 1 ] بتطبيق نموذج جود-تولمين، إذا تم اكتشاف عدد مماثل من أعمال شكسبير، فإنه يُقدّر أنيوكلمات11460{\displaystyle U^{\text{words}}\approx 11,460}سيتم العثور على كلمات فريدة. والهدف هو استخلاصيوكلمات{\displaystyle U^{\text{words}}}لت={\displaystyle t=\infty }يقدر ثيستيد وإفرون أنيوكلمات(ت)35000{\displaystyle U^{\text{words}}(t\to \infty )\approx 35{,}000}وهذا يعني أن شكسبير كان على الأرجح يعرف أكثر من ضعف عدد الكلمات التي استخدمها فعلياً في جميع كتاباته. [ 1 ]

انظر أيضاً

مراجع

  1. 1 2 3 4 5 إيفرون، برادلي؛ ثيستيد، رونالد (1976). "تقدير عدد الأنواع غير المرئية: كم عدد الكلمات التي عرفها شكسبير؟". Biometrika . 63 (3): 435–447 . doi : 10.2307/2335721 . JSTOR 2335721 . 
  2. 1 2 3 4 أورليتسكي، ألون؛ سوريش، أناندا ثيرثا؛ وو، ييهونغ (22-11-2016). "التنبؤ الأمثل بعدد الأنواع غير المرئية" . وقائع الأكاديمية الوطنية للعلوم . 113 (47): 13283-13288 . Bibcode : 2016PNAS..11313283O . doi : 10.1073/pnas.1607774113 . PMC 5127330. PMID 27830649 .  
  3. 1 2 3 4 5 6 7 8 أورليتسكي، ألون؛ سوريش، أناندا ثيرثا؛ وو، ييهونغ (23-11-2015). "تقدير عدد الأنواع غير المرئية: عصفور في اليد خير من لوغاريتم ن في الأدغال". arXiv : 1511.07428 [ math.ST ].
  4. جود، آي جيه؛ تولمين، جي إتش (1956). "عدد الأنواع الجديدة، وزيادة التغطية السكانية، عند زيادة حجم العينة" . بيومتريكا . 43 ( 1-2 ): 45-63 . doi : 10.1093/biomet/43.1-2.45 . ISSN 0006-3444 . 
  5. 1 2 بيبر، د. ب؛ ماريوت، ف. هـ. س؛ جاستون، ك. ج؛ هاريس، س. أ؛ سكوتلاند، ر. و (7 يوليو 2007). " التنبؤ بأعداد الأنواع غير المعروفة باستخدام منحنيات الاكتشاف" . وقائع الجمعية الملكية ب: العلوم البيولوجية . 274 (1618): 1651-1658 . doi : 10.1098/rspb.2007.0464 . PMC 2169286. PMID 17456460 .