CO2                 package:datasets                 R Documentation

_C_a_r_b_o_n _D_i_o_x_i_d_e _u_p_t_a_k_e _i_n _g_r_a_s_s _p_l_a_n_t_s

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

     The 'CO2' data frame has 84 rows and 5 columns of data from an
     experiment on the cold tolerance of the grass species _Echinochloa
     crus-galli_.

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

     CO2

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

     This data frame contains the following columns:

     _P_l_a_n_t an ordered factor with levels 'Qn1' < 'Qn2' < 'Qn3' < ... <
          'Mc1' giving a unique identifier for each plant.

     _T_y_p_e a factor with levels 'Quebec'  'Mississippi' giving the
          origin of the plant

     _T_r_e_a_t_m_e_n_t a factor with levels 'nonchilled'  'chilled'

     _c_o_n_c a numeric vector of ambient carbon dioxide concentrations
          (mL/L).

     _u_p_t_a_k_e a numeric vector of carbon dioxide uptake rates (umol/m^2
          sec).


_D_e_t_a_i_l_s:

     The CO2 uptake of six plants from Quebec and six plants from
     Mississippi was measured at several levels of ambient CO2
     concentration.  Half the plants of each type were chilled
     overnight before the experiment was conducted.

_S_o_u_r_c_e:

     Potvin, C., Lechowicz, M. J. and Tardif, S. (1990) "The
     statistical analysis of ecophysiological response curves obtained
     from experiments involving repeated measures", _Ecology_, *71*,
     1389-1400.

     Pinheiro, J. C. and Bates, D. M. (2000) _Mixed-effects Models in S
     and S-PLUS_, Springer.

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

     require(stats); require(graphics)
     coplot(uptake ~ conc | Plant, data = CO2, show.given = FALSE, type = "b")
     ## fit the data for the first plant
     fm1 <- nls(uptake ~ SSasymp(conc, Asym, lrc, c0),
        data = CO2, subset = Plant == 'Qn1')
     summary(fm1)
     ## fit each plant separately
     fmlist <- list()
     for (pp in levels(CO2$Plant)) {
       fmlist[[pp]] <- nls(uptake ~ SSasymp(conc, Asym, lrc, c0),
           data = CO2, subset = Plant == pp)
     }
     ## check the coefficients by plant
     sapply(fmlist, coef)

