Я пытаюсь соответствовать следующему.Я измерил количество затуханий в минуту (DPM) с течением времени, и этим я хочу оценить K согласно формуле уравнения = y ~ dpm_dat[1,1] + y*(1-exp(1-k*x)).Но моя припадок выглядит довольно странно.У кого-нибудь есть идея, что не так с кодом или как можно избежать такого зигзагообразного поведения кривой соответствия?Большое спасибо
Код:
dpm_dat <- read.csv("bla.csv", sep=";")
colnames(dpm_dat) <- c("DPM", "t")
fit <- nls(DPM ~ dpm_dat[1,1] + DPM*(1-exp(1-k*t)), data = dpm_dat, start = list("k"=1))
summary(fit)$parameters
library(ggplot2)
ggplot(dpm_dat, aes(x=t, y=DPM)) +
geom_point() +
stat_smooth(method = "nls", se = FALSE,
fullrange = TRUE, method.args = list(formula=y ~ dpm_dat[1,1] + y*(1-exp(1-k*x)),
start = list("k"= 0)))+ ggtitle(paste("k=",summary(fit)$parameters))
