library(readxl)
## Warning: package 'readxl' was built under R version 4.4.2
bird_Ohio <- read_excel("C:/Users/LENOVO/Downloads/bird_Ohio.xlsx")
## New names:
## • `` -> `...1`
env_Ohio <- read_excel("env_Ohio.xlsx")
## New names:
## • `` -> `...1`
library(vegan)
## Warning: package 'vegan' was built under R version 4.4.2
## Loading required package: permute
## Warning: package 'permute' was built under R version 4.4.2
## Loading required package: lattice
## This is vegan 2.6-8
bird_Ohio
## # A tibble: 210 × 49
##     ...1  ACFL  AMCR  AMGO  AMRE  AMRO  BAOR  BAWW  BGGN  BHCO  BLBW  BLJA  BRTH
##    <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
##  1     1     0     0     1     0     0     0     1     0     0     0     1     0
##  2     2     9     0     0     0     3     0     2     2     1     0     1     0
##  3     3     4     0     0     0     3     0     0     4     4     0     0     0
##  4     4     7     1     0     1     0     0     3     3     0     0     0     0
##  5     5     2     0     0     0     2     0     1     0     1     0     0     0
##  6     6     8     0     0     5     1     0     4     0     3     0     0     0
##  7     7     3     3     2     2     0     0     0     3     1     0     0     1
##  8     8    11     0     0     4     0     0     2     2     0     0     0     0
##  9     9    10     0     0     1     1     0     3     1     3     0     0     0
## 10    10     1     0     0     0     0     0     2     1     0     0     1     1
## # ℹ 200 more rows
## # ℹ 36 more variables: BTNW <dbl>, BWWA <dbl>, CACH <dbl>, CARW <dbl>,
## #   CEDW <dbl>, CERW <dbl>, CHSP <dbl>, COYE <dbl>, EAPH <dbl>, EATO <dbl>,
## #   EAWP <dbl>, ETTI <dbl>, FISP <dbl>, GCFL <dbl>, GRCA <dbl>, HOWA <dbl>,
## #   INBU <dbl>, KEWA <dbl>, LOWA <dbl>, NOCA <dbl>, NOPA <dbl>, OVEN <dbl>,
## #   PIWA <dbl>, PRAW <dbl>, RBGR <dbl>, REVI <dbl>, RWBL <dbl>, SCTA <dbl>,
## #   SUTA <dbl>, WBNU <dbl>, WEVI <dbl>, WEWA <dbl>, WOTH <dbl>, YBCH <dbl>, …
ncol(bird_Ohio)
## [1] 49
sp.rich<-specnumber(bird_Ohio, MARGIN=1)
as.data.frame(sp.rich)
##     sp.rich
## 1        11
## 2        16
## 3        15
## 4        17
## 5        16
## 6        17
## 7        23
## 8        19
## 9        18
## 10       15
## 11       18
## 12       19
## 13       18
## 14       18
## 15       17
## 16       20
## 17       13
## 18       13
## 19       16
## 20       17
## 21       22
## 22       14
## 23       21
## 24       19
## 25       14
## 26       15
## 27       16
## 28       16
## 29       18
## 30       13
## 31       19
## 32       21
## 33       18
## 34       17
## 35       15
## 36       22
## 37       18
## 38       22
## 39       11
## 40       16
## 41       29
## 42       23
## 43       13
## 44       15
## 45       18
## 46       15
## 47       14
## 48       10
## 49       21
## 50       17
## 51       12
## 52       23
## 53       17
## 54       17
## 55       25
## 56       17
## 57       13
## 58       14
## 59       11
## 60       21
## 61       15
## 62       25
## 63       15
## 64       12
## 65       15
## 66       20
## 67       18
## 68       18
## 69       19
## 70       12
## 71       12
## 72       14
## 73       17
## 74       16
## 75       18
## 76       18
## 77       19
## 78       14
## 79       16
## 80       14
## 81       13
## 82       17
## 83        7
## 84       16
## 85       17
## 86       23
## 87       17
## 88       15
## 89       15
## 90       15
## 91       15
## 92       18
## 93       12
## 94       10
## 95       16
## 96       20
## 97       23
## 98       19
## 99       18
## 100      23
## 101      22
## 102      25
## 103      22
## 104      14
## 105      22
## 106      15
## 107      16
## 108      18
## 109      19
## 110      17
## 111      13
## 112      16
## 113      16
## 114      16
## 115      15
## 116      15
## 117      23
## 118      24
## 119      21
## 120      13
## 121      16
## 122      18
## 123      20
## 124      17
## 125      19
## 126      21
## 127      15
## 128      18
## 129      14
## 130      19
## 131      15
## 132      14
## 133      18
## 134      16
## 135      20
## 136      27
## 137      23
## 138      18
## 139      20
## 140      19
## 141      26
## 142      15
## 143      10
## 144      21
## 145      26
## 146      17
## 147      14
## 148      13
## 149      12
## 150      18
## 151      29
## 152      24
## 153      15
## 154      24
## 155      28
## 156      15
## 157      18
## 158      16
## 159      17
## 160      19
## 161      16
## 162      21
## 163      16
## 164      29
## 165      20
## 166      21
## 167      22
## 168      23
## 169      25
## 170      15
## 171      12
## 172      19
## 173      18
## 174      15
## 175      15
## 176      18
## 177      15
## 178      14
## 179      18
## 180      13
## 181      14
## 182      21
## 183      14
## 184      17
## 185      18
## 186      22
## 187      14
## 188      18
## 189      23
## 190      24
## 191      13
## 192      20
## 193      18
## 194      21
## 195      13
## 196      23
## 197      19
## 198      15
## 199      19
## 200      18
## 201      18
## 202      16
## 203      19
## 204      10
## 205      16
## 206      16
## 207      11
## 208      24
## 209      23
## 210      20
sp.even<-specnumber(bird_Ohio,MARGIN=2)
as.data.frame(sp.even)
##      sp.even
## ...1     210
## ACFL     181
## AMCR      43
## AMGO      34
## AMRE      81
## AMRO      78
## BAOR      16
## BAWW     149
## BGGN     136
## BHCO     133
## BLBW      12
## BLJA     102
## BRTH       9
## BTNW      18
## BWWA      28
## CACH      54
## CARW      32
## CEDW      23
## CERW      61
## CHSP       5
## COYE      25
## EAPH      32
## EATO     141
## EAWP     127
## ETTI     147
## FISP       6
## GCFL      37
## GRCA      24
## HOWA     187
## INBU      52
## KEWA      62
## LOWA      35
## NOCA      70
## NOPA      12
## OVEN     196
## PIWA      10
## PRAW      26
## RBGR      24
## REVI     208
## RWBL       6
## SCTA     172
## SUTA      20
## WBNU     157
## WEVI      48
## WEWA     157
## WOTH     173
## YBCH      24
## YTVI      95
## YTWA      14
shannon<-diversity(bird_Ohio, index = "shannon")
shannon
##   [1] 1.8917399 2.5122269 2.5759511 2.5306341 2.5813695 2.5851922 2.8543118
##   [8] 2.6869090 2.6290101 2.4471082 2.4514742 2.6624436 2.4801533 2.4524231
##  [15] 2.4446459 2.5985159 2.1580650 2.0290131 2.3439848 2.3677896 2.6230152
##  [22] 2.1818266 2.3999167 2.3497389 2.1019228 2.1558323 2.1542864 2.2277483
##  [29] 2.3289428 1.8809617 2.4493401 2.3266446 2.2270348 1.9964440 2.0417370
##  [36] 2.4618897 2.1320689 2.2922720 1.5430953 2.0627115 2.6546098 2.3331170
##  [43] 1.7535317 1.8944254 1.9905828 1.8202840 1.7588983 1.2049611 2.1437307
##  [50] 2.0737966 1.6534827 2.1077805 1.9625262 1.7738141 2.0210011 1.8759969
##  [57] 1.4522335 1.5969316 1.4228140 2.2233828 1.6941131 2.2527933 1.4705667
##  [64] 1.4376671 1.6556983 2.0069460 1.7242396 1.6262868 1.9412633 1.2809809
##  [71] 1.3010760 1.5318794 1.5531853 1.5204937 1.7648837 1.6943305 1.7000499
##  [78] 1.3091997 1.5702975 1.3837769 1.3385135 1.4513903 0.9429183 1.4490505
##  [85] 1.7895111 1.9656594 1.5788314 1.5022623 1.4904512 1.4962931 1.4059042
##  [92] 1.6478061 1.2990656 1.0930459 1.5543977 1.7289643 1.6443609 1.5495614
##  [99] 1.4031313 1.7226377 1.8883235 1.9303606 1.4901640 1.5787674 1.5705318
## [106] 1.4178638 1.2837894 1.5926701 1.6081552 1.4082954 1.2234179 1.2528753
## [113] 1.3100790 1.4507992 1.2323738 1.3797947 1.7293521 1.4854109 1.4658616
## [120] 0.9669648 1.1795715 1.3672704 1.5350481 1.3090028 1.2126036 1.3047201
## [127] 1.2040135 1.3639687 1.1349332 1.4067959 1.2394246 1.3198667 1.1560033
## [134] 1.2668167 1.3633187 1.7184293 1.4705114 1.4134997 1.1657677 1.4382794
## [141] 1.6241226 0.9479850 0.7890530 1.4545973 1.4366209 1.3270892 1.1408832
## [148] 0.9003295 0.9349969 1.1903317 1.5749987 1.5003813 1.0810697 1.5488876
## [155] 1.5759803 0.9573217 1.2080605 1.1433944 1.0993195 1.1174376 1.1019272
## [162] 0.9943400 1.0646499 1.6364990 1.2202461 1.3055935 1.2551904 1.3908337
## [169] 1.5319849 0.9827985 0.9091475 1.3284004 1.0280144 1.0155201 0.8256294
## [176] 1.1407526 1.1095773 0.9370976 1.1957853 0.9132408 0.9268789 0.9070407
## [183] 0.9184634 1.0505849 1.1130730 1.1475092 0.8861070 1.1679846 1.3269686
## [190] 1.4025824 0.7316835 1.0732068 1.0880649 0.8289934 0.8676827 1.1044350
## [197] 1.0755235 1.0144104 1.1081660 1.0762302 1.0511711 0.9282779 1.2321592
## [204] 0.7472004 0.8967415 0.9696155 0.4879819 1.2601278 1.2607309 1.0734488
simpson<-diversity(bird_Ohio, index = "simpson")
simpson
##   [1] 0.7923875 0.9023669 0.9141051 0.8955078 0.9112426 0.9100346 0.9286265
##   [8] 0.9171429 0.9131944 0.8927116 0.8897290 0.9137893 0.8793388 0.8758573
##  [15] 0.8830959 0.8967347 0.8536155 0.8264046 0.8624852 0.8617998 0.8762865
##  [22] 0.8433163 0.8526139 0.8573088 0.8130987 0.8250548 0.8136574 0.8264701
##  [29] 0.8491358 0.7441406 0.8553590 0.8238062 0.8106576 0.7780612 0.7801904
##  [36] 0.8531268 0.7831946 0.8068698 0.6477631 0.7808963 0.8616564 0.7890625
##  [43] 0.6954194 0.7339693 0.7269136 0.6911844 0.6822566 0.5081633 0.7527571
##  [50] 0.7694515 0.6567901 0.7160000 0.7361333 0.6479362 0.6747189 0.6863983
##  [57] 0.5658574 0.6459054 0.5912465 0.7710506 0.6341785 0.7706531 0.5559896
##  [64] 0.5789628 0.6272199 0.7150879 0.6403025 0.6024793 0.7088021 0.4988000
##  [71] 0.4968140 0.5861804 0.5995339 0.5778835 0.6374853 0.6346974 0.6039282
##  [78] 0.4959013 0.5915848 0.5347080 0.5234094 0.5429308 0.4238292 0.5318678
##  [85] 0.6617778 0.6964237 0.5921019 0.5775047 0.5791837 0.5853587 0.5223585
##  [92] 0.6010774 0.5181661 0.4416889 0.5807846 0.6220482 0.5677222 0.5565515
##  [99] 0.5220452 0.5952296 0.6635526 0.6503184 0.5207111 0.6158734 0.5447686
## [106] 0.5379747 0.4695946 0.5760346 0.5821278 0.5175781 0.5051503 0.4630629
## [113] 0.4962000 0.5478637 0.4671492 0.5235294 0.6090305 0.5106574 0.5185606
## [120] 0.3701524 0.4405887 0.5040980 0.5602648 0.4785325 0.4391298 0.4572742
## [127] 0.4574581 0.4963772 0.4396221 0.5037577 0.4670187 0.4978497 0.4104902
## [134] 0.4898061 0.4876454 0.5769861 0.5051541 0.5225468 0.3938079 0.5078587
## [141] 0.5466634 0.3572150 0.3180077 0.4964437 0.4751863 0.4889245 0.4303415
## [148] 0.3341988 0.3763435 0.4262014 0.5111597 0.5117908 0.4075500 0.5420794
## [155] 0.5161432 0.3428280 0.4343764 0.4120049 0.4084227 0.3921670 0.4009981
## [162] 0.3344865 0.4009879 0.5320156 0.4317769 0.4590037 0.4250527 0.4708680
## [169] 0.5338560 0.3637872 0.3495779 0.4747466 0.3657520 0.3795146 0.2896172
## [176] 0.4219639 0.4203375 0.3414050 0.4325555 0.3431800 0.3372680 0.3009006
## [183] 0.3442333 0.3724143 0.3960015 0.3812779 0.3234457 0.4113292 0.4541780
## [190] 0.4803133 0.2576901 0.3719623 0.3935743 0.2681985 0.3192152 0.3765746
## [197] 0.3795178 0.3771101 0.3910828 0.3889234 0.3824131 0.3371236 0.4337904
## [204] 0.2852947 0.3188013 0.3529198 0.1674397 0.4239074 0.4304125 0.3842678
inv.simpson<-diversity(bird_Ohio, index = "invsimpson")
inv.simpson
##   [1]  4.816667 10.242424 11.642140  9.570093 11.266667 11.115385 14.010811
##   [8] 12.068966 11.520000  9.320675  9.068571 11.599483  8.287671  8.055249
##  [15]  8.554023  9.683794  6.831325  5.760522  7.271945  7.235880  8.083192
##  [22]  6.382284  6.784902  7.008138  5.350417  5.716075  5.366460  5.762696
##  [29]  6.628478  3.908397  6.913669  5.675570  5.281437  4.505747  4.549391
##  [36]  6.808596  4.612431  5.177854  2.838998  4.564050  7.228380  4.740741
##  [43]  3.283203  3.758965  3.661844  3.238179  3.147193  2.033195  4.044605
##  [50]  4.337483  2.913669  3.521127  3.789793  2.840394  3.074264  3.188759
##  [57]  2.303391  2.824104  2.446462  4.367778  2.733574  4.360207  2.252199
##  [64]  2.375087  2.682547  3.509854  2.780114  2.515593  3.434090  1.995211
##  [71]  1.987337  2.416512  2.497090  2.369014  2.758509  2.737457  2.524795
##  [78]  1.983739  2.448489  2.149188  2.098237  2.187852  1.735597  2.136149
##  [85]  2.956636  3.294065  2.451592  2.366890  2.376334  2.411723  2.093621
##  [92]  2.506752  2.075404  1.791116  2.385409  2.645840  2.313327  2.255053
##  [99]  2.092248  2.470536  2.972233  2.859744  2.086424  2.603309  2.196685
## [106]  2.164384  1.885350  2.358683  2.393076  2.072874  2.020815  1.862416
## [113]  1.984915  2.211722  1.876698  2.098765  2.557744  2.043558  2.077105
## [120]  1.587686  1.787593  2.016527  2.274096  1.917665  1.782944  1.842551
## [127]  1.843175  1.985613  1.784510  2.015144  1.876238  1.991436  1.696325
## [134]  1.960039  1.951773  2.363989  2.020831  2.094446  1.649642  2.031937
## [141]  2.205866  1.555730  1.466292  1.985875  1.905438  1.956658  1.755438
## [148]  1.501950  1.603447  1.742772  2.045658  2.048302  1.687906  2.183785
## [155]  2.066727  1.521672  1.767960  1.700695  1.690396  1.645189  1.669444
## [162]  1.502599  1.669415  2.136824  1.759872  1.848441  1.739290  1.889888
## [169]  2.145260  1.571801  1.537463  1.903843  1.576670  1.611642  1.407692
## [176]  1.729996  1.725142  1.518384  1.762287  1.522487  1.508906  1.430412
## [183]  1.524933  1.593408  1.655633  1.616235  1.478078  1.698742  1.832099
## [190]  1.924236  1.347146  1.592261  1.649007  1.366491  1.468893  1.604041
## [197]  1.611650  1.605420  1.642259  1.636456  1.619205  1.508577  1.766130
## [204]  1.399178  1.468000  1.545403  1.201114  1.735832  1.755657  1.624083
fish.alp<-fisher.alpha(bird_Ohio)
fish.alp
##   [1]  5.642066  7.896600  6.488217  7.565648  7.896600  7.275326 11.680602
##   [8]  8.578491  7.703286  7.612305  7.916266  8.839863  9.316717  9.454690
##  [15]  7.814550  9.354043  4.970140  4.970140  6.781181  7.415136 11.156386
##  [22]  5.110244  9.198544  7.705935  5.389613  5.679263  6.377738  6.055108
##  [29]  6.765902  4.927590  7.224522  8.920683  7.027652  6.430958  5.381526
##  [36]  8.310228  6.891754  8.984035  3.690822  5.538659 11.970603  9.586698
##  [43]  4.350582  5.112918  6.765902  5.254828  4.718079  3.192559  7.976087
##  [50]  5.576450  3.718729  9.354332  5.646225  6.096825 10.837223  5.798205
##  [57]  4.212510  4.392532  3.217794  7.064403  4.829737  8.991549  4.986195
##  [64]  3.546422  4.711590  6.648249  5.988234  6.114229  6.065899  3.560586
##  [71]  3.546422  4.192993  5.427384  5.030672  5.726751  5.687885  6.259378
##  [78]  4.238987  4.858106  4.121275  3.717707  5.349775  1.659647  4.890459
##  [85]  4.931276  7.342929  5.082987  4.281400  4.257961  4.223961  4.342856
##  [92]  5.369842  3.173985  2.544462  4.545684  6.048954  7.602821  5.778996
##  [99]  5.414712  7.383522  6.606168  7.981444  7.105981  3.601971  6.964669
## [106]  4.071694  4.557373  5.109555  5.469168  4.810390  3.335706  4.489371
## [113]  4.415941  4.311580  4.062475  3.967505  6.848485  7.622560  6.259936
## [120]  3.401558  4.396019  5.031600  5.661968  4.672620  5.506903  6.274631
## [127]  3.935634  4.959159  3.588003  5.297514  3.882807  3.497519  5.042383
## [134]  4.164308  5.639364  8.224130  6.763796  4.777415  5.769443  5.149475
## [141]  7.797693  3.904973  2.305868  5.871991  7.955229  4.394232  3.450861
## [148]  3.204387  2.846412  4.802582  8.995673  6.871785  3.756033  6.745033
## [155]  8.446768  3.800210  4.713438  4.046218  4.379899  5.104476  4.033430
## [162]  5.926942  4.014640  8.614307  5.345689  5.647828  6.075704  6.342263
## [169]  6.915739  3.675202  2.750911  4.849896  4.653569  3.632950  3.708700
## [176]  4.551192  3.574193  3.328464  4.513385  3.010270  3.314740  5.730437
## [183]  3.296916  4.218990  4.513385  5.923958  3.292541  4.471297  6.113701
## [190]  6.406272  3.014282  5.165024  4.453883  5.647828  2.957055  6.182347
## [197]  4.772041  3.494964  4.741473  4.409155  4.409155  3.815290  4.655622
## [204]  2.103496  3.815290  3.778652  2.415700  6.323830  5.950260  4.993876
Div.Ind<-cbind.data.frame(shannon, simpson, inv.simpson,fish.alp)
Div.Ind
##       shannon   simpson inv.simpson  fish.alp
## 1   1.8917399 0.7923875    4.816667  5.642066
## 2   2.5122269 0.9023669   10.242424  7.896600
## 3   2.5759511 0.9141051   11.642140  6.488217
## 4   2.5306341 0.8955078    9.570093  7.565648
## 5   2.5813695 0.9112426   11.266667  7.896600
## 6   2.5851922 0.9100346   11.115385  7.275326
## 7   2.8543118 0.9286265   14.010811 11.680602
## 8   2.6869090 0.9171429   12.068966  8.578491
## 9   2.6290101 0.9131944   11.520000  7.703286
## 10  2.4471082 0.8927116    9.320675  7.612305
## 11  2.4514742 0.8897290    9.068571  7.916266
## 12  2.6624436 0.9137893   11.599483  8.839863
## 13  2.4801533 0.8793388    8.287671  9.316717
## 14  2.4524231 0.8758573    8.055249  9.454690
## 15  2.4446459 0.8830959    8.554023  7.814550
## 16  2.5985159 0.8967347    9.683794  9.354043
## 17  2.1580650 0.8536155    6.831325  4.970140
## 18  2.0290131 0.8264046    5.760522  4.970140
## 19  2.3439848 0.8624852    7.271945  6.781181
## 20  2.3677896 0.8617998    7.235880  7.415136
## 21  2.6230152 0.8762865    8.083192 11.156386
## 22  2.1818266 0.8433163    6.382284  5.110244
## 23  2.3999167 0.8526139    6.784902  9.198544
## 24  2.3497389 0.8573088    7.008138  7.705935
## 25  2.1019228 0.8130987    5.350417  5.389613
## 26  2.1558323 0.8250548    5.716075  5.679263
## 27  2.1542864 0.8136574    5.366460  6.377738
## 28  2.2277483 0.8264701    5.762696  6.055108
## 29  2.3289428 0.8491358    6.628478  6.765902
## 30  1.8809617 0.7441406    3.908397  4.927590
## 31  2.4493401 0.8553590    6.913669  7.224522
## 32  2.3266446 0.8238062    5.675570  8.920683
## 33  2.2270348 0.8106576    5.281437  7.027652
## 34  1.9964440 0.7780612    4.505747  6.430958
## 35  2.0417370 0.7801904    4.549391  5.381526
## 36  2.4618897 0.8531268    6.808596  8.310228
## 37  2.1320689 0.7831946    4.612431  6.891754
## 38  2.2922720 0.8068698    5.177854  8.984035
## 39  1.5430953 0.6477631    2.838998  3.690822
## 40  2.0627115 0.7808963    4.564050  5.538659
## 41  2.6546098 0.8616564    7.228380 11.970603
## 42  2.3331170 0.7890625    4.740741  9.586698
## 43  1.7535317 0.6954194    3.283203  4.350582
## 44  1.8944254 0.7339693    3.758965  5.112918
## 45  1.9905828 0.7269136    3.661844  6.765902
## 46  1.8202840 0.6911844    3.238179  5.254828
## 47  1.7588983 0.6822566    3.147193  4.718079
## 48  1.2049611 0.5081633    2.033195  3.192559
## 49  2.1437307 0.7527571    4.044605  7.976087
## 50  2.0737966 0.7694515    4.337483  5.576450
## 51  1.6534827 0.6567901    2.913669  3.718729
## 52  2.1077805 0.7160000    3.521127  9.354332
## 53  1.9625262 0.7361333    3.789793  5.646225
## 54  1.7738141 0.6479362    2.840394  6.096825
## 55  2.0210011 0.6747189    3.074264 10.837223
## 56  1.8759969 0.6863983    3.188759  5.798205
## 57  1.4522335 0.5658574    2.303391  4.212510
## 58  1.5969316 0.6459054    2.824104  4.392532
## 59  1.4228140 0.5912465    2.446462  3.217794
## 60  2.2233828 0.7710506    4.367778  7.064403
## 61  1.6941131 0.6341785    2.733574  4.829737
## 62  2.2527933 0.7706531    4.360207  8.991549
## 63  1.4705667 0.5559896    2.252199  4.986195
## 64  1.4376671 0.5789628    2.375087  3.546422
## 65  1.6556983 0.6272199    2.682547  4.711590
## 66  2.0069460 0.7150879    3.509854  6.648249
## 67  1.7242396 0.6403025    2.780114  5.988234
## 68  1.6262868 0.6024793    2.515593  6.114229
## 69  1.9412633 0.7088021    3.434090  6.065899
## 70  1.2809809 0.4988000    1.995211  3.560586
## 71  1.3010760 0.4968140    1.987337  3.546422
## 72  1.5318794 0.5861804    2.416512  4.192993
## 73  1.5531853 0.5995339    2.497090  5.427384
## 74  1.5204937 0.5778835    2.369014  5.030672
## 75  1.7648837 0.6374853    2.758509  5.726751
## 76  1.6943305 0.6346974    2.737457  5.687885
## 77  1.7000499 0.6039282    2.524795  6.259378
## 78  1.3091997 0.4959013    1.983739  4.238987
## 79  1.5702975 0.5915848    2.448489  4.858106
## 80  1.3837769 0.5347080    2.149188  4.121275
## 81  1.3385135 0.5234094    2.098237  3.717707
## 82  1.4513903 0.5429308    2.187852  5.349775
## 83  0.9429183 0.4238292    1.735597  1.659647
## 84  1.4490505 0.5318678    2.136149  4.890459
## 85  1.7895111 0.6617778    2.956636  4.931276
## 86  1.9656594 0.6964237    3.294065  7.342929
## 87  1.5788314 0.5921019    2.451592  5.082987
## 88  1.5022623 0.5775047    2.366890  4.281400
## 89  1.4904512 0.5791837    2.376334  4.257961
## 90  1.4962931 0.5853587    2.411723  4.223961
## 91  1.4059042 0.5223585    2.093621  4.342856
## 92  1.6478061 0.6010774    2.506752  5.369842
## 93  1.2990656 0.5181661    2.075404  3.173985
## 94  1.0930459 0.4416889    1.791116  2.544462
## 95  1.5543977 0.5807846    2.385409  4.545684
## 96  1.7289643 0.6220482    2.645840  6.048954
## 97  1.6443609 0.5677222    2.313327  7.602821
## 98  1.5495614 0.5565515    2.255053  5.778996
## 99  1.4031313 0.5220452    2.092248  5.414712
## 100 1.7226377 0.5952296    2.470536  7.383522
## 101 1.8883235 0.6635526    2.972233  6.606168
## 102 1.9303606 0.6503184    2.859744  7.981444
## 103 1.4901640 0.5207111    2.086424  7.105981
## 104 1.5787674 0.6158734    2.603309  3.601971
## 105 1.5705318 0.5447686    2.196685  6.964669
## 106 1.4178638 0.5379747    2.164384  4.071694
## 107 1.2837894 0.4695946    1.885350  4.557373
## 108 1.5926701 0.5760346    2.358683  5.109555
## 109 1.6081552 0.5821278    2.393076  5.469168
## 110 1.4082954 0.5175781    2.072874  4.810390
## 111 1.2234179 0.5051503    2.020815  3.335706
## 112 1.2528753 0.4630629    1.862416  4.489371
## 113 1.3100790 0.4962000    1.984915  4.415941
## 114 1.4507992 0.5478637    2.211722  4.311580
## 115 1.2323738 0.4671492    1.876698  4.062475
## 116 1.3797947 0.5235294    2.098765  3.967505
## 117 1.7293521 0.6090305    2.557744  6.848485
## 118 1.4854109 0.5106574    2.043558  7.622560
## 119 1.4658616 0.5185606    2.077105  6.259936
## 120 0.9669648 0.3701524    1.587686  3.401558
## 121 1.1795715 0.4405887    1.787593  4.396019
## 122 1.3672704 0.5040980    2.016527  5.031600
## 123 1.5350481 0.5602648    2.274096  5.661968
## 124 1.3090028 0.4785325    1.917665  4.672620
## 125 1.2126036 0.4391298    1.782944  5.506903
## 126 1.3047201 0.4572742    1.842551  6.274631
## 127 1.2040135 0.4574581    1.843175  3.935634
## 128 1.3639687 0.4963772    1.985613  4.959159
## 129 1.1349332 0.4396221    1.784510  3.588003
## 130 1.4067959 0.5037577    2.015144  5.297514
## 131 1.2394246 0.4670187    1.876238  3.882807
## 132 1.3198667 0.4978497    1.991436  3.497519
## 133 1.1560033 0.4104902    1.696325  5.042383
## 134 1.2668167 0.4898061    1.960039  4.164308
## 135 1.3633187 0.4876454    1.951773  5.639364
## 136 1.7184293 0.5769861    2.363989  8.224130
## 137 1.4705114 0.5051541    2.020831  6.763796
## 138 1.4134997 0.5225468    2.094446  4.777415
## 139 1.1657677 0.3938079    1.649642  5.769443
## 140 1.4382794 0.5078587    2.031937  5.149475
## 141 1.6241226 0.5466634    2.205866  7.797693
## 142 0.9479850 0.3572150    1.555730  3.904973
## 143 0.7890530 0.3180077    1.466292  2.305868
## 144 1.4545973 0.4964437    1.985875  5.871991
## 145 1.4366209 0.4751863    1.905438  7.955229
## 146 1.3270892 0.4889245    1.956658  4.394232
## 147 1.1408832 0.4303415    1.755438  3.450861
## 148 0.9003295 0.3341988    1.501950  3.204387
## 149 0.9349969 0.3763435    1.603447  2.846412
## 150 1.1903317 0.4262014    1.742772  4.802582
## 151 1.5749987 0.5111597    2.045658  8.995673
## 152 1.5003813 0.5117908    2.048302  6.871785
## 153 1.0810697 0.4075500    1.687906  3.756033
## 154 1.5488876 0.5420794    2.183785  6.745033
## 155 1.5759803 0.5161432    2.066727  8.446768
## 156 0.9573217 0.3428280    1.521672  3.800210
## 157 1.2080605 0.4343764    1.767960  4.713438
## 158 1.1433944 0.4120049    1.700695  4.046218
## 159 1.0993195 0.4084227    1.690396  4.379899
## 160 1.1174376 0.3921670    1.645189  5.104476
## 161 1.1019272 0.4009981    1.669444  4.033430
## 162 0.9943400 0.3344865    1.502599  5.926942
## 163 1.0646499 0.4009879    1.669415  4.014640
## 164 1.6364990 0.5320156    2.136824  8.614307
## 165 1.2202461 0.4317769    1.759872  5.345689
## 166 1.3055935 0.4590037    1.848441  5.647828
## 167 1.2551904 0.4250527    1.739290  6.075704
## 168 1.3908337 0.4708680    1.889888  6.342263
## 169 1.5319849 0.5338560    2.145260  6.915739
## 170 0.9827985 0.3637872    1.571801  3.675202
## 171 0.9091475 0.3495779    1.537463  2.750911
## 172 1.3284004 0.4747466    1.903843  4.849896
## 173 1.0280144 0.3657520    1.576670  4.653569
## 174 1.0155201 0.3795146    1.611642  3.632950
## 175 0.8256294 0.2896172    1.407692  3.708700
## 176 1.1407526 0.4219639    1.729996  4.551192
## 177 1.1095773 0.4203375    1.725142  3.574193
## 178 0.9370976 0.3414050    1.518384  3.328464
## 179 1.1957853 0.4325555    1.762287  4.513385
## 180 0.9132408 0.3431800    1.522487  3.010270
## 181 0.9268789 0.3372680    1.508906  3.314740
## 182 0.9070407 0.3009006    1.430412  5.730437
## 183 0.9184634 0.3442333    1.524933  3.296916
## 184 1.0505849 0.3724143    1.593408  4.218990
## 185 1.1130730 0.3960015    1.655633  4.513385
## 186 1.1475092 0.3812779    1.616235  5.923958
## 187 0.8861070 0.3234457    1.478078  3.292541
## 188 1.1679846 0.4113292    1.698742  4.471297
## 189 1.3269686 0.4541780    1.832099  6.113701
## 190 1.4025824 0.4803133    1.924236  6.406272
## 191 0.7316835 0.2576901    1.347146  3.014282
## 192 1.0732068 0.3719623    1.592261  5.165024
## 193 1.0880649 0.3935743    1.649007  4.453883
## 194 0.8289934 0.2681985    1.366491  5.647828
## 195 0.8676827 0.3192152    1.468893  2.957055
## 196 1.1044350 0.3765746    1.604041  6.182347
## 197 1.0755235 0.3795178    1.611650  4.772041
## 198 1.0144104 0.3771101    1.605420  3.494964
## 199 1.1081660 0.3910828    1.642259  4.741473
## 200 1.0762302 0.3889234    1.636456  4.409155
## 201 1.0511711 0.3824131    1.619205  4.409155
## 202 0.9282779 0.3371236    1.508577  3.815290
## 203 1.2321592 0.4337904    1.766130  4.655622
## 204 0.7472004 0.2852947    1.399178  2.103496
## 205 0.8967415 0.3188013    1.468000  3.815290
## 206 0.9696155 0.3529198    1.545403  3.778652
## 207 0.4879819 0.1674397    1.201114  2.415700
## 208 1.2601278 0.4239074    1.735832  6.323830
## 209 1.2607309 0.4304125    1.755657  5.950260
## 210 1.0734488 0.3842678    1.624083  4.993876
summary(Div.Ind)
##     shannon         simpson        inv.simpson        fish.alp     
##  Min.   :0.488   Min.   :0.1674   Min.   : 1.201   Min.   : 1.660  
##  1st Qu.:1.171   1st Qu.:0.4272   1st Qu.: 1.746   1st Qu.: 4.214  
##  Median :1.453   Median :0.5277   Median : 2.117   Median : 5.112  
##  Mean   :1.559   Mean   :0.5678   Mean   : 3.192   Mean   : 5.521  
##  3rd Qu.:1.891   3rd Qu.:0.6944   3rd Qu.: 3.272   3rd Qu.: 6.638  
##  Max.   :2.854   Max.   :0.9286   Max.   :14.011   Max.   :11.971
env_Ohio
## # A tibble: 210 × 6
##     ...1 all_stem_den big_stem_bas can_covr_mea can_heig_mea Habitat  
##    <dbl>        <dbl>        <dbl>        <dbl>        <dbl> <chr>    
##  1     1       11017.         29.1         92.2         22.6 dry-oak  
##  2     2        3200.         25.9         92.7         22.9 dry-oak  
##  3     3        7975.         29.5         94.3         24.3 dry-oak  
##  4     4       15294.         27.0         97.9         18.9 dry-mesic
##  5     5       17181.         36.2         97.9         26   dry-mesic
##  6     6        7018.         27.7         91.7         26.8 wet-mesic
##  7     7       26318.         38.1         92.2         15.5 dry-oak  
##  8     8       10459.         26.9         94.3         29   dry-mesic
##  9     9        5366.         25.8         99.5         27.7 wet-mesic
## 10    10       21003.         13.9         99.5         21.2 dry-oak  
## # ℹ 200 more rows
Ohio.env.Div<-cbind.data.frame(env_Ohio, Div.Ind)
Ohio.env.Div
##     ...1 all_stem_den big_stem_bas can_covr_mea can_heig_mea   Habitat
## 1      1   11016.5982    29.104188     92.20000    22.600000   dry-oak
## 2      2    3199.9803    25.923036     92.72000    22.900000   dry-oak
## 3      3    7974.6286    29.470877     94.28000    24.300000   dry-oak
## 4      4   15294.4149    27.040870     97.92000    18.900000 dry-mesic
## 5      5   17180.6870    36.192879     97.92000    26.000000 dry-mesic
## 6      6    7017.9108    27.663756     91.68000    26.800000 wet-mesic
## 7      7   26317.7188    38.088936     92.20000    15.500000   dry-oak
## 8      8   10458.7869    26.917779     94.28000    29.000000 dry-mesic
## 9      9    5365.8287    25.780004     99.48000    27.700000 wet-mesic
## 10    10   21003.0879    13.873316     99.48000    21.200000   dry-oak
## 11    11    5949.9095    35.891250     97.40000    24.200000 dry-mesic
## 12    12   20158.9946    19.835275     92.72000    21.000000 dry-mesic
## 13    13   12813.3328    16.903017     89.08000    12.600000   dry-oak
## 14    14    8877.0193    38.192311     97.40000    24.650000   dry-oak
## 15    15    8100.6114    18.553685     97.92000    24.500000 wet-mesic
## 16    16   13376.1148    18.041125     94.28000    26.100000 dry-mesic
## 17    17    7160.8282    22.766701      2.76000    25.300000   dry-oak
## 18    18   27098.1501    23.189384     98.44000    19.300000   dry-oak
## 19    19   10749.4861    21.345035     97.40000    19.900000   dry-oak
## 20    20   19115.4746    24.594643     90.64000    20.650000   dry-oak
## 21    21    3626.1817    37.749615     95.84000    39.400000 dry-mesic
## 22    22   12482.1648    34.777772     90.64000    22.000000   dry-oak
## 23    23   12747.9358    23.063846     99.48000    20.600000   dry-oak
## 24    24    4766.6014    19.451376     97.92000    26.500000 dry-mesic
## 25    25    6854.9825    19.273113     94.28000    28.600000 dry-mesic
## 26    26    4015.1219    29.121760     95.32000    23.000000 dry-mesic
## 27    27   10748.1450    35.978796     91.16000    26.000000 dry-mesic
## 28    28   30600.9000    26.873688     93.24000    27.100000 dry-mesic
## 29    29   31500.6084    20.085862     91.16000    27.000000   dry-oak
## 30    30    4498.1482    23.795657     95.84000    25.700000 dry-mesic
## 31    31    8189.2017    45.056406     93.24000    30.100000 dry-mesic
## 32    32    7530.3359    26.487734     50.08000    25.900000   dry-oak
## 33    33   11090.0421    28.899434     96.88000    17.900000   dry-oak
## 34    34    5601.3067    39.079333     91.16000    27.500000 dry-mesic
## 35    35   15536.5986    48.183433     88.04000    28.100000 wet-mesic
## 36    36   20287.6597    25.212928     92.20000    23.500000 wet-mesic
## 37    37    3350.8915    36.182914     98.44000    22.100000   dry-oak
## 38    38    7313.0805    33.735482     93.24000    30.000000 dry-mesic
## 39    39    4881.4612    21.474137     96.36000    25.900000 wet-mesic
## 40    40     706.0402    19.170207     92.72000    30.700000 wet-mesic
## 41    41   25250.6116    53.206300     94.80000    25.500000 wet-mesic
## 42    42    8764.8418    39.151999     94.80000    30.700000 dry-mesic
## 43    43   18382.9805    24.533023     90.98667    18.966667   dry-oak
## 44    44    4028.9272    22.170667     96.36000    29.100000 wet-mesic
## 45    45   17514.1433    15.353375     98.96000    30.100000 wet-mesic
## 46    46    8889.4835    21.182268     98.96000    23.600000   dry-oak
## 47    47    3555.6827    28.542819     96.88000    30.750000 dry-mesic
## 48    48   14417.3464    42.396967     97.92000    20.300000 dry-mesic
## 49    49    4982.5157    31.442504     95.49333    27.933333 wet-mesic
## 50    50   10381.3196    31.795863     93.76000    29.300000 dry-mesic
## 51    51    8099.2703    40.314825     94.80000    25.100000 wet-mesic
## 52    52   13399.7020    56.411250     87.00000    32.000000 wet-mesic
## 53    53   13295.9652    29.939159     97.92000    23.500000   dry-oak
## 54    54   14722.7982    16.751400     95.32000    13.600000 wet-mesic
## 55    55   16610.4115    16.342020     97.92000    18.400000   dry-oak
## 56    56   13037.2939     9.689551     80.76000    15.200000 wet-mesic
## 57    57   14753.4850    44.522594     95.84000    29.000000 dry-mesic
## 58    58    9948.1499    21.641789     90.12000    22.250000   dry-oak
## 59    59   11920.3301    38.436105     92.20000    21.500000   dry-oak
## 60    60   17032.9052     9.486563     92.02667    10.433333 dry-mesic
## 61    61    5143.2087    32.615201     96.88000    20.100000   dry-oak
## 62    62    9632.5223    19.721298     91.68000    21.200000 wet-mesic
## 63    63   12353.4997    36.139811     96.36000    26.700000 dry-mesic
## 64    64   18680.9638    28.144823     94.80000    21.400000   dry-oak
## 65    65   28685.6761    20.426766     96.88000    19.000000 wet-mesic
## 66    66    4594.0226    31.083216     93.41333    33.716667 dry-mesic
## 67    67   11257.3095    31.371081     96.18667    27.600000 wet-mesic
## 68    68   10814.4893    25.454088     94.28000    25.500000 wet-mesic
## 69    69    2156.4604    27.053685     93.76000    26.800000 wet-mesic
## 70    70   12417.1617    42.563004     89.08000    26.050000 dry-mesic
## 71    71   10814.4893    28.819690     94.80000    24.900000 dry-mesic
## 72    72   22986.3911    23.350458     97.40000    20.900000 wet-mesic
## 73    73    2193.8529    33.285927     95.32000    26.000000 dry-mesic
## 74    74    5047.9127    46.833699     93.76000    23.500000   dry-oak
## 75    75    7398.9885    34.091804     91.68000    28.200000 wet-mesic
## 76    76   13834.2128    33.181279     87.00000    25.000000   dry-oak
## 77    77   10940.4721    35.020675     95.32000    23.050000 wet-mesic
## 78    78   10146.6293    15.161798     85.96000    12.400000 wet-mesic
## 79    79   13045.4720    23.485101     88.56000    23.133333   dry-oak
## 80    80    3673.2248    34.509633     97.40000    30.600000 wet-mesic
## 81    81    7181.7331    30.754454     82.84000    20.250000   dry-oak
## 82    82   23264.6262    26.316107     90.12000    20.500000   dry-oak
## 83    83   11640.7538    20.824820     51.64000    25.800000 dry-mesic
## 84    84   28181.7449    19.802207     90.12000    24.700000   dry-oak
## 85    85   28195.5502    23.732624     95.32000    23.300000   dry-oak
## 86    86   10562.5236    19.380590     94.80000    35.250000 wet-mesic
## 87    87   21151.3168    23.157354     97.40000    24.800000   dry-oak
## 88    88   18697.4514    23.077747     91.16000    27.700000   dry-oak
## 89    89   14468.1503    21.461322     94.80000    26.100000 dry-mesic
## 90    90    5830.4480    32.953344     89.60000    27.200000 wet-mesic
## 91    91   20518.7204    15.732301     97.40000    23.200000   dry-oak
## 92    92   26653.8574    32.208053     95.32000    25.300000   dry-oak
## 93    93    3822.7948    40.413335     95.32000    28.750000 dry-mesic
## 94    94   15410.3000    32.868810     87.86667    19.833333   dry-oak
## 95    95   13920.1208    33.785986     95.84000    23.150000   dry-oak
## 96    96    6028.7179    16.245712     97.40000    27.100000 dry-mesic
## 97    97   33686.9678    21.094222     75.04000    21.400000   dry-oak
## 98    98    9793.2153    17.508008     94.28000    18.900000 dry-mesic
## 99    99   14101.3249    27.096552     95.32000    29.150000 dry-mesic
## 100  100   16834.6352    29.316522     88.90667    31.933333 wet-mesic
## 101  101   21995.4101    21.477122     96.88000    24.500000 dry-mesic
## 102  102    7744.9090    17.898485     92.72000    26.700000 dry-mesic
## 103  103    9593.7887    16.635563     91.68000    25.050000   dry-oak
## 104  104   24285.9001    24.399258     93.76000    30.150000 wet-mesic
## 105  105   12615.4317    25.290861     95.14667    34.466667 dry-mesic
## 106  106   10258.0191    25.972316     94.80000    35.850000 wet-mesic
## 107  107   13179.7643    25.789784     92.72000    22.200000   dry-oak
## 108  108   17257.2602    53.510899     90.98667    32.366667   dry-oak
## 109  109    7756.0321    44.066147     92.20000    27.550000   dry-oak
## 110  110   11881.5964    45.079313     96.36000    27.250000   dry-oak
## 111  111    2386.1799    25.552304     97.40000    26.150000   dry-oak
## 112  112    4528.4411    35.258899     90.12000    36.450000 dry-mesic
## 113  113    6856.3237    28.086362     94.80000    22.650000 dry-mesic
## 114  114    7451.5275    27.317615     93.24000    25.400000   dry-oak
## 115  115   15007.7390    17.066930     93.24000    32.000000 dry-mesic
## 116  116   14483.2968    22.782755     92.72000    28.300000 dry-mesic
## 117  117   16589.9004    17.396458     91.68000    19.050000   dry-oak
## 118  118    8625.0537    30.562182     91.16000    24.850000 dry-mesic
## 119  119    7195.5384    31.069054     92.72000    20.350000   dry-oak
## 120  120   28248.0891    36.652811     93.24000    24.800000   dry-oak
## 121  121   21141.5349    27.194025     96.88000    21.250000   dry-oak
## 122  122   23623.0109    33.590943     89.08000    26.900000 dry-mesic
## 123  123    8387.0247    19.212298     87.86667    23.933333 wet-mesic
## 124  124    9283.9195    18.587863     90.29333    18.000000 dry-mesic
## 125  125    9900.0814    16.116506     91.33333    26.466667 wet-mesic
## 126  126    8503.3569    11.243813     89.25333    17.433333   dry-oak
## 127  127   18974.2141     4.049234     62.90667    10.300000   dry-oak
## 128  128   10092.9869    15.391874     82.66667    21.100000 wet-mesic
## 129  129   20349.2177    11.932597     96.53333    14.633333   dry-oak
## 130  130    9126.0401    20.767856     88.90667    21.633333 dry-mesic
## 131  131   12472.9081    20.210830     80.93333    14.666667 wet-mesic
## 132  132   12232.5657    14.047589     81.97333    18.633333 dry-mesic
## 133  133   14276.2703    35.419747     86.48000    24.533333   dry-oak
## 134  134    9470.3037    19.816457     74.00000    19.933333   dry-oak
## 135  135   16205.4308    28.896716     87.17333    24.266667 dry-mesic
## 136  136   24215.5856    16.347861     58.74667    21.816667 dry-mesic
## 137  137   35732.9076    13.770257     75.38667    19.266667   dry-oak
## 138  138   17616.9860     9.126119     27.20000     5.366667   dry-oak
## 139  139   40402.3467    16.039614     38.64000    16.100000   dry-oak
## 140  140   12669.1273    21.789901     92.02667    23.200000 wet-mesic
## 141  141   12628.4742    24.606410     40.72000    13.400000 wet-mesic
## 142  142    4672.6465    41.939241     89.60000    25.966667   dry-oak
## 143  143    4487.4191    19.608511     76.77333    15.466667   dry-oak
## 144  144   21179.3744    15.307731     86.48000    13.800000   dry-oak
## 145  145   21661.6911     5.284569     16.10667    19.433333   dry-oak
## 146  146   14327.0210    24.792838     90.64000    28.866667 wet-mesic
## 147  147   17115.3681    31.592483     81.28000    21.266667 wet-mesic
## 148  148    7468.4090    18.921078     89.25333    16.133333   dry-oak
## 149  149   11978.1804    19.483913     85.44000    22.300000   dry-oak
## 150  150   12349.7921    24.521517     64.64000    24.666667   dry-oak
## 151  151    7753.3498    16.900129     84.74667    17.066667 wet-mesic
## 152  152   13473.4616    22.683557     90.64000    17.500000 dry-mesic
## 153  153    6794.9750    16.587686     91.68000    14.666667 dry-mesic
## 154  154   18811.8642    28.354145     95.49333    23.733333 wet-mesic
## 155  155   17860.6422    29.345955     85.44000    27.133333 wet-mesic
## 156  156   19429.4453    14.661126     68.80000    16.700000 wet-mesic
## 157  157   13096.2228    19.571580     92.02667    17.700000   dry-oak
## 158  158   10633.1007    24.871309     87.52000    17.223333   dry-oak
## 159  159    9839.4956    19.421346     81.62667    17.233333 wet-mesic
## 160  160   13903.5019    25.149450     90.98667    25.000000   dry-oak
## 161  161    7019.0674    20.364822     87.52000    19.366667   dry-oak
## 162  162   12054.4910    10.341282     85.78667    19.666667 wet-mesic
## 163  163    9442.6931    18.427794     84.05333    17.933333   dry-oak
## 164  164   27387.6927    12.906460     75.04000    19.066667   dry-oak
## 165  165   20938.6631    28.784504     92.02667    21.583333 dry-mesic
## 166  166   14979.3656    28.462722     92.37333    22.900000 dry-mesic
## 167  167    4543.2718    17.290489     67.76000    26.066667 wet-mesic
## 168  168   13498.3899    13.315438     74.00000    15.866667 wet-mesic
## 169  169   21045.5292    19.118196     97.22667    23.583333 wet-mesic
## 170  170    2858.6616    16.429112     94.10667    20.600000   dry-oak
## 171  171    3457.3105    29.545064     90.29333    33.333333 wet-mesic
## 172  172   12586.9270    29.319863     95.49333    23.833333   dry-oak
## 173  173    9569.1229    12.142800     62.56000    12.733333   dry-oak
## 174  174    9627.2890    20.551505     78.16000    15.600000 dry-mesic
## 175  175   19614.6728     4.411145     61.86667    10.000000   dry-oak
## 176  176   13437.5415    21.081519     76.08000    20.900000   dry-oak
## 177  177   14094.6192    19.109340     90.98667    17.466667 wet-mesic
## 178  178   15720.1692    23.404051     88.90667    21.133333   dry-oak
## 179  179   11112.4726    26.803913     89.94667    24.066667   dry-oak
## 180  180   13826.0347    23.429406     90.29333    20.566667   dry-oak
## 181  181   14837.4735    29.854889     81.97333    17.733333 wet-mesic
## 182  182   17745.2040    10.091687     57.70667    13.600000   dry-oak
## 183  183   13238.3775    27.585043     74.69333    17.466667 dry-mesic
## 184  184   10748.0137    23.852671     86.13333    27.166667 wet-mesic
## 185  185   26095.1520     9.328745     63.25333     8.400000   dry-oak
## 186  186   18309.9837     6.276581     95.84000     9.800000 dry-mesic
## 187  187   10487.3697    20.320281     86.13333    15.266667   dry-oak
## 188  188   12817.8032    18.314022     84.74667    17.790000 dry-mesic
## 189  189   29962.9391    20.851176     80.58667    21.000000   dry-oak
## 190  190   22338.0700     5.870011     34.13333    18.400000   dry-oak
## 191  191   12999.4543    16.983580     77.46667    11.133333   dry-oak
## 192  192   10126.4872    12.635413     57.01333     8.966667   dry-oak
## 193  193    4120.0153    33.505456     90.29333    23.533333   dry-oak
## 194  194   23438.8619     3.375055     42.10667     7.400000   dry-oak
## 195  195    7526.3125    45.167039     88.21333    25.133333   dry-oak
## 196  196   25378.3826    12.674609     80.58667    18.843333 dry-mesic
## 197  197   12485.4255    40.605261     82.32000    23.156667 dry-mesic
## 198  198   11716.1172    29.021132     92.72000    21.166667 wet-mesic
## 199  199   18906.2129    23.600876     74.34667    19.816667 dry-mesic
## 200  200   12960.5894    25.533512     90.64000    21.410000   dry-oak
## 201  201    2509.9275    24.636020     80.93333    22.066667 dry-mesic
## 202  202   14193.1758    26.843018     86.48000    20.466667 dry-mesic
## 203  203   12231.6717    29.587419     94.80000    27.066667 wet-mesic
## 204  204    6738.3345    40.580331     87.86667    24.466667   dry-oak
## 205  205    9116.2050    17.105751     83.36000    11.500000   dry-oak
## 206  206    7415.6074    25.667489     96.18667    17.433333   dry-oak
## 207  207   16024.9895    19.691871     95.84000    24.000000   dry-oak
## 208  208   19835.8204     7.986256     70.53333    16.200000   dry-oak
## 209  209   16300.4110     8.486471     51.81333    13.533333 wet-mesic
## 210  210   28912.0038    22.976962     36.56000    20.866667 wet-mesic
##       shannon   simpson inv.simpson  fish.alp
## 1   1.8917399 0.7923875    4.816667  5.642066
## 2   2.5122269 0.9023669   10.242424  7.896600
## 3   2.5759511 0.9141051   11.642140  6.488217
## 4   2.5306341 0.8955078    9.570093  7.565648
## 5   2.5813695 0.9112426   11.266667  7.896600
## 6   2.5851922 0.9100346   11.115385  7.275326
## 7   2.8543118 0.9286265   14.010811 11.680602
## 8   2.6869090 0.9171429   12.068966  8.578491
## 9   2.6290101 0.9131944   11.520000  7.703286
## 10  2.4471082 0.8927116    9.320675  7.612305
## 11  2.4514742 0.8897290    9.068571  7.916266
## 12  2.6624436 0.9137893   11.599483  8.839863
## 13  2.4801533 0.8793388    8.287671  9.316717
## 14  2.4524231 0.8758573    8.055249  9.454690
## 15  2.4446459 0.8830959    8.554023  7.814550
## 16  2.5985159 0.8967347    9.683794  9.354043
## 17  2.1580650 0.8536155    6.831325  4.970140
## 18  2.0290131 0.8264046    5.760522  4.970140
## 19  2.3439848 0.8624852    7.271945  6.781181
## 20  2.3677896 0.8617998    7.235880  7.415136
## 21  2.6230152 0.8762865    8.083192 11.156386
## 22  2.1818266 0.8433163    6.382284  5.110244
## 23  2.3999167 0.8526139    6.784902  9.198544
## 24  2.3497389 0.8573088    7.008138  7.705935
## 25  2.1019228 0.8130987    5.350417  5.389613
## 26  2.1558323 0.8250548    5.716075  5.679263
## 27  2.1542864 0.8136574    5.366460  6.377738
## 28  2.2277483 0.8264701    5.762696  6.055108
## 29  2.3289428 0.8491358    6.628478  6.765902
## 30  1.8809617 0.7441406    3.908397  4.927590
## 31  2.4493401 0.8553590    6.913669  7.224522
## 32  2.3266446 0.8238062    5.675570  8.920683
## 33  2.2270348 0.8106576    5.281437  7.027652
## 34  1.9964440 0.7780612    4.505747  6.430958
## 35  2.0417370 0.7801904    4.549391  5.381526
## 36  2.4618897 0.8531268    6.808596  8.310228
## 37  2.1320689 0.7831946    4.612431  6.891754
## 38  2.2922720 0.8068698    5.177854  8.984035
## 39  1.5430953 0.6477631    2.838998  3.690822
## 40  2.0627115 0.7808963    4.564050  5.538659
## 41  2.6546098 0.8616564    7.228380 11.970603
## 42  2.3331170 0.7890625    4.740741  9.586698
## 43  1.7535317 0.6954194    3.283203  4.350582
## 44  1.8944254 0.7339693    3.758965  5.112918
## 45  1.9905828 0.7269136    3.661844  6.765902
## 46  1.8202840 0.6911844    3.238179  5.254828
## 47  1.7588983 0.6822566    3.147193  4.718079
## 48  1.2049611 0.5081633    2.033195  3.192559
## 49  2.1437307 0.7527571    4.044605  7.976087
## 50  2.0737966 0.7694515    4.337483  5.576450
## 51  1.6534827 0.6567901    2.913669  3.718729
## 52  2.1077805 0.7160000    3.521127  9.354332
## 53  1.9625262 0.7361333    3.789793  5.646225
## 54  1.7738141 0.6479362    2.840394  6.096825
## 55  2.0210011 0.6747189    3.074264 10.837223
## 56  1.8759969 0.6863983    3.188759  5.798205
## 57  1.4522335 0.5658574    2.303391  4.212510
## 58  1.5969316 0.6459054    2.824104  4.392532
## 59  1.4228140 0.5912465    2.446462  3.217794
## 60  2.2233828 0.7710506    4.367778  7.064403
## 61  1.6941131 0.6341785    2.733574  4.829737
## 62  2.2527933 0.7706531    4.360207  8.991549
## 63  1.4705667 0.5559896    2.252199  4.986195
## 64  1.4376671 0.5789628    2.375087  3.546422
## 65  1.6556983 0.6272199    2.682547  4.711590
## 66  2.0069460 0.7150879    3.509854  6.648249
## 67  1.7242396 0.6403025    2.780114  5.988234
## 68  1.6262868 0.6024793    2.515593  6.114229
## 69  1.9412633 0.7088021    3.434090  6.065899
## 70  1.2809809 0.4988000    1.995211  3.560586
## 71  1.3010760 0.4968140    1.987337  3.546422
## 72  1.5318794 0.5861804    2.416512  4.192993
## 73  1.5531853 0.5995339    2.497090  5.427384
## 74  1.5204937 0.5778835    2.369014  5.030672
## 75  1.7648837 0.6374853    2.758509  5.726751
## 76  1.6943305 0.6346974    2.737457  5.687885
## 77  1.7000499 0.6039282    2.524795  6.259378
## 78  1.3091997 0.4959013    1.983739  4.238987
## 79  1.5702975 0.5915848    2.448489  4.858106
## 80  1.3837769 0.5347080    2.149188  4.121275
## 81  1.3385135 0.5234094    2.098237  3.717707
## 82  1.4513903 0.5429308    2.187852  5.349775
## 83  0.9429183 0.4238292    1.735597  1.659647
## 84  1.4490505 0.5318678    2.136149  4.890459
## 85  1.7895111 0.6617778    2.956636  4.931276
## 86  1.9656594 0.6964237    3.294065  7.342929
## 87  1.5788314 0.5921019    2.451592  5.082987
## 88  1.5022623 0.5775047    2.366890  4.281400
## 89  1.4904512 0.5791837    2.376334  4.257961
## 90  1.4962931 0.5853587    2.411723  4.223961
## 91  1.4059042 0.5223585    2.093621  4.342856
## 92  1.6478061 0.6010774    2.506752  5.369842
## 93  1.2990656 0.5181661    2.075404  3.173985
## 94  1.0930459 0.4416889    1.791116  2.544462
## 95  1.5543977 0.5807846    2.385409  4.545684
## 96  1.7289643 0.6220482    2.645840  6.048954
## 97  1.6443609 0.5677222    2.313327  7.602821
## 98  1.5495614 0.5565515    2.255053  5.778996
## 99  1.4031313 0.5220452    2.092248  5.414712
## 100 1.7226377 0.5952296    2.470536  7.383522
## 101 1.8883235 0.6635526    2.972233  6.606168
## 102 1.9303606 0.6503184    2.859744  7.981444
## 103 1.4901640 0.5207111    2.086424  7.105981
## 104 1.5787674 0.6158734    2.603309  3.601971
## 105 1.5705318 0.5447686    2.196685  6.964669
## 106 1.4178638 0.5379747    2.164384  4.071694
## 107 1.2837894 0.4695946    1.885350  4.557373
## 108 1.5926701 0.5760346    2.358683  5.109555
## 109 1.6081552 0.5821278    2.393076  5.469168
## 110 1.4082954 0.5175781    2.072874  4.810390
## 111 1.2234179 0.5051503    2.020815  3.335706
## 112 1.2528753 0.4630629    1.862416  4.489371
## 113 1.3100790 0.4962000    1.984915  4.415941
## 114 1.4507992 0.5478637    2.211722  4.311580
## 115 1.2323738 0.4671492    1.876698  4.062475
## 116 1.3797947 0.5235294    2.098765  3.967505
## 117 1.7293521 0.6090305    2.557744  6.848485
## 118 1.4854109 0.5106574    2.043558  7.622560
## 119 1.4658616 0.5185606    2.077105  6.259936
## 120 0.9669648 0.3701524    1.587686  3.401558
## 121 1.1795715 0.4405887    1.787593  4.396019
## 122 1.3672704 0.5040980    2.016527  5.031600
## 123 1.5350481 0.5602648    2.274096  5.661968
## 124 1.3090028 0.4785325    1.917665  4.672620
## 125 1.2126036 0.4391298    1.782944  5.506903
## 126 1.3047201 0.4572742    1.842551  6.274631
## 127 1.2040135 0.4574581    1.843175  3.935634
## 128 1.3639687 0.4963772    1.985613  4.959159
## 129 1.1349332 0.4396221    1.784510  3.588003
## 130 1.4067959 0.5037577    2.015144  5.297514
## 131 1.2394246 0.4670187    1.876238  3.882807
## 132 1.3198667 0.4978497    1.991436  3.497519
## 133 1.1560033 0.4104902    1.696325  5.042383
## 134 1.2668167 0.4898061    1.960039  4.164308
## 135 1.3633187 0.4876454    1.951773  5.639364
## 136 1.7184293 0.5769861    2.363989  8.224130
## 137 1.4705114 0.5051541    2.020831  6.763796
## 138 1.4134997 0.5225468    2.094446  4.777415
## 139 1.1657677 0.3938079    1.649642  5.769443
## 140 1.4382794 0.5078587    2.031937  5.149475
## 141 1.6241226 0.5466634    2.205866  7.797693
## 142 0.9479850 0.3572150    1.555730  3.904973
## 143 0.7890530 0.3180077    1.466292  2.305868
## 144 1.4545973 0.4964437    1.985875  5.871991
## 145 1.4366209 0.4751863    1.905438  7.955229
## 146 1.3270892 0.4889245    1.956658  4.394232
## 147 1.1408832 0.4303415    1.755438  3.450861
## 148 0.9003295 0.3341988    1.501950  3.204387
## 149 0.9349969 0.3763435    1.603447  2.846412
## 150 1.1903317 0.4262014    1.742772  4.802582
## 151 1.5749987 0.5111597    2.045658  8.995673
## 152 1.5003813 0.5117908    2.048302  6.871785
## 153 1.0810697 0.4075500    1.687906  3.756033
## 154 1.5488876 0.5420794    2.183785  6.745033
## 155 1.5759803 0.5161432    2.066727  8.446768
## 156 0.9573217 0.3428280    1.521672  3.800210
## 157 1.2080605 0.4343764    1.767960  4.713438
## 158 1.1433944 0.4120049    1.700695  4.046218
## 159 1.0993195 0.4084227    1.690396  4.379899
## 160 1.1174376 0.3921670    1.645189  5.104476
## 161 1.1019272 0.4009981    1.669444  4.033430
## 162 0.9943400 0.3344865    1.502599  5.926942
## 163 1.0646499 0.4009879    1.669415  4.014640
## 164 1.6364990 0.5320156    2.136824  8.614307
## 165 1.2202461 0.4317769    1.759872  5.345689
## 166 1.3055935 0.4590037    1.848441  5.647828
## 167 1.2551904 0.4250527    1.739290  6.075704
## 168 1.3908337 0.4708680    1.889888  6.342263
## 169 1.5319849 0.5338560    2.145260  6.915739
## 170 0.9827985 0.3637872    1.571801  3.675202
## 171 0.9091475 0.3495779    1.537463  2.750911
## 172 1.3284004 0.4747466    1.903843  4.849896
## 173 1.0280144 0.3657520    1.576670  4.653569
## 174 1.0155201 0.3795146    1.611642  3.632950
## 175 0.8256294 0.2896172    1.407692  3.708700
## 176 1.1407526 0.4219639    1.729996  4.551192
## 177 1.1095773 0.4203375    1.725142  3.574193
## 178 0.9370976 0.3414050    1.518384  3.328464
## 179 1.1957853 0.4325555    1.762287  4.513385
## 180 0.9132408 0.3431800    1.522487  3.010270
## 181 0.9268789 0.3372680    1.508906  3.314740
## 182 0.9070407 0.3009006    1.430412  5.730437
## 183 0.9184634 0.3442333    1.524933  3.296916
## 184 1.0505849 0.3724143    1.593408  4.218990
## 185 1.1130730 0.3960015    1.655633  4.513385
## 186 1.1475092 0.3812779    1.616235  5.923958
## 187 0.8861070 0.3234457    1.478078  3.292541
## 188 1.1679846 0.4113292    1.698742  4.471297
## 189 1.3269686 0.4541780    1.832099  6.113701
## 190 1.4025824 0.4803133    1.924236  6.406272
## 191 0.7316835 0.2576901    1.347146  3.014282
## 192 1.0732068 0.3719623    1.592261  5.165024
## 193 1.0880649 0.3935743    1.649007  4.453883
## 194 0.8289934 0.2681985    1.366491  5.647828
## 195 0.8676827 0.3192152    1.468893  2.957055
## 196 1.1044350 0.3765746    1.604041  6.182347
## 197 1.0755235 0.3795178    1.611650  4.772041
## 198 1.0144104 0.3771101    1.605420  3.494964
## 199 1.1081660 0.3910828    1.642259  4.741473
## 200 1.0762302 0.3889234    1.636456  4.409155
## 201 1.0511711 0.3824131    1.619205  4.409155
## 202 0.9282779 0.3371236    1.508577  3.815290
## 203 1.2321592 0.4337904    1.766130  4.655622
## 204 0.7472004 0.2852947    1.399178  2.103496
## 205 0.8967415 0.3188013    1.468000  3.815290
## 206 0.9696155 0.3529198    1.545403  3.778652
## 207 0.4879819 0.1674397    1.201114  2.415700
## 208 1.2601278 0.4239074    1.735832  6.323830
## 209 1.2607309 0.4304125    1.755657  5.950260
## 210 1.0734488 0.3842678    1.624083  4.993876
library("rstatix")
## Warning: package 'rstatix' was built under R version 4.4.2
## 
## Attaching package: 'rstatix'
## The following object is masked from 'package:stats':
## 
##     filter
Ohio.env.Div %>%
  group_by(Habitat) %>%
  get_summary_stats(shannon, type = "mean_sd")
