pollution <- read.csv2("../R Scripts/avgpm25.csv")


colClasses = c("numeric", "character", "factor", "numeric", "numeric")

summary(pollution$pm25)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   3.383   8.549  10.047   9.836  11.356  18.441
pollution
##     ï..X      pm25  fips region  longitude latitude
## 1     NA  9.771185  1003   east  -87.74826 30.59278
## 2     NA  9.993817  1027   east  -85.84286 33.26581
## 3     NA 10.688618  1033   east  -87.72596 34.73148
## 4     NA 11.337424  1049   east  -85.79892 34.45913
## 5     NA 12.119764  1055   east  -86.03212 34.01860
## 6     NA 10.827805  1069   east  -85.35039 31.18973
## 7     NA 11.583928  1073   east  -86.82805 33.52787
## 8     NA 11.261996  1089   east  -86.58823 34.73079
## 9     NA  9.414423  1097   east  -88.13967 30.72226
## 10    NA 11.391494  1103   east  -86.91892 34.50702
## 11    NA 12.384795  1113   east  -85.10110 32.37600
## 12    NA 10.649500  1117   east  -86.69867 33.26912
## 13    NA 11.333821  1121   east  -86.17828 33.36850
## 14    NA 12.302436  1125   east  -87.51169 33.23560
## 15    NA 11.024508  1127   east  -87.28541 33.81989
## 16    NA 10.835515  5001   east  -91.42941 34.36000
## 17    NA 10.436235  5003   east  -91.78535 33.18542
## 18    NA 11.114743  5035   east  -90.27280 35.19770
## 19    NA 10.823157  5045   east  -92.37994 35.11948
## 20    NA 10.745523  5051   east  -93.07589 34.54531
## 21    NA 10.450351  5067   east  -91.22758 35.61326
## 22    NA 10.723959  5107   east  -90.76197 34.46837
## 23    NA 10.441737  5113   east  -94.25363 34.48549
## 24    NA 11.097645  5115   east  -93.08445 35.33048
## 25    NA 11.840086  5119   east  -92.29447 34.76654
## 26    NA 11.105003  5131   east  -94.35177 35.29392
## 27    NA 10.887024  5139   east  -92.61279 33.20718
## 28    NA 10.696524  5143   east  -94.19868 36.05138
## 29    NA 10.641453  5145   east  -91.73056 35.25011
## 30    NA  9.821778  9001   east  -73.35482 41.21290
## 31    NA  8.894819  9003   east  -72.71545 41.76997
## 32    NA  7.029871  9005   east  -73.21269 41.75426
## 33    NA 10.001609  9009   east  -72.94032 41.39076
## 34    NA  8.821875  9011   east  -72.09295 41.44611
## 35    NA 10.412954 10001   east  -75.56027 39.09271
## 36    NA 11.551756 10003   east  -75.61171 39.69561
## 37    NA 10.837718 10005   east  -75.34227 38.65142
## 38    NA 11.023823 11001   east  -77.01322 38.91361
## 39    NA  7.503461 12001   east  -82.37995 29.67644
## 40    NA  6.877201 12009   east  -80.68519 28.22008
## 41    NA  6.977536 12011   east  -80.20610 26.13479
## 42    NA  7.375229 12017   east  -82.46916 28.89464
## 43    NA  8.378378 12031   east  -81.63887 30.32087
## 44    NA  8.961842 12033   east  -87.27736 30.50204
## 45    NA  8.341545 12057   east  -82.41208 27.96635
## 46    NA  6.744966 12071   east  -81.84626 26.60254
## 47    NA  9.795967 12073   east  -84.26089 30.46293
## 48    NA  5.000000 12086   east  -80.30226 25.75643
## 49    NA  7.285535 12095   east  -81.40416 28.54713
## 50    NA  6.348540 12099   east  -80.17639 26.63774
## 51    NA  7.738196 12103   east  -82.72777 27.88965
## 52    NA  7.679519 12105   east  -81.75830 27.99244
## 53    NA  7.584909 12111   east  -80.36087 27.34825
## 54    NA  6.800815 12115   east  -82.39888 27.16493
## 55    NA  7.500355 12117   east  -81.31083 28.70991
## 56    NA  7.213773 12127   east  -81.13895 29.04904
## 57    NA 11.959208 13021   east  -83.66521 32.82969
## 58    NA 10.824470 13051   east  -81.10716 32.04756
## 59    NA 11.292279 13059   east  -83.38011 33.95475
## 60    NA 12.913903 13063   east  -84.36560 33.56618
## 61    NA 12.018416 13067   east  -84.55928 33.93554
## 62    NA 12.145414 13089   east  -84.25768 33.79094
## 63    NA 11.746439 13095   east  -84.16321 31.57287
## 64    NA 11.817821 13121   east  -84.41933 33.77960
## 65    NA  9.793812 13127   east  -81.49629 31.18763
## 66    NA 12.145278 13135   east  -84.05845 33.96706
## 67    NA 11.188334 13139   east  -83.84296 34.29680
## 68    NA 11.422325 13153   east  -83.66323 32.55806
## 69    NA  9.709880 13185   east  -83.27718 30.84853
## 70    NA 12.200660 13215   east  -84.94218 32.49006
## 71    NA 10.949059 13223   east  -84.83112 33.90278
## 72    NA 12.392875 13245   east  -82.02840 33.42949
## 73    NA 10.672707 13295   east  -85.29183 34.81869
## 74    NA 11.358464 13303   east  -82.78261 32.94936
## 75    NA 12.121381 13319   east  -83.22660 32.83589
## 76    NA  9.353740 17001   east  -91.26805 39.96545
## 77    NA 11.059005 17019   east  -88.21675 40.13682
## 78    NA 11.878008 17031   east  -87.76782 41.83765
## 79    NA 10.901309 17043   east  -88.06359 41.85603
## 80    NA 11.269838 17065   east  -88.53624 38.09678
## 81    NA 10.336829 17083   east  -90.34214 39.08361
## 82    NA 10.581457 17089   east  -88.33844 41.91823
## 83    NA  9.256495 17097   east  -87.98380 42.31332
## 84    NA 10.275219 17099   east  -88.90429 41.33216
## 85    NA 10.003862 17111   east  -88.35377 42.29640
## 86    NA 10.469588 17113   east  -88.88358 40.49525
## 87    NA 11.631926 17115   east  -88.95323 39.85314
## 88    NA 12.702123 17119   east  -90.01884 38.81251
## 89    NA 11.116070 17143   east  -89.66506 40.75249
## 90    NA 10.104734 17157   east  -89.79688 38.07921
## 91    NA  9.655263 17161   east  -90.50753 41.49143
## 92    NA 12.331060 17163   east  -90.00207 38.53268
## 93    NA 11.039065 17167   east  -89.64650 39.76385
## 94    NA 10.505179 17197   east  -88.02984 41.52743
## 95    NA 10.034725 17201   east  -89.08163 42.31112
## 96    NA 13.093119 18003   east  -85.10007 41.09028
## 97    NA 13.144659 18019   east  -85.72746 38.38648
## 98    NA 11.855302 18035   east  -85.39608 40.21189
## 99    NA 13.023923 18037   east  -86.90055 38.35009
## 100   NA 12.046176 18039   east  -85.89940 41.62609
## 101   NA 12.367141 18043   east  -85.86329 38.31693
## 102   NA 11.630029 18051   east  -87.55836 38.31770
## 103   NA 11.443252 18065   east  -85.40599 39.92353
## 104   NA 11.692452 18067   east  -86.12641 40.47544
## 105   NA 11.811781 18083   east  -87.41894 38.69742
## 106   NA 12.753287 18089   east  -87.39525 41.52960
## 107   NA 10.646946 18091   east  -86.77926 41.59922
## 108   NA 12.089511 18095   east  -85.71920 40.14594
## 109   NA 13.883217 18097   east  -86.14181 39.79279
## 110   NA 12.197354 18127   east  -87.08665 41.52351
## 111   NA 13.435254 18141   east  -86.24193 41.66478
## 112   NA 12.224040 18147   east  -87.01320 38.02437
## 113   NA 12.195055 18157   east  -86.88088 40.40948
## 114   NA 12.828878 18163   east  -87.55696 37.98989
## 115   NA 12.260916 18167   east  -87.39127 39.45689
## 116   NA 10.473696 19013   east  -92.34473 42.48759
## 117   NA 11.699045 19045   east  -90.42740 41.87117
## 118   NA 10.917263 19103   east  -91.56887 41.67491
## 119   NA 11.123664 19111   east  -91.41661 40.57683
## 120   NA 10.407355 19113   east  -91.63056 42.02999
## 121   NA  9.288892 19137   east  -95.15925 41.01310
## 122   NA 12.679211 19139   east  -91.08480 41.47246
## 123   NA  8.828504 19147   east  -94.68872 43.09030
## 124   NA  9.402015 19153   east  -93.63094 41.62794
## 125   NA 10.831346 19155   east  -95.68040 41.30853
## 126   NA 11.960760 19163   east  -90.57886 41.57213
## 127   NA  9.281105 19177   east  -91.94881 40.74022
## 128   NA  8.938962 20091   east  -94.75377 38.93218
## 129   NA  9.337790 20107   east  -94.80532 38.22383
## 130   NA  9.057115 20173   east  -97.36347 37.68133
## 131   NA  9.322162 20177   east  -95.71026 39.03915
## 132   NA  8.681833 20191   east  -97.43256 37.26461
## 133   NA  9.741182 20209   east  -94.69276 39.10371
## 134   NA 11.281826 21013   east  -83.70225 36.69139
## 135   NA 11.405537 21019   east  -82.66458 38.42536
## 136   NA 12.665325 21029   east  -85.67205 38.00518
## 137   NA 11.801188 21037   east  -84.43330 39.03652
## 138   NA  9.685132 21043   east  -83.04932 38.30785
## 139   NA 11.322465 21047   east  -87.48199 36.85444
## 140   NA 12.051753 21059   east  -87.10510 37.75065
## 141   NA 12.395794 21067   east  -84.49464 38.02963
## 142   NA 11.103714 21073   east  -84.86839 38.20365
## 143   NA 11.765029 21093   east  -85.91493 37.73892
## 144   NA 11.859195 21101   east  -87.56808 37.80262
## 145   NA 13.058814 21111   east  -85.70402 38.20924
## 146   NA 11.903645 21117   east  -84.53935 39.02674
## 147   NA 11.670030 21145   east  -88.65111 37.05971
## 148   NA 10.343433 21151   east  -84.27224 37.69968
## 149   NA 11.554662 21183   east  -86.87936 37.45651
## 150   NA 10.090150 21195   east  -82.41111 37.44779
## 151   NA 11.511463 21227   east  -86.42455 36.98691
## 152   NA 10.851174 22017   east  -93.81391 32.52756
## 153   NA  8.874310 22019   east  -93.29176 30.23744
## 154   NA  9.970397 22033   east  -91.11965 30.48012
## 155   NA  9.761818 22047   east  -91.28719 30.27517
## 156   NA  9.079453 22051   east  -90.15140 29.93524
## 157   NA  9.799970 22055   east  -92.03888 30.20739
## 158   NA  9.703961 22073   east  -92.11586 32.50938
## 159   NA  8.898718 22079   east  -92.47659 31.26208
## 160   NA 12.178354 22087   east  -89.89026 29.92209
## 161   NA  9.134486 22105   east  -90.46295 30.60740
## 162   NA  8.528276 22109   east  -90.73592 29.54117
## 163   NA 11.802678 22121   east  -91.26305 30.45209
## 164   NA  7.794286 23001   east  -70.21766 44.13507
## 165   NA  6.162443 23003   east  -68.26166 46.63736
## 166   NA  8.182713 23005   east  -70.32798 43.78448
## 167   NA  4.504539 23009   east  -68.41571 44.51274
## 168   NA  7.334426 23011   east  -69.78455 44.39996
## 169   NA  8.259664 23017   east  -70.70517 44.32452
## 170   NA  7.384147 23019   east  -68.72247 45.17193
## 171   NA  5.186250 23021   east  -69.33369 45.56518
## 172   NA 11.439402 24003   east  -76.57689 39.05786
## 173   NA 11.023192 24005   east  -76.61263 39.37221
## 174   NA 11.665486 24015   east  -75.95146 39.59263
## 175   NA  9.529465 24025   east  -76.29024 39.53467
## 176   NA 10.669508 24031   east  -77.12439 39.07261
## 177   NA 10.926625 24033   east  -76.88218 38.90633
## 178   NA 10.919743 24043   east  -77.77444 39.61470
## 179   NA 11.370975 24510   east  -76.61702 39.30796
## 180   NA  9.194519 25003   east  -73.20989 42.39613
## 181   NA  8.282754 25005   east  -71.10213 41.77833
## 182   NA  7.949325 25009   east  -70.97083 42.63547
## 183   NA  9.244576 25013   east  -72.57131 42.12756
## 184   NA  8.031137 25017   east  -71.27557 42.45909
## 185   NA  8.550007 25023   east  -70.81896 41.97888
## 186   NA  9.788672 25025   east  -71.07349 42.33495
## 187   NA  8.834979 25027   east  -71.84021 42.32964
## 188   NA  8.935905 26005   east  -85.90387 42.59273
## 189   NA  8.404685 26017   east  -83.91741 43.63666
## 190   NA  9.208338 26021   east  -86.42682 41.97124
## 191   NA 14.623500 26033   east  -84.49799 46.32227
## 192   NA  9.140998 26049   east  -83.70050 43.01188
## 193   NA  9.253624 26065   east  -84.46983 42.66775
## 194   NA 10.071916 26077   east  -85.56132 42.25964
## 195   NA 10.229916 26081   east  -85.61270 42.98864
## 196   NA  9.535325 26091   east  -84.06593 41.91763
## 197   NA  9.669258 26099   east  -82.95176 42.59114
## 198   NA  6.745091 26101   east  -86.13540 44.32379
## 199   NA  6.153301 26113   east  -85.14557 44.32747
## 200   NA 10.370516 26115   east  -83.48413 41.91905
## 201   NA  9.090500 26121   east  -86.21617 43.25468
## 202   NA 10.020057 26125   east  -83.30560 42.58823
## 203   NA  9.993553 26139   east  -86.04975 42.92997
## 204   NA  9.888886 26147   east  -82.59785 42.91695
## 205   NA 10.104034 26161   east  -83.75947 42.26033
## 206   NA 10.880669 26163   east  -83.19795 42.33724
## 207   NA  9.568676 27037   east  -93.12013 44.75880
## 208   NA  9.649893 27053   east  -93.36236 44.97133
## 209   NA  9.594459 27109   east  -92.44255 44.00748
## 210   NA 10.205015 27123   east  -93.10659 44.99092
## 211   NA  6.263498 27137   east  -92.40455 47.36915
## 212   NA  9.422119 27139   east  -93.51083 44.69634
## 213   NA  8.559982 27145   east  -94.49589 45.54086
## 214   NA 10.287751 28001   east  -91.35435 31.53554
## 215   NA 10.758182 28011   east  -90.83014 33.80328
## 216   NA 10.645153 28033   east  -89.99207 34.90434
## 217   NA 11.998575 28035   east  -89.27621 31.26302
## 218   NA  9.768253 28043   east  -89.80345 33.77504
## 219   NA  8.993276 28047   east  -89.07954 30.42790
## 220   NA 11.447412 28049   east  -90.30537 32.29433
## 221   NA  9.705737 28059   east  -88.63314 30.45571
## 222   NA 12.157404 28067   east  -89.17378 31.65227
## 223   NA 11.497955 28075   east  -88.67819 32.40008
## 224   NA 11.304253 28081   east  -88.68251 34.27690
## 225   NA 11.375990 29021   east  -94.82701 39.72254
## 226   NA  9.823971 29037   east  -94.37666 38.68937
## 227   NA 10.199264 29039   east  -93.88540 37.73050
## 228   NA  9.609904 29047   east  -94.47403 39.25367
## 229   NA 10.135240 29077   east  -93.31609 37.21630
## 230   NA 11.173592 29095   east  -94.46324 39.02401
## 231   NA 10.935242 29099   east  -90.49927 38.28584
## 232   NA 12.523549 29510   east  -90.24281 38.62772
## 233   NA  9.376956 31055   east  -96.05094 41.26403
## 234   NA  7.367368 31079   east  -98.41773 40.89416
## 235   NA  8.503896 31109   east  -96.68139 40.79670
## 236   NA 10.006475 31153   east  -96.03110 41.14287
## 237   NA  8.759718 31177   east  -96.19888 41.51628
## 238   NA  5.711976 33001   east  -71.43608 43.51520
## 239   NA 10.208299 33005   east  -72.24208 42.91180
## 240   NA  6.834998 33009   east  -71.89463 43.89607
## 241   NA  8.349729 33011   east  -71.58274 42.89558
## 242   NA  8.515929 33013   east  -71.63725 43.27770
## 243   NA  7.598258 33015   east  -71.07762 42.96287
## 244   NA  8.827018 34001   east  -74.62157 39.44772
## 245   NA  9.035905 34003   east  -74.05904 40.93046
## 246   NA  9.876226 34007   east  -75.02926 39.86921
## 247   NA  9.117745 34015   east  -75.13620 39.75278
## 248   NA 10.656258 34017   east  -74.06900 40.74035
## 249   NA  9.024159 34021   east  -74.71165 40.25291
## 250   NA  8.632513 34023   east  -74.38392 40.49404
## 251   NA  8.073528 34027   east  -74.50044 40.86717
## 252   NA  8.677116 34029   east  -74.21090 39.91882
## 253   NA  9.538970 34031   east  -74.21683 40.94649
## 254   NA 10.432357 34039   east  -74.30165 40.65775
## 255   NA  9.551347 34041   east  -75.03070 40.80302
## 256   NA  8.233488 36001   east  -73.84908 42.66210
## 257   NA 11.367129 36005   east  -73.87321 40.84676
## 258   NA  7.771279 36013   east  -79.34572 42.24070
## 259   NA 10.034822 36029   east  -78.79630 42.85991
## 260   NA  4.625279 36031   east  -73.73759 44.16219
## 261   NA 10.884358 36047   east  -73.95225 40.64981
## 262   NA  8.251590 36055   east  -77.63032 43.15996
## 263   NA  9.704766 36059   east  -73.60254 40.71595
## 264   NA 11.784900 36061   east  -73.97353 40.76966
## 265   NA  8.904244 36063   east  -78.84810 43.14609
## 266   NA  7.926522 36067   east  -76.17912 43.04843
## 267   NA  8.591969 36071   east  -74.27948 41.40612
## 268   NA 10.114641 36081   east  -73.82038 40.71429
## 269   NA  9.936354 36085   east  -74.14092 40.58220
## 270   NA  7.500272 36101   east  -77.36315 42.28362
## 271   NA  8.924974 36103   east  -73.02623 40.83331
## 272   NA  9.611338 36119   east  -73.79401 41.05717
## 273   NA 11.272163 37001   east  -79.41392 36.07604
## 274   NA  9.319477 37021   east  -82.53146 35.59526
## 275   NA 10.383642 37033   east  -79.32001 36.40641
## 276   NA 11.806276 37035   east  -81.24757 35.69561
## 277   NA  9.707310 37037   east  -79.28286 35.70231
## 278   NA 11.376835 37051   east  -78.90725 35.05621
## 279   NA 12.433921 37057   east  -80.19781 35.81681
## 280   NA  9.750804 37061   east  -77.98801 34.92412
## 281   NA 10.954336 37063   east  -78.90036 35.99834
## 282   NA 10.062751 37065   east  -77.64031 35.91624
## 283   NA 10.971491 37067   east  -80.24460 36.11420
## 284   NA 11.308123 37071   east  -81.16645 35.28455
## 285   NA 11.022960 37081   east  -79.84698 36.05264
## 286   NA 11.045066 37087   east  -82.96053 35.52338
## 287   NA  9.445623 37099   east  -83.18386 35.31347
## 288   NA  9.809391 37107   east  -77.63059 35.25056
## 289   NA 11.033742 37111   east  -82.03840 35.68694
## 290   NA  9.444095 37117   east  -77.12383 35.83833
## 291   NA 11.315136 37119   east  -80.83086 35.23179
## 292   NA 10.194477 37121   east  -82.14284 35.97769
## 293   NA 10.506809 37123   east  -79.90331 35.34010
## 294   NA  9.971608 37147   east  -77.39108 35.58991
## 295   NA 10.710838 37155   east  -79.10653 34.65837
## 296   NA 11.681061 37159   east  -80.51550 35.61734
## 297   NA 10.150782 37173   east  -83.44653 35.42182
## 298   NA 10.747831 37183   east  -78.66575 35.79651
## 299   NA  9.434592 37189   east  -81.70011 36.21793
## 300   NA 11.059944 37191   east  -77.99534 35.36471
## 301   NA  8.598159 38017   east  -97.08444 46.89958
## 302   NA  9.719606 39009   east  -82.07508 39.35691
## 303   NA 13.424673 39017   east  -84.51319 39.42786
## 304   NA 12.776898 39023   east  -83.82636 39.92131
## 305   NA 11.575359 39025   east  -84.17742 39.05957
## 306   NA 12.593666 39035   east  -81.66385 41.46030
## 307   NA 11.909242 39049   east  -82.99653 39.99051
## 308   NA 11.875767 39057   east  -83.94940 39.71256
## 309   NA 13.770987 39061   east  -84.49379 39.17953
## 310   NA 12.807829 39081   east  -80.72585 40.36070
## 311   NA 12.228452 39087   east  -82.54751 38.53040
## 312   NA 10.415077 39093   east  -82.13582 41.37027
## 313   NA 11.556389 39095   east  -83.60847 41.64927
## 314   NA 12.434767 39099   east  -80.70650 41.05256
## 315   NA 11.149050 39103   east  -81.87205 41.11258
## 316   NA 13.217910 39113   east  -84.21859 39.75026
## 317   NA 11.494418 39133   east  -81.25372 41.17219
## 318   NA 11.685097 39135   east  -84.63719 39.74386
## 319   NA 11.594706 39145   east  -82.94307 38.78893
## 320   NA 13.224641 39151   east  -81.38360 40.81941
## 321   NA 13.082987 39153   east  -81.52064 41.09560
## 322   NA 11.864689 39165   east  -84.22176 39.44015
## 323   NA 12.232856 40001   east  -94.64069 35.89642
## 324   NA  8.950644 40015   east  -98.34474 35.13654
## 325   NA 10.716408 40097   east  -95.21124 36.31482
## 326   NA 10.549134 40101   east  -95.36713 35.68659
## 327   NA  9.630629 40109   east  -97.49653 35.50967
## 328   NA  9.467113 40115   east  -94.83983 36.86137
## 329   NA 10.217874 40121   east  -95.71430 34.93971
## 330   NA 10.245121 40135   east  -94.78890 35.48435
## 331   NA 11.147216 40143   east  -95.93900 36.12476
## 332   NA 11.355195 42001   east  -77.18254 39.85438
## 333   NA 13.136062 42003   east  -79.95970 40.44057
## 334   NA 12.897254 42007   east  -80.31267 40.69541
## 335   NA 11.132102 42011   east  -75.91424 40.37774
## 336   NA 11.342894 42017   east  -75.02784 40.25494
## 337   NA 13.195630 42021   east  -78.78923 40.43547
## 338   NA  9.923461 42027   east  -77.83446 40.88728
## 339   NA 13.654429 42029   east  -75.68330 39.98881
## 340   NA 13.919460 42041   east  -77.13368 40.20120
## 341   NA 12.420925 42043   east  -76.80466 40.32449
## 342   NA 13.308193 42045   east  -75.35064 39.91298
## 343   NA 10.511073 42049   east  -80.05698 42.05203
## 344   NA  9.677219 42069   east  -75.63221 41.43348
## 345   NA 12.702044 42071   east  -76.28122 40.06767
## 346   NA 11.020534 42085   east  -80.31313 41.29443
## 347   NA 10.774001 42091   east  -75.31743 40.17060
## 348   NA 12.709182 42095   east  -75.32207 40.70286
## 349   NA 11.500995 42101   east  -75.14479 39.99801
## 350   NA 11.838172 42125   east  -80.13855 40.18261
## 351   NA 13.281318 42129   east  -79.56187 40.32611
## 352   NA 12.459742 42133   east  -76.74218 39.92730
## 353   NA  6.615918 44003   east  -71.47953 41.70412
## 354   NA  8.690118 44007   east  -71.45616 41.85285
## 355   NA  9.201633 45019   east  -79.99138 32.80964
## 356   NA 10.297114 45025   east  -80.15629 34.65417
## 357   NA 10.740024 45037   east  -81.90446 33.77124
## 358   NA 10.563050 45041   east  -79.74030 34.06507
## 359   NA 11.558594 45045   east  -82.35311 34.86127
## 360   NA 11.422803 45063   east  -81.23232 33.94922
## 361   NA  8.600291 45073   east  -83.02142 34.71504
## 362   NA 11.203930 45079   east  -80.97433 34.03540
## 363   NA  8.319356 46011   east  -96.79849 44.34542
## 364   NA  8.136039 46013   east  -98.38199 45.53867
## 365   NA  8.631982 46029   east  -97.15738 44.93933
## 366   NA  9.101323 46099   east  -96.74220 43.60144
## 367   NA 11.782392 47009   east  -83.96924 35.75428
## 368   NA 10.815339 47037   east  -86.76576 36.15660
## 369   NA 10.104499 47045   east  -89.35898 36.06046
## 370   NA 11.480995 47065   east  -85.23010 35.09552
## 371   NA 12.155269 47093   east  -83.95593 35.97311
## 372   NA  9.043446 47099   east  -87.37581 35.20526
## 373   NA 12.445652 47105   east  -84.29551 35.74490
## 374   NA 11.721120 47107   east  -84.59591 35.41638
## 375   NA  9.721455 47113   east  -88.82738 35.62865
## 376   NA  9.755862 47119   east  -87.06484 35.61640
## 377   NA 10.607539 47125   east  -87.37043 36.51062
## 378   NA 10.459144 47141   east  -85.49528 36.15783
## 379   NA 12.035383 47145   east  -84.52673 35.87976
## 380   NA 10.658054 47157   east  -89.92690 35.13338
## 381   NA 11.053976 47163   east  -82.38360 36.53036
## 382   NA 10.431050 47165   east  -86.48751 36.43644
## 383   NA  8.911126 48029   east  -98.50894 29.45323
## 384   NA 11.277835 48037   east  -94.25020 33.43807
## 385   NA 10.711332 48061   east  -97.55744 26.08299
## 386   NA 10.102954 48113   east  -96.78883 32.80358
## 387   NA  9.747652 48139   east  -96.79844 32.39143
## 388   NA 11.727956 48201   east  -95.36306 29.78778
## 389   NA  9.873713 48203   east  -94.35308 32.54868
## 390   NA 10.204761 48215   east  -98.15980 26.24229
## 391   NA  9.866045 48245   east  -94.06177 29.98423
## 392   NA  9.722218 48355   east  -97.45894 27.74715
## 393   NA  9.835390 48361   east  -93.84206 30.11460
## 394   NA 10.304592 48439   east  -97.26583 32.76383
## 395   NA  9.214648 48453   east  -97.76981 30.32181
## 396   NA  6.902222 50003   east  -73.10796 43.02942
## 397   NA  6.296837 50007   east  -73.12727 44.47568
## 398   NA  9.782239 50021   east  -73.05306 43.59508
## 399   NA  9.386446 51003   east  -78.53528 38.02418
## 400   NA 10.881701 51013   east  -77.10826 38.88034
## 401   NA  9.487579 51036   east  -77.07161 37.34487
## 402   NA 10.236205 51041   east  -77.54590 37.41994
## 403   NA 10.353459 51059   east  -77.24957 38.84085
## 404   NA 11.075066 51069   east  -78.24462 39.19275
## 405   NA  9.767818 51087   east  -77.47699 37.58475
## 406   NA 10.190762 51107   east  -77.56882 39.06463
## 407   NA  9.844875 51139   east  -78.48847 38.61346
## 408   NA 10.801909 51165   east  -78.84438 38.46651
## 409   NA 10.221530 51520   east  -82.17619 36.61112
## 410   NA  9.327619 51680   east  -79.17020 37.40367
## 411   NA 10.257793 51710   east  -76.25990 36.88575
## 412   NA 10.314598 51770   east  -79.95571 37.27690
## 413   NA  9.928187 51775   east  -80.05584 37.28690
## 414   NA  9.802806 51810   east  -76.08718 36.83450
## 415   NA 12.876411 54003   east  -77.98393 39.46299
## 416   NA 13.339803 54009   east  -80.59096 40.29642
## 417   NA 13.043723 54011   east  -82.32324 38.41292
## 418   NA 12.405182 54029   east  -80.57548 40.50472
## 419   NA 11.783240 54033   east  -80.34195 39.29704
## 420   NA 12.420487 54039   east  -81.60569 38.34893
## 421   NA 12.843694 54049   east  -80.19794 39.49985
## 422   NA 13.131861 54051   east  -80.71289 39.90745
## 423   NA 11.556837 54061   east  -79.99478 39.63479
## 424   NA 12.423321 54069   east  -80.67612 40.07877
## 425   NA 10.114421 54081   east  -81.22827 37.77134
## 426   NA 12.840365 54107   east  -81.53170 39.25943
## 427   NA  5.715239 55003   east  -90.72189 46.33752
## 428   NA 11.366288 55009   east  -88.03277 44.49349
## 429   NA 11.098744 55025   east  -89.39367 43.07150
## 430   NA  6.308198 55041   east  -88.77048 45.60437
## 431   NA 11.219191 55043   east  -90.66734 42.87649
## 432   NA 10.440923 55059   east  -87.95243 42.56867
## 433   NA 10.058572 55063   east  -91.18913 43.85941
## 434   NA  9.417947 55071   east  -87.75858 44.10232
## 435   NA 11.575934 55079   east  -87.96071 43.03411
## 436   NA 10.248115 55087   east  -88.41002 44.33653
## 437   NA 10.030949 55089   east  -87.94472 43.32925
## 438   NA  9.749491 55109   east  -92.50697 45.02847
## 439   NA  9.208060 55111   east  -89.90087 43.44451
## 440   NA  8.483095 55119   east  -90.49127 45.19183
## 441   NA  6.111594 55125   east  -89.53184 46.02458
## 442   NA 11.975980 55133   east  -88.25589 43.03635
## 443   NA  6.058860  2020   west -149.76210 61.19190
## 444   NA 11.101467  2090   west -147.56838 64.81859
## 445   NA  7.308113  2110   west -134.51158 58.35142
## 446   NA  7.147626  2170   west -149.48109 61.76274
## 447   NA  6.929844  4003   west -109.90432 31.75027
## 448   NA  6.132351  4005   west -111.51106 35.77144
## 449   NA  8.228339  4013   west -112.08791 33.49451
## 450   NA  5.328475  4019   west -111.08862 32.17841
## 451   NA 10.502862  4021   west -111.49811 32.96567
## 452   NA 11.326499  4023   west -110.90573 31.48190
## 453   NA  5.165413  4025   west -112.41471 34.65033
## 454   NA  9.465902  6001   west -122.09615 37.71666
## 455   NA  9.842297  6007   west -121.64631 39.64835
## 456   NA  6.768743  6009   west -120.55724 38.16757
## 457   NA  6.944759  6011   west -122.15239 39.18749
## 458   NA  8.545173  6013   west -122.06089 37.94242
## 459   NA 16.194524  6019   west -119.90347 36.63837
## 460   NA  7.196966  6023   west -123.98826 40.74740
## 461   NA  9.417947  6025   west -115.48759 32.96310
## 462   NA  6.431331  6027   west -117.69333 36.70904
## 463   NA 15.803780  6029   west -118.68334 35.29602
## 464   NA 18.440731  6031   west -119.81127 36.15514
## 465   NA  3.494351  6033   west -122.75179 39.02326
## 466   NA 16.661801  6037   west -118.23422 34.08851
## 467   NA  8.151318  6045   west -123.43022 39.36310
## 468   NA 15.015733  6047   west -120.67412 37.24578
## 469   NA  6.455655  6053   west -121.52895 36.45073
## 470   NA  5.280999  6057   west -120.81427 39.27078
## 471   NA 14.506272  6059   west -117.86166 33.73152
## 472   NA  7.969183  6061   west -120.94762 38.97158
## 473   NA  9.678393  6063   west -120.85120 40.02015
## 474   NA 17.429049  6065   west -116.80361 33.78331
## 475   NA  9.679342  6067   west -121.38927 38.56111
## 476   NA  5.995774  6069   west -121.28884 36.74575
## 477   NA  9.843627  6073   west -117.06299 32.90462
## 478   NA 10.741087  6075   west -122.43739 37.75988
## 479   NA 11.560022  6077   west -121.28131 37.94596
## 480   NA  7.308466  6079   west -120.58746 35.37391
## 481   NA  8.999636  6081   west -122.33193 37.53104
## 482   NA  8.415964  6083   west -120.07143 34.64938
## 483   NA  9.841446  6085   west -121.91282 37.30651
## 484   NA  6.196792  6087   west -121.97376 37.00235
## 485   NA  5.560546  6089   west -122.12120 40.68870
## 486   NA  5.702414  6093   west -122.50282 41.61494
## 487   NA 10.099460  6095   west -122.05183 38.23593
## 488   NA  6.925403  6097   west -122.76862 38.44009
## 489   NA 16.251904  6099   west -120.95880 37.61380
## 490   NA  7.309657  6101   west -121.65816 39.07619
## 491   NA 16.183583  6107   west -119.16609 36.23465
## 492   NA  9.617136  6111   west -119.04533 34.27705
## 493   NA  7.705003  6113   west -121.77571 38.63415
## 494   NA  8.469519  8001   west -104.86995 39.86672
## 495   NA  6.654664  8005   west -104.84145 39.64190
## 496   NA  6.771765  8013   west -105.19515 40.05809
## 497   NA  7.733877  8031   west -104.96549 39.72629
## 498   NA  5.732005  8035   west -104.88872 39.46117
## 499   NA  4.186090  8039   west -104.29216 39.29364
## 500   NA  5.971900  8041   west -104.74883 38.86544
## 501   NA  6.614878  8069   west -105.21301 40.53410
## 502   NA  9.328776  8077   west -108.50910 39.09523
## 503   NA  6.275940  8083   west -108.57584 37.35369
## 504   NA  7.680745  8123   west -104.72821 40.34668
## 505   NA  6.824419 15001   west -155.39316 19.68389
## 506   NA  5.783519 15003   west -158.03586 21.42613
## 507   NA  5.693945 15009   west -156.61584 20.90103
## 508   NA 10.329405 16001   west -116.28007 43.56030
## 509   NA  4.917140 16005   west -112.34627 42.77779
## 510   NA 10.792448 16009   west -116.63764 47.24004
## 511   NA  8.846641 16027   west -116.67363 43.61814
## 512   NA 10.376327 16041   west -111.85559 42.15321
## 513   NA 12.437898 16059   west -113.77824 44.94882
## 514   NA 13.346910 16079   west -116.04049 47.41788
## 515   NA  8.185202 30029   west -114.32128 48.26788
## 516   NA  6.875453 30031   west -111.18679 45.66387
## 517   NA  8.030033 30049   west -112.15393 46.81387
## 518   NA 11.464039 30053   west -115.35935 48.55862
## 519   NA  8.294273 30063   west -113.96861 46.93997
## 520   NA  7.187243 30081   west -114.11724 46.21017
## 521   NA  5.135916 30083   west -104.45825 47.77934
## 522   NA  6.170566 30089   west -115.08446 47.63311
## 523   NA  8.858588 30093   west -112.56254 45.96093
## 524   NA  6.672426 32003   west -115.11790 36.13922
## 525   NA  7.203337 32031   west -119.76363 39.70417
## 526   NA  5.425800 35001   west -106.61247 35.10105
## 527   NA  5.689097 35005   west -104.44499 33.35589
## 528   NA  7.838288 35013   west -106.78167 32.27412
## 529   NA  4.793644 35017   west -108.24675 32.70857
## 530   NA  5.678738 35025   west -103.32776 32.72844
## 531   NA  4.601408 35045   west -108.30406 36.57760
## 532   NA  4.195688 35049   west -106.00713 35.58883
## 533   NA  4.460193 38007   west -103.35967 46.99011
## 534   NA  7.057676 38015   west -100.64104 46.89068
## 535   NA  6.121487 38057   west -101.68637 47.28470
## 536   NA  4.978397 41017   west -121.33720 43.99606
## 537   NA 11.422585 41025   west -119.03972 43.46143
## 538   NA  6.544022 41029   west -122.78566 42.38983
## 539   NA  7.857464 41033   west -123.43236 42.36216
## 540   NA 11.201131 41035   west -121.65163 42.66203
## 541   NA  9.606104 41037   west -120.57698 42.71338
## 542   NA  7.981166 41039   west -123.08081 43.98925
## 543   NA  9.559339 41043   west -122.73143 44.54483
## 544   NA  7.138978 41051   west -122.61323 45.52247
## 545   NA  7.512537 41059   west -118.81731 45.63730
## 546   NA  6.722891 41061   west -118.03131 45.32501
## 547   NA  8.617499 41067   west -122.91514 45.49680
## 548   NA  4.324736 46033   west -103.51523 43.69311
## 549   NA  4.175901 46071   west -101.66929 43.76851
## 550   NA  6.105285 46103   west -103.10550 44.04519
## 551   NA  5.349604 48043   west -103.40222 29.99413
## 552   NA  7.412242 48135   west -102.43482 31.86721
## 553   NA 10.305085 48141   west -106.35181 31.77016
## 554   NA  7.740238 48303   west -101.85207 33.57008
## 555   NA  5.998857 48375   west -101.84221 35.25705
## 556   NA  7.974543 49003   west -112.48032 41.59579
## 557   NA  9.757741 49005   west -111.84212 41.75108
## 558   NA  9.631552 49011   west -111.94075 41.00435
## 559   NA  9.632232 49035   west -111.91557 40.65580
## 560   NA  6.624555 49045   west -112.69003 40.41542
## 561   NA  8.866858 49049   west -111.73671 40.23660
## 562   NA  9.195895 49057   west -111.97287 41.23218
## 563   NA  7.930779 53011   west -122.54783 45.69219
## 564   NA  7.939444 53033   west -122.19685 47.54874
## 565   NA  8.843770 53053   west -122.38357 47.17157
## 566   NA  7.821492 53061   west -122.09418 47.97684
## 567   NA  9.462253 53077   west -120.51294 46.48245
## 568   NA  5.272252 56005   west -105.49616 44.08580
## 569   NA  3.382626 56009   west -105.50272 42.85738
## 570   NA  8.505596 56013   west -108.68112 43.09821
## 571   NA  4.132739 56021   west -104.76208 41.19620
## 572   NA  4.955570 56029   west -108.99902 44.56896
## 573   NA  6.549239 56033   west -106.96813 44.78363
## 574   NA  5.632587 56035   west -109.99115 42.75095
## 575   NA  6.349710 56037   west -109.16820 41.64071
## 576   NA  4.565808 56039   west -110.67353 43.61212
boxplot(pollution$pm25, col = "blue")
abline(h=12)

