Pupation and Hatching rate in a feeding experiment with four varieties

Description

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.

Usage

data(Feeding)

Format

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

Rep

a factor with 32 levels indexing the 32 replications

Variety

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

Total

the total number of animals in each experimental unit

Pupating

number of individuals pupating in each unit, the others died

Hatching

number of individuals hatching from the pupae

Examples

data(Feeding)

# Larval mortality:

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

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

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(BSagri)
Loading required package: gamlss
Loading required package: splines
Loading required package: gamlss.data
Loading required package: gamlss.dist
Loading required package: MASS
Loading required package: nlme
Loading required package: parallel
 **********   GAMLSS Version 4.4-0  ********** 
For more on GAMLSS look at http://www.gamlss.org/
Type gamlssNews() to see new features/changes/bug fixes.

Loading required package: multcomp
Loading required package: mvtnorm
Loading required package: survival

Attaching package: 'survival'

The following object is masked from 'package:gamlss.data':

    leukemia

Loading required package: TH.data

Attaching package: 'TH.data'

The following object is masked from 'package:MASS':

    geyser

Loading required package: MCPAN
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/BSagri/Feeding.Rd_%03d_medium.png", width=480, height=480)
> ### Name: Feeding
> ### Title: Pupation and Hatching rate in a feeding experiment with four
> ###   varieties
> ### Aliases: Feeding
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(Feeding)
> 
> # Larval mortality:
> 
> Feeding$Lmort <- Feeding$Total - Feeding$Pupating
> 
> fit1<-glm(cbind(Pupating,Lmort)~Variety,data=Feeding, family=quasibinomial)
> anova(fit1, test="F")
Analysis of Deviance Table

Model: quasibinomial, link: logit

Response: cbind(Pupating, Lmort)

Terms added sequentially (first to last)


        Df Deviance Resid. Df Resid. Dev      F Pr(>F)
NULL                       31     59.623              
Variety  3   3.3623        28     56.260 0.6478 0.5909
> 
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>