## # A tibble: 3 × 5
##   Habitat   variable     n  mean    sd
##   <chr>     <fct>    <dbl> <dbl> <dbl>
## 1 dry-mesic shannon     60  1.67 0.527
## 2 dry-oak   shannon     95  1.46 0.507
## 3 wet-mesic shannon     55  1.61 0.452
Ohio.env.Div %>%
  group_by(Habitat) %>%
  get_summary_stats(simpson, type = "mean_sd")
## # A tibble: 3 × 5
##   Habitat   variable     n  mean    sd
##   <chr>     <fct>    <dbl> <dbl> <dbl>
## 1 dry-mesic simpson     60 0.606 0.18 
## 2 dry-oak   simpson     95 0.537 0.185
## 3 wet-mesic simpson     55 0.58  0.157
Ohio.env.Div %>%
  group_by(Habitat) %>%
  get_summary_stats(fish.alp, type = "mean_sd")
## # A tibble: 3 × 5
##   Habitat   variable     n  mean    sd
##   <chr>     <fct>    <dbl> <dbl> <dbl>
## 1 dry-mesic fish.alp    60  5.77  1.90
## 2 dry-oak   fish.alp    95  5.22  1.84
## 3 wet-mesic fish.alp    55  5.76  1.87
Ohio.env.Div %>%
  group_by(Habitat) %>%
  get_summary_stats(inv.simpson, type = "mean_sd")