hist(pollution$pm25, col="green")
abline(v = 12, lwd =2)
abline(v = median(pollution$pm25), col = "magenta", lwd = 4)
abline(v = mean(pollution$pm25), col = "blue", lwd = 4)

hist(pollution$pm25, col="green", breaks = 100)
rug(pollution$pm25)

barplot(table(pollution$region), col = "wheat", main = "Number of Counties in Each Region")

boxplot(pm25 ~ region, data=pollution, col = "red")

par(mfrow = c(2,1), mar = c(4,4,2,1))
hist(subset(pollution, region=="east")$pm25, col="green")
hist(subset(pollution, region=="west")$pm25, col="green")

Scatterplot

with(pollution, plot(latitude, pm25))
abline(h=12, lwd=4,lty=5)

par(mfrow = c(1,2), mar = c(5,4,2,1))

with(subset(pollution, region =="west"),
            
            plot(latitude, pm25, main = "Oeste"))


with(subset(pollution, region =="east"),
            
            plot(latitude, pm25, main = "Este"))

library(datasets)
data("cars")
with(cars, plot(speed, dist))

library(lattice)
state <- data.frame(state.x77, region=state.region)
xyplot(Life.Exp~Income | region, data = state, layout=c(4,1))

ggplot2

library(ggplot2)
data(mpg)
qplot(displ, hwy, data=mpg)

