Uncategorized

2-way ANOVA and Tukey HSD

#1 Read in .csv of ammonium data
 require(dplyr)
 setwd('C:/thesis')
 NH4 <- read.csv('ammonium.csv')

#2 calculate daily flux for summer using 
averages in light and dark
NH4$dailyflux <- 
((NH4$ave_H_dark*10)+(NH4$ave_H_light*14))

dat1
#3. Run a 2 way anova with interaction

ANOVA <- aov(dailyflux ~ Sed * Org, data = NH4_day2)
summary(ANOVA)
 
 summ

#4 We have significant differences for the organism 
treatment we still want to find out where the 
difference occurs so we run a post hoc test
 
 TukeyHSD(ANOVA, "Org", ordered = TRUE)
 plot(TukeyHSD(ANOVA, "Org"))

posthic

fig
Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s