## # A tibble: 3 × 5
##   Habitat   variable        n  mean    sd
##   <chr>     <fct>       <dbl> <dbl> <dbl>
## 1 dry-mesic inv.simpson    60  3.66  2.81
## 2 dry-oak   inv.simpson    95  2.98  2.44
## 3 wet-mesic inv.simpson    55  3.04  2.15
library("ggplot2")
## Warning: package 'ggplot2' was built under R version 4.4.2
Shanon.habitat<-ggplot(Ohio.env.Div, aes(x = Habitat, y = shannon, fill = Habitat)) +
  geom_boxplot() + 
  stat_summary(fun = mean, geom = "point", shape =21, size = 3, colour = "black", fill="yellow")

Simp.habitat<-ggplot(Ohio.env.Div, aes(x = Habitat, y = simpson, fill = Habitat)) +
  geom_boxplot() +
  stat_summary(fun = mean, geom = "point", shape =21, size = 3, colour = "black", fill="yellow")

InvSimp.habitat<-ggplot(Ohio.env.Div, aes(x = Habitat, y = inv.simpson, fill = Habitat)) +
  geom_boxplot()  + 
  stat_summary(fun = mean, geom = "point", shape =21, size = 3, colour = "black", fill="yellow")

Fish.habitat<-ggplot(Ohio.env.Div, aes(x = Habitat, y = fish.alp, fill = Habitat)) +
  geom_boxplot()  + 
  stat_summary(fun = mean, geom = "point", shape =21, size = 3, colour = "black", fill="yellow")