library(datasets)
hist(airquality$Ozone)

library(datasets)
with(airquality, plot(Wind, Ozone))

library(datasets)
airquality <- transform(airquality, Month=factor(Month))
boxplot(Ozone~Month, airquality, xlab="Month", ylab="Ozone (ppb)")

par("lty")
## [1] "solid"
par("col")
## [1] "black"
par("pch")
## [1] 1
par("bg")
## [1] "white"
par("mar")
## [1] 5.1 4.1 4.1 2.1
par("mfrow")
## [1] 1 1
library(datasets)
with(airquality, plot(Wind, Ozone))
title(main = "Ozone and Wind in New York City")

library(datasets)
with(airquality, plot(Wind, Ozone, main="Ozone and Wind in New York City"))

with(subset(airquality, Month == 5), points(Wind, Ozone, col = "blue"))

library(datasets)
with(airquality, plot(Wind, Ozone, main="Ozono y viento en la ciudad de Nueva York", type="n"))

with(subset(airquality, Month == 5), points(Wind, Ozone, col = "blue"))
with(subset(airquality, Month != 5), points(Wind, Ozone, col = "red"))

legend("topright", pch=1, col=c("blue","red"), legend=c("May", "Other Months"))

Regresión Lineal

library(datasets)
with(airquality, plot(Wind, Ozone, main="Ozono y viento en la ciudad de Nueva York", pch=20))

