Skip to content

MLE GPP_fun='satlightq10temp' : issue with GPP (GPP=0, no initialization) #405

Open
@sophieguillon

Description

@sophieguillon

Hello

I am new in using streamMetabolizer package, for both lake and small stream metabolism calculations.

I am using MLE models right now, and I have a problem with GPP_fun='satlightq10temp' option.
It leads to GPP = 0, and it seems that there is no alpha and Pmax parameters in initialization and in output, as when using satlight only.
The temperature dependence is ok for ER (only ER20 parameter)

Below the code and datafile I use for this test.

I thank you in advance for your answer

Sophie Guillon

datasub_github.txt

datasub=read.csv('datasub_github.txt',header=TRUE,sep="",dec = ".")

classic_specs <- specs( mm_name('mle', GPP_fun='linlight', ER_fun='constant') )
mm_classic <- streamMetabolizer::metab(classic_specs, data=datasub)

mm_classic@metab_daily

mle_specs <- specs(mm_name('mle', GPP_fun='satlightq10temp', ER_fun='q10temp'))
mm_classicTEMP <- streamMetabolizer::metab(mle_specs, data= datasub)

mm_classicTEMP@metab_daily
mm_classicTEMP@specs

mle_specs <- specs(mm_name('mle', GPP_fun='satlight', ER_fun='constant'))
mm_saturator <- streamMetabolizer::metab(mle_specs, data= datasub)

mm_saturator@specs

Session information

─ Session info ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
setting value
version R version 4.0.3 (2020-10-10)
os Ubuntu 18.04.5 LTS
system x86_64, linux-gnu
ui RStudio
language fr_FR
collate fr_FR.UTF-8
ctype fr_FR.UTF-8
tz GMT
date 2021-01-26

