Time course of 2,4,5-trichlorophenoxyacetic acid, and the corresponding 2,4,5-trichlorophenol and 2,4,5-trichloroanisole as recovered in diethylether extracts.

mccall81_245T

Format

A dataframe containing the following variables.

name

the name of the compound observed. Note that T245 is used as an acronym for 2,4,5-T. T245 is a legitimate object name in R, which is necessary for specifying models using mkinmod.

time

a numeric vector containing sampling times in days after treatment

value

a numeric vector containing concentrations in percent of applied radioactivity

soil

a factor containing the name of the soil

Source

McCall P, Vrona SA, Kelley SS (1981) Fate of uniformly carbon-14 ring labeled 2,4,5-Trichlorophenoxyacetic acid and 2,4-dichlorophenoxyacetic acid. J Agric Chem 29, 100-107 http://dx.doi.org/10.1021/jf00103a026

Examples

SFO_SFO_SFO <- mkinmod(T245 = list(type = "SFO", to = "phenol"), phenol = list(type = "SFO", to = "anisole"), anisole = list(type = "SFO"))
#> Successfully compiled differential equation model from auto-generated C code.
fit.1 <- mkinfit(SFO_SFO_SFO, subset(mccall81_245T, soil == "Commerce"), quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
summary(fit.1)$bpar
#> Warning: Could not calculate correlation; no covariance matrix
#> Estimate se_notrans t value Pr(>t) Lower Upper #> T245_0 1.038550e+02 NA NA NA NA NA #> k_T245_sink 1.636106e-02 NA NA NA NA NA #> k_T245_phenol 2.700936e-02 NA NA NA NA NA #> k_phenol_sink 1.788604e-10 NA NA NA NA NA #> k_phenol_anisole 4.050581e-01 NA NA NA NA NA #> k_anisole_sink 6.678742e-03 NA NA NA NA NA #> sigma 2.514628e+00 NA NA NA NA NA
endpoints(fit.1)
#> $ff #> T245_sink T245_phenol phenol_sink phenol_anisole anisole_sink #> 3.772401e-01 6.227599e-01 4.415672e-10 1.000000e+00 1.000000e+00 #> #> $SFORB #> logical(0) #> #> $distimes #> DT50 DT90 #> T245 15.982025 53.09114 #> phenol 1.711229 5.68458 #> anisole 103.784092 344.76329 #>
# k_phenol_sink is really small, therefore fix it to zero fit.2 <- mkinfit(SFO_SFO_SFO, subset(mccall81_245T, soil == "Commerce"), parms.ini = c(k_phenol_sink = 0), fixed_parms = "k_phenol_sink", quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
summary(fit.2)$bpar
#> Warning: Could not calculate correlation; no covariance matrix
#> Estimate se_notrans t value Pr(>t) Lower Upper #> T245_0 1.038550e+02 NA NA NA NA NA #> k_T245_sink 1.636106e-02 NA NA NA NA NA #> k_T245_phenol 2.700936e-02 NA NA NA NA NA #> k_phenol_anisole 4.050581e-01 NA NA NA NA NA #> k_anisole_sink 6.678742e-03 NA NA NA NA NA #> sigma 2.514628e+00 NA NA NA NA NA
endpoints(fit.1)
#> $ff #> T245_sink T245_phenol phenol_sink phenol_anisole anisole_sink #> 3.772401e-01 6.227599e-01 4.415672e-10 1.000000e+00 1.000000e+00 #> #> $SFORB #> logical(0) #> #> $distimes #> DT50 DT90 #> T245 15.982025 53.09114 #> phenol 1.711229 5.68458 #> anisole 103.784092 344.76329 #>
plot_sep(fit.2)