#Set directory and defining data source
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
t3.aov1 <- aov(t3~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=deerproject1)
summary(t3.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.0267 0.02667 0.159 0.729
## Deer:Honeysuckle 1 0.0002 0.00022 0.001 0.974
## Residuals 2 0.3354 0.16770
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.2784 0.27839 7.014 0.0244 *
## Honeysuckle 1 0.0000 0.00000 0.000 0.9949
## Deer:Honeysuckle 1 0.1459 0.14585 3.675 0.0842 .
## Residuals 10 0.3969 0.03969
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Modeling response variables-Nmin,t1.t2,t3,NetN. Setting others as explanatory and random variables
#Site as random effect
dproject=read.csv("dpdata.csv",header=TRUE)
print(dproject)
## ï..No Site Deer Honeysuckle t1 t2 t3 Nmin NetN
## 1 1 Kramer Con NH 2.99 2.25 2.17 1.39 1.32
## 2 2 Kramer Con H 3.62 2.31 1.75 1.24 1.32
## 3 3 Kramer Excl NH 3.98 1.77 1.74 1.15 1.17
## 4 4 Kramer Excl H 3.02 2.07 1.65 0.85 1.16
## 5 5 Rhinehart Con NH 3.47 2.67 1.82 0.84 0.88
## 6 6 Rhinehart Con H 4.46 2.79 1.83 0.96 0.99
## 7 7 Rhinehart Excl NH 5.18 2.31 0.97 0.67 0.76
## 8 8 Rhinehart Excl H 4.73 2.76 1.67 0.88 0.98
## 9 9 Bishop Con NH 2.66 1.83 2.10 1.30 1.34
## 10 10 Bishop Con H 3.02 2.31 1.74 1.57 1.64
## 11 11 Bishop Excl NH 4.10 1.95 NA 1.06 1.10
## 12 12 Bishop Excl H 3.59 2.07 1.76 1.86 1.63
## 13 13 Western Con NH 3.62 2.34 2.00 1.36 1.40
## 14 14 Western Con H 3.38 2.31 1.99 1.16 1.22
## 15 15 Western Excl NH 4.01 2.85 1.86 1.05 0.82
## 16 16 Western Excl H 3.17 2.40 2.06 1.02 0.70
## 17 17 College Con NH 3.29 2.34 1.92 1.04 1.13
## 18 18 College Con H 4.13 2.88 NA 2.79 2.00
## 19 19 College Excl NH 4.10 3.24 1.68 1.41 1.21
## 20 20 College Excl H 4.34 3.42 1.60 1.45 1.37
Nmin.aov1 <- aov(Nmin~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=dproject)
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=dproject)
summary(t1.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 4 3.516 0.8789
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 1.5568 1.5568 11.340 0.00560 **
## Honeysuckle 1 0.0002 0.0002 0.001 0.97171
## Deer:Honeysuckle 1 1.3005 1.3005 9.473 0.00958 **
## Residuals 12 1.6475 0.1373
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
NetN.aov1 <- aov(NetN~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=dproject)
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=dproject)
summary(t2.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 4 2.379 0.5947
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.0328 0.03281 0.353 0.563
## Honeysuckle 1 0.1566 0.15664 1.685 0.219
## Deer:Honeysuckle 1 0.0162 0.01625 0.175 0.683
## Residuals 12 1.1153 0.09294
t3.aov1 <- aov(t3~ Deer + Honeysuckle + Deer*Honeysuckle + Error(Site),data=dproject)
summary(t3.aov1)
##
## Error: Site
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.0267 0.02667 0.159 0.729
## Deer:Honeysuckle 1 0.0002 0.00023 0.001 0.974
## Residuals 2 0.3354 0.16770
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## Deer 1 0.2784 0.27839 7.014 0.0244 *
## Honeysuckle 1 0.0000 0.00000 0.000 0.9949
## Deer:Honeysuckle 1 0.1459 0.14585 3.675 0.0842 .
## Residuals 10 0.3969 0.03969
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1