Hello, I get an error message:Error in checkFunc(Func2, times, y, rho) : The number of derivatives returned by func() (22) must equal the length of the initial conditions vector (2) I try to optimise system of differential equations with 2 extra variables derived from the data.frame. I didn't manage to install XLConnect package, so I don't know if this is the source of the problem. Loading required package: XLConnectWarning message:In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package called ‘XLConnect’ I have missing data, so time frame is limited. I haven't got values for solution of the first equation, but I hope that it isn't problem, I've got the same message when putting random numbers as x values. The code is as follow:time=c(16,17,18,19,20,21,22,23,24,25,26)
#x=c(20.2,18.9,16.5)y=c(7.63,9.22,4.86,4.78,0.38,6.13,3.91,38.41,2.58,36.95,1.73) cE=c(15.05,38.01,41.09,31.41,3.54,15.92,24.01,25.29,14.82,43.93,2.45)cP=c(0.47,0.43,4.8,1.07,0.38,0.3,0.14,0.29,0.9,2.51,1.94) #df<-data.frame(time,y,cE,cP)#dfrequire(FME)require(XLConnect) #Initial values of the parametersparms=c(k1=500, k2=4500, k3=200,k4=2.42, k5=0.26,k6=12.2,k7=0.004,k8=55, k9=24,k10=8) #definition of the parameters functionderivs<-function(time,y,pars){ with(as.list(c(pars,y)),{ cE=c(15.05,38.01,41.09,31.41,3.54,15.92,24.01,25.29,14.82,43.93,2.45) cP=c(0.47,0.43,4.8,1.07,0.38,0.3,0.14,0.29,0.9,2.51,1.94) dx=(k1+(k2*cE^k10)/(k3^k10+cE^k10))/(1+cP/k6)-k4*((1+k5*cP)/(1+k7*cE))*x; #dRP_LH/dt dy=(1/k8)*k4*((1+k5*cP)/(1+k7*cE))*x-k9*y #dL/dt list(c(dx,dy)) })}initial<-c(x=x[1],y=y[1]) model_cost<-function(pars){ out<-ode(y=initial,time=time,func=derivs,parms=pars) cost<-modCost(model=out,obs=df,x="time") return(cost) model_cost(parms)$model# model fitting model_fit<-modFit(f=model_cost,p=parms)}model_cost(parms) Thank you in advance for any help. PS, I've tried to put cE and cP as initial conditions or/and as parameters, but it didn't solve the problem. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.