clementines               package:ade4               R Documentation

_F_r_u_i_t _P_r_o_d_u_c_t_i_o_n

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

     The 'clementines' is a data set containing the fruit production of
     20 clementine trees during 15 years.

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

     data(clementines)

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

     A data frame with 15 rows and 20 columns

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

     Tisn-Agostini, D. (1988) _Description par analyse en composantes
     principales de l'volution de  a production du clmentinier en
     association avec 12 types de porte-greffe_. Rapport technique, DEA
     Analyse et modlisation des systmes biologiques, Universit Lyon
     1.

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

     data(clementines)
     op <- par(no.readonly = TRUE) 
     par(mfrow = c(5,4)) ; par(mar = c(2,2,1,1))
     for (i in 1:20) {
         w0 <- 1:15 ; plot(w0, clementines[,i], type = "b")
         abline(lm(clementines[,i] ~ w0))
     }
     par(op)

     pca1 <- dudi.pca(clementines, scan = FALSE)
     s.corcircle(pca1$co, clab = 0.75)

     barplot(pca1$li[,1])

     op <- par(no.readonly = TRUE) 
     par(mfrow = c(5,4)) ; par(mar = c(2,2,1,1))
     clem0 <- pca1$tab
     croi <- 1:15
     alter <- c(rep(c(1,-1),7),1)
     for (i in 1:20) {
         y <- clem0[,i]
         plot(w0, y, type = "b", ylim = c(-2,2))
         z <- predict(lm(clem0[,i] ~ croi*alter))
         points(w0, z, pch = 20, cex = 2)
         for(j in 1:15) segments (j,y[j],j,z[j])
     }
     par(op)
     par(mfrow = c(1,1))

