From ba0ffaa9579322d0dca4973b4750432fccd3e83b Mon Sep 17 00:00:00 2001
From: Johannes Ranke As gmkin is an R package, you need to start R and load the gmkin package before you can run gmkin. The latter can be achieved by entering the command into the R console. This will also load the packages that gmkin depends on, most notably gWidgetsWWW2 and mkin. Loading the package only has to be done once after you have started R. Before you start gmkin, you should make sure that R is using the working directory that you would like to keep your gmkin project file(s) in. If you use the standard R GUI application on windows, you can change the working directory from the File menu (‘File’ -> ‘Change dir…’). Once you are sure that the working directory is what you want it to be, gmkin can be started by entering the R command This will cause the default browser to start up or, if it is already running, to pop up and open a new tab for displaying the gmkin user interface. In the R console, you should see a message that the httpd help server is started, if it wasn’t already started before. In the browser, you should see something like the screenshot below. The statusbar at the bottom of the gmkin window shows, among others, the working directory that gmkin uses. If the browser only shows “Loading ExtJS…” and nothing else happens for about 10 to 15 seconds, please enter the A gmkin project workspace contains datasets, kinetic models for fitting, and so-called fits, i.e. the results of fitting models to data. The project area also shows the current working directory, where project workspace files are saved using the file extension Once a project has been saved by the user, the project explorer to the left will show it in the project list. The current state of a project should repeatedly be saved during the work in order to avoid loosing data. This can be achieved by selecting the ‘Project’ tab in the center and pressing the ‘Save project to project file’ button. More conveniently, the keyboard shortcut Shift-F12 (was Ctrl-X in gmkin < 0.6.6) can be used to save the current status of the project. In the project file management area, datasets and models can be imported from one of the projects in the project list, once it has been selected in the droplist labelled ‘Import from:’. Once the project has been set up, you may want to minimize the project explorer, especially if you are limited in vertical screen space. When you select one of the datasets in the dataset explorer to the left, some summary information about the dataset is shown in the center, and the data itself is loaded into the data editor to the right. When you have added information about the units, or edited the data to the right, you should hit the button ‘Keep changes’. This will add a new entry in the dataset explorer. In the dataset editor to the right, you can override original data (for example in order to follow FOCUS recommendations for time zero samples or values below the limit of detection) by entering numbers in the override column. For entering new data manually, there are two possibilities. You can either change the Dataset title of the current dataset and edit the data in the Data editor to the right. Or, if the new data should have a different structure, e.g. different sampling times, observed variables and replicates, click on “New dataset”, edit the dataset title, sampling times, number of replicates and the list of observed variables, and press the button ‘Generate grid for entering kinetic data’, in order to prepare the Data editor to the right. For entering new data manually, there are two possibilities. You can either change the Dataset title of the current dataset and edit the data in the Data editor to the right. Or, if the new data should have a different structure, e.g. different sampling times, observed variables and replicates, click on “New dataset”, edit the dataset title, sampling times, number of replicates and the list of observed variables, and press the button ‘Generate grid for entering kinetic data’, in order to prepare the Data editor to the right. For sampling times and short names of the relevant compounds, a comma separated list must be entered, with a space after each comma. An example of filling out the respective fields is shown below. After that, the actual measurements can be entered into the Data editor to the right, in the column ‘value’. If everything is OK, press “Keep changes” to save the dataset in the current workspace. Note that you need to save the project file (see above) in order to be able to use the dataset that you created in a future gmkin session. In case you want to work with a larger dataset that is already available as a computer file e.g. in a spreadsheet application, you can export these data as a tab separated or comma separated text file and import it using the “Browse” and “Upload” buttons in the dataset editor. As an example, we can create a text file from one of the datasets shipped with the mkin package using the following R command: This produces a text file with comma separated values in the current working directory of R. Loading this text file into gmkin using the “Browse” and “Upload” buttons results in an an import configuration area showing the top lines of the imported file, and giving the possibility to chose the import options matching the file format. In the import configuration area, the following options can be specified. In the field “Comment lines”, the number of lines in the beginning of the file that should be ignored can be specified. The checkbox on the next line should be checked if the first line of the file contains the column names, i.e. the names of the observed variables when the data are in wide format. As the “Decimal” separator, comma “,” or period “.” can be selected. In the next line, it can be specified if the data are in wide or in long format. If in wide format, the only option left to specify is the title of the column containing the sampling times. If the data is in long format, the column headings specifying the columns containing the observed variables (default is “name”), the sampling times (default is “time”), the observed values (default is “value”) and, if present in the data, the relative errors (default is “err”) can be adapted. The default settings appearing if the long format is selected are shown below. In our example we have data in the wide format, and after adapting the “Separator” to a comma, we can press the button “Import using options specified below”, and the data should be imported. If successful, the dataset editor should show the sampling times and the names of the observed variables, and the Data editor should show the imported data in a grid for further editing or specifying overrides. After adapting the dataset title and possibly the units, press “Keep changes” to save the dataset in the current workspace. Again, you need to save the project file in order to be able to use the dataset that you created in a future gmkin session. The following screenshot shows the model editor for the model ‘Z.2a.ff’ from the project workspace ‘FOCUS_2006_Z’ provided with the gmkin package. In the first line the name of the model can be edited. You can also specify “min” or “max” for minimum or maximum use of formation fractions. Maximum use of formation fractions means that the differential equations in the degradation model are formulated using formation fractions. When you specify “min”, then formation fractions are only used for the parent compound when you use the FOMC, DFOP or the HS model for it. Pressing “Add observed variable” adds a line in the array of state variable specifications below. The observed variables to be added are usually transformation products (usually termed metabolites), but can also be the parent compound in a different compartment (e.g. “parent_sediment”). Pressing “Add observed variable” adds a line in the array of state variable specifications below. The observed variables to be added are usually transformation products (usually termed metabolites), but can also be the parent compound in a different compartment (e.g. “parent_sediment”). Only observed variable names that occur in previously defined datasets or models can be selected. For any observed variable other than the first one, only the SFO or the SFORB model can be selected. For each observed variables, a comma separated list of target variables can be specified. In addition, a pathway to the sink compartment can be selected. If too many observed variables have been added, complete lines can be removed from the model definition by pressing the button “Remove observed variable”. If the model definition is supposedly correct, press “Keep changes” to make it possible to select it for fitting in the listing of models to the left. If a dataset and a kinetic model are selected, the button “Configure fit” below in the Explorer window “Configuration” becomes active. Pressing it opens the “Configuration” area in the center and loads the dataset into the Data editor to the right, for viewing only. After pressing the button ‘Plot unoptimised’, the data and the kinetic model with default starting parameters are plotted in the ‘Plot’ area to the right. In the right part of the fit configuration area, the file format can be chosen, the legend can be turned off, and the observed variables for which the data and the model fit should be plotted can be selected as shown below. On systems running the Windows operating system, the enhanced metafile (emf) format can be additionally chosen. Changing the file format for plotting will also change the extension of the proposed filename for saving the plot. In many cases the starting parameters and the fit options do not need to be modified and the model fitting process can simply be started by pressing the “Run fit” button. In the R console, the progressive reduction in the model cost can be monitored and will be displayed in the following way: If plotting of the fitting progress was selected, a new separate graphics window should either pop up, or a graphics window previously started for this purpose will be reused. If your screen size allows for it, you can arrange the R plotting window and the R console in a way that you can see everything at the same time: Once a fit has successfully been completed, the most important results are shown in several tables in the “Result” area as shown below. The detailed summary can be accessed below these tables. The complete summary can be saved into a text file by specifying a suitable file name and pressing the button “Save summary”. Whenever a new fit has been configured or a run of a fit has been completed, the plotting area is updated with the abovementioned plot of the data and the current model solution. In addition, a confidence interval plot is shown below this conventional plot. In case a fit has been run and confidence intervals were successfully calculated for the fit ( i.e. if the model was not overparameterised and no other problems occurred), the confidence intervals are graphically displayed as bars as shown below. In addition, a confidence interval plot is shown below this conventional plot. In case a fit has been run and confidence intervals were successfully calculated for the fit (i.e. if the model was not overparameterised and no other problems occurred), the confidence intervals are graphically displayed as bars as shown below. Windows and Linux users running R 3.1.0 or later can make use of my drat package repository on github. Before installing gmkin, you should update your R packages already installed, e.g. by starting R, and pasting the command or similar, depending on your R installation. For installing gmkin or upgrading to the latest released version, please install the drat package add my repository and install as usual The latest changes to gmkin are recorded in the NEWS file, more details can be found in the commit history. Without further preparation, you can start the GUI by starting R and pasting the following commands into the R terminal (R console on Windows): You may also want to adapt the browser that R starts (using e.g. Users of the Installing gWidgetsWWW2 in this way yields a lot of warnings concerning overly long path names (see Notes on gWidgetsWWW2 above). Using You can put the commands into an R script (e.g. under Data taken from FOCUS (2006), Appendix 7, p. 350.Manual for gmkin
+ Manual for gmkin
Johannes Ranke
- 2020-02-07
-
+ 2021-02-12
+ Source: vignettes/gmkin_manual.Rmd
gmkin_manual.Rmd
Starting gmkin
+
gmkin()
gmkin()
command again in the R console. If this still does not bring up the gmkin GUI, please close the browser and try again..gmkinws
.
Entering data directly
-write.table(schaefer07_complex_case, sep = ",", dec = ".",
- row.names = FALSE, quote = FALSE,
- file = "schaefer07.csv")
+
write.table(schaefer07_complex_case, sep = ",", dec = ".",
+ row.names = FALSE, quote = FALSE,
+ file = "schaefer07.csv")
@@ -226,7 +236,7 @@
Plot options
Fitting the model
> Ordinary least squares optimisation
-Sum of squared residuals at call 1: 2388.077
-Sum of squared residuals at call 3: 2388.077
-Sum of squared residuals at call 4: 247.1962
-Sum of squared residuals at call 7: 200.6791
-Sum of squared residuals at call 10: 197.7231
-Sum of squared residuals at call 11: 197.0872
-Sum of squared residuals at call 14: 196.535
-Sum of squared residuals at call 15: 196.535
-Sum of squared residuals at call 16: 196.535
-Sum of squared residuals at call 17: 196.5334
-Sum of squared residuals at call 20: 196.5334
-Sum of squared residuals at call 25: 196.5334
-Negative log-likelihood at call 31: 26.64668
-Optimisation successfully terminated.
> Ordinary least squares optimisation
+1: 2388.077
+ Sum of squared residuals at call 3: 2388.077
+ Sum of squared residuals at call 4: 247.1962
+ Sum of squared residuals at call 7: 200.6791
+ Sum of squared residuals at call 10: 197.7231
+ Sum of squared residuals at call 11: 197.0872
+ Sum of squared residuals at call 14: 196.535
+ Sum of squared residuals at call 15: 196.535
+ Sum of squared residuals at call 16: 196.535
+ Sum of squared residuals at call 17: 196.5334
+ Sum of squared residuals at call 20: 196.5334
+ Sum of squared residuals at call 25: 196.5334
+ Sum of squared residuals at call -likelihood at call 31: 26.64668
+ Negative log Optimisation successfully terminated.
Confidence interval plots
All vignettes
-
-
+
+
@@ -134,7 +144,7 @@
diff --git a/docs/authors.html b/docs/authors.html
index b5c6fb1..9cd757e 100644
--- a/docs/authors.html
+++ b/docs/authors.html
@@ -10,23 +10,27 @@
-
+
-
+
-
+
+
+
+
+
-
-
+
+
-
+
-
-
+
+
@@ -53,7 +57,7 @@
-
+
@@ -143,7 +152,7 @@
diff --git a/docs/bootstrap-toc.css b/docs/bootstrap-toc.css
new file mode 100644
index 0000000..5a85941
--- /dev/null
+++ b/docs/bootstrap-toc.css
@@ -0,0 +1,60 @@
+/*!
+ * Bootstrap Table of Contents v0.4.1 (http://afeld.github.io/bootstrap-toc/)
+ * Copyright 2015 Aidan Feldman
+ * Licensed under MIT (https://github.com/afeld/bootstrap-toc/blob/gh-pages/LICENSE.md) */
+
+/* modified from https://github.com/twbs/bootstrap/blob/94b4076dd2efba9af71f0b18d4ee4b163aa9e0dd/docs/assets/css/src/docs.css#L548-L601 */
+
+/* All levels of nav */
+nav[data-toggle='toc'] .nav > li > a {
+ display: block;
+ padding: 4px 20px;
+ font-size: 13px;
+ font-weight: 500;
+ color: #767676;
+}
+nav[data-toggle='toc'] .nav > li > a:hover,
+nav[data-toggle='toc'] .nav > li > a:focus {
+ padding-left: 19px;
+ color: #563d7c;
+ text-decoration: none;
+ background-color: transparent;
+ border-left: 1px solid #563d7c;
+}
+nav[data-toggle='toc'] .nav > .active > a,
+nav[data-toggle='toc'] .nav > .active:hover > a,
+nav[data-toggle='toc'] .nav > .active:focus > a {
+ padding-left: 18px;
+ font-weight: bold;
+ color: #563d7c;
+ background-color: transparent;
+ border-left: 2px solid #563d7c;
+}
+
+/* Nav: second level (shown on .active) */
+nav[data-toggle='toc'] .nav .nav {
+ display: none; /* Hide by default, but at >768px, show it */
+ padding-bottom: 10px;
+}
+nav[data-toggle='toc'] .nav .nav > li > a {
+ padding-top: 1px;
+ padding-bottom: 1px;
+ padding-left: 30px;
+ font-size: 12px;
+ font-weight: normal;
+}
+nav[data-toggle='toc'] .nav .nav > li > a:hover,
+nav[data-toggle='toc'] .nav .nav > li > a:focus {
+ padding-left: 29px;
+}
+nav[data-toggle='toc'] .nav .nav > .active > a,
+nav[data-toggle='toc'] .nav .nav > .active:hover > a,
+nav[data-toggle='toc'] .nav .nav > .active:focus > a {
+ padding-left: 28px;
+ font-weight: 500;
+}
+
+/* from https://github.com/twbs/bootstrap/blob/e38f066d8c203c3e032da0ff23cd2d6098ee2dd6/docs/assets/css/src/docs.css#L631-L634 */
+nav[data-toggle='toc'] .nav > .active > ul {
+ display: block;
+}
diff --git a/docs/bootstrap-toc.js b/docs/bootstrap-toc.js
new file mode 100644
index 0000000..1cdd573
--- /dev/null
+++ b/docs/bootstrap-toc.js
@@ -0,0 +1,159 @@
+/*!
+ * Bootstrap Table of Contents v0.4.1 (http://afeld.github.io/bootstrap-toc/)
+ * Copyright 2015 Aidan Feldman
+ * Licensed under MIT (https://github.com/afeld/bootstrap-toc/blob/gh-pages/LICENSE.md) */
+(function() {
+ 'use strict';
+
+ window.Toc = {
+ helpers: {
+ // return all matching elements in the set, or their descendants
+ findOrFilter: function($el, selector) {
+ // http://danielnouri.org/notes/2011/03/14/a-jquery-find-that-also-finds-the-root-element/
+ // http://stackoverflow.com/a/12731439/358804
+ var $descendants = $el.find(selector);
+ return $el.filter(selector).add($descendants).filter(':not([data-toc-skip])');
+ },
+
+ generateUniqueIdBase: function(el) {
+ var text = $(el).text();
+ var anchor = text.trim().toLowerCase().replace(/[^A-Za-z0-9]+/g, '-');
+ return anchor || el.tagName.toLowerCase();
+ },
+
+ generateUniqueId: function(el) {
+ var anchorBase = this.generateUniqueIdBase(el);
+ for (var i = 0; ; i++) {
+ var anchor = anchorBase;
+ if (i > 0) {
+ // add suffix
+ anchor += '-' + i;
+ }
+ // check if ID already exists
+ if (!document.getElementById(anchor)) {
+ return anchor;
+ }
+ }
+ },
+
+ generateAnchor: function(el) {
+ if (el.id) {
+ return el.id;
+ } else {
+ var anchor = this.generateUniqueId(el);
+ el.id = anchor;
+ return anchor;
+ }
+ },
+
+ createNavList: function() {
+ return $(' ');
+ },
+
+ createChildNavList: function($parent) {
+ var $childList = this.createNavList();
+ $parent.append($childList);
+ return $childList;
+ },
+
+ generateNavEl: function(anchor, text) {
+ var $a = $('');
+ $a.attr('href', '#' + anchor);
+ $a.text(text);
+ var $li = $('');
+ $li.append($a);
+ return $li;
+ },
+
+ generateNavItem: function(headingEl) {
+ var anchor = this.generateAnchor(headingEl);
+ var $heading = $(headingEl);
+ var text = $heading.data('toc-text') || $heading.text();
+ return this.generateNavEl(anchor, text);
+ },
+
+ // Find the first heading level (``, then `
`, etc.) that has more than one element. Defaults to 1 (for `
`).
+ getTopLevel: function($scope) {
+ for (var i = 1; i <= 6; i++) {
+ var $headings = this.findOrFilter($scope, 'h' + i);
+ if ($headings.length > 1) {
+ return i;
+ }
+ }
+
+ return 1;
+ },
+
+ // returns the elements for the top level, and the next below it
+ getHeadings: function($scope, topLevel) {
+ var topSelector = 'h' + topLevel;
+
+ var secondaryLevel = topLevel + 1;
+ var secondarySelector = 'h' + secondaryLevel;
+
+ return this.findOrFilter($scope, topSelector + ',' + secondarySelector);
+ },
+
+ getNavLevel: function(el) {
+ return parseInt(el.tagName.charAt(1), 10);
+ },
+
+ populateNav: function($topContext, topLevel, $headings) {
+ var $context = $topContext;
+ var $prevNav;
+
+ var helpers = this;
+ $headings.each(function(i, el) {
+ var $newNav = helpers.generateNavItem(el);
+ var navLevel = helpers.getNavLevel(el);
+
+ // determine the proper $context
+ if (navLevel === topLevel) {
+ // use top level
+ $context = $topContext;
+ } else if ($prevNav && $context === $topContext) {
+ // create a new level of the tree and switch to it
+ $context = helpers.createChildNavList($prevNav);
+ } // else use the current $context
+
+ $context.append($newNav);
+
+ $prevNav = $newNav;
+ });
+ },
+
+ parseOps: function(arg) {
+ var opts;
+ if (arg.jquery) {
+ opts = {
+ $nav: arg
+ };
+ } else {
+ opts = arg;
+ }
+ opts.$scope = opts.$scope || $(document.body);
+ return opts;
+ }
+ },
+
+ // accepts a jQuery object, or an options object
+ init: function(opts) {
+ opts = this.helpers.parseOps(opts);
+
+ // ensure that the data attribute is in place for styling
+ opts.$nav.attr('data-toggle', 'toc');
+
+ var $topContext = this.helpers.createChildNavList(opts.$nav);
+ var topLevel = this.helpers.getTopLevel(opts.$scope);
+ var $headings = this.helpers.getHeadings(opts.$scope, topLevel);
+ this.helpers.populateNav($topContext, topLevel, $headings);
+ }
+ };
+
+ $(function() {
+ $('nav[data-toggle="toc"]').each(function(i, el) {
+ var $nav = $(el);
+ Toc.init($nav);
+ });
+ });
+})();
diff --git a/docs/index.html b/docs/index.html
index 6936d50..229d914 100644
--- a/docs/index.html
+++ b/docs/index.html
@@ -6,21 +6,21 @@
Installing gmkin and gWidgetsWWW2
update.packages()
"drat") install.packages(
drat::addRepo("jranke")
"gmkin") install.packages(
Usage
+ require(gmkin) gmkin()
options(browser="/usr/bin/firefox")
on linux, or setting the default browser on Windows).
Installation using the devtools package
devtools
package can also install gWidgetsWWW2 and gmkin directly from the respective github repositories:require(devtools)
-install_github("jverzani/gWidgetsWWW2", quick = TRUE)
-install_github("jranke/gmkin", quick = TRUE)
+ require(devtools)"jverzani/gWidgetsWWW2", quick = TRUE)
+ install_github("jranke/gmkin", quick = TRUE) install_github(
quick = TRUE
skips docs, multiple-architecture builds, demos, and vignettes, to make installation as fast and painless as possible.
Simplifying the start of gmkin under Windows
+ require(gmkin)
+ gmkin()
+ while (TRUE) {10)
+ Sys.sleep( }
C:\Users\<your username>\AppData\Roaming\gmkin_start.R
), and create a Windows shortcut, e.g. on the Desktop, with the target "C:\Program Files\R\R-3.x.y\bin\R.exe" --slave -f c:\Users\<your username>\AppData\Roaming\gmkin_start.R
, where R.x.y is the version of R that you would like to use with gmkin. This will start a non-interactive R console that runs the gmkin GUI. In the properties dialogue of the shortcut, you can also set the working directory in which it starts gmkin. Thanks to Jonas Klein of the Umweltbundesamt for the trick with using Sys.sleep in a while loop and for suggesting the AppData directory for the start script.FOCUS_2006_Z
+ FOCUS_2006_Z
Format
@@ -140,19 +148,16 @@
http://esdac.jrc.ec.europa.eu/projects/degradation-kinetics
- @@ -163,7 +168,7 @@ diff --git a/docs/reference/Rplot001.png b/docs/reference/Rplot001.png new file mode 100644 index 0000000..17a3580 Binary files /dev/null and b/docs/reference/Rplot001.png differ diff --git a/docs/reference/UBA_model_gallery.html b/docs/reference/UBA_model_gallery.html index 1e8abd8..4e7cce1 100644 --- a/docs/reference/UBA_model_gallery.html +++ b/docs/reference/UBA_model_gallery.html @@ -10,23 +10,27 @@ - + - + - + + + + + - - + + - + - - + + @@ -42,7 +46,6 @@ GUI to see their definitions. The models use the topology of the models defined by Michael Klein (2010) for DegKinManager, with the exception of the kinetic sorption model." /> - @@ -60,7 +63,7 @@ - +
data("UBA_model_gallery")+
data("UBA_model_gallery")
+data(UBA_model_gallery)
- @@ -170,7 +175,7 @@ diff --git a/docs/reference/gmkin.html b/docs/reference/gmkin.html index fa33831..c4226b3 100644 --- a/docs/reference/gmkin.html +++ b/docs/reference/gmkin.html @@ -6,27 +6,31 @@ -
gWidgetsWWW2
toolkit. — gmkin • gmkingWidgetsWWW2
toolkit.R/gmkin.R
gmkin.Rd
gmkin(script_name, show.log = FALSE)+
gmkin(script_name, show.log = FALSE)
In the R console, you should see a message that the httpd help server is started, if it wasn’t already started before.
In the browser, you should see something like the screenshot below.
The statusbar at the bottom of the gmkin window shows, among others, the working directory that gmkin uses.
@@ -465,18 +480,18 @@ $(document).ready(function () {A gmkin project workspace contains datasets, kinetic models for fitting, and so-called fits, i.e. the results of fitting models to data. The project area also shows the current working directory, where project workspace files are saved using the file extension .gmkinws
.
Once a project has been saved by the user, the project explorer to the left will show it in the project list.
The current state of a project should repeatedly be saved during the work in order to avoid loosing data. This can be achieved by selecting the ‘Project’ tab in the center and pressing the ‘Save project to project file’ button. More conveniently, the keyboard shortcut Shift-F12 (was Ctrl-X in gmkin < 0.6.6) can be used to save the current status of the project.
In the project file management area, datasets and models can be imported from one of the projects in the project list, once it has been selected in the droplist labelled ‘Import from:’.
Once the project has been set up, you may want to minimize the project explorer, especially if you are limited in vertical screen space.
When you select one of the datasets in the dataset explorer to the left, some summary information about the dataset is shown in the center, and the data itself is loaded into the data editor to the right.
When you have added information about the units, or edited the data to the right, you should hit the button ‘Keep changes’. This will add a new entry in the dataset explorer.
@@ -493,10 +508,10 @@ $(document).ready(function () {However, entering data in these fields is a prerequisite for entering a new dataset directly in the dataset editor.
For entering new data manually, there are two possibilities. You can either change the Dataset title of the current dataset and edit the data in the Data editor to the right. Or, if the new data should have a different structure, e.g. different sampling times, observed variables and replicates, click on “New dataset”, edit the dataset title, sampling times, number of replicates and the list of observed variables, and press the button ‘Generate grid for entering kinetic data’, in order to prepare the Data editor to the right.
+For entering new data manually, there are two possibilities. You can either change the Dataset title of the current dataset and edit the data in the Data editor to the right. Or, if the new data should have a different structure, e.g. different sampling times, observed variables and replicates, click on “New dataset”, edit the dataset title, sampling times, number of replicates and the list of observed variables, and press the button ‘Generate grid for entering kinetic data’, in order to prepare the Data editor to the right.
For sampling times and short names of the relevant compounds, a comma separated list must be entered, with a space after each comma. An example of filling out the respective fields is shown below.
After that, the actual measurements can be entered into the Data editor to the right, in the column ‘value’.
@@ -512,12 +527,12 @@ $(document).ready(function () {This produces a text file with comma separated values in the current working directory of R.
Loading this text file into gmkin using the “Browse” and “Upload” buttons
results in an an import configuration area showing the top lines of the imported file, and giving the possibility to chose the import options matching the file format.
In the import configuration area, the following options can be specified. In the field “Comment lines”, the number of lines in the beginning of the file that should be ignored can be specified.
@@ -526,12 +541,12 @@ $(document).ready(function () {As the “Decimal” separator, comma “,” or period “.” can be selected.
In the next line, it can be specified if the data are in wide or in long format. If in wide format, the only option left to specify is the title of the column containing the sampling times. If the data is in long format, the column headings specifying the columns containing the observed variables (default is “name”), the sampling times (default is “time”), the observed values (default is “value”) and, if present in the data, the relative errors (default is “err”) can be adapted. The default settings appearing if the long format is selected are shown below.
In our example we have data in the wide format, and after adapting the “Separator” to a comma, we can press the button “Import using options specified below”, and the data should be imported. If successful, the dataset editor should show the sampling times and the names of the observed variables, and the Data editor should show the imported data in a grid for further editing or specifying overrides.
After adapting the dataset title and possibly the units, press “Keep changes” to save the dataset in the current workspace. Again, you need to save the project file in order to be able to use the dataset that you created in a future gmkin session.
@@ -541,11 +556,11 @@ $(document).ready(function () {The following screenshot shows the model editor for the model ‘Z.2a.ff’ from the project workspace ‘FOCUS_2006_Z’ provided with the gmkin package.
In the first line the name of the model can be edited. You can also specify “min” or “max” for minimum or maximum use of formation fractions. Maximum use of formation fractions means that the differential equations in the degradation model are formulated using formation fractions. When you specify “min”, then formation fractions are only used for the parent compound when you use the FOMC, DFOP or the HS model for it.
-Pressing “Add observed variable” adds a line in the array of state variable specifications below. The observed variables to be added are usually transformation products (usually termed metabolites), but can also be the parent compound in a different compartment (e.g. “parent_sediment”).
+Pressing “Add observed variable” adds a line in the array of state variable specifications below. The observed variables to be added are usually transformation products (usually termed metabolites), but can also be the parent compound in a different compartment (e.g. “parent_sediment”).
Only observed variable names that occur in previously defined datasets or models can be selected. For any observed variable other than the first one, only the SFO or the SFORB model can be selected. For each observed variables, a comma separated list of target variables can be specified. In addition, a pathway to the sink compartment can be selected. If too many observed variables have been added, complete lines can be removed from the model definition by pressing the button “Remove observed variable”.
If the model definition is supposedly correct, press “Keep changes” to make it possible to select it for fitting in the listing of models to the left.
If a dataset and a kinetic model are selected, the button “Configure fit” below in the Explorer window “Configuration” becomes active. Pressing it opens the “Configuration” area in the center and loads the dataset into the Data editor to the right, for viewing only.
After pressing the button ‘Plot unoptimised’, the data and the kinetic model with default starting parameters are plotted in the ‘Plot’ area to the right.
In the right part of the fit configuration area, the file format can be chosen, the legend can be turned off, and the observed variables for which the data and the model fit should be plotted can be selected as shown below.
On systems running the Windows operating system, the enhanced metafile (emf) format can be additionally chosen. Changing the file format for plotting will also change the extension of the proposed filename for saving the plot.
@@ -604,7 +619,7 @@ Optimisation successfully terminated.If plotting of the fitting progress was selected, a new separate graphics window should either pop up, or a graphics window previously started for this purpose will be reused.
If your screen size allows for it, you can arrange the R plotting window and the R console in a way that you can see everything at the same time:
Once a fit has successfully been completed, the most important results are shown in several tables in the “Result” area as shown below.
The detailed summary can be accessed below these tables.
The complete summary can be saved into a text file by specifying a suitable file name and pressing the button “Save summary”.
Whenever a new fit has been configured or a run of a fit has been completed, the plotting area is updated with the abovementioned plot of the data and the current model solution.
-In addition, a confidence interval plot is shown below this conventional plot. In case a fit has been run and confidence intervals were successfully calculated for the fit ( i.e. if the model was not overparameterised and no other problems occurred), the confidence intervals are graphically displayed as bars as shown below.
+In addition, a confidence interval plot is shown below this conventional plot. In case a fit has been run and confidence intervals were successfully calculated for the fit (i.e. if the model was not overparameterised and no other problems occurred), the confidence intervals are graphically displayed as bars as shown below.