วิธีการคำนวณช่วงความมั่นใจ 95% สำหรับสมการที่ไม่ใช่เชิงเส้น?


10

ฉันมีสมการในการทำนายน้ำหนักของ manatees จากอายุของพวกเขาในวัน (dias ในโปรตุเกส):

R <- function(a, b, c, dias) c + a*(1 - exp(-b*dias))

ฉันสร้างแบบจำลองมันใน R โดยใช้ nls () และได้กราฟิกนี้:

ป้อนคำอธิบายรูปภาพที่นี่

ตอนนี้ฉันต้องการคำนวณช่วงความมั่นใจ 95% และพล็อตมันในกราฟิก ฉันใช้ขีด จำกัด ล่างและสูงกว่าสำหรับแต่ละตัวแปร a, b และ c ดังนี้:

lower a = a - 1.96*(standard error of a)
higher a = a + 1.96*(standard error of a)
(the same for b and c)

จากนั้นฉันพล็อตบรรทัดล่างโดยใช้ a, b, c และบรรทัดที่สูงกว่า แต่ฉันไม่แน่ใจว่าเป็นวิธีที่ถูกต้องหรือไม่ มันทำให้ฉันกราฟิกนี้:

ป้อนคำอธิบายรูปภาพที่นี่

นี่เป็นวิธีที่จะทำหรือฉันทำผิดหรือเปล่า?

คำตอบ:


14
  1. QA ในเว็บไซต์นี้จะอธิบายคณิตศาสตร์เพื่อสร้างแถบความมั่นใจรอบ ๆ เส้นโค้งที่สร้างขึ้นจากการถดถอยแบบไม่เชิงเส้น: รูปร่างของความเชื่อมั่นและการทำนายช่วงเวลาสำหรับการถดถอยแบบไม่เชิงเส้น

  2. ถ้าคุณอ่านต่อไปก็จะช่วยในการแยกแยะความแตกต่างความเชื่อมั่นช่วงสำหรับพารามิเตอร์จากความเชื่อมั่นของวงดนตรีสำหรับเส้นโค้ง

  3. ดูกราฟของคุณแน่นอนว่าคุณมีข้อมูลจากสัตว์สี่ตัววัดแต่ละตัวในหลายวัน ถ้าเป็นเช่นนั้นการปรับข้อมูลทั้งหมดให้พร้อมกันนั้นเป็นการละเมิดข้อสันนิษฐานข้อหนึ่งของการถดถอยว่าจุดข้อมูลแต่ละจุดนั้นมีความเป็นอิสระ คุณอาจพิจารณาปรับให้เหมาะกับการติดตามสัตว์แต่ละตัวหรือใช้โมเดลผสมเพื่อให้พอดีกับพวกมันทั้งหมดในคราวเดียว


5
+1 ประเด็นที่สามมีความสำคัญ: ความพยายามใด ๆ ในการคำนวณ CIs หรือกลุ่มความเชื่อมั่นที่สมมติว่าข้อมูลเหล่านี้เป็นอิสระจะสร้างช่วงเวลาที่ไม่เพียงพอ (นั่นคือสั้น ) มันอาจช่วยให้ OP เน้นย้ำว่าวิธีที่เสนอในคำถามนั้นผิดอย่างแน่นอน: มันละเลยความสัมพันธ์ (ที่แข็งแกร่ง) ในการประมาณค่าพารามิเตอร์และรวมการจำกัดความเชื่อมั่นอย่างไม่ถูกต้อง ผลลัพธ์สุทธิดูเหมือนจริง ๆ แล้วสมเหตุสมผล (โดยอุบัติเหตุและโชค) แต่สำหรับตาของฉันมันยังไม่ได้อนุรักษ์อย่างเพียงพอ
whuber
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.