我有一个类似于以下内容的数据集:
time <- c(0, 10, 20, 30, 40, 50, 60, 70 , 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190)
temp <- c(30, 28, 24, 20, 17, 13, 10, 9, 8, 7.5, 7, 6.5, 6, 5.8, 5.5, 5.3, 5.1, 5, 5, 5)
Ts = 5 # ambient temperature of surroundings
T0 = 30 # start temperature
plot(temp, time)
使用牛顿冷却定律,其中 Temp(time) ~ Ts + (T0)exp^-c*time 我试图在 R 中拟合一个指数模型来计算c
(即冷却速率)的值。
从这里的其他问题来看,我觉得 usingnls()
是合适的,但我不确定如何设置模型,因为我不知道c
. 任何意见,将不胜感激。