rm(list=ls()) #clears memory
setwd("/rstuff") #sets path

#miceants
impap<-round(rnorm(8, mean=14.5, sd=4),digits=0)
impaa<-round(rnorm(8, mean=30, sd=4),digits=0)
imaap<-round(rnorm(8, mean=20, sd=4),digits=0)
imaaa<-round(rnorm(8, mean=27, sd=4),digits=0)

mpap<-c(impap, rep(NA,20-length(impap)))
mpaa<-c(impaa, rep(NA,20-length(impaa)))
maap<-c(imaap, rep(NA,20-length(imaap)))
maaa<-c(imaaa, rep(NA,20-length(imaaa)))

#brown snake
dd<-rnorm(20, mean=0.41, sd=.08)
ld<-rnorm(20, mean=0.3, sd=.08)
dw<-rnorm(20, mean=0.19, sd=.08)
lw<-rnorm(20, mean=0.25, sd=.08)


#rotifers
inf0<-round(rnorm(10, mean=1450, sd=125),digits=0)
if0<-round(rnorm(10, mean=1400, sd=125),digits=0)
inf100<-round(rnorm(10, mean=1400, sd=120),digits=0)
if100<-round(rnorm(10, mean=1750, sd=130),digits=0)
inf200<-round(rnorm(10,mean=2200, sd=125),digits=0)
if200<-round(rnorm(10,mean=1100, sd=120),digits=0)

nf0<-c(inf0, rep(NA,20-length(inf0)))
f0<-c(if0, rep(NA,20-length(if0)))
nf100<-c(inf100, rep(NA,20-length(inf100)))
f100<-c(if100, rep(NA,20-length(if100)))
nf200<-c(inf200, rep(NA,20-length(inf200)))
f200<-c(if200, rep(NA,20-length(if200)))

#cladocerans
icnf0<-round(rnorm(10, mean=80, sd=8),digits=0)
icf0<-round(rnorm(10,mean=40, sd=8),digits=0)
icnf100<-round(rnorm(10, mean=58, sd=8),digits=0)
icf100<-round(rnorm(10, mean=19, sd=8),digits=0)
icnf200<-round(rnorm(10, mean=69, sd=8),digits=0)
icf200<-round(rnorm(10, mean=34, sd=8),digits=0)

cnf0<-c(icnf0, rep(NA,20-length(icnf0)))
cf0<-c(icf0, rep(NA,20-length(icf0)))
cnf100<-c(icnf100, rep(NA,20-length(icnf100)))
cf100<-c(icf100, rep(NA,20-length(icf100)))
cnf200<-c(icnf200, rep(NA,20-length(icnf200)))
cf200<-c(icf200, rep(NA,20-length(icf200)))

twofac<-cbind(mpap,mpaa,maap,maaa,dd,ld,dw,lw,nf0,f0,nf100,f100,nf200,f200,cnf0,cf0,cnf100,cf100,cnf200,cf200)

write.csv(twofac, file = "2facex.csv")