model <- lm(Ozone ~ Wind, airquality)
abline(model, lwd=2)

summary(model)
## 
## Call:
## lm(formula = Ozone ~ Wind, data = airquality)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -51.572 -18.854  -4.868  15.234  90.000 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  96.8729     7.2387   13.38  < 2e-16 ***
## Wind         -5.5509     0.6904   -8.04 9.27e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 26.47 on 114 degrees of freedom
##   (37 observations deleted due to missingness)
## Multiple R-squared:  0.3619, Adjusted R-squared:  0.3563 
## F-statistic: 64.64 on 1 and 114 DF,  p-value: 9.272e-13
par(mfrow=c(1,2))
with(airquality, {
  plot(Wind, Ozone, main="Ozone and Wind")
  plot(Solar.R, Ozone, main="Ozone and Solar Radiation")
})

par(mfrow=c(1,3), mar=c(4,4,2,1), oma=c(0,0,2,0))
with(airquality,{
  plot(Wind,Ozone,main="Ozone and Wind")
  plot(Solar.R, Ozone, main = "Ozone and Solar Radiation")
  plot(Temp, Ozone, main = "Ozone and Temperature")
  mtext("Ozone and Weather in New York City", outer=TRUE)
})

Base Plotting Demonstration

Histograma aleatorio