library("gridExtra")
## Warning: package 'gridExtra' was built under R version 4.4.2
grid.arrange(Shanon.habitat, Simp.habitat, InvSimp.habitat, Fish.habitat,
             nrow=2, ncol=2) 

grid.arrange(Shanon.habitat, Simp.habitat, InvSimp.habitat, Fish.habitat, nrow=2, ncol=2) 

shannon_aov <- aov(Ohio.env.Div$shannon ~ Ohio.env.Div$Habitat)

summary(shannon_aov)
##                       Df Sum Sq Mean Sq F value Pr(>F)  
## Ohio.env.Div$Habitat   2   1.66  0.8308   3.335 0.0375 *
## Residuals            207  51.57  0.2491                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
simpson_aov <- aov(Ohio.env.Div$simpson ~ Ohio.env.Div$Habitat)

summary(simpson_aov)
##                       Df Sum Sq Mean Sq F value Pr(>F)  
## Ohio.env.Div$Habitat   2  0.189 0.09426   3.015 0.0512 .
## Residuals            207  6.471 0.03126                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
inv.simpson_aov <- aov(Ohio.env.Div$inv.simpson ~ Ohio.env.Div$Habitat)

summary(inv.simpson_aov)
##                       Df Sum Sq Mean Sq F value Pr(>F)
## Ohio.env.Div$Habitat   2   18.2   9.109   1.481   0.23
## Residuals            207 1273.6   6.153
fish.alp_aov <- aov(Ohio.env.Div$fish.alp ~ Ohio.env.Div$Habitat)

