Untitled

Have categorical x vars and continuous y vars, will run ANOVA on random data sets

First run ANOVA on bare bones random data

ANOVA

Summary Output

aovSUM <- function(dataFrame=data){
     data <- cbind(as.factor(rep(c("Control","Toxic"),each=5)),c(rgamma(5,shape=5,scale=5),rgamma(5,shape=5,scale=10)))
     data
  x <- data[,1]
  y <- data[,2]
  aovMOD <- aov(y~x)
  myOUT <- summary(aovMOD)
  return(myOUT)
}

aovSUM()

##             Df Sum Sq Mean Sq F value Pr(>F)  
## x            1   1420  1419.6   4.106 0.0773 .
## Residuals    8   2766   345.7                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

##With dummy data

randomData <- (cbind(as.factor(rep(c("Maine","New Jersey"),each=5)),c(rgamma(8,shape=5,scale=5),rgamma(2,shape=5,scale=10))))

aovSUM(dataFrame = randomData)

##             Df Sum Sq Mean Sq F value Pr(>F)
## x            1    3.6    3.57   0.028  0.872
## Residuals    8 1023.9  127.99

ANOVA

Plot Output

aovPLT <- function(dataFrame=data){
     data <- cbind(as.factor(rep(c("Control","Toxic"),each=5)),c(rgamma(5,shape=5,scale=5),rgamma(5,shape=5,scale=10)))
  x <- data[,1]
  y <- data[,2]
  aovMOD <- aov(y~x)
  boxplot(y~x,col=c("darkolivegreen1","darkgreen"), ylab="Fecundity")
}

aovPLT()

##With dummy data


randomData <- (cbind(as.factor(rep(c("Maine","New Jersey"),each=5)),c(rgamma(8,shape=5,scale=5),rgamma(2,shape=5,scale=10))))

aovPLT(dataFrame = randomData)

## Error in plot.new(): unable to open file 'C:\Users\Lauren\AppData\Local\Temp\RtmpKMGnTR\file37e043f477ee' for writing

Read in dataframe

Dummy data modeling expected results for cope selection experiment

These are my randomDataBling

These data were modeled looking at variables: Temp HAB and POP, two way interaction

terms and three way interaction terms. Model assumed fecundity would be normally distributed, sd=15, variables explain 40% variation in fecundity, and a relationship factor of 0.5

randomDataBling <- read.csv(file="copeDD", header=TRUE, sep=" ")

ANOVA

Summary Output

##With randomDataBling

randomDataBling <- (cbind(as.factor(rep(c("Maine","New Jersey"),each=5)),c(rgamma(8,shape=5,scale=5),rgamma(2,shape=5,scale=10))))

aovSUM(dataFrame = randomDataBling)

##             Df Sum Sq Mean Sq F value   Pr(>F)    
## x            1 1523.7  1523.7   28.81 0.000671 ***
## Residuals    8  423.1    52.9                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

ANOVA

Plot Output

##With randomDataBling


randomDataBling <- (cbind(as.factor(rep(c("Maine","New Jersey"),each=5)),c(rgamma(8,shape=5,scale=5),rgamma(2,shape=5,scale=10))))

aovPLT(dataFrame = randomDataBling)