إنتروبيا العينة

إنتروبيا العينة (SampEn؛ أو بشكل أدق إنتروبيا K₂ أو إنتروبيا ارتباط تاكنز-جراسبرجر-بروكاسيا ) هي تعديل للإنتروبيا التقريبية (ApEn؛ أو بشكل أدق "إنتروبيا بروكاسيا-كوهين")، وتُستخدم لتقييم تعقيد الإشارات الفيزيولوجية وغيرها من إشارات السلاسل الزمنية ، لتشخيص الحالات المرضية على سبيل المثال. [ 1 ] تتميز SampEn بميزتين عن ApEn: استقلالها عن طول البيانات وسهولة تطبيقها نسبيًا. كما يوجد فرق حسابي بسيط: في ApEn، تتضمن المقارنة بين متجه النموذج (انظر أدناه) وبقية المتجهات أيضًا مقارنة مع نفسه. وهذا يضمن أن الاحتمالاتجأنام(ر){\displaystyle C_{i}'^{m}(r)}لا تساوي القيم صفرًا أبدًا. وبالتالي، من الممكن دائمًا حساب لوغاريتم الاحتمالات. ولأن مقارنات القوالب مع نفسها تُقلل من قيم ApEn، تُفسَّر الإشارات على أنها أكثر انتظامًا مما هي عليه في الواقع. لا تُؤخذ هذه المقارنات الذاتية في الحسبان في SampEn. مع ذلك، ولأن SampEn يستخدم تكاملات الارتباط مباشرةً، فهو ليس مقياسًا حقيقيًا للمعلومات، بل هو تقريب. تتوفر الأسس والاختلافات بينه وبين ApEn، بالإضافة إلى دليل إرشادي خطوة بخطوة لتطبيقه، على الرابط التالي: [ 2 ]

إن SampEn مطابق تمامًا لـ "إنتروبيا الارتباط" K₂ لغراسبرغر وبروكاسيا [ 3 باستثناء أن الأخير يقترح ضرورة مراعاة حدود معينة للحصول على نتيجة ثابتة عند تغيير المتغيرات. أما SampEn فلا يأخذ في الاعتبار مثل هذه الحدود أو خصائص الثبات.

يوجد أيضًا إصدار متعدد المقاييس من SampEn، اقترحه كوستا وآخرون. [ 4 ] يمكن استخدام SampEn في البحوث الطبية الحيوية والميكانيكية الحيوية، على سبيل المثال لتقييم التحكم في وضعية الجسم. [ 5 ] [ 6 ]

تعريف

مثل الإنتروبيا التقريبية (ApEn)، تُعد إنتروبيا العينة ( SampEn ) مقياسًا للتعقيد . [ 1 ] لكنها لا تشمل الأنماط المتشابهة ذاتيًا كما تفعل ApEn. بالنسبة لبُعد تضمين مُحددم{\displaystyle m}، تسامحر{\displaystyle r} وعدد نقاط البياناتشمال{\displaystyle N}، SampEn هو اللوغاريتم الطبيعي السالب لاحتمالية أن تكون مجموعتان من نقاط البيانات المتزامنة بطولم{\displaystyle m} مسافة<ر{\displaystyle <r}ثم مجموعتان من نقاط البيانات المتزامنة بطولم+1{\displaystyle m+1} كما أن للمسافة أهمية.<ر{\displaystyle <r}ونحن نمثله من خلالSأمصهـن(م،ر،شمال){\displaystyle SampEn(m,r,N)}(أو عن طريقSأمصهـن(م،ر،τ،شمال){\displaystyle SampEn(m,r,\tau ,N)}بما في ذلك وقت أخذ العيناتτ{\displaystyle \tau }).

لنفترض الآن أن لدينا مجموعة بيانات متسلسلة زمنية بطولشمال={x1،x2،x3،...،xشمال}{\displaystyle N={\{x_{1},x_{2},x_{3},...,x_{N}\}}}بفترة زمنية ثابتةτ{\displaystyle \tau }نُعرّف متجه قالب بطولم{\displaystyle m}بحيثXم(أنا)={xأنا،xأنا+1،xأنا+2،...،xأنا+م-1}{\displaystyle X_{m}(i)={\{x_{i},x_{i+1},x_{i+2},...,x_{i+m-1}\}}}ودالة المسافةد[Xم(أنا)،Xم(ج)]{\displaystyle d[X_{m}(i),X_{m}(j)]}(i≠j) هي مسافة تشيبيشيف (ولكن يمكن أن تكون أي دالة مسافة، بما في ذلك المسافة الإقليدية ). نُعرّف إنتروبيا العينة على أنها

Sأمصهـن=-lnأب{\displaystyle SampEn=-\ln {A \over B}}

أين

أ{\displaystyle A}= عدد أزواج متجهات القوالب التي تحتوي علىد[Xم+1(أنا)،Xم+1(ج)]<ر{\displaystyle d[X_{m+1}(i),X_{m+1}(j)]<r}

ب{\displaystyle B}= عدد أزواج متجهات القوالب التي تحتوي علىد[Xم(أنا)،Xم(ج)]<ر{\displaystyle d[X_{m}(i),X_{m}(j)]<r}

يتضح من التعريف أنأ{\displaystyle A}ستكون قيمتها دائماً أصغر من أو تساويب{\displaystyle B}. لذلك،Sأمصهـن(م،ر،τ){\displaystyle SampEn(m,r,\tau )}ستكون القيمة دائمًا إما صفرًا أو قيمة موجبة. قيمة أصغر لـSأمصهـن{\displaystyle SampEn}يشير ذلك أيضاً إلى مزيد من التشابه الذاتي في مجموعة البيانات أو إلى ضوضاء أقل.

بشكل عام، نأخذ قيمةم{\displaystyle m}يكون2{\displaystyle 2}وقيمةر{\displaystyle r}يكون0.2×sتد{\displaystyle 0.2\times std}حيث يرمز std إلى الانحراف المعياري، والذي يجب حسابه على مجموعة بيانات كبيرة جدًا. على سبيل المثال، قيمة r البالغة 6 مللي ثانية مناسبة لحسابات إنتروبيا العينة لفترات معدل ضربات القلب، لأن هذا يتوافق مع0.2×sتد{\displaystyle 0.2\times std}لعدد سكان كبير جداً.

أخذ العينات متعدد المقاييس

التعريف المذكور أعلاه هو حالة خاصة من sampEn متعدد المقاييس معدلتا=1{\displaystyle \delta =1}، أيندلتا{\displaystyle \delta }يُطلق على هذا اسم مُعامل التخطي. في SampEn متعدد المقاييس، تُعرَّف متجهات القالب بفاصل زمني مُحدد بين عناصرها، ويُحدد ذلك بواسطة قيمةدلتا{\displaystyle \delta }ويتم تعريف متجه القالب المعدل على النحو التالي: Xم،دلتا(أنا)=xأنا،xأنا+دلتا،xأنا+2×دلتا،...،xأنا+(م-1)×دلتا{\displaystyle X_{m,\delta }(i)={x_{i},x_{i+\delta },x_{i+2\times \delta },...,x_{i+(m-1)\times \delta }}} ويمكن كتابة sampEn على النحو التالي Sأمصهـن(م،ر،دلتا)=-lnأدلتابدلتا{\displaystyle SampEn\left(m,r,\delta \right)=-\ln {A_{\delta } \over B_{\delta }}} ونحن نحسبأدلتا{\displaystyle A_{\delta }}وبدلتا{\displaystyle B_{\delta }}كما في السابق.

تطبيق

يمكن تطبيق مفهوم إنتروبيا العينة بسهولة في العديد من لغات البرمجة المختلفة. فيما يلي مثال مكتوب بلغة بايثون.

من itertools استورد combinations ومن math استورد logدالة construct_templates ( timeseries_data : list , m : int = 2 ): num_windows = len ( timeseries_data ) - m + 1 return [ timeseries_data [ x : x + m ] for x in range ( 0 , num_windows )]دالة get_matches ( القوالب : قائمة ، r : عدد عشري ) -> عدد صحيح : تُرجع طول ( القائمة ( تصفية ( دالة لامدا x : تطابق ( x [ 0 x [ 1 r تركيبات ( القوالب ، 2 )))) )دالة is_match ( template_1 : قائمة , template_2 : قائمة , r : عدد عشري ) -> منطقي : تُرجع all ([ abs ( x - y ) < r for ( x , y ) in zip ( template_1 , template_2 )])دالة ` sample_entropy` تأخذ قائمة من البيانات الزمنية ، وحجم نافذة من نوع عدد صحيح ، وقيمة عشوائية من نوع عدد عشري ، وتعيد : `B = get_matches ( construct_templates ( timeseries_data , window_size ), r ) `، و`A = get_matches ( construct_templates ( timeseries_data , window_size + 1 ), r )`. ثم تُرجع الدالة ` -log ( A / B )`.

مثال مكافئ في بايثون العددية.

استيراد numpyدالة بناء القوالب ( بيانات السلاسل الزمنية ، م ): عدد النوافذ = طول ( بيانات السلاسل الزمنية ) - م + 1 إرجاع مصفوفة numpy ([ بيانات السلاسل الزمنية [ س : س + م ] لـ س في النطاق ( 0 ، عدد النوافذ )])دالة get_matches ( القوالب ، r ) -> عدد صحيح : تُرجع طول ( قائمة ( تصفية ( دالة لامدا x : is_match ( x [ 0 x [ 1 r تركيبات ( قوالب )))) )دالة التوليفات ( x ) : idx = numpy.stack ( numpy.triu_indices ( len ( x ) , k = 1 ) , axis = -1 ) return x [ idx ]دالة is_match ( template_1 , template_2 , r ) -> bool : return numpy . all ([ abs ( x - y ) < r for ( x , y ) in zip ( template_1 , template_2 )])دالة ` sample_entropy` تأخذ بيانات السلاسل الزمنية ، وحجم النافذة ، ومعامل ` r` كمعاملات . تقوم هذه الدالة بحساب عدد التطابقات في كل تكرار باستخدام ` get_matches`، ثم تُنشئ قالبًا باستخدام ` construct_templates` ، وتُعيد قيمة ` A` . بعد ذلك ، تُعيد الدالة ` -numpy.log ( A / B ) ` .

يمكن العثور على مثال مكتوب بلغات أخرى:

انظر أيضاً

مراجع

  1. 1 2 ريتشمان، جيه إس؛ مورمان، جيه آر (2000). "تحليل السلاسل الزمنية الفسيولوجية باستخدام الإنتروبيا التقريبية وإنتروبيا العينة" . المجلة الأمريكية لعلم وظائف الأعضاء. القلب والدورة الدموية . 278 (6): H2039–49. doi : 10.1152/ajpheart.2000.278.6.H2039 . PMID 10843903 . 
  2. ديلجادو-بونال، ألفونسو؛ مارشاك، ألكسندر (يونيو 2019). "الإنتروبيا التقريبية وإنتروبيا العينة: دليل شامل" . إنتروبي . 21 ( 6): 541. Bibcode : 2019Entrp..21..541D . doi : 10.3390/e21060541 . PMC 7515030. PMID 33267255 .  
  3. غراسبرغر، بيتر؛ بروكاسيا، إيتامار (1983). "تقدير إنتروبيا كولموغوروف من إشارة فوضوية". مجلة Physical Review A. 28 ( 4): 2591(R). Bibcode : 1983PhRvA..28.2591G . doi : 10.1103/PhysRevA.28.2591 .
  4. كوستا، مادالينا؛ غولدبيرغر، آري؛ بينغ، سي.-كيه. (2005). "تحليل إنتروبيا متعدد المقاييس للإشارات البيولوجية". مجلة Physical Review E. 71 ( 2) 021906. Bibcode : 2005PhRvE..71b1906C . doi : 10.1103/PhysRevE.71.021906 . PMID 15783351 . 
  5. بلازكيفيتش، ميشالينا؛ كيدزيوريك، جوستينا؛ هاداموس، آنا (مارس 2021). "تأثير المدخلات البصرية والتلاعب بمنطقة الدعم على التحكم في وضعية الجسم لدى الأفراد بعد كسر فقري ناتج عن هشاشة العظام" . إنتروبي . 23 ( 3): 375. Bibcode : 2021Entrp..23..375B . doi : 10.3390/e23030375 . PMC 8004071. PMID 33804770 .  
  6. ^ هداموس، آنا؛ بيالوسزيفسكي، داريوش؛ بلاسكيويتش، ميشالينا؛ كوالسكا، ألكسندرا ج.؛ أوربانياك، إديتا؛ ويدرا، كميل ت.؛ وياديرنا، كارولينا؛ بوراتينسكي، رافائيل؛ كوبزا، أنيسكا؛ Marczyński, Wojciech (فبراير 2021). "تقييم فعالية إعادة التأهيل بعد جراحة استبدال الركبة بالكامل باستخدام عينة من الإنتروبيا والمقاييس الكلاسيكية لتوازن الجسم" . إنتروبيا . 23 (2): 164. بيب كود : 2021Entrp..23..164H . دوى : 10.3390/e23020164 . بمك 7911395 . بميد 33573057 .