summary(fish.alp_aov)
##                       Df Sum Sq Mean Sq F value Pr(>F)
## Ohio.env.Div$Habitat   2   15.4   7.687   2.206  0.113
## Residuals            207  721.4   3.485
distance_matrix<-vegdist(Ohio.env.Div[,10:10], method="bray", binary=FALSE)
adonis2(distance_matrix ~ Habitat, data=Ohio.env.Div)
## Permutation test for adonis under reduced model
## Permutation: free
## Number of permutations: 999
## 
## adonis2(formula = distance_matrix ~ Habitat, data = Ohio.env.Div)
##           Df SumOfSqs      R2      F Pr(>F)  
## Model      2   0.1221 0.02272 2.4058  0.063 .
## Residual 207   5.2533 0.97728                
## Total    209   5.3754 1.00000                
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
shannon_Tukey<-TukeyHSD(shannon_aov, conf.level=.95, ordered = TRUE)
shannon_Tukey
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
##     factor levels have been ordered
## 
## Fit: aov(formula = Ohio.env.Div$shannon ~ Ohio.env.Div$Habitat)
## 
## $`Ohio.env.Div$Habitat`
##                           diff          lwr       upr     p adj
## wet-mesic-dry-oak   0.14147731 -0.058161073 0.3411157 0.2180524
## dry-mesic-dry-oak   0.20188576  0.007587024 0.3961845 0.0396316
## dry-mesic-wet-mesic 0.06040845 -0.159546381 0.2803633 0.7935133
simpson_Tukey<-TukeyHSD(simpson_aov, conf.level=.95, ordered = TRUE)
simpson_Tukey
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
##     factor levels have been ordered
## 
## Fit: aov(formula = Ohio.env.Div$simpson ~ Ohio.env.Div$Habitat)
## 
## $`Ohio.env.Div$Habitat`
##                           diff           lwr       upr     p adj
## wet-mesic-dry-oak   0.04373077 -0.0269916799 0.1144532 0.3125264
## dry-mesic-dry-oak   0.06936909  0.0005382234 0.1382000 0.0477573
## dry-mesic-wet-mesic 0.02563832 -0.0522812927 0.1035579 0.7177126
inv.simpson_Tukey<-TukeyHSD(inv.simpson_aov, conf.level=.95, ordered = TRUE)
inv.simpson_Tukey 
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
##     factor levels have been ordered
## 
## Fit: aov(formula = Ohio.env.Div$inv.simpson ~ Ohio.env.Div$Habitat)
## 
## $`Ohio.env.Div$Habitat`
##                           diff        lwr      upr     p adj
## wet-mesic-dry-oak   0.06042921 -0.9317123 1.052571 0.9886676
## dry-mesic-dry-oak   0.67188001 -0.2937251 1.637485 0.2301703
## dry-mesic-wet-mesic 0.61145080 -0.4816572 1.704559 0.3853924
fish.alp_Tukey<-TukeyHSD(fish.alp_aov, conf.level=.95, ordered = TRUE)
fish.alp_Tukey
##   Tukey multiple comparisons of means
##     95% family-wise confidence level
##     factor levels have been ordered
## 
## Fit: aov(formula = Ohio.env.Div$fish.alp ~ Ohio.env.Div$Habitat)
## 
## $`Ohio.env.Div$Habitat`
##                            diff        lwr       upr     p adj
## wet-mesic-dry-oak   0.538910024 -0.2077869 1.2856069 0.2061985
## dry-mesic-dry-oak   0.547843742 -0.1788816 1.2745691 0.1789658
## dry-mesic-wet-mesic 0.008933717 -0.8137517 0.8316192 0.9996378
shannon_Tukey_plot <- as.data.frame(shannon_Tukey$`Ohio.env.Div$Habitat`)

