สักงอย ดอท คอม
ผมออกแบบมันและทำมาหลายเวอร์ชั่นมากยังไม่ได้ตีพิมพ์สักที เฮ้อออออ
จากรูป เปรียบเทียบผลระหว่างโมเดลที่ผมเสนอไปในปี2011 และอีก 10 ปีถัดมาเป็นเวอร์ชั่นโควิด 2021 นี้ เหตุผลที่ต้องปรับปรุงเปลี่ยนแปลงเพราะเริ่มข้อมูลใหม่ๆที่เจอเกี่ยวกับปฏิกิริยาของยาที่มีต่อเชื้อมากขึ้น และเหตุผลที่สำคัญเลยโมเดลที่ผมเสนอไปตอนแรกมันดันไปบอกว่าถ้าแบ่งยากินของครั้งต่อวันหรือทุก 12 ชมมันจะกำจัดเชื้อได้เร็วมากซึ่งมันขัดกับที่เข้าเจอกันจริงๆ โมเดลในเวอร์ชั่นใหม่นี้ก็ได้ปรับปรุงเพื่อที่จะอธิบายผลที่เจอได้มากขึ้น
เห็นคนถามคำถามนี้ในเวบpantip
คำถามลักษณะแบบนี้เราสามารถเขียนcode ในภาษาWolfram แก้ปัญหาได้ง่ายๆเลยครับ ด้วยการใช้คำสั่งอย่าง RegionDifference เพื่อดูความแตกต่างกันของพื้นที่ โดยจากโจทย์เราก็เพียงเอา พื้นที่ของสี่เหลี่ยมเป็นตัวตั้งแล้วลบออกด้วยพื้นที่วงกลมสำหรับพื้นที่ A ส่วนพื้นที่ B ก็ทำตรงกันข้ามคือเอาพื้นที่วงกลมตั้งลบออกด้วยพื้นที่สี่เหลี่ยม คำตอบที่ได้จาก Wolfram หาร4 ก็จะได้คำตอบสำหรับ พื้นที่ 1/4ของพื้นที่ที่ต่างกัน code ตัวอย่างด้านล่างนี้ผมใช้ Manipulate ครอบอีกทีเพื่อดูค่าที่รัศมีวงกลมหรือด้านของสี่เหลี่ยมต่างๆ
ggbiplot เป็น package อันหนึ่งที่ช่วยให้วาดกราฟ PCA ได้ง่ายขึ้น แต่มันถูกเขียนมาหลายปีแล้ว และก็ไม่ค่อยได้อัพเดทมาสักระยะ เวลาติดตั้งกับ R รุ่นใหม่ๆ( >4 ) ก็อาจจะมีปัญหาอยู่บ้าง วิธีการลงก็ทำตามนี้นะครับ
Sys.setenv(R_REMOTES_NO_ERRORS_FROM_WARNINGS="true")
devtools::install_github("vqv/ggbiplot")หวังว่าคงช่วยได้
พอเรารู้วิธีการดึงข้อมูลมาแล้วเราก็สามารถนำไปประยุกต์ในเรื่องต่างๆได้เช่น ลองทำเป็น API แบบง่ายๆสำหรับดูแนวโน้วการระบาด ในตัวอย่าง API ที่ผมลองทำนี้ก็เอาค่า new cases ย้อนหลัง 3 วันมาคูณกันแล้วหารากที่ 3 เพื่อดูค่าเฉลี่ยแบบ exponential แล้วก็ลองหาความชันของค่านี้ดูแบบ 7 วันย้อนหลังว่าแนวโน้มการควบคุมดีขึ้นบ้างไหม ซึ่งแน่นอนเราอยากให้ค่านี้เป็นลบ
หลังจากที่ลง opencl sdk แล้วก็ให้เพิ่มไฟล์ local ในโฟลเดอร์ make ของ cmdstan โดยที่ใน local พิมพ์ตามนี้ครับ
STAN_OPENCL=true
OPENCL_DEVICE_ID=0
OPENCL_PLATFORM_ID=0
CC = g++
CXXFLAGS += -fpermissive
LDFLAGS_OPENCL= -L"$(AMDAPPSDKROOT)lib\x86_64" -lOpenCLจากนั้นก็ให้ make ตามปกติครับ