Dear R Community,
I am trying to optimize a water quality model that I am using. Based on
conversations with others more familiar with what I am doing I plan to
implement DEOptim to do this. The water quality model is interfaced through a
GUI. I have the input file necessary to alter parameters and run the model as a
text file.
To do the optimization I have figured out the general procedure but I need some
help on the specific methods and commands that may be helpful.
Here is an example of a section of the input file:
.
2.50000 Phytoplankton Maximum Growth Rate Constant @20 øC for Group 1 (1/day)
1.08000 Phytoplankton Growth Temperature Coefficient for Group 1
0.200000 Phytoplankton Respiration Rate Constant @20 øC for Group 1 (1/day)
1.04500 Phytoplankton Respiration Temperature Coefficient for Group 1
3.000000E-02 Phytoplankton Death Rate Constant (Non-Zoo Predation) for Group 1
(1/day)
0.000000 Phytoplankton Zooplankton Grazing Rate Constant for Group 1 (1/day)
2.500000E-02 Phytoplankton Half-Saturation Constant for N Uptake for Group 1
(mg N/L)
.
The first numbers of each line are the parameter values that I am interested in.
I have familiarized myself with DEOptim but need more guidance on writing the
function to that will run the water quality model and feed DEOptim.
In general I need to write a function to: (1) open and write over the model
input file, inserting parameter values from DEOptim to the specific lines where
they belong in the input file, (2) run the input file in the water quality
model, (3) read the output files, (4) select values of interest in the output,
(5) use values of interest in the equation to be optimized, then the result of
the equation will be what is minimized in DEOptim. Steps 2-5 I am reasonably
familiar with. I am not as familiar with how I would work with what I describe
below for step (1).
(1) (a)open the input file, (b)find the lines in the input file to insert
parameter values into(values used in DEOptim optimization between lower and
upper), (c) insert those parameter values into the input file in the proper
spot, and then run this altered input file(step 2). Is this the right way to
start?
#function
optimfxn<- function (x) {
#(1)(a) would I do something like this?
txt<-readLines("inp_File.inp")
#From here I need help finishing step 1. I have not worked much with text files
in R. I have read a lot of help files but am not getting anywhere.
#Do I have to use textConnection() to write to the text file?
#How do I incorporate the parameter values selected in each iteration of
DEOptim into the input file?
Here is an example of what I would be running in DEOptim:
optimization<-DEoptim( fn=optimfxn, lower=c(phyt_g_rt_coef=1.01000,
phyt_mx_g_rt=0.98000),upper=c(phyt_g_rt_coef=1.20000, phyt_mx_g_rt=4.44000)
Any help you could provide would be very much appreciated. Thank you in
advance.
Galen
_________________________________________________________________
Hotmail has tools for the New Busy. Search, chat and e-mail from your inbox.
N:WL:en-US:WM_HMP:042010_1
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