─ Packages ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────
package * version date lib source
assertthat 0.2.1 2019-03-21 [1] CRAN (R 4.0.3)
callr 3.5.1 2020-10-13 [1] CRAN (R 4.0.3)
cli 2.2.0 2020-11-20 [1] CRAN (R 4.0.3)
codetools 0.2-16 2018-12-24 [4] CRAN (R 4.0.0)
colorspace 2.0-0 2020-11-11 [1] CRAN (R 4.0.3)
crayon 1.3.4 2017-09-16 [1] CRAN (R 4.0.3)
crosstalk 1.1.1 2021-01-12 [1] CRAN (R 4.0.3)
curl 4.3 2019-12-02 [1] CRAN (R 4.0.3)
data.table 1.13.6 2020-12-30 [1] CRAN (R 4.0.3)
desc 1.2.0 2018-05-01 [1] CRAN (R 4.0.3)
deSolve 1.28 2020-03-08 [1] CRAN (R 4.0.3)
devtools 2.3.2 2020-09-18 [1] CRAN (R 4.0.3)
digest 0.6.27 2020-10-24 [1] CRAN (R 4.0.3)
dplyr * 1.0.3 2021-01-15 [1] CRAN (R 4.0.3)
ellipsis 0.3.1 2020-05-15 [1] CRAN (R 4.0.3)
fansi 0.4.2 2021-01-15 [1] CRAN (R 4.0.3)
farver 2.0.3 2020-01-16 [1] CRAN (R 4.0.3)
fs 1.5.0 2020-07-31 [1] CRAN (R 4.0.3)
generics 0.1.0 2020-10-31 [1] CRAN (R 4.0.3)
ggplot2 * 3.3.3 2020-12-30 [1] CRAN (R 4.0.3)
glue 1.4.2 2020-08-27 [1] CRAN (R 4.0.3)
gridExtra 2.3 2017-09-09 [1] CRAN (R 4.0.3)
gtable 0.3.0 2019-03-25 [1] CRAN (R 4.0.3)
htmltools 0.5.1 2021-01-12 [1] CRAN (R 4.0.3)
htmlwidgets 1.5.3 2020-12-10 [1] CRAN (R 4.0.3)
httr 1.4.2 2020-07-20 [1] CRAN (R 4.0.3)
inline 0.3.17 2020-12-01 [1] CRAN (R 4.0.3)
jsonlite 1.7.2 2020-12-09 [1] CRAN (R 4.0.3)
labeling 0.4.2 2020-10-20 [1] CRAN (R 4.0.3)
LakeMetabolizer * 1.5.3 2021-01-19 [1] Github (GLEON/LakeMetabolizer@5e42781)
lazyeval 0.2.2 2019-03-15 [1] CRAN (R 4.0.3)
lifecycle 0.2.0 2020-03-06 [1] CRAN (R 4.0.3)
loo 2.4.1 2020-12-09 [1] CRAN (R 4.0.3)
lubridate 1.7.9.2 2020-11-13 [1] CRAN (R 4.0.3)
magrittr 2.0.1 2020-11-17 [1] CRAN (R 4.0.3)
matrixStats 0.57.0 2020-09-25 [1] CRAN (R 4.0.3)
memoise 1.1.0 2017-04-21 [1] CRAN (R 4.0.3)
munsell 0.5.0 2018-06-12 [1] CRAN (R 4.0.3)
pillar 1.4.7 2020-11-20 [1] CRAN (R 4.0.3)
pkgbuild 1.2.0 2020-12-15 [1] CRAN (R 4.0.3)
pkgconfig 2.0.3 2019-09-22 [1] CRAN (R 4.0.3)
pkgload 1.1.0 2020-05-29 [1] CRAN (R 4.0.3)
plotly * 4.9.3 2021-01-10 [1] CRAN (R 4.0.3)
plyr 1.8.6 2020-03-03 [1] CRAN (R 4.0.3)
prettyunits 1.1.1 2020-01-24 [1] CRAN (R 4.0.3)
processx 3.4.5 2020-11-30 [1] CRAN (R 4.0.3)
ps 1.5.0 2020-12-05 [1] CRAN (R 4.0.3)
purrr 0.3.4 2020-04-17 [1] CRAN (R 4.0.3)
R6 2.5.0 2020-10-28 [1] CRAN (R 4.0.3)
Rcpp 1.0.6 2021-01-15 [1] CRAN (R 4.0.3)
RcppParallel 5.0.2 2020-06-24 [1] CRAN (R 4.0.3)
remotes 2.2.0 2020-07-21 [1] CRAN (R 4.0.3)
rLakeAnalyzer * 1.11.4.1 2019-06-09 [1] CRAN (R 4.0.3)
rlang 0.4.10 2020-12-30 [1] CRAN (R 4.0.3)
rprojroot 2.0.2 2020-11-15 [1] CRAN (R 4.0.3)
rstan 2.21.2 2020-07-27 [1] CRAN (R 4.0.3)
rstudioapi 0.13 2020-11-12 [1] CRAN (R 4.0.3)
scales 1.1.1 2020-05-11 [1] CRAN (R 4.0.3)
sessioninfo 1.1.1 2018-11-05 [1] CRAN (R 4.0.3)
StanHeaders 2.21.0-7 2020-12-17 [1] CRAN (R 4.0.3)
streamMetabolizer * 0.11.4 2021-01-19 [1] Github (4a7ba9c)
testthat 3.0.1 2020-12-17 [1] CRAN (R 4.0.3)
tibble 3.0.5 2021-01-15 [1] CRAN (R 4.0.3)
tidyr * 1.1.2 2020-08-27 [1] CRAN (R 4.0.3)
tidyselect 1.1.0 2020-05-11 [1] CRAN (R 4.0.3)
unitted 0.2.9 2021-01-19 [1] Github (appling/unitted@d1f1172)
usethis 2.0.0 2020-12-10 [1] CRAN (R 4.0.3)
V8 3.4.0 2020-11-04 [1] CRAN (R 4.0.3)
vctrs 0.3.6 2020-12-17 [1] CRAN (R 4.0.3)
viridisLite 0.3.0 2018-02-01 [1] CRAN (R 4.0.3)
withr 2.4.0 2021-01-16 [1] CRAN (R 4.0.3)
yaml 2.2.1 2020-02-01 [1] CRAN (R 4.0.3)

[1] /home/sguillon/R/x86_64-pc-linux-gnu-library/4.0
[2] /usr/local/lib/R/site-library
[3] /usr/lib/R/site-library
[4] /usr/lib/R/library

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions