ข้อมูล

เกณฑ์แบบจำลองที่ขับเคลื่อนด้วยกิจกรรม

เกณฑ์แบบจำลองที่ขับเคลื่อนด้วยกิจกรรม



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-แกมมา)/(2-แกมมา))*((1-(เอปซิลอน**(2-แกมมา)))/(1-(เอปซิลอน**(1-แกมมา))))

ดังนั้นตามสมการที่ 4 ค่าธรณีประตูคือ2.020396043515874แต่นั่นไม่ใช่สิ่งที่รูปที่ 4 (B) แสดง ดังนั้นฉันจึงสับสนมาก

ฉันจะขอบคุณมากถ้ามีคนสามารถล้างสิ่งนี้ให้ฉันได้


ฉันคิดว่าคุณเพิ่งคำนวณผิด อาจจะเป็นวินาทีที่สอง $$ ของการกระจายกิจกรรม การคำนวณทั้งหมดเป็นเพียงพีชคณิตง่ายๆ และปริพันธ์ของกำลังที่แน่นอน สำหรับการกระจายกฎอำนาจ $F(a)propto a^{-gamma}$ กับ $ain[epsilon,1]$ เพื่อให้มี $int^1_epsilon F(a) :da=1$ เช่น การทำให้เป็นมาตรฐานที่ถูกต้อง เรามี: $$ 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 egamma-1 $$ กรณี $n=gamma-1$ มีลอการิทึมและไม่เกี่ยวข้องในกรณีของเรา สำหรับ $gamma=2.1$, $epsilon=.001$, $n=1$ and $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}{+sqrt{}}=0.36330095782279 $$ ซึ่งอยู่ระหว่าง $0.2$ ถึง $0.5$ ฉันไม่เก่งในการอธิบายสิ่งต่าง ๆ ดังนั้นโปรดแจ้งให้เราทราบหากมีข้อสงสัย! ขอให้โชคดี!


(ยาวไปหน่อยสำหรับความคิดเห็น อาจจะมีประโยชน์)

สมมติว่าNSควรจะแสดงถึงระดับของโหนด และใช้สูตรจาก Wikipedia สำหรับการแจกแจง Pareto ฉันมี

แปลเป็นรหัส 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*(แกมมา-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(gamma=2.1)/(pmean(gamma=2.1)+prms(gamma=2.1))

ฉันได้รับ 0.588 (ยังไม่ถูกต้อง แต่ใกล้กว่า… ??)