Lepi: Insect counts of 12 Species
In BSagri: Safety Assessment in Agricultural Field Trials

Description Usage Format Source Examples

Simulated data, inpired by a real field investigating the potential impact of genetically modified crop on several insect species belonging to the same order. The trial was designed as a randomized complete block design with 8 blocks (Block), and a total of 24 plots. In each block, three treatments (Treatment) were randomized: a conventional variety treated with insecticides (Ins), a genetically modified variety (GM) without insecticide treatment, and the near-isogenic variety (Iso) the to genetically modified variety, without insecticide treatment. Individuals were counted (after classification to the species level) in two different dates in each year of the trial, where the the second date was of higher importance for assessment of impacts of GM variety on non-target species. In total 12 Species were observed during the trial.

1	data(Lepi)

A data frame with 144 observations on the following 17 variables.

Year: a numeric vector, the year 1, 2, 3
Date: a numeric vector, 1 and 2 separating the 2 sampling date in each year
Block: a numeric vector, with values 1-8, indicator variable for the 8 blocks
Treatment: a factor with three levels identifying the varieties: GM is the genetically modified variety, Ins the conventional variety with insecticide treatment and Iso the near isognic line without insecticide treatment
Plot: a factor with 24 levels, identifying the individual plots
Sp1: counts of taxon 1
Sp2: counts of taxon 2
Sp3: counts of taxon 3
Sp4: counts of taxon 4
Sp5: counts of taxon 5
Sp6: counts of taxon 6
Sp7: counts of taxon 7
Sp8: counts of taxon 8
Sp9: counts of taxon 9
Sp10: counts of taxon 10
Sp11: counts of taxon 11
Sp12: counts of taxon 12

Simulated data.

data(Lepi)

str(Lepi)

summary(Lepi)

SPEC<-names(Lepi)[-(1:5)]

# Occurrence

occur<-lapply(X=Lepi[,SPEC], FUN=function(x){length(which(x>0))})

unlist(occur)

# Species with reasonable occurence in the whole data:

SPEC2<-SPEC[c(1,2,3,6,8,9,11)]

pairs(Lepi[,SPEC2])

# 


layout(matrix(1:2, ncol=1 ))
par(mar=c(2,8,2,1))

boxplot(Sp2 ~ Treatment*Year, data=Lepi, main="Species 2",
 las=1, horizontal=TRUE, col=c("red","white","white"))

boxplot(Sp3 ~ Treatment*Year, data=Lepi, main="Species 3",
 las=1, horizontal=TRUE, col=c("red","white","white"))