shannon_Tukey_plot$comparison <- rownames(shannon_Tukey_plot)

shannon_Tukey_plot
##                           diff          lwr       upr      p adj
## wet-mesic-dry-oak   0.14147731 -0.058161073 0.3411157 0.21805244
## dry-mesic-dry-oak   0.20188576  0.007587024 0.3961845 0.03963162
## dry-mesic-wet-mesic 0.06040845 -0.159546381 0.2803633 0.79351326
##                              comparison
## wet-mesic-dry-oak     wet-mesic-dry-oak
## dry-mesic-dry-oak     dry-mesic-dry-oak
## dry-mesic-wet-mesic dry-mesic-wet-mesic
simpson_Tukey_plot <- as.data.frame(simpson_Tukey$`Ohio.env.Div$Habitat`)

simpson_Tukey_plot$comparison <- rownames(simpson_Tukey_plot)

simpson_Tukey_plot
##                           diff           lwr       upr      p adj
## wet-mesic-dry-oak   0.04373077 -0.0269916799 0.1144532 0.31252636
## dry-mesic-dry-oak   0.06936909  0.0005382234 0.1382000 0.04775727
## dry-mesic-wet-mesic 0.02563832 -0.0522812927 0.1035579 0.71771258
##                              comparison
## wet-mesic-dry-oak     wet-mesic-dry-oak
## dry-mesic-dry-oak     dry-mesic-dry-oak
## dry-mesic-wet-mesic dry-mesic-wet-mesic
inv.simpson_Tukey_plot <- as.data.frame(inv.simpson_Tukey$`Ohio.env.Div$Habitat`)

