Feeding                package:BSagri                R Documentation

_P_u_p_a_t_i_o_n _a_n_d _H_a_t_c_h_i_n_g _r_a_t_e _i_n _a _f_e_e_d_i_n_g _e_x_p_e_r_i_m_e_n_t _w_i_t_h _f_o_u_r _v_a_r_i_e_t_i_e_s

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

     Larvae of a non-target organism were fed with plant material
     derived from a novel variety(Novum), material from three standard
     varieties (NStandard: the standard variety most similar to Novum,
     and two additional standard varieties S1 and S2). Objective was to
     assess the impact of Novum on the pupation and hatching rate of an
     animal that potentially feeds on plant material compared to
     accepted standard varieties.

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

     data(Feeding)

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

     A data frame with 32 observations on the following 5 variables.

     '_R_e_p' a factor with 32 levels indexing the 32 replications 

     '_V_a_r_i_e_t_y' a factor with 4 levels: 'S1' and 'S2' are two standard
          varieties, 'Novum' is a novel variety, and 'NStandard' is the
          standard variety most similar to 'Novum' 

     '_T_o_t_a_l' the total number of animals in each experimental unit

     '_P_u_p_a_t_i_n_g' number of individuals pupating in each unit, the others
          died

     '_H_a_t_c_h_i_n_g' number of individuals hatching from the pupae

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

     ...

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

     data(Feeding)

     # Larval mortality:

     Feeding$Lmort <- Feeding$Total - Feeding$Pupating

     # Pupae mortality

     Feeding$Pmort <- Feeding$Pupating - Feeding$Hatching

     # Total mortality

     Feeding$Tmort <- Feeding$Total - Feeding$Hatching

     fit1<-glm(cbind(Pupating,Lmort)~Variety,data=Feeding, family=quasibinomial)
     anova(fit1, test="F")

     fit2<-glm(cbind(Hatching,Pmort)~Variety,data=Feeding, family=quasibinomial)
     anova(fit2, test="F")

     fit3<-glm(cbind(Total,Tmort)~Variety,data=Feeding, family=quasibinomial)
     anova(fit3, test="F")