x <- rnorm(100)
hist(x)

boxplot(x)

summary(x)
##     Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
## -2.22742 -0.63403 -0.01828  0.01784  0.61539  2.09843
y <- rnorm(100)

plot(x,y)

z <- rnorm(100)
plot(x, z)

example(points)
## 
## points> require(stats) # for rnorm
## 
## points> plot(-4:4, -4:4, type = "n")  # setting up coord. system

## 
## points> points(rnorm(200), rnorm(200), col = "red")
## 
## points> points(rnorm(100)/2, rnorm(100)/2, col = "blue", cex = 1.5)
## 
## points> op <- par(bg = "light blue")
## 
## points> x <- seq(0, 2*pi, len = 51)
## 
## points> ## something "between type='b' and type='o'":
## points> plot(x, sin(x), type = "o", pch = 21, bg = par("bg"), col = "blue", cex = .6,
## points+  main = 'plot(..., type="o", pch=21, bg=par("bg"))')

## 
## points> par(op)
## 
## points> ## Not run: 
## points> ##D ## The figure was produced by calls like
## points> ##D png("pch.png", height = 0.7, width = 7, res = 100, units = "in")
## points> ##D par(mar = rep(0,4))
## points> ##D plot(c(-1, 26), 0:1, type = "n", axes = FALSE)
## points> ##D text(0:25, 0.6, 0:25, cex = 0.5)
## points> ##D points(0:25, rep(0.3, 26), pch = 0:25, bg = "grey")
## points> ## End(Not run)
## points> 
## points> ##-------- Showing all the extra & some char graphics symbols ---------
## points> pchShow <-
## points+   function(extras = c("*",".", "o","O","0","+","-","|","%","#"),
## points+            cex = 3, ## good for both .Device=="postscript" and "x11"
## points+            col = "red3", bg = "gold", coltext = "brown", cextext = 1.2,
## points+            main = paste("plot symbols :  points (...  pch = *, cex =",
## points+                         cex,")"))
## points+   {
## points+     nex <- length(extras)
## points+     np  <- 26 + nex
## points+     ipch <- 0:(np-1)
## points+     k <- floor(sqrt(np))
## points+     dd <- c(-1,1)/2
## points+     rx <- dd + range(ix <- ipch %/% k)
## points+     ry <- dd + range(iy <- 3 + (k-1)- ipch %% k)
## points+     pch <- as.list(ipch) # list with integers & strings
## points+     if(nex > 0) pch[26+ 1:nex] <- as.list(extras)
## points+     plot(rx, ry, type = "n", axes  =  FALSE, xlab = "", ylab = "", main = main)
## points+     abline(v = ix, h = iy, col = "lightgray", lty = "dotted")
## points+     for(i in 1:np) {
## points+       pc <- pch[[i]]
## points+       ## 'col' symbols with a 'bg'-colored interior (where available) :
## points+       points(ix[i], iy[i], pch = pc, col = col, bg = bg, cex = cex)
## points+       if(cextext > 0)
## points+           text(ix[i] - 0.3, iy[i], pc, col = coltext, cex = cextext)
## points+     }
## points+   }
## 
## points> pchShow()

