这可能是一个容易回答的问题,最近令我感到沮丧。我的最终目标是将观察长度数据叠加到模型的长度图上。我有一个按时间运行的模型 - 写如下:
times <-seq(1, nrow(Mussel_Water_Data),1)
我可以使用730
代替nrow(Mussel_Water_Data)
,因为我在数据集中有多少观察值来运行模型(每日水温等)。但是,我的长度数据只有455
个观察值。因此,存在差异。
我从数据集中提取了长度数据:
Mussel_Daily_LengthAS=data.frame(Mussel_Daily_Length$AS_length)
完成模型:
sol <- ode(state, times, Deb, param)
然后我试图绘制它:
plot(sol, which = "L", obs= Mussel_Daily_Length$AS_length)
我收到以下错误消息:
Error in obs[, 1] : incorrect number of dimensions
我的问题是:即使我有不同数量的观察结果,我怎样才能绘制两个?
谢谢!