SCOC             package:marelacTeaching             R Documentation

_a _S_e_d_i_m_e_n_t _C_o_m_m_u_n_i_t_y _O_x_y_g_e_n _C_o_n_s_u_m_p_t_i_o_n (_S_C_O_C) _d_a_t_a_s_e_t

_D_e_s_c_r_i_p_t_i_o_n:

     This literature dataset, compiled by Andersson et al. (2004)
     contains 584 measurements of sediment community oxygen consumption
     rates, as a function of water depth, and performed in deep-water
     sediments, either by in situ incubations or via modelling of
     oxygen microprofiles.

     It is used in the book to demonstrate how one can obtain
     order-of-magnitude estimates of model parameters (i.c. sediment
     oxygen consumption rate, a measure of deposition flux) by
     performing log-log regression with water depth.

_U_s_a_g_e:

     SCOC

_F_o_r_m_a_t:

     a dataframe with 584 rows, and with following columns: 
      Depth.m, the water depth at which the measurement was performed. 
      SCOC.mmol/m2/d, the oxygen consumption rate of the sediment,
     [mmolO2/m2/d]

_A_u_t_h_o_r(_s):

     Karline Soetaert <k.soetaert@nioo.knaw.nl>

_R_e_f_e_r_e_n_c_e_s:

     Andersson, H., Wijsman, J., Herman, P., Middelburg, J., Soetaert,
     K., Heip, C., 2004. Respiration patterns in the deep ocean.
     Geophysical Research Letters 31, LO3304.

_S_e_e _A_l_s_o:

     'Zoogrowth', a dataset containing zooplankton maximal growth rates 
      see the paper of Andersson et al. for a description of the
     original literature sources of this dataset

_E_x_a_m_p_l_e_s:

     plot(SCOC[,1],SCOC[,2],log="xy",xlab="water depth, m",ylab="" ,
          main="SCOC, mmol O2/m2/d",pch=16,xaxt="n",yaxt="n",cex.main=1)

     axis(1,at=c(0.5,5,50,500,5000),labels=c("0.5","5","50","500","5000"))
     axis(2,at=c(0.1,1,10,100),labels=c("0.1","1","10","100"))

     ll <- lm(log(SCOC[,2])~ log(SCOC[,1]))
     rr <- summary(ll)$r.squared
     A  <- exp(coef(ll)[1])
     B  <- (coef(ll)[2])
     curve(A*x^B,add=TRUE,lwd=2)
     AA <- round(A*100)/100
     BB <- round(B*100)/100
     expr <- substitute(y==A*x^B,list(A=AA,B=BB))
     text(1,.1,expr,adj=0)
     expr2 <- substitute(r^2==rr,list(rr=round(rr*100)/100))
     text(1,0.04,expr2, adj=0)

