Profj <- read.table("C:/Users/User/Desktop/Mahidol-Oxford Bangkok/ProfJosh1.csv",sep = ",", strip.white = T, header = T)
summary(Profj)
## Treatment Weeks Aggr. Rep X13C
## Control :60 Min. :1 <53 :27 Min. :1 Min. :-27.73
## Labelled:60 1st Qu.:2 >2000:30 1st Qu.:1 1st Qu.:-23.04
## Median :3 >250 :30 Median :2 Median :-22.36
## Mean :3 >53 :33 Mean :2 Mean :-21.95
## 3rd Qu.:4 3rd Qu.:3 3rd Qu.:-21.63
## Max. :5 Max. :3 Max. : 13.37
##
## X15N N. C. CN
## Min. : 7.29 Min. :0.0500 Min. : 0.650 Min. :11.20
## 1st Qu.: 8.72 1st Qu.:0.0800 1st Qu.: 0.990 1st Qu.:12.70
## Median : 10.31 Median :0.1300 Median : 1.630 Median :13.00
## Mean : 31.31 Mean :0.1763 Mean : 2.522 Mean :13.76
## 3rd Qu.: 43.97 3rd Qu.:0.2700 3rd Qu.: 3.770 3rd Qu.:13.80
## Max. :255.05 Max. :0.6000 Max. :11.180 Max. :23.50
## NA's :3 NA's :3 NA's :3 NA's :3
AnovaModel.1 <- aov(X15N~Treatment, data=Profj)
summary(AnovaModel.1)
## Df Sum Sq Mean Sq F value Pr(>F)
## Treatment 1 62534 62534 99.22 <2e-16 ***
## Residuals 115 72483 630
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 3 observations deleted due to missingness
AnovaModel.2 <- aov(X13C~Treatment, data=Profj)
summary(AnovaModel.2)
## Df Sum Sq Mean Sq F value Pr(>F)
## Treatment 1 167.2 167.23 14.83 0.000192 ***
## Residuals 118 1330.8 11.28
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1