inv.simpson_Tukey_plot$comparison <- rownames(inv.simpson_Tukey_plot)

inv.simpson_Tukey_plot
##                           diff        lwr      upr     p adj
## wet-mesic-dry-oak   0.06042921 -0.9317123 1.052571 0.9886676
## dry-mesic-dry-oak   0.67188001 -0.2937251 1.637485 0.2301703
## dry-mesic-wet-mesic 0.61145080 -0.4816572 1.704559 0.3853924
##                              comparison
## wet-mesic-dry-oak     wet-mesic-dry-oak
## dry-mesic-dry-oak     dry-mesic-dry-oak
## dry-mesic-wet-mesic dry-mesic-wet-mesic
fish.alp_Tukey_plot <- as.data.frame(fish.alp_Tukey$`Ohio.env.Div$Habitat`)

fish.alp_Tukey_plot$comparison <- rownames(fish.alp_Tukey_plot)

fish.alp_Tukey_plot
##                            diff        lwr       upr     p adj
## wet-mesic-dry-oak   0.538910024 -0.2077869 1.2856069 0.2061985
## dry-mesic-dry-oak   0.547843742 -0.1788816 1.2745691 0.1789658
## dry-mesic-wet-mesic 0.008933717 -0.8137517 0.8316192 0.9996378
##                              comparison
## wet-mesic-dry-oak     wet-mesic-dry-oak
## dry-mesic-dry-oak     dry-mesic-dry-oak
## dry-mesic-wet-mesic dry-mesic-wet-mesic
shannon_Tukey_plot_result<-ggplot(shannon_Tukey_plot, aes(x = comparison, y = diff)) +
  geom_point() +
  geom_errorbar(aes(ymin = lwr, ymax = upr), width = 0.2) +
  labs(
    title = "Tukey HSD Test Results for Shannon Diversity",
    x = "Habitat Comparison",
    y = "Pairwise Difference in Mean"
  )

