如何将Weibull函数拟合到时间浓度曲线？

Question

我有以下时间集中数据

Time <- c(0,0.5,1,1.5,2,3,4,6,8,12,12.5,13,14,15,16,18,20,24,36,48,60,72,96)
Concentration <- c(0,30,156,181,257,198,94,40,33,28,26,25,24,21,19,17,15,14,13,11,9,7,2)

我想将Weibull函数拟合到该数据并优化以下四个参数：

function(psi,x) {
max <- psi[1]
scale <-psi[2]
shape <- psi[3]
tlag <- psi[4]
f <- max * (1-exp(-((x-tlag)^shape)/scale))
return(f)
}

我尝试使用fitdistr包，但似乎只保留分发。

感谢您的帮助！