## 
## points> pchShow(c("o","O","0"), cex = 2.5)

## 
## points> pchShow(NULL, cex = 4, cextext = 0, main = NULL)

## 
## points> ## No test: 
## points> ##D ## ------------ test code for various pch specifications -------------
## points> ##D # Try this in various font families (including Hershey)
## points> ##D # and locales.  Use sign = -1 asserts we want Latin-1.
## points> ##D # Standard cases in a MBCS locale will not plot the top half.
## points> ##D TestChars <- function(sign = 1, font = 1, ...)
## points> ##D {
## points> ##D    MB <- l10n_info()$MBCS
## points> ##D    r <- if(font == 5) { sign <- 1; c(32:126, 160:254)
## points> ##D        } else if(MB) 32:126 else 32:255
## points> ##D    if (sign == -1) r <- c(32:126, 160:255)
## points> ##D    par(pty = "s")
## points> ##D    plot(c(-1,16), c(-1,16), type = "n", xlab = "", ylab = "",
## points> ##D         xaxs = "i", yaxs = "i",
## points> ##D         main = sprintf("sign = %d, font = %d", sign, font))
## points> ##D    grid(17, 17, lty = 1) ; mtext(paste("MBCS:", MB))
## points> ##D    for(i in r) try(points(i%%16, i%/%16, pch = sign*i, font = font,...))
## points> ##D }
## points> ##D TestChars()
## points> ##D try(TestChars(sign = -1))
## points> ##D TestChars(font = 5)  # Euro might be at 160 (0+10*16).
## points> ##D                      # macOS has apple at 240 (0+15*16).
## points> ##D try(TestChars(-1, font = 2))  # bold
## points> ## End(No test)
## points> 
## points>