class: center, middle, inverse, title-slide # Technical tests with xaringan ### Thomas Petzoldt ### 2021-06-24 --- ### IMPORTANT NOTE This slide set is for testing only. It is made with the **R** package **xaringan**. It works with **Chrome** and recent versions of **Edge**, wheres scrolling may have problems with Firefox. The source code of these slides is found [here](https://github.com/dynamic-R/hacking-limnology/tree/main/docs). More about **xaringan** can be found in the wiki: https://github.com/yihui/xaringan/wiki . The name of the repo points to an upcoming workshop about open data in aquatic sciences https://aquaticdatasciopensci.github.io/ --- ## A Scrollable Code .scrollable[ ``` Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Duis autem vel eum iriure dolor in hendrerit in vulputate velit esse molestie consequat, vel illum dolore eu feugiat nulla facilisis at vero eros et accumsan et iusto odio dignissim qui blandit praesent luptatum zzril delenit augue duis dolore te feugait nulla facilisi. Lorem ipsum dolor sit amet, consectetuer adipiscing elit, sed diam nonummy nibh euismod tincidunt ut laoreet dolore magna aliquam erat volutpat. Ut wisi enim ad minim veniam, quis nostrud exerci tation ullamcorper suscipit lobortis nisl ut aliquip ex ea commodo consequat. Duis autem vel eum iriure dolor in hendrerit in vulputate velit esse molestie consequat, vel illum dolore eu feugiat nulla facilisis at vero eros et accumsan et iusto odio dignissim qui blandit praesent luptatum zzril delenit augue duis dolore te feugait nulla facilisi. Nam liber tempor cum soluta nobis eleifend option congue nihil imperdiet doming id quod mazim placerat facer possim assum. Lorem ipsum dolor sit amet, consectetuer adipiscing elit, sed diam nonummy nibh euismod tincidunt ut laoreet dolore magna aliquam erat volutpat. Ut wisi enim ad minim veniam, quis nostrud exerci tation ullamcorper suscipit lobortis nisl ut aliquip ex ea commodo consequat. ``` ] --- ## Scrollable with Images .scrollable[ ```r library("deSolve") library("rootSolve") chemostat <- function(time, init, parms) { with(as.list(c(init, parms)), { mu <- mumax * P/(kp + P) # Monod equation dAlg <- mu * Alg - D * Alg dP <- D *(P0 - P) - 1/Y * mu * Alg list(c(dAlg, dP), mu=mu) }) } parms <- c( mumax = 0.5, # 1/d kp = 0.01, # half saturation constant (mg/L) Y = 41, # yield coefficient (stoichiometric C:P ratio) D = 0.1, # 1/d P0 = 0.05 # P in inflow (mg/L) ) times <- seq(0, 40, 0.1) # (d) init <- c(Alg=0.01, P=0.05) # Phytoplankton C and Phosphorus P (mg/L) ## ============================================================================= ## Dynamic simulation ## ============================================================================= out <- ode(init, times, chemostat, parms) plot(out) ## ============================================================================= ## Steady state solution ## ============================================================================= state <- data.frame( D = seq(0, 0.6, length.out = 100), X = 0, S = 0 ) for (i in 1:nrow(state)) { parms["D"] <- state$D[i] times <- c(0, Inf) out <- runsteady(init, times, chemostat, parms) state[i, 2:3] <- out$y } par(mfrow = c(3, 1)) ``` <!-- --> ```r plot(S ~ D, data = state, type = "l") plot(X ~ D, data = state, type = "l") plot(S * X ~ D, data = state, type = "l") ``` <!-- --> ] --- class: page-font-large * Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua. * At vero eos et accusam et justo duo dolores et ea rebum. Stet clita kasd gubergren, no sea takimata sanctus est Lorem ipsum dolor sit amet. Lorem ipsum dolor sit amet, consetetur sadipscing elitr, sed diam nonumy eirmod tempor invidunt ut labore et dolore magna aliquyam erat, sed diam voluptua.