แบบจำลองคณิตศาสตร์การลดลงของเชื้อมาลาเรียระหว่างการรักษาด้วยยาอาร์ทิซูเนต

ผมออกแบบมันและทำมาหลายเวอร์ชั่นมากยังไม่ได้ตีพิมพ์สักที เฮ้อออออ

จากรูป เปรียบเทียบผลระหว่างโมเดลที่ผมเสนอไปในปี2011 และอีก 10 ปีถัดมาเป็นเวอร์ชั่นโควิด 2021 นี้ เหตุผลที่ต้องปรับปรุงเปลี่ยนแปลงเพราะเริ่มข้อมูลใหม่ๆที่เจอเกี่ยวกับปฏิกิริยาของยาที่มีต่อเชื้อมากขึ้น และเหตุผลที่สำคัญเลยโมเดลที่ผมเสนอไปตอนแรกมันดันไปบอกว่าถ้าแบ่งยากินของครั้งต่อวันหรือทุก 12 ชมมันจะกำจัดเชื้อได้เร็วมากซึ่งมันขัดกับที่เข้าเจอกันจริงๆ โมเดลในเวอร์ชั่นใหม่นี้ก็ได้ปรับปรุงเพื่อที่จะอธิบายผลที่เจอได้มากขึ้น

It doesn’t matter how beautiful your theory is, it doesn’t matter how smart you are. If it doesn’t agree with experiment, it’s wrong.” Richard P. Feynman

หาพื้นที่แรเงา

เห็นคนถามคำถามนี้ในเวบpantip

ช่วยหาพื้นที่ส่วนที่แรเงาหน่อยครับ – Pantip

คำถามลักษณะแบบนี้เราสามารถเขียนcode ในภาษาWolfram แก้ปัญหาได้ง่ายๆเลยครับ ด้วยการใช้คำสั่งอย่าง RegionDifference เพื่อดูความแตกต่างกันของพื้นที่ โดยจากโจทย์เราก็เพียงเอา พื้นที่ของสี่เหลี่ยมเป็นตัวตั้งแล้วลบออกด้วยพื้นที่วงกลมสำหรับพื้นที่ A ส่วนพื้นที่ B ก็ทำตรงกันข้ามคือเอาพื้นที่วงกลมตั้งลบออกด้วยพื้นที่สี่เหลี่ยม คำตอบที่ได้จาก Wolfram หาร4 ก็จะได้คำตอบสำหรับ พื้นที่ 1/4ของพื้นที่ที่ต่างกัน code ตัวอย่างด้านล่างนี้ผมใช้ Manipulate ครอบอีกทีเพื่อดูค่าที่รัศมีวงกลมหรือด้านของสี่เหลี่ยมต่างๆ

ggbiplot

ggbiplot เป็น package อันหนึ่งที่ช่วยให้วาดกราฟ PCA ได้ง่ายขึ้น แต่มันถูกเขียนมาหลายปีแล้ว และก็ไม่ค่อยได้อัพเดทมาสักระยะ เวลาติดตั้งกับ R รุ่นใหม่ๆ( >4 ) ก็อาจจะมีปัญหาอยู่บ้าง วิธีการลงก็ทำตามนี้นะครับ

Sys.setenv(R_REMOTES_NO_ERRORS_FROM_WARNINGS="true")
devtools::install_github("vqv/ggbiplot")

หวังว่าคงช่วยได้

ลองดึงข้อมูลผู้ติดเชื้อโควิดของไทยจากเว็บมาplotด้วยภาษาWolfram

พอเรารู้วิธีการดึงข้อมูลมาแล้วเราก็สามารถนำไปประยุกต์ในเรื่องต่างๆได้เช่น ลองทำเป็น API แบบง่ายๆสำหรับดูแนวโน้วการระบาด ในตัวอย่าง API ที่ผมลองทำนี้ก็เอาค่า new cases ย้อนหลัง 3 วันมาคูณกันแล้วหารากที่ 3 เพื่อดูค่าเฉลี่ยแบบ exponential แล้วก็ลองหาความชันของค่านี้ดูแบบ 7 วันย้อนหลังว่าแนวโน้มการควบคุมดีขึ้นบ้างไหม ซึ่งแน่นอนเราอยากให้ค่านี้เป็นลบ

cmdstan กับ opencl

หลังจากที่ลง 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 ตามปกติครับ

จาก Stan Math Library: OpenCL CPU/GPU Support (mc-stan.org)