Malgorzata Wieteska
2017-Feb-20 15:26 UTC
[R] optimisation of the system of differential equations
I try to do optimisation of the system of the differential equations. I've managed to make my exemplary code working, however after changing equations on the exact model and data on the real one, I get an error messages and code doesn't work. I have a partial data for L (concentrations were assessed only for 3 days, the cycle lasted 24. My data for external variables (concentrations) of 2 other variables cE and cP have more datapoints (cP was assessed on each day) I couldn't fit more parameters than number of equations (excluded parameters for the initial values, which weren't assessed, so I try to get them as a result ?of parameterisation). The most sensitive parameters were labelled as k1 and k2, the values for other were derived from the literature and embedded in the system of equations. I've attached my code. Thanks in advance for any hint why this code doesn't work or information were I can find examples for optimisation of the similar system of the equations. Kind Regards,Malgorzata
Jeff Newmiller
2017-Feb-20 15:48 UTC
[R] optimisation of the system of differential equations
Go back and read the Posting Guide. 1) Missing code. Only certain types of attachments get through, and the recommendation is to include everything in the body of the email. Whatever you did, it did not get through to us. 2) Please make it clear what results you obtained and what results you expected. This is not the "do your research for you" mailing list, so please keep your question focused on how R does things, not on whether your theory is correct. -- Sent from my phone. Please excuse my brevity. On February 20, 2017 7:26:43 AM PST, Malgorzata Wieteska via R-help <r-help at r-project.org> wrote:> >I try to do optimisation of the system of the differential equations. >I've managed to make my exemplary code working, however after changing >equations on the exact model and data on the real one, I get an error >messages and code doesn't work. >I have a partial data for L (concentrations were assessed only for 3 >days, the cycle lasted 24. My data for external variables >(concentrations) of 2 other variables cE and cP have more datapoints >(cP was assessed on each day) >I couldn't fit more parameters than number of equations (excluded >parameters for the initial values, which weren't assessed, so I try to >get them as a result ?of parameterisation). The most sensitive >parameters were labelled as k1 and k2, the values for other were >derived from the literature and embedded in the system of equations. >I've attached my code. Thanks in advance for any hint why this code >doesn't work or information were I can find examples for optimisation >of the similar system of the equations. >Kind Regards,Malgorzata >______________________________________________ >R-help at 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.
Hi Malgorzata, There is currently another thread on the list about attaching R code to emails. The problem appears to be that if you attach a file with the extension ".R" it is sent in a way that is blocked by the mail server. If you want to attach an R code file to your message, change the extension to ".txt". You can also just copy it into your message, although if it is very long it may be hard to read. Jim On Tue, Feb 21, 2017 at 2:26 AM, Malgorzata Wieteska via R-help <r-help at r-project.org> wrote:> > I try to do optimisation of the system of the differential equations. I've managed to make my exemplary code working, however after changing equations on the exact model and data on the real one, I get an error messages and code doesn't work. > I have a partial data for L (concentrations were assessed only for 3 days, the cycle lasted 24. My data for external variables (concentrations) of 2 other variables cE and cP have more datapoints (cP was assessed on each day) > I couldn't fit more parameters than number of equations (excluded parameters for the initial values, which weren't assessed, so I try to get them as a result of parameterisation). The most sensitive parameters were labelled as k1 and k2, the values for other were derived from the literature and embedded in the system of equations. > I've attached my code. Thanks in advance for any hint why this code doesn't work or information were I can find examples for optimisation of the similar system of the equations. > Kind Regards,Malgorzata > ______________________________________________ > R-help at 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.