Orange package:datasets R Documentation _G_r_o_w_t_h _o_f _O_r_a_n_g_e _T_r_e_e_s _D_e_s_c_r_i_p_t_i_o_n: The 'Orange' data frame has 35 rows and 3 columns of records of the growth of orange trees. _U_s_a_g_e: Orange _F_o_r_m_a_t: This data frame contains the following columns: _T_r_e_e an ordered factor indicating the tree on which the measurement is made. The ordering is according to increasing maximum diameter. _a_g_e a numeric vector giving the age of the tree (days since 1968/12/31) _c_i_r_c_u_m_f_e_r_e_n_c_e a numeric vector of trunk circumferences (mm). This is probably "circumference at breast height", a standard measurement in forestry. _S_o_u_r_c_e: Draper, N. R. and Smith, H. (1998), _Applied Regression Analysis (3rd ed)_, Wiley (exercise 24.N). 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(circumference ~ age | Tree, data = Orange, show.given = FALSE) fm1 <- nls(circumference ~ SSlogis(age, Asym, xmid, scal), data = Orange, subset = Tree == 3) plot(circumference ~ age, data = Orange, subset = Tree == 3, xlab = "Tree age (days since 1968/12/31)", ylab = "Tree circumference (mm)", las = 1, main = "Orange tree data and fitted model (Tree 3 only)") age <- seq(0, 1600, length.out = 101) lines(age, predict(fm1, list(age = age)))