Эта кривая:
cu=function(t){((-(c-t)*I/(a-t)*(1+exp((b-t)*I))+1))}
представляет затухание в переменной отклика как функцию t и еще 4 параметров, один из которых я не знаю («I»). Кривая должна проходить как можно ближе к некоторым точкам (см. Ниже). Можно ли будет использовать эту информацию для оптимизации I? Я смотрел на другие связанные вопросы, но я не знаю, как адаптировать мой случай к функции оптимизации, такой как «optim». Как мне ввести «t» и «I» в функцию оптимизации или использовать известные наблюдения ??
Здесь приведены значения для a, b, c, разумные начальные значения для I и значения дляответ при различных значениях t
a=39.9
b=42.5
c=47
c(2.9,3.1)
y(41.5)=230
y(42.5)=3.4
y(44)=1.27
y(47)=1
Этот фрагмент кода позволяет визуализировать кривую и точки, когда I = 3:
cu=function(t){((-(c-t)*I/(a-t)*(1+exp((b-t)*I))+1))}
I=3
a=39.9
b=42.5
c=47
curve(cu,from=a,to=c,col="blue",n=100000,
ylim=c((1),(2000)),
xlim=c(33,50),ylab="response",xlab="predictor", add=F)
points(41.5,(230),cex=1,pch=21,bg="blue")
points(44,(76.4/60),cex=1,pch=21,bg="blue")
points(b,(3.4),pch=21,cex=1,bg="blue")
points(c,(1),pch=21,bg="blue")#