beavers                package:stats                R Documentation

_B_o_d_y _T_e_m_p_e_r_a_t_u_r_e _S_e_r_i_e_s _o_f _T_w_o _B_e_a_v_e_r_s

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

     Reynolds (1994) describes a small part of a study of the long-term
     temperature dynamics of beaver _Castor canadensis_ in
     north-central Wisconsin.  Body temperature was measured by
     telemetry every 10 minutes for four females, but data from a one
     period of less than a day for each of two animals is used there.

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

     data(beavers)

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

     The 'beaver1' data frame has 114 rows and 4 columns on body
     temperature measurements at 10 minute intervals.

     The 'beaver2' data frame has 100 rows and 4 columns on body
     temperature measurements at 10 minute intervals.

     The variables are as follows:

     _d_a_y Day of observation (in days since the beginning of 1990),
          December 12-13 ('beaver1') and November 3-4 ('beaver2').

     _t_i_m_e Time of observation, in the form '0330' for 3:30am

     _t_e_m_p Measured body temperature in degrees Celsius.

     _a_c_t_i_v Indicator of activity outside the retreat.

_N_o_t_e:

     The observation at 22:20 is missing in 'beaver1'.

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

     P. S. Reynolds (1994) Time-series analyses of beaver body
     temperatures.  Chapter 11 of Lange, N., Ryan, L., Billard, L.,
     Brillinger, D., Conquest, L. and Greenhouse, J. eds (1994) _Case
     Studies in Biometry._ New York: John Wiley and Sons.

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

     data(beavers)

     (yl <- range(beaver1$temp, beaver2$temp))

     beaver.plot <- function(bdat, ...) {
       nam <- deparse(substitute(bdat))  
       attach(bdat)
       # Hours since start of day:
       hours <- time %/% 100 + 24*(day - day[1]) + (time %% 100)/60
       plot (hours, temp, type = "l", ...,
             main = paste(nam, "body temperature"))
       abline(h = 37.5, col = "gray", lty = 2)
       is.act <- activ == 1
       points(hours[is.act], temp[is.act], col = 2, cex = .8)
     }
     op <- par(mfrow = c(2,1), mar = c(3,3,4,2), mgp = .9* 2:0)
      beaver.plot(beaver1, ylim = yl)
      beaver.plot(beaver2, ylim = yl)
     par(op)

