LOWESS ใน Mathematica

ตั้งแต่มี RLink มานี่การทำอะไรที่เกี่ยวกับทางสถิติใน Mathematica มันดูง่ายไปหมด แต่ก่อนเวลาผมทำพวก local regression ผมต้อง export ข้อมูลจาก Mathematica ไปทำใน R หรือ Stata แต่เดี่ยวนี้ไม่ต้องแล้ว ก็ทำมัน ใน Mathematica นี่แหละด้วยการใช้ RLink เรียกคำสั่ง lowess (stats) จาก R อีกที

อันนี่เป็นตัวอย่างที่ผมใช้ทำ LOWESS (http://en.wikipedia.org/wiki/Local_regression)ครับ ผมมีข้อมูลที่เป็นคู่อันดับ (x,y) โดย x ของผมก็คือเวลา และ y คือข้อมูลจากการวัด วิธีการที่ทำก็คือสร้าง function ที่ชื่อ Lowess โดยให้มันอ่านค่า List ซึ่งเป็นคู่อันดับของข้อมูลแล้วแปลงเป็น Data Frame เพื่อเอาไปใช้กับคำสั่ง lowess ใน R ครับ

lowess

ตัวอย่างที่ได้ครับ

gam1