Doğrusal olmayan denklem için% 95 güven aralığı nasıl hesaplanır?

Manatların ağırlığını günlerden (gün olarak, dias, portekizce) tahmin etmek için bir denklemim var:

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

Ben nls () kullanarak, R modelledim ve bu grafik var:

Şimdi% 95 güven aralığını hesaplamak ve grafikte çizmek istiyorum. Her a, b ve c değişkeni için alt ve üst limitleri şu şekilde kullandım:

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

sonra daha düşük a, b, c kullanarak daha düşük bir çizgi ve daha yüksek a, b, c kullanarak daha yüksek bir çizgi çizerim. Ama bunu yapmanın doğru yolu olup olmadığından emin değilim. Bana şu grafiği veriyor:

Bunu yapmanın yolu mu yoksa yanlış mı yapıyorum?

1. Bu sitedeki bu KG, doğrusal olmayan regresyon tarafından oluşturulan eğriler etrafında güven bantları oluşturma matematiğini açıklar: Doğrusal olmayan regresyon için güven şekli ve tahmin aralıkları

2. Daha fazla okursanız , parametrelerin güven aralıklarını eğri için güven bantlarından ayırt etmeye yardımcı olacaktır .

3. Grafiğinize baktığımızda, her gün birçok günde ölçülen dört hayvandan verileriniz olduğu anlaşılıyor. Eğer öyleyse, tüm verilerin bir defada yerleştirilmesi regresyon varsayımlarından birini ihlal eder - her bir veri noktasının bağımsız olması (veya her artığın bağımsız "hatası" olması). Her bir hayvanın izini ayrı ayrı takmayı düşünebilir veya hepsini bir araya getirmek için karışık bir model kullanabilirsiniz.