Report for Module 04

Using knitr to generate report for module 06 using module 04 data. The goal of this report is to:

getwd()
## [1] "C:/GEOG 5680/module06"
temp = read.csv("Temperature.csv")
head(temp)
##          Sample     Date    DateNr dDay1 dDay2 dDay3 Station Area X31UE_ED50
## 1 DANT.19900110 19900110   10/1/90     7     9     9    DANT   WZ   681379.6
## 2 DANT.19900206 19900206    6/2/90    34    36    36    DANT   WZ   681379.6
## 3 DANT.19900308 19900308    8/3/90    64    66    66    DANT   WZ   681379.6
## 4 DANT.19900404 19900404    4/4/90    91    93    93    DANT   WZ   681379.6
## 5 DANT.19900509 19900509    9/5/90   126   128   128    DANT   WZ   681379.6
## 6 DANT.19900620 19900620 6/20/1990   168   170   170    DANT   WZ   681379.6
##   X31UN_ED50 Year Month Season Salinity Temperature CHLFa
## 1    5920571 1990     1 winter    29.19         4.0   1.3
## 2    5920571 1990     2 winter    27.37         6.0    NA
## 3    5920571 1990     3 spring    24.99         7.3  21.1
## 4    5920571 1990     4 spring    28.79         8.2  25.0
## 5    5920571 1990     5 spring    33.28        17.4  10.2
## 6    5920571 1990     6 summer    32.69        18.1   6.2
mean(temp$Salinity, na.rm = TRUE)
## [1] 29.70221
mean(temp$Temperature, na.rm = TRUE)
## [1] 12.20743
sd(temp$Salinity, na.rm = TRUE)
## [1] 5.421593
sd(temp$Temperature, na.rm = TRUE)
## [1] 5.412521
tapply(temp$Temperature, temp$Station, mean, na.rm = TRUE)
##     DANT     DREI       G6     GROO     HAMM     HANS     HUIB     LODS 
## 12.05908 12.77160 10.66570 12.25853 12.50000 13.54894 11.85353 12.61192 
##     MARS      N02      N10      N20      N70      R03      R50      R70 
## 12.39607 11.37433 12.45878 12.19753 12.14789 12.80000 13.94911 13.98022 
##     SOEL     T004     T010     T100     T135     T175     T235     VLIS 
## 13.21203 11.28425 12.37517 11.94766 11.75512 11.55431 11.32355 12.79292 
##      W02      W20      W70     WISS     ZIJP     ZUID 
## 10.51829 11.87937 12.18243 12.45090 12.55904 11.83928
tapply(temp$Salinity, temp$Station, mean, na.rm = TRUE)
##     DANT     DREI       G6     GROO     HAMM     HANS     HUIB     LODS 
## 28.94468 29.83813 30.59610 14.57571 31.35738 18.29124 29.52631 30.73000 
##     MARS      N02      N10      N20      N70      R03      R50      R70 
## 28.03457 28.78695 30.45223 31.93275 34.94147 30.54978 33.49888 33.67944 
##     SOEL     T004     T010     T100     T135     T175     T235     VLIS 
## 18.09477 32.08114 32.72095 34.56746 34.67667 34.75135 34.87812 29.55476 
##      W02      W20      W70     WISS     ZIJP     ZUID 
## 32.05992 33.31554 35.02006 31.75164 29.39230 28.83997
tapply(temp$Temperature, temp$Year, mean, na.rm = TRUE)
##     1990     1991     1992     1993     1994     1995     1996     1997 
## 11.99146 11.00885 11.55799 11.45385 12.07944 12.32080 10.72728 12.33687 
##     1998     1999     2000     2001     2002     2003     2004     2005 
## 12.13386 13.10577 12.61539 12.74162 13.05500 12.69017 12.56788 12.56112
tapply(temp$Salinity, temp$Year, mean, na.rm = TRUE)
##     1990     1991     1992     1993     1994     1995     1996     1997 
## 30.89873 30.88625 30.27116 29.52839 29.07707 28.94954 30.66670 30.54525 
##     1998     1999     2000     2001     2002     2003     2004     2005 
## 29.55748 29.03291 28.99534 28.22390 28.89323 29.60946 30.45629 30.35355
aggregate(Temperature ~ Year + Month, data = temp, FUN = mean, na.rm = TRUE)
##     Year Month Temperature
## 1   1990     1    6.313889
## 2   1991     1    5.210000
## 3   1992     1    5.858140
## 4   1993     1    4.965217
## 5   1994     1    5.691111
## 6   1995     1    5.778000
## 7   1996     1    3.645238
## 8   1997     1    2.146857
## 9   1998     1    6.685926
## 10  1999     1    5.717778
## 11  2000     1    5.622963
## 12  2001     1    4.754375
## 13  2002     1    4.427778
## 14  2003     1    4.848485
## 15  2004     1    5.558636
## 16  2005     1    5.877692
## 17  1990     2    6.475000
## 18  1991     2    3.576923
## 19  1992     2    5.040625
## 20  1993     2    5.231111
## 21  1994     2    4.814706
## 22  1995     2    6.400000
## 23  1996     2    1.751351
## 24  1997     2    2.793600
## 25  1998     2    5.023438
## 26  1999     2    4.280000
## 27  2000     2    5.635517
## 28  2001     2    5.717826
## 29  2002     2    7.634286
## 30  2003     2    3.864545
## 31  2004     2    5.739286
## 32  2005     2    4.670455
## 33  1990     3    8.192857
## 34  1991     3    5.286207
## 35  1992     3    6.525806
## 36  1993     3    5.521622
## 37  1994     3    5.815686
## 38  1995     3    6.800000
## 39  1996     3    2.950000
## 40  1997     3    6.747551
## 41  1998     3    6.965862
## 42  1999     3    6.684898
## 43  2000     3    6.984773
## 44  2001     3    5.268108
## 45  2002     3    7.487778
## 46  2003     3    6.042927
## 47  2004     3    5.784407
## 48  2005     3    5.051000
## 49  1990     4    9.153333
## 50  1991     4    8.365714
## 51  1992     4    8.030769
## 52  1993     4    8.584615
## 53  1994     4    7.742857
## 54  1995     4    8.386364
## 55  1996     4    6.392683
## 56  1997     4    8.507045
## 57  1998     4    8.965862
## 58  1999     4   10.036364
## 59  2000     4    8.857907
## 60  2001     4    8.117222
## 61  2002     4    9.589038
## 62  2003     4    9.197736
## 63  2004     4    9.695370
## 64  2005     4    9.044444
## 65  1990     5   14.357576
## 66  1991     5   10.185000
## 67  1992     5   12.074324
## 68  1993     5   13.170270
## 69  1994     5   11.987097
## 70  1995     5   11.572549
## 71  1996     5   10.101429
## 72  1997     5   12.223214
## 73  1998     5   13.181455
## 74  1999     5   13.028776
## 75  2000     5   13.255000
## 76  2001     5   12.031867
## 77  2002     5   12.416441
## 78  2003     5   13.620857
## 79  2004     5   12.478364
## 80  2005     5   11.883061
## 81  1990     6   16.277273
## 82  1991     6   13.059259
## 83  1992     6   16.887838
## 84  1993     6   16.024390
## 85  1994     6   14.637255
## 86  1995     6   14.632609
## 87  1996     6   14.996364
## 88  1997     6   15.570200
## 89  1998     6   16.557660
## 90  1999     6   16.019091
## 91  2000     6   15.611111
## 92  2001     6   15.330877
## 93  2002     6   16.343036
## 94  2003     6   18.874643
## 95  2004     6   15.784259
## 96  2005     6   15.061846
## 97  1990     7   17.423333
## 98  1991     7   17.724324
## 99  1992     7   18.745000
## 100 1993     7   17.142857
## 101 1994     7   19.314035
## 102 1995     7   18.335849
## 103 1996     7   16.509091
## 104 1997     7   18.483594
## 105 1998     7   17.108049
## 106 1999     7   18.442623
## 107 2000     7   16.645918
## 108 2001     7   19.071818
## 109 2002     7   17.872241
## 110 2003     7   19.644000
## 111 2004     7   17.292632
## 112 2005     7   18.799111
## 113 1990     8   18.976667
## 114 1991     8   19.550000
## 115 1992     8   18.905000
## 116 1993     8   17.376316
## 117 1994     8   19.903065
## 118 1995     8   20.385000
## 119 1996     8   18.390741
## 120 1997     8   20.919815
## 121 1998     8   18.749200
## 122 1999     8   19.174603
## 123 2000     8   18.756102
## 124 2001     8   19.292273
## 125 2002     8   20.403500
## 126 2003     8   21.038269
## 127 2004     8   19.780185
## 128 2005     8   17.692105
## 129 1990     9   16.831250
## 130 1991     9   17.896000
## 131 1992     9   15.958621
## 132 1993     9   15.660526
## 133 1994     9   15.753846
## 134 1995     9   17.243902
## 135 1996     9   15.953191
## 136 1997     9   17.910833
## 137 1998     9   15.998444
## 138 1999     9   18.999268
## 139 2000     9   17.424048
## 140 2001     9   15.645429
## 141 2002     9   18.297895
## 142 2003     9   17.982400
## 143 2004     9   16.174694
## 144 2005     9   17.883455
## 145 1990    10   13.851724
## 146 1991    10   13.647727
## 147 1992    10   13.095135
## 148 1993    10   13.648485
## 149 1994    10   12.884375
## 150 1995    10   15.133333
## 151 1996    10   13.725000
## 152 1997    10   13.920000
## 153 1998    10   11.752941
## 154 1999    10   13.080385
## 155 2000    10   13.697297
## 156 2001    10   14.866410
## 157 2002    10   13.262703
## 158 2003    10   12.246765
## 159 2004    10   12.790000
## 160 2005    10   14.498125
## 161 1990    11   10.345161
## 162 1991    11    8.400000
## 163 1992    11    8.900833
## 164 1993    11    8.580000
## 165 1994    11   11.212195
## 166 1995    11   11.479487
## 167 1996    11    9.386410
## 168 1997    11   10.031200
## 169 1998    11    9.037097
## 170 1999    11    7.207857
## 171 2000    11    9.189630
## 172 2001    11   10.069130
## 173 2002    11   10.565357
## 174 2003    11    9.401000
## 175 2004    11   10.301786
## 176 2005    11   11.562759
## 177 1990    12    6.766667
## 178 1991    12    6.556757
## 179 1992    12    7.576977
## 180 1993    12    5.532258
## 181 1994    12    9.012821
## 182 1995    12    5.417143
## 183 1996    12    6.581818
## 184 1997    12    6.399259
## 185 1998    12    6.080741
## 186 1999    12    5.637200
## 187 2000    12    7.550952
## 188 2001    12    7.821786
## 189 2002    12    5.246552
## 190 2003    12    6.464762
## 191 2004    12    7.444643
## 192 2005    12    7.111481
aggregate(Salinity ~ Year + Month, data = temp, FUN =mean, na.rm = TRUE)
##     Year Month Salinity
## 1   1990     1 30.27000
## 2   1991     1 29.77267
## 3   1992     1 29.99535
## 4   1993     1 28.19087
## 5   1994     1 27.31000
## 6   1995     1 27.97180
## 7   1996     1 30.66524
## 8   1997     1 31.24865
## 9   1998     1 30.24886
## 10  1999     1 28.94469
## 11  2000     1 28.49000
## 12  2001     1 29.25000
## 13  2002     1 28.19519
## 14  2003     1 28.65606
## 15  2004     1 28.64591
## 16  2005     1 30.72360
## 17  1990     2 27.94625
## 18  1991     2 30.51923
## 19  1992     2 30.47500
## 20  1993     2 27.25280
## 21  1994     2 28.97088
## 22  1995     2 25.15593
## 23  1996     2 30.01189
## 24  1997     2 28.73720
## 25  1998     2 30.55347
## 26  1999     2 26.74842
## 27  2000     2 28.62345
## 28  2001     2 26.99696
## 29  2002     2 26.19952
## 30  2003     2 26.75606
## 31  2004     2 25.41143
## 32  2005     2 29.20818
## 33  1990     3 30.22250
## 34  1991     3 30.27966
## 35  1992     3 30.16500
## 36  1993     3 28.32056
## 37  1994     3 27.95980
## 38  1995     3 28.45118
## 39  1996     3 30.35333
## 40  1997     3 30.33224
## 41  1998     3 26.51448
## 42  1999     3 27.82719
## 43  2000     3 26.74886
## 44  2001     3 27.01111
## 45  2002     3 27.08886
## 46  2003     3 27.52780
## 47  2004     3 29.26576
## 48  2005     3 29.38240
## 49  1990     4 30.80133
## 50  1991     4 30.84543
## 51  1992     4 29.69103
## 52  1993     4 29.49949
## 53  1994     4 27.25381
## 54  1995     4 28.28159
## 55  1996     4 29.99561
## 56  1997     4 29.89523
## 57  1998     4 30.20926
## 58  1999     4 26.46606
## 59  2000     4 27.78628
## 60  2001     4 25.79417
## 61  2002     4 29.27327
## 62  2003     4 29.42528
## 63  2004     4 30.29222
## 64  2005     4 29.66933
## 65  1990     5 30.71281
## 66  1991     5 29.63042
## 67  1992     5 29.55811
## 68  1993     5 30.34541
## 69  1994     5 29.02492
## 70  1995     5 29.05375
## 71  1996     5 30.96757
## 72  1997     5 29.87750
## 73  1998     5 30.07436
## 74  1999     5 28.90551
## 75  2000     5 29.11696
## 76  2001     5 27.69520
## 77  2002     5 28.68169
## 78  2003     5 27.61743
## 79  2004     5 30.51545
## 80  2005     5 29.39306
## 81  1990     6 31.01000
## 82  1991     6 31.08407
## 83  1992     6 30.35611
## 84  1993     6 30.54561
## 85  1994     6 28.83816
## 86  1995     6 28.32043
## 87  1996     6 30.93800
## 88  1997     6 30.86431
## 89  1998     6 29.16064
## 90  1999     6 29.79455
## 91  2000     6 29.80463
## 92  2001     6 28.86246
## 93  2002     6 29.54439
## 94  2003     6 27.59964
## 95  2004     6 30.69167
## 96  2005     6 30.86169
## 97  1990     7 31.43033
## 98  1991     7 30.83892
## 99  1992     7 30.59800
## 100 1993     7 30.94595
## 101 1994     7 30.38172
## 102 1995     7 29.41472
## 103 1996     7 31.25352
## 104 1997     7 30.05484
## 105 1998     7 29.86098
## 106 1999     7 30.31197
## 107 2000     7 29.52940
## 108 2001     7 29.14891
## 109 2002     7 30.07845
## 110 2003     7 31.02720
## 111 2004     7 31.00930
## 112 2005     7 30.15822
## 113 1990     8 31.99467
## 114 1991     8 31.43059
## 115 1992     8 30.80000
## 116 1993     8 30.53184
## 117 1994     8 30.16051
## 118 1995     8 29.85034
## 119 1996     8 31.21519
## 120 1997     8 30.58132
## 121 1998     8 30.79800
## 122 1999     8 30.69317
## 123 2000     8 29.97169
## 124 2001     8 30.04818
## 125 2002     8 29.98167
## 126 2003     8 31.83269
## 127 2004     8 31.41574
## 128 2005     8 30.64439
## 129 1990     9 30.41125
## 130 1991     9 31.36160
## 131 1992     9 31.25241
## 132 1993     9 28.51892
## 133 1994     9 30.61487
## 134 1995     9 29.36300
## 135 1996     9 30.43723
## 136 1997     9 30.03083
## 137 1998     9 30.19549
## 138 1999     9 30.03122
## 139 2000     9 29.35476
## 140 2001     9 27.70486
## 141 2002     9 29.12079
## 142 2003     9 31.38060
## 143 2004     9 31.03040
## 144 2005     9 31.44145
## 145 1990    10 31.91690
## 146 1991    10 31.58044
## 147 1992    10 31.80270
## 148 1993    10 30.09182
## 149 1994    10 30.50844
## 150 1995    10 30.29039
## 151 1996    10 31.50750
## 152 1997    10 31.62538
## 153 1998    10 26.48765
## 154 1999    10 29.21885
## 155 2000    10 29.77703
## 156 2001    10 28.18974
## 157 2002    10 27.88711
## 158 2003    10 30.37000
## 159 2004    10 31.61171
## 160 2005    10 30.68656
## 161 1990    11 31.53484
## 162 1991    11 31.11848
## 163 1992    11 28.24167
## 164 1993    11 29.79150
## 165 1994    11 28.81049
## 166 1995    11 30.25692
## 167 1996    11 30.36205
## 168 1997    11 31.85209
## 169 1998    11 27.67930
## 170 1999    11 28.35071
## 171 2000    11 29.72741
## 172 2001    11 27.49217
## 173 2002    11 28.56679
## 174 2003    11 31.03300
## 175 2004    11 31.16250
## 176 2005    11 31.25586
## 177 1990    12 31.10160
## 178 1991    12 31.42028
## 179 1992    12 29.92356
## 180 1993    12 28.58935
## 181 1994    12 28.97487
## 182 1995    12 29.41657
## 183 1996    12 29.53676
## 184 1997    12 31.20455
## 185 1998    12 29.48020
## 186 1999    12 28.57640
## 187 2000    12 27.88381
## 188 2001    12 28.03893
## 189 2002    12 29.41345
## 190 2003    12 29.75091
## 191 2004    12 30.50429
## 192 2005    12 30.53444
#using table() to determine # of observations
#at each station
table(temp$Station)
## 
## DANT DREI   G6 GROO HAMM HANS HUIB LODS MARS  N02  N10  N20  N70  R03  R50  R70 
##  300  293  278  296  295  309  296  294  296  402  665  266  268  161  106  106 
## SOEL T004 T010 T100 T135 T175 T235 VLIS  W02  W20  W70 WISS ZIJP ZUID 
##  295  339  261  258  259  258  258  421  272  191  190  296  296  303
#each year
table(temp$Year)
## 
## 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 
##  367  392  438  436  590  590  583  636  608  570  563  568  545  550  540  552
#each station, each year
table(temp$Station, temp$Year)
##       
##        1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003
##   DANT   12   12   18   19   17   13   19   20   22   21   22   21   21   21
##   DREI   13   13   13   13   19   21   20   20   21   20   20   20   20   20
##   G6     12   16   20   17   21   24   23   28   23   21   13   12   12   12
##   GROO   12   12   18   19   17   12   18   19   21   21   21   21   21   21
##   HAMM   13   13   13   13   20   21   20   21   20   20   21   20   20   20
##   HANS   13   16   13   13   21   24   19   19   18   18   29   29   20   19
##   HUIB   12   11   18   19   17   12   19   19   21   21   21   22   21   21
##   LODS   13   13   13   13   20   21   20   21   20   20   20   20   20   20
##   MARS   12   12   18   18   16   12   19   19   22   21   21   21   21   21
##   N02    12   17   20   17   27   29   31   40   37   29   24   27   23   23
##   N10    47   45   40   46   53   56   43   52   45   43   36   37   30   31
##   N20    12   12   11   11   18   18   18   18   18   18   19   19   19   19
##   N70    12   11   11   11   18   18   18   18   18   19   18   20   20   19
##   R03     0    5    7    7   13   17   17   23   20    9    8    7    7    7
##   R50     0    5    7    7   10    7    7    7    7    7    7    7    7    7
##   R70     0    5    7    7   10    7    7    7    7    7    7    7    7    7
##   SOEL   13   13   13   12   20   22   20   21   20   20   20   20   20   20
##   T004   12   14   20   16   26   27   27   33   28   20   19   19   19   20
##   T010   12   10   11   10   18   18   18   18   18   18   18   18   18   19
##   T100   12   10   10   10   18   18   18   18   18   18   18   18   18   18
##   T135   12   10   10   10   18   18   18   18   18   18   17   19   18   19
##   T175   12   10   10   10   18   18   18   18   18   18   18   18   18   18
##   T235   12   10   10   10   18   18   18   18   18   18   18   18   18   18
##   VLIS   14   18   21   24   35   36   22   23   22   22   30   31   30   32
##   W02    12   17   20   17   21   23   23   28   23   16   12   12   12   12
##   W20    12   12   11   11   12   12   12   12   12   13   12   12   12   12
##   W70    12   12   11   11   12   12   12   12   12   12   12   12   12   12
##   WISS   13   13   13   13   20   21   20   22   20   20   21   20   20   20
##   ZIJP   13   13   13   13   20   21   20   23   20   20   20   20   20   20
##   ZUID   11   12   18   19   17   14   19   21   21   22   21   21   21   22
##       
##        2004 2005
##   DANT   21   21
##   DREI   20   20
##   G6     12   12
##   GROO   21   22
##   HAMM   20   20
##   HANS   19   19
##   HUIB   21   21
##   LODS   20   20
##   MARS   21   22
##   N02    22   24
##   N10    29   32
##   N20    17   19
##   N70    18   19
##   R03     7    7
##   R50     7    7
##   R70     7    7
##   SOEL   20   21
##   T004   19   20
##   T010   18   19
##   T100   18   18
##   T135   18   18
##   T175   18   18
##   T235   18   18
##   VLIS   31   30
##   W02    12   12
##   W20    12   12
##   W70    12   12
##   WISS   20   20
##   ZIJP   20   20
##   ZUID   22   22