Uncategorized

2-way ANOVA and Tukey HSD

#1 Read in .csv of ammonium data
 require(dplyr)
 setwd('C:/Users/kkuntz/Documents/Personal_Doc/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 for sediment and organism

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 )

Twitter picture

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

Facebook photo

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

Google+ photo

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

Connecting to %s