shannon_Tukey_plot_result

shannon_Tukey_plot_result_meandiff<-shannon_Tukey_plot_result + 
  theme(axis.text.x = element_text(angle = 45, hjust = 1)) + 
  geom_hline(yintercept=0, linetype="dashed", color = "red")

shannon_Tukey_plot_result_meandiff

simpson_Tukey_plot_result<-ggplot(simpson_Tukey_plot, aes(x = comparison, y = diff)) +
  geom_point() +
  geom_errorbar(aes(ymin = lwr, ymax = upr), width = 0.2) +
  labs(
    title = "Tukey HSD Test Results for Simpson Diversity",
    x = "Habitat Comparison",
    y = "Pairwise Difference in Mean"
  )

simpson_Tukey_plot_result

simpson_Tukey_plot_result_meandiff<-simpson_Tukey_plot_result + 
  theme(axis.text.x = element_text(angle = 45, hjust = 1)) + 
  geom_hline(yintercept=0, linetype="dashed", color = "red")

simpson_Tukey_plot_result_meandiff

inv.simpson_Tukey_plot_result<-ggplot(inv.simpson_Tukey_plot, aes(x = comparison, y = diff)) +
  geom_point() +
  geom_errorbar(aes(ymin = lwr, ymax = upr), width = 0.2) +
  labs(
    title = "Tukey HSD Test Results for inv.Simpson Diversity",
    x = "Habitat Comparison",
    y = "Pairwise Difference in Mean"
  )

inv.simpson_Tukey_plot_result

inv.simpson_Tukey_plot_result_meandiff<-inv.simpson_Tukey_plot_result + 
  theme(axis.text.x = element_text(angle = 45, hjust = 1)) + 
  geom_hline(yintercept=0, linetype="dashed", color = "red")

inv.simpson_Tukey_plot_result_meandiff

fish.alp_Tukey_plot_result<-ggplot(fish.alp_Tukey_plot, aes(x = comparison, y = diff)) +
  geom_point() +
  geom_errorbar(aes(ymin = lwr, ymax = upr), width = 0.2) +
  labs(
    title = "Tukey HSD Test Results for fish alpha Diversity",
    x = "Habitat Comparison",
    y = "Pairwise Difference in Mean"
  )

fish.alp_Tukey_plot_result

fish.alp_Tukey_plot_result_meandiff<-fish.alp_Tukey_plot_result + 
  theme(axis.text.x = element_text(angle = 45, hjust = 1)) + 
  geom_hline(yintercept=0, linetype="dashed", color = "red")

fish.alp_Tukey_plot_result_meandiff

grid.arrange(shannon_Tukey_plot_result_meandiff,Shanon.habitat,nrow=1, ncol=2)

grid.arrange(simpson_Tukey_plot_result_meandiff,Simp.habitat,nrow=1, ncol=2)

grid.arrange(inv.simpson_Tukey_plot_result_meandiff,InvSimp.habitat,nrow=1, ncol=2)

grid.arrange(fish.alp_Tukey_plot_result_meandiff,Fish.habitat,nrow=1, ncol=2)