معلومة

عتبة النموذج المدفوعة بالنشاط

عتبة النموذج المدفوعة بالنشاط



We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

لذلك أحاول معرفة عتبة لنموذج مدفوع بالنشاط. هذا موصوف في هذه الورقة: https://www.nature.com/articles/srep00469

هذا هو أنني لا أرى كيف حصل الباحثونR0بين 0.2 و 0.5 (الصفحة 6). إذا كانوا يستخدمونجاما = -2.1وإبسيلون = 10 ^ -3ويتم تحديد العتبة في المعادلة 4 ، تعطيني حساباتي عتبة2.020396043515874.

الآن إذا كان توزيع قانون السلطة (النشاط) عبارة عن رقم عشوائي ، أي هذا: http://mathworld.wolfram.com/RandomNumber.html

ينتج عن هذا متوسط ​​نشاط0.005489691793543477باستخدام هذه المعادلة:

((1-gamma) / (2-gamma)) * ((1- (epsilon ** (2-gamma))) / (1- (epsilon ** (1-gamma))))

إذن بعد المعادلة 4 ، تكون العتبة2.020396043515874، ولكن هذا ليس ما يظهره الشكل 4 (ب) لذلك أنا في حيرة من أمري.

سأكون ممتنًا حقًا إذا تمكن شخص ما من مسح هذا لي.


أعتقد أنك أخطأت في بعض الحسابات ، ربما للحظة الثانية $من توزيع النشاط. كل الحسابات هي مجرد جبر بسيط وتكاملات محددة للقوى. لتوزيع قانون السلطة $ F (a) propto a ^ {- gamma} $ with $ a in [ epsilon، 1] $ ، من أجل الحصول على $ int ^ 1_ epsilon F (a) : da = 1 $ ie التسوية الصحيحة لدينا: $$ F (a) = frac { gamma-1} { epsilon ^ {1- gamma} -1} a ^ {- gamma} $$ لحظات التوزيعات $= int ^ 1_ epsilon a ^ nF (a) : da $ هي: $$ = frac { gamma-1} { gamma-1-n} frac { epsilon ^ {1 + n- gamma} -1} { epsilon ^ {1- gamma} -1} quad if quad n ne gamma-1 $$ الحالة $ n = gamma-1 $ لها لوغاريتم وهي ليست ذات صلة في حالتنا. بالنسبة إلى $ gamma = 2.1 $ ، $ epsilon = .001 $ ، $ n = 1 $ و $ n = 2 $ نحصل على: $$ = frac { gamma-1} { gamma-2} frac { epsilon ^ {2- gamma} -1} { epsilon ^ {1- gamma} -1} = 0.00548969179354348 = frac { gamma-1} { gamma-3} frac { epsilon ^ {3- gamma} -1} { epsilon ^ {1- gamma} -1} = 0.00061164650162922 $$ بهذه الأرقام أحصل على: $$ R_0 = frac {2}{+ الجذر التربيعي {}} = 0.36330095782279 $$ والذي يتراوح بين 0.2 دولار و 0.5 دولار. أنا لست الأفضل في شرح الأشياء ، لذا يرجى إعلامي إذا كانت هناك أي شكوك! حظا طيبا وفقك الله!


(فترة طويلة جدًا للتعليق ، وربما تكون مفيدة)

افترض أنأمن المفترض أن تمثل درجة العقدة ، وباستخدام الصيغ من ويكيبيديا لتوزيع باريتو ، لدي

  • يعني (الذي أفترضه هو) =جاما * إبسيلون / (جاما -1)

    و

  • التباين (يساوي-^2، أو= يعني ^ 2 + فار) =إبسيلون ^ 2 * جاما / ((جاما -1) ^ 2 * (جاما 2))

الترجمة إلى كود R (بدلاً من الجبر):

pmean <- function (gamma = 2.8، epsilon = 0.001) {gamma * epsilon / (gamma-1)} pvar <- function (gamma = 2.8، epsilon = 0.001) {epsilon ^ 2 * gamma / ((gamma-1) ^ 2 * (gamma-2))} prms <- function (gamma = 2.8، epsilon = 0.001) {sqrt (pvar (gamma = gamma، epsilon = epsilon) + pmean (gamma = gamma، epsilon = epsilon) ^ 2) }

عندما أستخدم قيمةجاما = 2.1في شرح الشكل وأدخله في المعادلة (4):

2 * pmean (جاما = 2.1) / (pmean (جاما = 2.1) + prms (جاما = 2.1))

أحصل على 0.588 (ما زلت غير صحيح ، لكن أقرب ... ؟؟)


شاهد الفيديو: Сборка модели 38 t. Часть третья. (أغسطس 2022).