setwd("C:/Users/Paul/Documents/Rwork")
deerproject1=read.csv("deerExcldata2.csv",header=TRUE)
print(deerproject1)
## ï..sample.. Site Deer Honeysuckle t1 t2 t3 Nmin NetN
## 1 1 Kramer Con NH 1.55 1.16 2.17 1.39 1.32
## 2 2 Kramer Con H 1.86 1.18 1.75 1.24 1.32
## 3 5 Rhinehart Con NH 1.71 1.31 1.82 0.84 0.88
## 4 6 Rhinehart Con H 2.29 1.43 1.83 0.96 0.99
## 5 9 Bishop Con NH 1.31 0.90 2.10 1.30 1.34
## 6 10 Bishop Con H 1.52 1.16 1.74 1.57 1.64
## 7 13 Western Con NH 1.78 1.15 2.00 1.36 1.40
## 8 14 Western Con H 1.73 1.18 1.99 1.16 1.22
## 9 17 College Con NH 1.62 1.14 1.92 1.04 1.13
## 10 18 College Con H 2.08 1.45 NA 2.79 2.00
## 11 3 Kramer Excl NH 2.03 0.90 1.74 1.15 1.17
## 12 4 Kramer Excl H 1.28 0.88 1.65 0.85 1.16
## 13 7 Rhinehart Excl NH 2.58 1.15 0.97 0.67 0.76
## 14 8 Rhinehart Excl H 2.32 1.35 1.67 0.88 0.98
## 15 11 Bishop Excl NH 1.96 0.93 NA 1.06 1.10
## 16 12 Bishop Excl H 1.73 0.99 1.76 1.86 1.63
## 17 15 Western Excl NH 2.02 1.43 1.86 1.05 0.82
## 18 16 Western Excl H 1.56 1.18 2.06 1.02 0.70
## 19 19 College Excl NH 2.10 1.65 1.68 1.41 1.21
## 20 20 College Excl H 2.20 1.72 1.60 1.45 1.37
Nmin.aov1 <- aov(Nmin~ Deer + Honeysuckle +Deer*Honeysuckle + Error(Site),data=deerproject1)
summary(Nmin.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 4 1.627 0.4067
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.2531 0.25312 1.716 0.215
## Honeysuckle 1 0.3150 0.31501 2.135 0.170
## Deer:Honeysuckle 1 0.0572 0.05725 0.388 0.545
## Residuals 12 1.7704 0.14753
t1.aov1 <- aov(t1~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=deerproject1)
summary(t1.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 4 0.9831 0.2458
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.2714 0.2714 6.224 0.0282 *
## Honeysuckle 1 0.0004 0.0004 0.009 0.9248
## Deer:Honeysuckle 1 0.4836 0.4836 11.089 0.0060 **
## Residuals 12 0.5233 0.0436
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
NetN.aov1 <- aov(NetN~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=deerproject1)
summary(NetN.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 4 0.8832 0.2208
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.2738 0.27378 5.505 0.0370 *
## Honeysuckle 1 0.1767 0.17672 3.553 0.0839 .
## Deer:Honeysuckle 1 0.0051 0.00512 0.103 0.7538
## Residuals 12 0.5968 0.04973
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
t2.aov1 <- aov(t2~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=deerproject1)
summary(t2.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 4 0.6705 0.1676
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.0007 0.00072 0.024 0.879
## Honeysuckle 1 0.0320 0.03200 1.070 0.321
## Deer:Honeysuckle 1 0.0231 0.02312 0.773 0.397
## Residuals 12 0.3590 0.02991