This is to show the different ways to check your data and also go through a couple more subsetting examples

# load packages
library(tidyselect)

# set working directory
setwd("C:/Users/siobh/OneDrive - The University Of British Columbia/TA/Biol 403 - Parfrey 2020/403 lab scripts")


# read in data frame
seagrass = readRDS("Rank6_top_split_taxa.rds")

this is the dataframe that is called “top_split_taxa” in lab 5 and I grouped it at Rank6

summary(seagrass)
##    SampleID         BarcodeSequence    LinkerPrimerSequence   swab_id         
##  Length:2384        Length:2384        Mode:logical         Length:2384       
##  Class :character   Class :character   NA's:2384            Class :character  
##  Mode  :character   Mode  :character                        Mode  :character  
##                                                                               
##                                                                               
##                                                                               
##                                                                               
##  barcode_plate      barcode_well       project_name       person_responsible
##  Length:2384        Length:2384        Length:2384        Length:2384       
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##      date              Run_id          sequence_file         region         
##  Length:2384        Length:2384        Length:2384        Length:2384       
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##      site              latitude       longitude      sample_method     
##  Length:2384        Min.   :51.67   Min.   :-128.5   Length:2384       
##  Class :character   1st Qu.:51.68   1st Qu.:-128.4   Class :character  
##  Mode  :character   Median :51.81   Median :-128.2   Mode  :character  
##                     Mean   :51.84   Mean   :-128.3                     
##                     3rd Qu.:51.93   3rd Qu.:-128.1                     
##                     Max.   :52.06   Max.   :-128.1                     
##                                                                        
##  host_species        host_type         survey_type        dna_extraction_plate
##  Length:2384        Length:2384        Length:2384        Length:2384         
##  Class :character   Class :character   Class :character   Class :character    
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character    
##                                                                               
##                                                                               
##                                                                               
##                                                                               
##  dna_extraction_well    growth          sample_growth       quadrat_id       
##  Length:2384         Length:2384        Length:2384        Length:2384       
##  Class :character    Class :character   Class :character   Class :character  
##  Mode  :character    Mode  :character   Mode  :character   Mode  :character  
##                                                                              
##                                                                              
##                                                                              
##                                                                              
##  meso_quadrat_density_shoots_m2 meso_quadrat_biomass_g_m2 meso_quadrats_epimass
##  Min.   : 8.00                  Min.   : 2.50             Min.   : 0.100       
##  1st Qu.:20.00                  1st Qu.:15.80             1st Qu.: 0.700       
##  Median :24.00                  Median :22.60             Median : 1.500       
##  Mean   :31.65                  Mean   :22.11             Mean   : 3.878       
##  3rd Qu.:40.00                  3rd Qu.:28.20             3rd Qu.:10.500       
##  Max.   :80.00                  Max.   :42.40             Max.   :10.500       
##  NA's   :736                    NA's   :736               NA's   :1520         
##  meso_shoot_id      meso_shoot_length meso_shoot_width meso_number_of_blades
##  Length:2384        Min.   : 41.9     Min.   :0.3000   Min.   :4.000        
##  Class :character   1st Qu.: 89.6     1st Qu.:0.4000   1st Qu.:5.000        
##  Mode  :character   Median :107.2     Median :0.4000   Median :5.000        
##                     Mean   :107.2     Mean   :0.5516   Mean   :5.505        
##                     3rd Qu.:124.0     3rd Qu.:0.6000   3rd Qu.:6.000        
##                     Max.   :198.5     Max.   :1.1000   Max.   :8.000        
##                     NA's   :896       NA's   :896      NA's   :896          
##  meso_total_dry_wt meso_blade_dry_wt meso_microepiphyte_wt
##  Min.   :0.500     Min.   :0.100     Min.   : 3.60        
##  1st Qu.:1.200     1st Qu.:0.200     1st Qu.:17.70        
##  Median :1.700     Median :0.300     Median :23.20        
##  Mean   :1.798     Mean   :0.329     Mean   :25.64        
##  3rd Qu.:2.200     3rd Qu.:0.400     3rd Qu.:35.40        
##  Max.   :5.200     Max.   :1.000     Max.   :52.70        
##  NA's   :896       NA's   :896       NA's   :896          
##  meso_macroepiphyte_wt meso_smithora_wt meso_bryozoans        broken         
##  Min.   :0.0000        Min.   :0.000    Length:2384        Length:2384       
##  1st Qu.:0.0000        1st Qu.:0.000    Class :character   Class :character  
##  Median :0.0000        Median :0.000    Mode  :character   Mode  :character  
##  Mean   :0.2204        Mean   :0.228                                         
##  3rd Qu.:0.0000        3rd Qu.:0.000                                         
##  Max.   :4.2000        Max.   :4.100                                         
##  NA's   :896           NA's   :896                                           
##  smithora_shoot     smithora_region      ocean_id           mean_depth    
##  Length:2384        Length:2384        Length:2384        Min.   : 0.500  
##  Class :character   Class :character   Class :character   1st Qu.: 4.000  
##  Mode  :character   Mode  :character   Mode  :character   Median : 5.000  
##                                                           Mean   : 5.117  
##                                                           3rd Qu.: 6.000  
##                                                           Max.   :10.000  
##                                                                           
##    ctd_depth     conductivity    temperature       pressure    
##  Min.   :0.50   Min.   :33.85   Min.   :10.28   Min.   :10.64  
##  1st Qu.:4.08   1st Qu.:35.09   1st Qu.:11.54   1st Qu.:14.25  
##  Median :4.37   Median :37.64   Median :13.93   Median :14.54  
##  Mean   :4.33   Mean   :37.14   Mean   :13.42   Mean   :14.39  
##  3rd Qu.:4.61   3rd Qu.:38.84   3rd Qu.:14.81   3rd Qu.:14.78  
##  Max.   :9.49   Max.   :39.60   Max.   :15.93   Max.   :19.71  
##                                                                
##       par         fluorometry_chlorophyll   turbidity      dissolved_oxygen
##  Min.   : 80.88   Min.   :-0.31           Min.   : 0.600   Min.   : 79.82  
##  1st Qu.:141.25   1st Qu.: 0.92           1st Qu.: 0.790   1st Qu.:108.37  
##  Median :246.00   Median : 1.26           Median : 1.260   Median :124.20  
##  Mean   :243.00   Mean   : 1.27           Mean   : 3.509   Mean   :117.57  
##  3rd Qu.:259.93   3rd Qu.: 1.53           3rd Qu.: 1.470   3rd Qu.:132.04  
##  Max.   :608.45   Max.   : 2.37           Max.   :16.820   Max.   :137.81  
##                                                                            
##     salinity     speed_of_sound dissolved_oxygen_concentration     chl_a       
##  Min.   :29.80   Min.   :1486   Min.   : 7.37                  Min.   :0.0074  
##  1st Qu.:30.82   1st Qu.:1490   1st Qu.: 9.50                  1st Qu.:0.0253  
##  Median :31.13   Median :1499   Median :10.34                  Median :0.0788  
##  Mean   :31.01   Mean   :1497   Mean   :10.07                  Mean   :0.0668  
##  3rd Qu.:31.25   3rd Qu.:1502   3rd Qu.:11.13                  3rd Qu.:0.0825  
##  Max.   :31.46   Max.   :1505   Max.   :11.24                  Max.   :0.2918  
##                                                                NA's   :720     
##      phaeo        nitrate_nitrite_ugL phosphate_ugL     silicate_ugL   
##  Min.   :0.6352   Min.   : 0.4240     Min.   :0.2321   Min.   : 2.296  
##  1st Qu.:0.8792   1st Qu.: 0.5497     1st Qu.:0.3311   1st Qu.: 4.444  
##  Median :2.3793   Median : 5.8840     Median :0.6839   Median :15.939  
##  Mean   :2.2349   Mean   : 4.7549     Mean   :0.6123   Mean   :12.076  
##  3rd Qu.:3.0255   3rd Qu.: 8.4623     3rd Qu.:0.8835   3rd Qu.:18.826  
##  Max.   :3.9090   Max.   :14.2314     Max.   :1.3045   Max.   :25.904  
##  NA's   :720      NA's   :720         NA's   :720      NA's   :720     
##  mg.region.mean.zm.density.plot mg.region.se.zm.density.plot
##  Min.   : 7.43                  Min.   :0.4100              
##  1st Qu.:11.04                  1st Qu.:0.6500              
##  Median :11.04                  Median :0.7700              
##  Mean   :11.82                  Mean   :0.7286              
##  3rd Qu.:13.64                  3rd Qu.:0.7700              
##  Max.   :15.67                  Max.   :1.0300              
##  NA's   :560                    NA's   :560                 
##  mg.region.mean.zm.density.msq mg.site.mean.zm.density.plot
##  Min.   :118.9                 Min.   : 5.58               
##  1st Qu.:176.7                 1st Qu.: 8.67               
##  Median :176.7                 Median :12.17               
##  Mean   :189.1                 Mean   :12.14               
##  3rd Qu.:218.2                 3rd Qu.:15.04               
##  Max.   :250.7                 Max.   :19.13               
##  NA's   :560                   NA's   :1152                
##  mg.site.se.mean.zm.density.plot mg.site.zm.shoot.biomass.site.g
##  Min.   :0.4900                  Min.   :1.100                  
##  1st Qu.:0.5100                  1st Qu.:1.530                  
##  Median :0.9100                  Median :1.760                  
##  Mean   :0.9825                  Mean   :1.707                  
##  3rd Qu.:1.5200                  3rd Qu.:1.900                  
##  Max.   :1.6600                  Max.   :2.110                  
##  NA's   :1152                    NA's   :560                    
##  mg.site.sd.zm.shoot.biomass.g mg.region.mean.zm.shoot.biomass.g
##  Min.   :0.2900                Min.   :1.350                    
##  1st Qu.:0.4300                1st Qu.:1.700                    
##  Median :0.8700                Median :1.830                    
##  Mean   :0.8125                Mean   :1.941                    
##  3rd Qu.:1.2500                3rd Qu.:2.380                    
##  Max.   :1.3000                Max.   :2.380                    
##  NA's   :560                   NA's   :560                      
##  mg.region.sd.zm.shoot.biomass.g mg.site.mean.zm.shoot.biomass.g.msq
##  Min.   :0.5700                  Min.   : 98.56                     
##  1st Qu.:0.7400                  1st Qu.:220.83                     
##  Median :0.9000                  Median :299.22                     
##  Mean   :0.8357                  Mean   :336.82                     
##  3rd Qu.:0.9000                  3rd Qu.:496.33                     
##  Max.   :1.0300                  Max.   :645.92                     
##  NA's   :560                     NA's   :1152                       
##  mg.site.se.zm.shoot.biomass.g.msq mg.region.mean.zm.shoot.biomass.g.msq
##  Min.   :11.74                     Min.   :160.3                        
##  1st Qu.:29.67                     1st Qu.:371.7                        
##  Median :35.68                     Median :420.8                        
##  Mean   :49.09                     Mean   :370.0                        
##  3rd Qu.:88.86                     3rd Qu.:420.8                        
##  Max.   :92.38                     Max.   :459.7                        
##  NA's   :1152                      NA's   :560                          
##  mg.region.se.zm.shoot.biomass.g.msq mg.site.mean.shoot.length.cm
##  Min.   :16.48                       Min.   : 92.81              
##  1st Qu.:41.19                       1st Qu.:100.48              
##  Median :41.19                       Median :109.07              
##  Mean   :39.85                       Mean   :114.22              
##  3rd Qu.:49.01                       3rd Qu.:121.68              
##  Max.   :49.28                       Max.   :152.08              
##  NA's   :560                         NA's   :560                 
##  mg.site.se.shoot.length.cm mg.region.mean.shoot.length.cm
##  Min.   : 5.45              Min.   : 97.24                
##  1st Qu.: 9.33              1st Qu.:103.78                
##  Median :12.57              Median :112.70                
##  Mean   :14.16              Mean   :119.16                
##  3rd Qu.:17.89              3rd Qu.:139.45                
##  Max.   :30.64              Max.   :139.45                
##  NA's   :560                NA's   :560                   
##  mg.region.se.shoot.length.cm mg.site.mean.shoot.epi.biomass.g
##  Min.   :6.860                Min.   :0.0300                  
##  1st Qu.:7.280                1st Qu.:0.0600                  
##  Median :7.960                Median :0.1700                  
##  Mean   :7.673                Mean   :0.4152                  
##  3rd Qu.:7.960                3rd Qu.:0.3100                  
##  Max.   :8.180                Max.   :2.1900                  
##  NA's   :560                  NA's   :1152                    
##  mg.site.se.shoot.epi.biomass.g mg.region.mean.shoot.epi.biomass.g
##  Min.   :0.0000                 Min.   :0.0400                    
##  1st Qu.:0.0300                 1st Qu.:0.1900                    
##  Median :0.0600                 Median :0.2100                    
##  Mean   :0.0836                 Mean   :0.5688                    
##  3rd Qu.:0.1200                 3rd Qu.:1.1700                    
##  Max.   :0.2600                 Max.   :1.1700                    
##  NA's   :1152                   NA's   :560                       
##  mg.region.se.shoot.epi.biomass.g mg.site.mean.meso.biomass.g.per.zm.g
##  Min.   :0.0100                   Min.   :0.0100                      
##  1st Qu.:0.0500                   1st Qu.:0.0200                      
##  Median :0.0700                   Median :0.0400                      
##  Mean   :0.1037                   Mean   :0.0531                      
##  3rd Qu.:0.1900                   3rd Qu.:0.0700                      
##  Max.   :0.1900                   Max.   :0.1400                      
##  NA's   :560                      NA's   :1504                        
##  mg.site.se.meso.biomass.g.per.zm.g mg.region.mean.meso.biomass.g.per.zm.g
##  Min.   :0.0100                     Min.   :0.0200                        
##  1st Qu.:0.0100                     1st Qu.:0.0300                        
##  Median :0.0200                     Median :0.0400                        
##  Mean   :0.0195                     Mean   :0.0711                        
##  3rd Qu.:0.0300                     3rd Qu.:0.1300                        
##  Max.   :0.0300                     Max.   :0.1300                        
##  NA's   :1504                       NA's   :560                           
##  mg.region.se.meso.biomass.g.per.zm.g  bed_area_m2     Description       
##  Min.   :0.0100                       Min.   :  7400   Length:2384       
##  1st Qu.:0.0100                       1st Qu.: 18500   Class :character  
##  Median :0.0200                       Median :191076   Mode  :character  
##  Mean   :0.0204                       Mean   :211586                     
##  3rd Qu.:0.0300                       3rd Qu.:340000                     
##  Max.   :0.0300                       Max.   :569520                     
##  NA's   :560                                                             
##  total_abundance     Rank1              Rank2              Rank3          
##  Min.   :  3141   Length:2384        Length:2384        Length:2384       
##  1st Qu.: 16263   Class :character   Class :character   Class :character  
##  Median : 23035   Mode  :character   Mode  :character   Mode  :character  
##  Mean   : 27267                                                           
##  3rd Qu.: 34284                                                           
##  Max.   :103567                                                           
##                                                                           
##     Rank4              Rank5              Rank6              Rank7          
##  Length:2384        Length:2384        Length:2384        Length:2384       
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##    abundance       relative_abundance 
##  Min.   :    0.0   Min.   :0.0000000  
##  1st Qu.:    7.0   1st Qu.:0.0002844  
##  Median :  100.0   Median :0.0047309  
##  Mean   : 1704.2   Mean   :0.0625000  
##  3rd Qu.:  989.2   3rd Qu.:0.0463993  
##  Max.   :75443.0   Max.   :0.8370517  
## 

Now, lets look at the data a bit bette

table(seagrass$region)
## 
##   choked    goose mcmullin  triquet 
##      944      448      576      416

So we see that we have 4 regions here. If we only want chocked, we use the code we used in lab and check the data again

choked = subset(seagrass, seagrass$region =="choked")

table(choked$region)
## 
## choked 
##    944

Now that we only have choked now and that the number of observations in choked (944) is the same for both dataframes

Now, let’s say we wanted to know what rows we would keep before subsetting, we would use which()

which(seagrass$region == "choked")
##   [1]    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15
##  [16]   16   17   18   19   20   21   22   23   24   25   26   27   28   29   30
##  [31]   31   32   33   34   35   36   37   38   39   40   41   42   43   44   45
##  [46]   46   47   48   49   50   51   52   53   54   55   56   57   58   59   60
##  [61]   61   62   63   64   65   66   67   68   69   70   71   72   73   74   75
##  [76]   76   77   78   79   80   81   82   83   84   85   86   87   88   89   90
##  [91]   91   92   93   94   95   96   97   98   99  100  101  102  103  104  105
## [106]  106  107  108  109  110  111  112  113  114  115  116  117  118  119  120
## [121]  121  122  123  124  125  126  127  128  129  130  131  132  133  134  135
## [136]  136  137  138  139  140  141  142  143  144  145  146  147  148  149  150
## [151]  151  152  153  154  155  156  157  158  159  160  161  162  163  164  165
## [166]  166  167  168  169  170  171  172  173  174  175  176  177  178  179  180
## [181]  181  182  183  184  185  186  187  188  189  190  191  192  193  194  195
## [196]  196  197  198  199  200  201  202  203  204  205  206  207  208  209  210
## [211]  211  212  213  214  215  216  217  218  219  220  221  222  223  224  225
## [226]  226  227  228  229  230  231  232  233  234  235  236  237  238  239  240
## [241]  241  242  243  244  245  246  247  248  249  250  251  252  253  254  255
## [256]  256  257  258  259  260  261  262  263  264  265  266  267  268  269  270
## [271]  271  272  273  274  275  276  277  278  279  280  281  282  283  284  285
## [286]  286  287  288  289  290  291  292  293  294  295  296  297  298  299  300
## [301]  301  302  303  304  305  306  307  308  309  310  311  312  313  314  315
## [316]  316  317  318  319  320  321  322  323  324  325  326  327  328  329  330
## [331]  331  332  333  334  335  336  337  338  339  340  341  342  343  344  345
## [346]  346  347  348  349  350  351  352  353  354  355  356  357  358  359  360
## [361]  361  362  363  364  365  366  367  368  369  370  371  372  373  374  375
## [376]  376  377  378  379  380  381  382  383  384  385  386  387  388  389  390
## [391]  391  392  393  394  395  396  397  398  399  400  401  402  403  404  405
## [406]  406  407  408  409  410  411  412  413  414  415  416  417  418  419  420
## [421]  421  422  423  424  425  426  427  428  429  430  431  432  433  434  435
## [436]  436  437  438  439  440  441  442  443  444  445  446  447  448  449  450
## [451]  451  452  453  454  455  456  457  458  459  460  461  462  463  464  465
## [466]  466  467  468  469  470  471  472  473  474  475  476  477  478  479  480
## [481]  481  482  483  484  485  486  487  488  489  490  491  492  493  494  495
## [496]  496  497  498  499  500  501  502  503  504  505  506  507  508  509  510
## [511]  511  512  513  514  515  516  517  518  519  520  521  522  523  524  525
## [526]  526  527  528  529  530  531  532  533  534  535  536  537  538  539  540
## [541]  541  542  543  544  545  546  547  548  549  550  551  552  553  554  555
## [556]  556  557  558  559  560  561  562  563  564  565  566  567  568  569  570
## [571] 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590
## [586] 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605
## [601] 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620
## [616] 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635
## [631] 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650
## [646] 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665
## [661] 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680
## [676] 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695
## [691] 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710
## [706] 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725
## [721] 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740
## [736] 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755
## [751] 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770
## [766] 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785
## [781] 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800
## [796] 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815
## [811] 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830
## [826] 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845
## [841] 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860
## [856] 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875
## [871] 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890
## [886] 2236 2237 2238 2239 2240 2241 2242 2243 2244 2245 2246 2247 2248 2249 2250
## [901] 2251 2252 2253 2254 2255 2256 2257 2258 2259 2260 2261 2262 2263 2264 2265
## [916] 2266 2267 2268 2269 2270 2271 2272 2273 2341 2342 2343 2344 2345 2346 2347
## [931] 2348 2349 2350 2351 2352 2353 2354 2355 2356 2357 2358 2359 2360 2361

this is more useful if you are trying to remove specific numeric data

which(seagrass$relative_abundance>0.8)
## [1]  763 2278 2279 2292 2325 2369

This tell use that in the seagrass dataset, the only rows where the relative abundance is greater than 0.8 are the rows 763 2278 2279 2292 2325 and 2369

We could then only keep these rows if that’s something we were interested in

# subset
seagrass_08 = subset(seagrass, seagrass$relative_abundance>0.8)

#we can check with a table
table(seagrass_08) #see that it's a 6x6 column
## < table of extent 6 x 6 x 0 x 6 x 2 x 6 x 1 x 1 x 2 x 1 x 1 x 3 x 4 x 4 x 4 x 2 x 2 x 2 x 1 x 2 x 6 x 2 x 3 x 5 x 3 x 4 x 0 x 4 x 4 x 1 x 3 x 4 x 2 x 4 x 1 x 1 x 1 x 2 x 1 x 1 x 4 x 3 x 4 x 4 x 4 x 4 x 4 x 4 x 4 x 4 x 4 x 4 x 4 x 1 x 1 x 1 x 1 x 1 x 2 x 2 x 2 x 3 x 3 x 3 x 3 x 2 x 2 x 3 x 3 x 2 x 2 x 3 x 3 x 2 x 2 x 3 x 3 x 2 x 2 x 1 x 1 x 2 x 1 x 4 x 6 x 6 x 2 x 1 x 1 x 1 x 1 x 1 x 1 x 6 x 6 >
# we can check with summary to make sure the values ma sense
summary(seagrass_08$relative_abundance)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##  0.8012  0.8160  0.8206  0.8196  0.8231  0.8371

We can also use other ways to check the data

head(seagrass_08$relative_abundance)
## [1] 0.8184569 0.8227821 0.8231488 0.8151271 0.8370517 0.8012140
list(seagrass_08$relative_abundance)
## [[1]]
## [1] 0.8184569 0.8227821 0.8231488 0.8151271 0.8370517 0.8012140

Do you see a difference between these two ways of checking? What about with the entire seagrass dataset? The output with the full dataset it at the end of the page. For now, let’s do one more exmaple

Subset choked but only keep relative abundance greater or equal to 0.5

choked_05 = subset(seagrass, seagrass$region=="choked"&seagrass$relative_abundance>=0.5)

let’s check that it worked

#see the minimum value
min(choked_05$relative_abundance)
## [1] 0.5024889
#see the maximum value
max(choked_05$relative_abundance)
## [1] 0.787434
#or just use summary to see everything at once
summary(choked_05)
##    SampleID         BarcodeSequence    LinkerPrimerSequence   swab_id         
##  Length:44          Length:44          Mode:logical         Length:44         
##  Class :character   Class :character   NA's:44              Class :character  
##  Mode  :character   Mode  :character                        Mode  :character  
##                                                                               
##                                                                               
##                                                                               
##                                                                               
##  barcode_plate      barcode_well       project_name       person_responsible
##  Length:44          Length:44          Length:44          Length:44         
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##      date              Run_id          sequence_file         region         
##  Length:44          Length:44          Length:44          Length:44         
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##      site              latitude       longitude      sample_method     
##  Length:44          Min.   :51.67   Min.   :-128.1   Length:44         
##  Class :character   1st Qu.:51.67   1st Qu.:-128.1   Class :character  
##  Mode  :character   Median :51.67   Median :-128.1   Mode  :character  
##                     Mean   :51.67   Mean   :-128.1                     
##                     3rd Qu.:51.68   3rd Qu.:-128.1                     
##                     Max.   :51.68   Max.   :-128.1                     
##                                                                        
##  host_species        host_type         survey_type        dna_extraction_plate
##  Length:44          Length:44          Length:44          Length:44           
##  Class :character   Class :character   Class :character   Class :character    
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character    
##                                                                               
##                                                                               
##                                                                               
##                                                                               
##  dna_extraction_well    growth          sample_growth       quadrat_id       
##  Length:44           Length:44          Length:44          Length:44         
##  Class :character    Class :character   Class :character   Class :character  
##  Mode  :character    Mode  :character   Mode  :character   Mode  :character  
##                                                                              
##                                                                              
##                                                                              
##                                                                              
##  meso_quadrat_density_shoots_m2 meso_quadrat_biomass_g_m2 meso_quadrats_epimass
##  Min.   :12.00                  Min.   : 9.60             Min.   : 0.400       
##  1st Qu.:16.00                  1st Qu.:15.40             1st Qu.: 1.000       
##  Median :20.00                  Median :26.60             Median : 1.200       
##  Mean   :21.28                  Mean   :22.26             Mean   : 3.564       
##  3rd Qu.:24.00                  3rd Qu.:28.20             3rd Qu.: 3.200       
##  Max.   :40.00                  Max.   :35.80             Max.   :10.500       
##  NA's   :19                     NA's   :19                NA's   :19           
##  meso_shoot_id      meso_shoot_length meso_shoot_width meso_number_of_blades
##  Length:44          Min.   : 77.0     Min.   :0.3000   Min.   :4.0          
##  Class :character   1st Qu.:106.5     1st Qu.:0.5000   1st Qu.:5.0          
##  Mode  :character   Median :113.5     Median :0.7000   Median :5.0          
##                     Mean   :116.3     Mean   :0.6909   Mean   :5.5          
##                     3rd Qu.:130.8     3rd Qu.:0.8750   3rd Qu.:6.0          
##                     Max.   :147.1     Max.   :1.1000   Max.   :8.0          
##                     NA's   :22        NA's   :22       NA's   :22           
##  meso_total_dry_wt meso_blade_dry_wt meso_microepiphyte_wt
##  Min.   :1.300     Min.   :0.30      Min.   :10.40        
##  1st Qu.:1.900     1st Qu.:0.40      1st Qu.:16.80        
##  Median :2.600     Median :0.40      Median :19.60        
##  Mean   :2.459     Mean   :0.45      Mean   :19.70        
##  3rd Qu.:2.875     3rd Qu.:0.50      3rd Qu.:23.85        
##  Max.   :5.200     Max.   :1.00      Max.   :26.20        
##  NA's   :22        NA's   :22        NA's   :22           
##  meso_macroepiphyte_wt meso_smithora_wt meso_bryozoans        broken         
##  Min.   :0.0000        Min.   :0.0000   Length:44          Length:44         
##  1st Qu.:0.0000        1st Qu.:0.0000   Class :character   Class :character  
##  Median :0.0000        Median :0.0000   Mode  :character   Mode  :character  
##  Mean   :0.4636        Mean   :0.5182                                        
##  3rd Qu.:0.0000        3rd Qu.:0.2750                                        
##  Max.   :4.2000        Max.   :4.1000                                        
##  NA's   :22            NA's   :22                                            
##  smithora_shoot     smithora_region      ocean_id           mean_depth   
##  Length:44          Length:44          Length:44          Min.   : 0.50  
##  Class :character   Class :character   Class :character   1st Qu.: 4.00  
##  Mode  :character   Mode  :character   Mode  :character   Median : 5.00  
##                                                           Mean   : 4.75  
##                                                           3rd Qu.: 5.00  
##                                                           Max.   :10.00  
##                                                                          
##    ctd_depth      conductivity    temperature       pressure    
##  Min.   :0.500   Min.   :33.85   Min.   :10.28   Min.   :10.64  
##  1st Qu.:4.520   1st Qu.:35.00   1st Qu.:11.18   1st Qu.:14.69  
##  Median :4.570   Median :35.02   Median :11.54   Median :14.74  
##  Mean   :4.956   Mean   :34.98   Mean   :11.28   Mean   :14.84  
##  3rd Qu.:5.045   3rd Qu.:35.04   3rd Qu.:11.54   3rd Qu.:14.77  
##  Max.   :9.490   Max.   :35.72   Max.   :13.31   Max.   :19.71  
##                                                                 
##       par         fluorometry_chlorophyll   turbidity      dissolved_oxygen
##  Min.   : 80.88   Min.   :-0.310          Min.   : 0.600   Min.   : 79.82  
##  1st Qu.:138.20   1st Qu.: 1.532          1st Qu.: 0.705   1st Qu.: 95.55  
##  Median :141.25   Median : 2.370          Median :10.420   Median :108.37  
##  Mean   :176.35   Mean   : 1.939          Mean   : 8.894   Mean   :101.34  
##  3rd Qu.:154.68   3rd Qu.: 2.370          3rd Qu.:16.820   3rd Qu.:108.37  
##  Max.   :608.45   Max.   : 2.370          Max.   :16.820   Max.   :124.31  
##                                                                            
##     salinity     speed_of_sound dissolved_oxygen_concentration
##  Min.   :29.80   Min.   :1486   Min.   : 7.370                
##  1st Qu.:30.56   1st Qu.:1489   1st Qu.: 8.620                
##  Median :30.56   Median :1490   Median : 9.720                
##  Mean   :30.74   Mean   :1489   Mean   : 9.121                
##  3rd Qu.:30.82   3rd Qu.:1490   3rd Qu.: 9.720                
##  Max.   :31.43   Max.   :1495   Max.   :10.790                
##                                                               
##      chl_a             phaeo       nitrate_nitrite_ugL phosphate_ugL   
##  Min.   :0.04490   Min.   :2.159   Min.   : 1.497      Min.   :0.2697  
##  1st Qu.:0.04490   1st Qu.:2.360   1st Qu.: 5.884      1st Qu.:0.6839  
##  Median :0.06185   Median :2.379   Median : 5.884      Median :0.6839  
##  Mean   :0.07328   Mean   :2.668   Mean   : 7.579      Mean   :0.8072  
##  3rd Qu.:0.08250   3rd Qu.:2.679   3rd Qu.: 8.462      3rd Qu.:0.8835  
##  Max.   :0.29180   Max.   :3.909   Max.   :14.231      Max.   :1.3045  
##                                                                        
##   silicate_ugL   mg.region.mean.zm.density.plot mg.region.se.zm.density.plot
##  Min.   : 6.29   Min.   :11.04                  Min.   :0.77                
##  1st Qu.:15.94   1st Qu.:11.04                  1st Qu.:0.77                
##  Median :15.94   Median :11.04                  Median :0.77                
##  Mean   :17.52   Mean   :11.04                  Mean   :0.77                
##  3rd Qu.:18.83   3rd Qu.:11.04                  3rd Qu.:0.77                
##  Max.   :25.90   Max.   :11.04                  Max.   :0.77                
##                  NA's   :11                     NA's   :11                  
##  mg.region.mean.zm.density.msq mg.site.mean.zm.density.plot
##  Min.   :176.7                 Min.   :11.89               
##  1st Qu.:176.7                 1st Qu.:11.89               
##  Median :176.7                 Median :11.89               
##  Mean   :176.7                 Mean   :11.89               
##  3rd Qu.:176.7                 3rd Qu.:11.89               
##  Max.   :176.7                 Max.   :11.89               
##  NA's   :11                    NA's   :37                  
##  mg.site.se.mean.zm.density.plot mg.site.zm.shoot.biomass.site.g
##  Min.   :1.66                    Min.   :1.760                  
##  1st Qu.:1.66                    1st Qu.:1.850                  
##  Median :1.66                    Median :1.900                  
##  Mean   :1.66                    Mean   :1.864                  
##  3rd Qu.:1.66                    3rd Qu.:1.900                  
##  Max.   :1.66                    Max.   :1.900                  
##  NA's   :37                      NA's   :11                     
##  mg.site.sd.zm.shoot.biomass.g mg.region.mean.zm.shoot.biomass.g
##  Min.   :0.430                 Min.   :2.38                     
##  1st Qu.:1.040                 1st Qu.:2.38                     
##  Median :1.300                 Median :2.38                     
##  Mean   :1.068                 Mean   :2.38                     
##  3rd Qu.:1.300                 3rd Qu.:2.38                     
##  Max.   :1.300                 Max.   :2.38                     
##  NA's   :11                    NA's   :11                       
##  mg.region.sd.zm.shoot.biomass.g mg.site.mean.zm.shoot.biomass.g.msq
##  Min.   :0.9                     Min.   :352.5                      
##  1st Qu.:0.9                     1st Qu.:352.5                      
##  Median :0.9                     Median :352.5                      
##  Mean   :0.9                     Mean   :352.5                      
##  3rd Qu.:0.9                     3rd Qu.:352.5                      
##  Max.   :0.9                     Max.   :352.5                      
##  NA's   :11                      NA's   :37                         
##  mg.site.se.zm.shoot.biomass.g.msq mg.region.mean.zm.shoot.biomass.g.msq
##  Min.   :55.49                     Min.   :420.8                        
##  1st Qu.:55.49                     1st Qu.:420.8                        
##  Median :55.49                     Median :420.8                        
##  Mean   :55.49                     Mean   :420.8                        
##  3rd Qu.:55.49                     3rd Qu.:420.8                        
##  Max.   :55.49                     Max.   :420.8                        
##  NA's   :37                        NA's   :11                           
##  mg.region.se.zm.shoot.biomass.g.msq mg.site.mean.shoot.length.cm
##  Min.   :41.19                       Min.   :100.5               
##  1st Qu.:41.19                       1st Qu.:121.7               
##  Median :41.19                       Median :121.7               
##  Mean   :41.19                       Mean   :122.7               
##  3rd Qu.:41.19                       3rd Qu.:121.7               
##  Max.   :41.19                       Max.   :152.1               
##  NA's   :11                          NA's   :11                  
##  mg.site.se.shoot.length.cm mg.region.mean.shoot.length.cm
##  Min.   : 5.88              Min.   :139.4                 
##  1st Qu.:17.89              1st Qu.:139.4                 
##  Median :17.89              Median :139.4                 
##  Mean   :17.66              Mean   :139.4                 
##  3rd Qu.:17.89              3rd Qu.:139.4                 
##  Max.   :30.64              Max.   :139.4                 
##  NA's   :11                 NA's   :11                    
##  mg.region.se.shoot.length.cm mg.site.mean.shoot.epi.biomass.g
##  Min.   :7.96                 Min.   :2.19                    
##  1st Qu.:7.96                 1st Qu.:2.19                    
##  Median :7.96                 Median :2.19                    
##  Mean   :7.96                 Mean   :2.19                    
##  3rd Qu.:7.96                 3rd Qu.:2.19                    
##  Max.   :7.96                 Max.   :2.19                    
##  NA's   :11                   NA's   :37                      
##  mg.site.se.shoot.epi.biomass.g mg.region.mean.shoot.epi.biomass.g
##  Min.   :0.26                   Min.   :1.17                      
##  1st Qu.:0.26                   1st Qu.:1.17                      
##  Median :0.26                   Median :1.17                      
##  Mean   :0.26                   Mean   :1.17                      
##  3rd Qu.:0.26                   3rd Qu.:1.17                      
##  Max.   :0.26                   Max.   :1.17                      
##  NA's   :37                     NA's   :11                        
##  mg.region.se.shoot.epi.biomass.g mg.site.mean.meso.biomass.g.per.zm.g
##  Min.   :0.19                     Min.   :0.14                        
##  1st Qu.:0.19                     1st Qu.:0.14                        
##  Median :0.19                     Median :0.14                        
##  Mean   :0.19                     Mean   :0.14                        
##  3rd Qu.:0.19                     3rd Qu.:0.14                        
##  Max.   :0.19                     Max.   :0.14                        
##  NA's   :11                       NA's   :37                          
##  mg.site.se.meso.biomass.g.per.zm.g mg.region.mean.meso.biomass.g.per.zm.g
##  Min.   :0.03                       Min.   :0.13                          
##  1st Qu.:0.03                       1st Qu.:0.13                          
##  Median :0.03                       Median :0.13                          
##  Mean   :0.03                       Mean   :0.13                          
##  3rd Qu.:0.03                       3rd Qu.:0.13                          
##  Max.   :0.03                       Max.   :0.13                          
##  NA's   :37                         NA's   :11                            
##  mg.region.se.meso.biomass.g.per.zm.g  bed_area_m2     Description       
##  Min.   :0.03                         Min.   :340000   Length:44         
##  1st Qu.:0.03                         1st Qu.:340000   Class :character  
##  Median :0.03                         Median :340000   Mode  :character  
##  Mean   :0.03                         Mean   :340000                     
##  3rd Qu.:0.03                         3rd Qu.:340000                     
##  Max.   :0.03                         Max.   :340000                     
##  NA's   :11                                                              
##  total_abundance    Rank1              Rank2              Rank3          
##  Min.   : 3141   Length:44          Length:44          Length:44         
##  1st Qu.:14559   Class :character   Class :character   Class :character  
##  Median :24798   Mode  :character   Mode  :character   Mode  :character  
##  Mean   :29074                                                           
##  3rd Qu.:39049                                                           
##  Max.   :95105                                                           
##                                                                          
##     Rank4              Rank5              Rank6              Rank7          
##  Length:44          Length:44          Length:44          Length:44         
##  Class :character   Class :character   Class :character   Class :character  
##  Mode  :character   Mode  :character   Mode  :character   Mode  :character  
##                                                                             
##                                                                             
##                                                                             
##                                                                             
##    abundance     relative_abundance
##  Min.   : 1803   Min.   :0.5025    
##  1st Qu.: 8440   1st Qu.:0.5422    
##  Median :14434   Median :0.5848    
##  Mean   :17431   Mean   :0.5887    
##  3rd Qu.:22264   3rd Qu.:0.6200    
##  Max.   :52479   Max.   :0.7874    
## 
head(seagrass$relative_abundance)
## [1] 0.004622496 0.017402338 0.000000000 0.000000000 0.017221064 0.037342518
list(seagrass$relative_abundance)
## [[1]]
##    [1] 4.622496e-03 1.740234e-02 0.000000e+00 0.000000e+00 1.722106e-02
##    [6] 3.734252e-02 4.250884e-02 1.060455e-02 1.237197e-01 3.625487e-03
##   [11] 4.513732e-02 5.827971e-02 1.894317e-02 5.365721e-02 2.528777e-02
##   [16] 0.000000e+00 3.824757e-03 6.171766e-03 1.303894e-03 8.257997e-03
##   [21] 2.433936e-02 1.990612e-02 2.399166e-02 1.310848e-01 7.823366e-04
##   [26] 1.416898e-02 1.106572e-01 7.927677e-02 8.431850e-03 3.755216e-02
##   [31] 0.000000e+00 2.475452e-03 3.053057e-03 1.732816e-03 1.006684e-02
##   [36] 2.888027e-02 9.819292e-03 1.163462e-02 1.393679e-01 1.155211e-03
##   [41] 1.939104e-02 1.592541e-02 1.694034e-01 3.548148e-02 3.721429e-02
##   [46] 5.730878e-04 2.254145e-03 1.872087e-03 5.348819e-04 3.018262e-03
##   [51] 8.749140e-03 1.352487e-02 4.126232e-03 3.396500e-02 2.980057e-03
##   [56] 1.230228e-02 5.757622e-02 5.167342e-01 5.501643e-03 4.026897e-02
##   [61] 1.250813e-04 2.426577e-03 9.506179e-03 4.002602e-04 1.255816e-02
##   [66] 1.490969e-02 9.431130e-03 1.751138e-02 4.855656e-02 5.203382e-03
##   [71] 1.460950e-02 3.617351e-02 1.589533e-01 1.412668e-01 1.693601e-02
##   [76] 0.000000e+00 0.000000e+00 3.098435e-03 0.000000e+00 0.000000e+00
##   [81] 6.395487e-03 4.846270e-03 1.747835e-03 1.195678e-02 1.191706e-03
##   [86] 3.972352e-03 1.481687e-02 5.195043e-01 2.705172e-02 4.250417e-03
##   [91] 0.000000e+00 1.786584e-02 3.602376e-03 2.872164e-03 1.538312e-02
##   [96] 6.362574e-02 7.633142e-02 1.869341e-02 4.337455e-02 1.168338e-03
##  [101] 4.829131e-02 1.260831e-02 9.249343e-04 1.217019e-03 1.429267e-01
##  [106] 1.290622e-04 1.626184e-03 6.969361e-04 8.259983e-04 3.007150e-02
##  [111] 7.841821e-02 9.826798e-02 5.446426e-03 1.419685e-02 0.000000e+00
##  [116] 1.164141e-02 1.910121e-03 1.290622e-04 1.290622e-04 1.449111e-01
##  [121] 0.000000e+00 6.354219e-03 9.627604e-04 8.472292e-04 1.474949e-02
##  [126] 6.681557e-02 6.851003e-02 3.646936e-02 1.725267e-02 5.391458e-04
##  [131] 1.987138e-02 3.234875e-03 1.540417e-04 0.000000e+00 1.762622e-01
##  [136] 0.000000e+00 1.033247e-02 2.328076e-03 6.016375e-04 1.216354e-02
##  [141] 7.523085e-02 1.307384e-01 2.288838e-02 1.593031e-02 4.446886e-04
##  [146] 1.459625e-02 2.301917e-03 1.569489e-04 3.400560e-04 1.420388e-01
##  [151] 0.000000e+00 5.191257e-03 3.708041e-03 6.635441e-04 1.494926e-02
##  [156] 8.290398e-02 1.127244e-01 1.280250e-02 1.912568e-02 3.903201e-04
##  [161] 1.327088e-02 6.245121e-04 2.732240e-04 2.341920e-04 1.023419e-01
##  [166] 0.000000e+00 6.014787e-03 3.922687e-03 6.894420e-04 1.595226e-02
##  [171] 9.782945e-02 1.035352e-01 2.037420e-02 1.854361e-02 8.796329e-04
##  [176] 1.673680e-02 3.898914e-03 2.852864e-04 5.230250e-04 1.257162e-01
##  [181] 0.000000e+00 0.000000e+00 4.119698e-03 0.000000e+00 2.761556e-03
##  [186] 1.041242e-02 1.072932e-02 1.222328e-03 1.780977e-01 1.557336e-02
##  [191] 1.000498e-02 6.374213e-02 2.607633e-01 9.688080e-03 1.860655e-02
##  [196] 0.000000e+00 4.185852e-04 7.800906e-03 0.000000e+00 9.475246e-03
##  [201] 1.742836e-02 1.681951e-02 1.217702e-03 1.407968e-01 2.914875e-02
##  [206] 2.682750e-02 1.434606e-01 3.013813e-02 7.648693e-03 2.648503e-02
##  [211] 0.000000e+00 8.381047e-03 3.999668e-03 9.957680e-05 4.978840e-03
##  [216] 3.364036e-02 4.729898e-03 6.804415e-04 1.141648e-01 9.360219e-03
##  [221] 6.356319e-03 2.837939e-02 1.521533e-01 4.022903e-02 8.281470e-03
##  [226] 0.000000e+00 0.000000e+00 1.153566e-03 0.000000e+00 3.238080e-04
##  [231] 6.476160e-04 1.639278e-03 0.000000e+00 3.703554e-03 1.042257e-02
##  [236] 6.071400e-04 3.068081e-02 2.978426e-01 4.058326e-01 5.059500e-04
##  [241] 0.000000e+00 0.000000e+00 1.054305e-02 1.394583e-04 0.000000e+00
##  [246] 6.136167e-03 9.622626e-03 0.000000e+00 3.690068e-02 3.497615e-02
##  [251] 2.426575e-03 8.334031e-02 2.912727e-01 4.808524e-02 7.809667e-03
##  [256] 0.000000e+00 0.000000e+00 0.000000e+00 2.314761e-04 3.240666e-03
##  [261] 3.911946e-02 2.060137e-03 9.953473e-04 1.193954e-01 4.513784e-03
##  [266] 5.317007e-02 1.893475e-02 2.728409e-01 1.944168e-01 1.326358e-02
##  [271] 0.000000e+00 2.135611e-04 0.000000e+00 6.406834e-04 0.000000e+00
##  [276] 2.135611e-03 2.669514e-04 3.203417e-04 1.793914e-02 2.349172e-03
##  [281] 1.708489e-03 3.828083e-02 4.366257e-01 2.437800e-01 3.416978e-03
##  [286] 0.000000e+00 0.000000e+00 0.000000e+00 1.237815e-03 2.217752e-03
##  [291] 1.340967e-02 3.146114e-03 2.991387e-03 7.591934e-02 1.284233e-02
##  [296] 1.170767e-02 7.014286e-02 3.771726e-01 2.506576e-02 1.748414e-02
##  [301] 1.129995e-04 1.581992e-04 1.105135e-02 0.000000e+00 3.028385e-03
##  [306] 2.018170e-02 2.259989e-03 1.468993e-03 1.419273e-01 8.994757e-03
##  [311] 1.145815e-02 2.870186e-03 5.613813e-02 4.098490e-01 9.604954e-03
##  [316] 9.469473e-05 0.000000e+00 3.314315e-04 9.469473e-05 1.444095e-03
##  [321] 3.693094e-03 9.469473e-05 1.112663e-03 3.205417e-02 2.367368e-03
##  [326] 3.195947e-03 5.847399e-03 7.874340e-01 1.422788e-02 1.704505e-03
##  [331] 0.000000e+00 2.558805e-02 3.973299e-03 0.000000e+00 1.589320e-02
##  [336] 6.436745e-02 2.876669e-02 7.104259e-02 1.684679e-02 5.085823e-03
##  [341] 9.917355e-02 1.112524e-02 0.000000e+00 0.000000e+00 5.578512e-02
##  [346] 0.000000e+00 5.933807e-03 1.969410e-03 5.882654e-04 2.654867e-02
##  [351] 1.178833e-01 3.895340e-02 1.422068e-02 2.984807e-02 2.225178e-03
##  [356] 1.755589e-01 1.276280e-02 1.278838e-04 0.000000e+00 8.852115e-02
##  [361] 0.000000e+00 2.915695e-03 1.194046e-03 0.000000e+00 1.782739e-02
##  [366] 3.040653e-02 6.675553e-02 2.196490e-02 6.145174e-02 2.554704e-03
##  [371] 2.826835e-02 4.470732e-03 1.666111e-04 2.776852e-04 1.591136e-01
##  [376] 0.000000e+00 1.277248e-03 0.000000e+00 0.000000e+00 8.429837e-03
##  [381] 2.669448e-01 3.576294e-03 1.702997e-03 2.741826e-02 1.106948e-03
##  [386] 2.978542e-01 1.592302e-02 9.366485e-04 8.514986e-04 2.733311e-02
##  [391] 0.000000e+00 4.731958e-03 7.326902e-04 1.037978e-03 2.292710e-02
##  [396] 1.719074e-01 1.114300e-02 7.174258e-03 1.694346e-02 8.120650e-03
##  [401] 1.973684e-01 4.240444e-02 4.884601e-04 1.160093e-03 1.828673e-02
##  [406] 0.000000e+00 1.910220e-03 0.000000e+00 0.000000e+00 2.929004e-02
##  [411] 1.862464e-01 4.457179e-03 7.322509e-03 2.833493e-02 7.959249e-03
##  [416] 1.219357e-01 7.004139e-03 1.273480e-03 2.865330e-03 2.737982e-02
##  [421] 0.000000e+00 5.126745e-03 4.746986e-04 0.000000e+00 2.088674e-02
##  [426] 1.609228e-01 2.895661e-02 1.519035e-03 1.402260e-01 7.690117e-03
##  [431] 4.784962e-02 6.171081e-03 0.000000e+00 0.000000e+00 7.016045e-02
##  [436] 0.000000e+00 0.000000e+00 1.837823e-03 0.000000e+00 2.017123e-03
##  [441] 2.693980e-02 9.502891e-03 5.244520e-03 3.541172e-03 1.197723e-01
##  [446] 1.631628e-02 1.434847e-01 3.585997e-03 2.106773e-03 2.241248e-03
##  [451] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 4.554985e-03
##  [456] 5.928711e-03 4.554985e-03 1.662931e-03 9.399176e-03 2.012147e-01
##  [461] 4.699588e-03 1.171282e-01 1.865375e-02 4.410382e-03 1.735232e-03
##  [466] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
##  [471] 6.955336e-03 1.295602e-03 0.000000e+00 3.409478e-04 2.857143e-02
##  [476] 0.000000e+00 4.364132e-02 1.264235e-01 2.783498e-01 1.193317e-02
##  [481] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 7.680492e-03
##  [486] 1.056068e-02 6.336406e-03 0.000000e+00 1.365207e-01 1.899002e-01
##  [491] 1.728111e-03 1.152074e-03 1.071429e-01 9.581413e-02 0.000000e+00
##  [496] 0.000000e+00 2.345100e-02 1.086152e-03 0.000000e+00 9.429770e-03
##  [501] 1.729943e-01 1.476179e-02 8.837324e-03 2.146137e-01 3.460874e-02
##  [506] 3.120217e-02 3.174525e-02 1.974821e-04 4.887682e-03 2.666008e-03
##  [511] 0.000000e+00 1.394322e-02 7.322700e-03 5.015548e-04 4.985455e-02
##  [516] 1.055271e-01 4.393620e-02 1.534758e-02 1.093390e-02 6.319591e-03
##  [521] 8.205437e-02 2.848831e-02 6.018658e-04 1.203732e-03 1.594944e-02
##  [526] 0.000000e+00 1.175986e-03 1.117187e-03 0.000000e+00 5.674134e-02
##  [531] 1.831011e-01 4.292350e-02 1.217146e-02 1.264185e-02 4.351150e-03
##  [536] 9.378491e-02 1.399424e-02 0.000000e+00 0.000000e+00 1.816899e-02
##  [541] 2.047418e-04 5.732771e-04 2.047418e-03 3.685353e-04 6.351091e-02
##  [546] 8.791614e-02 4.971131e-02 6.920274e-03 5.732771e-03 1.920478e-02
##  [551] 3.861431e-02 4.668114e-03 4.504320e-04 6.961222e-04 8.394415e-03
##  [556] 0.000000e+00 2.429441e-01 7.378712e-04 0.000000e+00 3.308123e-02
##  [561] 2.759024e-01 1.666359e-02 1.906167e-03 1.085901e-01 3.689356e-03
##  [566] 4.660887e-02 4.857652e-03 0.000000e+00 2.459571e-04 1.401955e-02
##  [571] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
##  [576] 5.436912e-03 6.563480e-03 0.000000e+00 7.674373e-01 2.840909e-03
##  [581] 7.200235e-03 8.571708e-03 2.889890e-03 3.575627e-03 1.126567e-03
##  [586] 0.000000e+00 5.656749e-04 1.131350e-04 7.070936e-04 7.580043e-03
##  [591] 2.302297e-02 5.091074e-03 0.000000e+00 6.056963e-01 1.951578e-03
##  [596] 6.762643e-02 2.027944e-02 7.373572e-02 5.291888e-02 1.131350e-04
##  [601] 3.684760e-03 5.499643e-03 0.000000e+00 0.000000e+00 3.316284e-02
##  [606] 5.961612e-02 3.277787e-02 3.849750e-04 2.180058e-01 2.964307e-02
##  [611] 3.690260e-02 2.067866e-02 6.698565e-02 1.319914e-03 7.534510e-03
##  [616] 0.000000e+00 0.000000e+00 1.078477e-03 0.000000e+00 1.797462e-04
##  [621] 1.150376e-03 1.150376e-03 5.032894e-04 3.321710e-02 2.156954e-03
##  [626] 1.689614e-03 5.191070e-02 3.589532e-01 3.451127e-03 2.839990e-03
##  [631] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 3.343160e-04
##  [636] 1.701972e-03 8.266723e-03 3.647084e-04 2.656293e-02 2.765705e-03
##  [641] 1.884327e-03 1.016017e-01 3.297572e-02 0.000000e+00 7.598091e-04
##  [646] 0.000000e+00 1.547378e-03 2.532073e-04 6.752194e-04 3.300135e-02
##  [651] 4.023183e-02 1.567072e-02 6.555255e-03 1.547378e-02 1.434841e-03
##  [656] 6.065721e-02 5.064146e-04 1.969390e-04 1.969390e-04 4.501463e-04
##  [661] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.443664e-02
##  [666] 5.399473e-02 2.546093e-02 3.072871e-03 1.097454e-02 0.000000e+00
##  [671] 1.068188e-01 0.000000e+00 0.000000e+00 0.000000e+00 1.902253e-03
##  [676] 0.000000e+00 4.298995e-04 2.149498e-04 0.000000e+00 3.777742e-02
##  [681] 8.361545e-02 4.938471e-02 2.686872e-04 4.836369e-03 5.911118e-04
##  [686] 2.315009e-01 5.373744e-04 2.686872e-04 0.000000e+00 8.114353e-03
##  [691] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.956659e-02
##  [696] 9.062994e-02 8.001404e-02 1.052816e-03 3.509388e-03 0.000000e+00
##  [701] 1.786278e-01 0.000000e+00 0.000000e+00 0.000000e+00 1.579224e-03
##  [706] 0.000000e+00 0.000000e+00 0.000000e+00 4.154919e-04 2.955928e-02
##  [711] 5.514171e-02 3.139932e-02 3.502003e-03 1.341445e-02 5.045259e-04
##  [716] 1.184746e-01 5.638819e-04 3.561359e-04 0.000000e+00 2.967799e-03
##  [721] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
##  [726] 2.828779e-02 1.110162e-02 0.000000e+00 1.026900e-01 6.084543e-03
##  [731] 1.174210e-02 5.924424e-02 2.113578e-02 0.000000e+00 4.910333e-03
##  [736] 0.000000e+00 6.199262e-04 0.000000e+00 5.018450e-04 6.464945e-03
##  [741] 4.841328e-02 7.380074e-03 0.000000e+00 6.219926e-01 9.033210e-03
##  [746] 5.177860e-02 7.321033e-03 5.904059e-03 1.771218e-04 6.110701e-03
##  [751] 1.871900e-04 4.648551e-03 0.000000e+00 0.000000e+00 1.341528e-03
##  [756] 5.272517e-03 4.398964e-03 0.000000e+00 2.127726e-02 4.024584e-03
##  [761] 3.899791e-03 7.328487e-02 8.184569e-01 9.047515e-04 0.000000e+00
##  [766] 3.348599e-03 2.314473e-03 0.000000e+00 7.386615e-04 1.920520e-03
##  [771] 7.928301e-03 3.481558e-02 8.371498e-04 1.351258e-01 1.629980e-02
##  [776] 7.091151e-03 3.517014e-01 1.993894e-01 3.397843e-03 0.000000e+00
##  [781] 4.569563e-03 7.343941e-04 0.000000e+00 0.000000e+00 3.459812e-02
##  [786] 5.948592e-02 1.550388e-02 0.000000e+00 1.819665e-01 5.401877e-02
##  [791] 6.128111e-02 7.278662e-02 2.301102e-02 1.468788e-02 5.467156e-03
##  [796] 5.501458e-03 0.000000e+00 2.200583e-04 1.485394e-03 2.475656e-03
##  [801] 4.181108e-03 5.837047e-02 2.750729e-04 6.024096e-02 1.747813e-01
##  [806] 1.358860e-02 8.554767e-02 3.223854e-02 3.575948e-03 6.051604e-04
##  [811] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 5.649981e-03
##  [816] 2.362719e-02 4.776802e-02 7.471050e-04 8.274188e-02 1.120657e-03
##  [821] 8.587038e-02 0.000000e+00 0.000000e+00 0.000000e+00 5.136347e-04
##  [826] 0.000000e+00 6.024500e-04 0.000000e+00 2.900685e-04 2.084031e-02
##  [831] 4.607627e-02 2.300466e-02 5.131981e-04 6.113751e-02 2.208983e-02
##  [836] 1.497646e-01 2.186670e-03 1.338778e-04 1.115648e-04 5.065042e-03
##  [841] 0.000000e+00 1.252244e-04 0.000000e+00 3.965438e-04 5.134199e-03
##  [846] 1.381642e-02 4.670869e-02 1.460951e-03 4.641650e-02 6.428184e-03
##  [851] 6.313395e-02 6.678633e-04 1.043536e-04 0.000000e+00 4.800267e-04
##  [856] 0.000000e+00 2.665671e-04 0.000000e+00 2.665671e-04 8.210268e-03
##  [861] 1.716692e-02 2.884257e-02 0.000000e+00 8.807379e-02 3.358746e-03
##  [866] 9.617743e-02 1.066269e-03 0.000000e+00 0.000000e+00 1.066269e-03
##  [871] 0.000000e+00 8.826991e-04 0.000000e+00 0.000000e+00 3.592095e-02
##  [876] 2.782954e-02 3.842193e-02 0.000000e+00 1.195077e-01 2.353864e-03
##  [881] 2.251128e-01 1.152413e-03 1.716359e-04 0.000000e+00 2.133189e-03
##  [886] 0.000000e+00 0.000000e+00 1.000400e-04 6.002401e-04 1.845738e-02
##  [891] 3.258804e-02 3.461385e-02 1.575630e-03 3.081232e-02 1.510604e-02
##  [896] 1.346539e-01 2.425970e-03 1.500600e-04 0.000000e+00 3.001200e-04
##  [901] 0.000000e+00 1.833757e-02 6.692545e-04 3.078570e-03 5.621737e-03
##  [906] 5.528042e-02 2.074689e-02 6.705930e-02 1.278276e-01 1.065453e-01
##  [911] 1.673136e-02 7.361799e-03 4.631241e-02 1.367956e-01 9.235711e-03
##  [916] 1.423027e-03 5.821475e-03 4.010349e-03 4.398448e-03 2.069858e-03
##  [921] 3.906856e-02 1.060802e-02 3.169470e-02 2.589909e-01 4.592497e-02
##  [926] 1.112549e-02 1.655886e-02 2.470893e-02 6.752911e-02 7.373868e-03
##  [931] 7.899175e-03 8.165439e-03 4.348984e-03 0.000000e+00 5.857815e-03
##  [936] 8.609213e-03 8.875477e-03 6.576728e-02 8.822224e-02 1.087246e-01
##  [941] 3.922961e-02 6.541227e-02 5.192154e-02 2.378628e-02 1.952605e-03
##  [946] 0.000000e+00 9.760095e-03 4.470273e-04 2.533155e-03 2.086127e-03
##  [951] 1.661451e-02 1.467740e-02 4.023245e-02 1.203993e-01 4.075399e-02
##  [956] 1.236775e-02 3.515124e-01 3.062137e-02 5.960364e-03 1.117568e-02
##  [961] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
##  [966] 0.000000e+00 1.183336e-02 2.771924e-02 6.808235e-03 2.382882e-02
##  [971] 3.566218e-03 2.125142e-01 2.193224e-01 0.000000e+00 0.000000e+00
##  [976] 0.000000e+00 3.448564e-02 5.344021e-03 0.000000e+00 0.000000e+00
##  [981] 6.596526e-03 6.028724e-02 3.022712e-02 8.625585e-02 5.435872e-02
##  [986] 1.235805e-02 2.321309e-01 8.633935e-02 0.000000e+00 1.336005e-03
##  [991] 0.000000e+00 0.000000e+00 3.012048e-03 2.454261e-03 7.809014e-02
##  [996] 8.355645e-02 2.510040e-02 2.788934e-02 2.677376e-03 5.131638e-03
## [1001] 5.232039e-02 0.000000e+00 0.000000e+00 0.000000e+00 4.462294e-03
## [1006] 0.000000e+00 1.199087e-03 2.862337e-03 1.005686e-03 5.190887e-02
## [1011] 8.532859e-02 3.589525e-02 6.896685e-02 1.249371e-02 1.469849e-03
## [1016] 5.067110e-02 1.547209e-04 1.934012e-04 0.000000e+00 1.783159e-02
## [1021] 0.000000e+00 7.216743e-04 2.465720e-03 1.683907e-03 7.613664e-02
## [1026] 7.427231e-02 2.357469e-02 1.894395e-02 1.076497e-02 1.022372e-03
## [1031] 2.964879e-02 0.000000e+00 6.615348e-04 2.405581e-04 1.040414e-02
## [1036] 0.000000e+00 3.882384e-04 2.411165e-03 1.900325e-03 9.025522e-02
## [1041] 7.023029e-02 2.952655e-02 1.215799e-02 3.821083e-03 7.560432e-04
## [1046] 6.385500e-02 2.043360e-04 6.947424e-04 1.103414e-03 3.535013e-03
## [1051] 0.000000e+00 0.000000e+00 4.893124e-03 0.000000e+00 3.966006e-02
## [1056] 5.124903e-02 2.742725e-02 1.386814e-01 7.854751e-03 1.545197e-03
## [1061] 9.271182e-02 0.000000e+00 1.030131e-03 0.000000e+00 1.712593e-02
## [1066] 0.000000e+00 1.605950e-03 9.297608e-04 1.690474e-04 7.429634e-02
## [1071] 6.909813e-02 3.892317e-02 2.157890e-01 6.677373e-03 1.436903e-03
## [1076] 6.326600e-02 0.000000e+00 4.226185e-04 0.000000e+00 1.010058e-02
## [1081] 0.000000e+00 1.311726e-01 1.367219e-03 6.433971e-04 2.469037e-02
## [1086] 2.348400e-02 2.066913e-02 3.056136e-03 2.697442e-01 2.243848e-02
## [1091] 4.246421e-02 3.860383e-03 9.650957e-04 6.514396e-03 1.367219e-03
## [1096] 3.115629e-04 3.653854e-01 3.466137e-03 0.000000e+00 1.943373e-02
## [1101] 9.409199e-02 4.295673e-02 7.789072e-04 1.794992e-01 1.686334e-02
## [1106] 1.799276e-02 3.232465e-03 2.726175e-04 2.920902e-03 3.582973e-03
## [1111] 0.000000e+00 5.835353e-03 0.000000e+00 0.000000e+00 1.714696e-02
## [1116] 7.316635e-02 5.188976e-02 3.959063e-02 1.434599e-01 9.516115e-03
## [1121] 4.084747e-02 6.436844e-02 1.091660e-01 3.546099e-02 1.705719e-03
## [1126] 0.000000e+00 3.463537e-03 0.000000e+00 0.000000e+00 2.212815e-02
## [1131] 3.944583e-02 1.269963e-02 9.813354e-03 1.287281e-01 3.771407e-02
## [1136] 3.155667e-02 5.406966e-02 1.472003e-01 2.212815e-02 7.696748e-04
## [1141] 2.277616e-03 2.910287e-03 5.187903e-03 0.000000e+00 5.061369e-04
## [1146] 1.170442e-02 7.022650e-03 1.474124e-02 8.889029e-02 8.756169e-02
## [1151] 2.619259e-02 3.473365e-02 2.290270e-01 1.581678e-03 7.085917e-03
## [1156] 0.000000e+00 1.053001e-03 0.000000e+00 0.000000e+00 1.677002e-02
## [1161] 1.260481e-01 4.278304e-02 2.031902e-02 2.620803e-02 2.301002e-03
## [1166] 2.283452e-01 1.560002e-03 5.850006e-04 0.000000e+00 9.165009e-03
## [1171] 0.000000e+00 9.691027e-04 2.793296e-03 9.691027e-04 5.535287e-02
## [1176] 5.751910e-02 5.672101e-02 1.550564e-02 3.876411e-03 5.814616e-03
## [1181] 8.636415e-02 3.078326e-03 5.700604e-04 2.280242e-04 3.363357e-03
## [1186] 0.000000e+00 2.016847e-03 8.304662e-04 0.000000e+00 2.190849e-02
## [1191] 4.745521e-02 5.172618e-02 7.205283e-02 3.480049e-03 3.045043e-03
## [1196] 1.171748e-01 0.000000e+00 3.559141e-04 0.000000e+00 3.598687e-03
## [1201] 3.318339e-04 7.005383e-04 1.659170e-03 8.111496e-04 4.955387e-02
## [1206] 3.683357e-02 2.507190e-02 3.288843e-02 2.138485e-03 1.253595e-03
## [1211] 5.379397e-02 0.000000e+00 4.793157e-04 2.580931e-04 3.428951e-03
## [1216] 0.000000e+00 5.002223e-04 2.834593e-03 3.890618e-04 7.047577e-02
## [1221] 5.002223e-02 3.973988e-02 2.256558e-02 4.668742e-03 6.669631e-04
## [1226] 5.874833e-02 0.000000e+00 1.111605e-03 2.223210e-04 4.168519e-03
## [1231] 0.000000e+00 9.962565e-04 1.358532e-03 6.339814e-04 4.754861e-02
## [1236] 5.708852e-02 5.412994e-02 5.252989e-02 1.098901e-02 4.558628e-03
## [1241] 7.260597e-02 3.622751e-04 1.449100e-03 1.811375e-04 9.087067e-03
## [1246] 0.000000e+00 1.939595e-03 1.330008e-03 6.095871e-04 3.047936e-02
## [1251] 3.724023e-02 5.458576e-02 1.060127e-01 5.929620e-03 8.312552e-04
## [1256] 8.667221e-02 1.274591e-03 0.000000e+00 0.000000e+00 4.821280e-03
## [1261] 0.000000e+00 6.539420e-04 0.000000e+00 0.000000e+00 3.269710e-03
## [1266] 8.419504e-03 1.275187e-02 6.294192e-03 1.768505e-01 4.904565e-03
## [1271] 5.497200e-02 2.787428e-02 4.621327e-01 8.092533e-03 2.705685e-02
## [1276] 0.000000e+00 1.396177e-03 2.888643e-04 0.000000e+00 3.504887e-02
## [1281] 2.946416e-02 7.703048e-03 3.658948e-03 2.774541e-01 1.068798e-02
## [1286] 5.700255e-02 1.863175e-02 1.889172e-01 2.089452e-02 1.839103e-02
## [1291] 0.000000e+00 0.000000e+00 5.984440e-04 8.065985e-04 7.831811e-03
## [1296] 7.186532e-02 9.627143e-03 3.590664e-03 2.746078e-01 3.200375e-03
## [1301] 1.731585e-01 9.418989e-03 1.177113e-01 1.079801e-02 2.594125e-02
## [1306] 0.000000e+00 1.402446e-03 1.458544e-03 0.000000e+00 3.253674e-03
## [1311] 5.441490e-03 2.973185e-02 1.065859e-03 1.436105e-02 2.827331e-02
## [1316] 1.043420e-02 2.260743e-02 6.102883e-01 4.712218e-03 4.319533e-03
## [1321] 0.000000e+00 1.269573e-03 0.000000e+00 0.000000e+00 5.645366e-02
## [1326] 6.432501e-02 4.299619e-02 9.479475e-03 4.909014e-03 4.231909e-04
## [1331] 9.250952e-02 0.000000e+00 3.385527e-04 0.000000e+00 1.294964e-01
## [1336] 0.000000e+00 5.030181e-04 7.693218e-04 2.958930e-03 3.506332e-02
## [1341] 1.296011e-02 2.683750e-02 1.352231e-02 3.964966e-03 0.000000e+00
## [1346] 4.426559e-02 0.000000e+00 3.550716e-04 0.000000e+00 2.174814e-02
## [1351] 0.000000e+00 2.595963e-03 1.050053e-03 1.195893e-03 1.003384e-01
## [1356] 3.377669e-02 2.910979e-02 2.721386e-02 3.704352e-03 4.375219e-04
## [1361] 8.657100e-02 2.916813e-04 1.458406e-04 1.750088e-04 2.727220e-02
## [1366] 0.000000e+00 3.881234e-04 9.703086e-04 2.231710e-03 1.144964e-01
## [1371] 2.192897e-02 3.386377e-02 9.314962e-03 4.269358e-03 4.851543e-04
## [1376] 9.169416e-02 0.000000e+00 0.000000e+00 0.000000e+00 6.695129e-03
## [1381] 0.000000e+00 4.978905e-04 3.930715e-04 0.000000e+00 7.035979e-02
## [1386] 7.811640e-02 4.247792e-02 2.188098e-02 1.185504e-01 6.813239e-04
## [1391] 1.878882e-02 0.000000e+00 1.572286e-04 0.000000e+00 4.457430e-02
## [1396] 0.000000e+00 1.651664e-04 3.716244e-04 0.000000e+00 2.857379e-02
## [1401] 6.152449e-02 4.777438e-02 3.621273e-02 6.730531e-03 0.000000e+00
## [1406] 7.705013e-02 0.000000e+00 0.000000e+00 0.000000e+00 4.988025e-02
## [1411] 0.000000e+00 1.674321e-03 0.000000e+00 0.000000e+00 3.348641e-03
## [1416] 1.731726e-02 1.148106e-02 4.496747e-03 8.452928e-02 3.539992e-03
## [1421] 2.741102e-02 2.478951e-01 1.020379e-01 1.349024e-02 0.000000e+00
## [1426] 0.000000e+00 4.564239e-03 0.000000e+00 2.467156e-04 6.723000e-03
## [1431] 4.305187e-02 2.627521e-02 2.602850e-02 1.440819e-01 5.551101e-03
## [1436] 2.528835e-02 4.681428e-02 1.569728e-01 3.108617e-02 1.233578e-03
## [1441] 0.000000e+00 3.508237e-03 0.000000e+00 0.000000e+00 6.233476e-02
## [1446] 2.725239e-02 2.826927e-02 1.794794e-02 2.567623e-02 2.135448e-03
## [1451] 5.979256e-02 4.865772e-02 1.431767e-01 2.237136e-03 1.967663e-02
## [1456] 0.000000e+00 1.539320e-01 0.000000e+00 0.000000e+00 1.479042e-02
## [1461] 9.683541e-02 2.433443e-02 1.674387e-03 4.329966e-01 2.701345e-02
## [1466] 3.415750e-02 1.395323e-03 0.000000e+00 0.000000e+00 0.000000e+00
## [1471] 0.000000e+00 4.247872e-01 6.215932e-04 0.000000e+00 1.611361e-02
## [1476] 1.085397e-01 6.282873e-02 4.446782e-03 2.997992e-02 1.797839e-02
## [1481] 5.326575e-02 4.590227e-03 8.606675e-04 2.199484e-03 0.000000e+00
## [1486] 2.399637e-03 2.983332e-03 3.664310e-03 3.567028e-04 1.112264e-02
## [1491] 2.662300e-02 1.462481e-02 1.177119e-02 3.599455e-02 2.426552e-01
## [1496] 2.004021e-02 1.636293e-01 1.491666e-03 9.339127e-03 5.350542e-03
## [1501] 0.000000e+00 4.320821e-04 1.188226e-03 0.000000e+00 6.151769e-02
## [1506] 5.482042e-02 3.985957e-02 7.410208e-02 2.916554e-03 8.101539e-04
## [1511] 1.165001e-01 2.538482e-03 2.160410e-04 0.000000e+00 7.021334e-04
## [1516] 2.079002e-04 6.237006e-04 9.147609e-04 4.989605e-04 5.525988e-02
## [1521] 3.650728e-02 3.958420e-02 1.305613e-02 1.355509e-02 6.237006e-04
## [1526] 1.054470e-01 0.000000e+00 4.573805e-04 5.405405e-04 2.494802e-04
## [1531] 0.000000e+00 5.908070e-04 0.000000e+00 0.000000e+00 4.572847e-02
## [1536] 5.518138e-02 4.383788e-02 7.798653e-02 7.326007e-03 0.000000e+00
## [1541] 1.827957e-01 5.908070e-04 5.908070e-04 0.000000e+00 0.000000e+00
## [1546] 0.000000e+00 7.403386e-04 0.000000e+00 0.000000e+00 1.673165e-02
## [1551] 3.301910e-02 4.940526e-02 9.283846e-02 1.441192e-02 2.911998e-03
## [1556] 9.866245e-02 8.390504e-04 0.000000e+00 0.000000e+00 0.000000e+00
## [1561] 0.000000e+00 1.372621e-03 0.000000e+00 0.000000e+00 1.354319e-02
## [1566] 1.584004e-01 2.122987e-02 1.871340e-02 2.310578e-02 8.510249e-03
## [1571] 2.357247e-01 7.320644e-03 0.000000e+00 2.745242e-04 2.882504e-03
## [1576] 0.000000e+00 0.000000e+00 1.544306e-04 1.544306e-04 2.656207e-02
## [1581] 3.548507e-01 2.053927e-02 1.636965e-03 0.000000e+00 1.111900e-03
## [1586] 1.769775e-02 3.706335e-04 0.000000e+00 0.000000e+00 1.221546e-01
## [1591] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 3.410306e-02
## [1596] 1.446084e-01 3.559974e-02 1.375526e-02 0.000000e+00 0.000000e+00
## [1601] 2.045471e-02 0.000000e+00 0.000000e+00 0.000000e+00 9.475447e-02
## [1606] 0.000000e+00 3.227815e-04 1.406405e-03 0.000000e+00 2.307887e-02
## [1611] 1.553270e-01 5.028474e-02 5.510341e-03 9.222327e-05 7.608420e-04
## [1616] 2.616835e-02 4.611164e-04 0.000000e+00 0.000000e+00 9.222327e-02
## [1621] 0.000000e+00 2.320589e-04 6.381621e-04 0.000000e+00 1.055868e-02
## [1626] 1.048326e-01 2.651273e-02 4.066833e-02 2.320589e-04 1.218309e-03
## [1631] 3.666531e-02 5.453385e-03 0.000000e+00 0.000000e+00 1.324476e-01
## [1636] 0.000000e+00 0.000000e+00 1.943918e-04 1.133952e-04 3.084350e-02
## [1641] 1.088270e-01 6.265896e-02 1.292705e-02 0.000000e+00 1.133952e-04
## [1646] 1.381802e-02 0.000000e+00 2.591891e-04 0.000000e+00 1.323484e-01
## [1651] 0.000000e+00 7.928867e-04 7.362519e-04 0.000000e+00 3.256499e-02
## [1656] 1.640709e-01 4.411848e-02 4.887580e-02 5.097129e-04 1.076060e-03
## [1661] 1.466840e-02 4.870590e-03 0.000000e+00 0.000000e+00 1.251628e-01
## [1666] 0.000000e+00 5.878767e-04 5.443303e-04 0.000000e+00 2.233931e-02
## [1671] 2.104163e-01 2.286187e-02 1.267201e-02 2.830517e-04 1.545898e-03
## [1676] 1.378244e-02 5.443303e-03 0.000000e+00 0.000000e+00 1.708762e-01
## [1681] 0.000000e+00 0.000000e+00 8.292905e-04 0.000000e+00 3.416677e-02
## [1686] 1.069785e-01 1.226106e-01 1.384915e-02 0.000000e+00 5.390389e-04
## [1691] 1.513455e-02 0.000000e+00 0.000000e+00 0.000000e+00 9.984658e-02
## [1696] 0.000000e+00 1.891985e-04 1.341589e-03 4.299966e-04 9.752322e-03
## [1701] 9.418645e-02 2.793258e-02 2.923977e-03 1.737186e-03 1.031992e-04
## [1706] 2.729618e-02 6.879945e-05 8.599931e-05 2.235982e-04 1.187650e-01
## [1711] 7.008622e-02 5.371501e-02 9.178533e-02 1.091737e-01 6.780975e-04
## [1716] 1.453066e-03 7.817495e-02 0.000000e+00 0.000000e+00 1.259324e-03
## [1721] 4.843553e-04 0.000000e+00 0.000000e+00 0.000000e+00 2.906132e-04
## [1726] 5.635918e-02 6.847190e-02 1.037589e-01 1.651234e-01 5.505781e-04
## [1731] 9.009460e-04 4.669903e-02 0.000000e+00 0.000000e+00 2.002102e-04
## [1736] 5.005256e-04 0.000000e+00 0.000000e+00 2.502628e-04 0.000000e+00
## [1741] 7.782226e-02 5.732586e-02 7.333867e-02 8.887110e-02 0.000000e+00
## [1746] 2.642114e-03 1.446757e-01 0.000000e+00 0.000000e+00 6.405124e-04
## [1751] 3.202562e-04 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1756] 4.836436e-02 7.145288e-02 8.272979e-02 1.995820e-01 5.832888e-04
## [1761] 4.034414e-03 5.619015e-02 0.000000e+00 2.916444e-04 4.860740e-04
## [1766] 6.173140e-03 0.000000e+00 0.000000e+00 0.000000e+00 1.555437e-03
## [1771] 5.266752e-02 6.316688e-02 1.199317e-01 1.349552e-01 0.000000e+00
## [1776] 6.828852e-04 1.660265e-01 0.000000e+00 0.000000e+00 0.000000e+00
## [1781] 4.268032e-04 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1786] 3.700583e-02 3.286286e-02 6.972825e-02 5.722913e-02 2.808792e-04
## [1791] 0.000000e+00 2.534232e-01 0.000000e+00 4.915385e-04 7.021979e-04
## [1796] 1.615055e-03 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1801] 6.649889e-02 5.943278e-02 1.195431e-01 9.224664e-02 6.775724e-04
## [1806] 2.032717e-03 1.552609e-01 4.839803e-04 0.000000e+00 0.000000e+00
## [1811] 8.711645e-04 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1816] 4.025012e-02 2.902776e-02 1.621007e-01 4.707155e-02 1.430300e-03
## [1821] 1.778707e-03 1.132871e-01 1.466975e-04 1.466975e-04 0.000000e+00
## [1826] 9.718708e-04 0.000000e+00 1.466975e-04 0.000000e+00 3.777460e-03
## [1831] 4.963643e-02 2.349790e-02 1.012629e-01 2.537313e-02 1.033295e-03
## [1836] 9.567547e-04 9.950249e-02 0.000000e+00 8.802143e-04 1.530808e-04
## [1841] 2.449292e-03 0.000000e+00 0.000000e+00 0.000000e+00 1.561424e-02
## [1846] 4.929510e-02 2.562638e-02 1.746809e-01 3.961554e-02 1.509034e-03
## [1851] 1.359490e-03 8.825809e-02 8.428838e-04 0.000000e+00 0.000000e+00
## [1856] 0.000000e+00 0.000000e+00 1.359490e-04 0.000000e+00 8.020991e-04
## [1861] 4.238473e-02 5.537038e-02 1.083224e-01 1.327060e-01 1.261763e-04
## [1866] 9.778666e-04 1.058830e-01 3.469849e-04 3.680143e-04 3.154408e-04
## [1871] 7.149992e-04 0.000000e+00 0.000000e+00 1.156616e-04 5.677935e-04
## [1876] 4.639356e-02 3.493087e-02 7.829113e-02 8.094128e-02 0.000000e+00
## [1881] 6.066605e-04 1.410326e-01 0.000000e+00 5.428015e-04 8.940260e-04
## [1886] 7.343785e-04 0.000000e+00 0.000000e+00 0.000000e+00 4.789425e-04
## [1891] 8.715896e-02 8.089291e-02 4.181715e-02 6.714242e-02 3.481137e-04
## [1896] 2.828423e-03 5.726470e-02 0.000000e+00 2.610852e-04 1.740568e-04
## [1901] 1.566511e-03 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1906] 7.014289e-02 5.407752e-02 3.402542e-02 8.305045e-02 3.947265e-04
## [1911] 2.881503e-03 2.992027e-02 0.000000e+00 0.000000e+00 3.157812e-04
## [1916] 9.473435e-04 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1921] 8.004522e-02 7.131148e-02 3.414358e-02 8.145845e-02 5.935557e-04
## [1926] 3.193895e-03 5.551159e-02 1.130582e-04 4.522329e-04 3.109101e-04
## [1931] 1.215376e-03 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [1936] 6.725265e-02 8.737640e-02 2.908691e-02 8.575724e-02 7.517493e-04
## [1941] 3.006997e-03 5.967733e-02 0.000000e+00 2.139594e-03 3.469612e-04
## [1946] 1.272191e-03 2.313075e-04 0.000000e+00 0.000000e+00 0.000000e+00
## [1951] 7.137918e-02 8.055903e-02 2.787790e-02 5.137183e-02 9.121000e-04
## [1956] 3.795513e-03 7.596911e-02 1.618242e-04 2.103715e-03 1.515263e-03
## [1961] 2.089003e-03 1.618242e-04 5.884516e-05 2.353807e-04 5.884516e-05
## [1966] 9.098369e-02 9.131191e-02 3.800834e-02 6.925526e-02 4.266912e-04
## [1971] 6.269078e-03 5.287688e-02 0.000000e+00 9.518495e-04 2.297568e-04
## [1976] 4.266912e-03 0.000000e+00 0.000000e+00 0.000000e+00 9.518495e-04
## [1981] 6.560430e-02 1.952263e-02 6.077539e-02 3.949365e-02 1.241722e-03
## [1986] 3.897627e-03 4.473648e-02 1.172737e-03 9.312914e-04 0.000000e+00
## [1991] 1.103753e-03 0.000000e+00 0.000000e+00 0.000000e+00 6.208609e-04
## [1996] 3.449282e-02 2.433757e-02 3.355901e-02 3.256683e-02 2.918174e-04
## [2001] 1.925995e-03 6.752656e-02 0.000000e+00 3.501809e-04 2.334540e-04
## [2006] 2.101086e-03 0.000000e+00 3.501809e-04 0.000000e+00 1.050543e-03
## [2011] 9.669112e-02 4.136682e-02 7.752458e-02 5.289027e-02 2.151827e-03
## [2016] 7.701895e-03 4.090824e-02 1.058276e-03 3.880344e-04 2.469310e-04
## [2021] 1.575655e-03 2.939655e-04 0.000000e+00 1.646207e-04 1.881379e-04
## [2026] 1.281413e-01 3.982735e-02 1.070310e-01 5.629883e-02 4.504444e-03
## [2031] 9.103011e-03 3.680198e-02 3.536325e-03 7.260895e-04 4.033830e-04
## [2036] 1.976577e-03 2.823681e-04 0.000000e+00 1.075688e-04 3.764908e-04
## [2041] 1.052812e-01 4.093507e-02 1.121265e-01 6.852175e-02 4.415238e-03
## [2046] 7.290276e-03 3.877879e-02 1.574426e-03 5.476264e-04 0.000000e+00
## [2051] 1.471746e-03 0.000000e+00 0.000000e+00 1.369066e-04 0.000000e+00
## [2056] 9.886209e-02 1.381311e-02 5.545526e-02 4.594227e-02 5.760532e-03
## [2061] 1.407679e-02 1.797124e-02 4.462384e-04 1.622685e-04 3.042535e-04
## [2066] 6.977546e-03 6.085069e-04 0.000000e+00 1.217014e-04 0.000000e+00
## [2071] 7.439798e-02 2.718946e-02 5.558088e-02 6.262693e-02 3.108551e-03
## [2076] 8.123679e-03 4.339537e-02 7.874995e-04 1.409210e-03 8.289468e-04
## [2081] 1.119078e-03 0.000000e+00 0.000000e+00 2.072367e-04 2.486840e-04
## [2086] 4.817964e-02 3.102972e-02 4.636291e-02 9.367052e-02 2.180074e-03
## [2091] 1.075503e-02 9.795800e-02 3.633457e-04 2.906765e-04 7.266914e-04
## [2096] 3.124773e-03 0.000000e+00 0.000000e+00 0.000000e+00 7.266914e-04
## [2101] 8.163207e-02 4.119882e-02 5.656686e-02 6.748327e-02 1.275944e-03
## [2106] 7.003516e-03 8.381536e-02 0.000000e+00 4.253147e-04 8.789838e-04
## [2111] 2.296700e-03 1.134173e-04 0.000000e+00 1.984802e-04 4.536690e-04
## [2116] 9.102289e-02 3.458512e-02 9.206009e-02 8.040057e-02 1.072961e-03
## [2121] 3.791130e-03 4.163090e-02 0.000000e+00 0.000000e+00 0.000000e+00
## [2126] 1.537911e-03 6.437768e-04 0.000000e+00 2.145923e-04 1.680973e-03
## [2131] 7.650339e-02 4.742008e-02 6.169871e-02 8.239520e-02 4.193408e-03
## [2136] 1.334676e-02 5.691912e-02 1.052109e-04 1.367742e-03 1.052109e-04
## [2141] 3.081178e-03 1.803616e-04 7.515068e-05 2.254520e-04 3.006027e-04
## [2146] 7.757679e-02 5.279400e-02 8.787922e-02 9.837155e-02 2.516261e-03
## [2151] 1.153682e-02 5.284148e-02 0.000000e+00 2.848597e-04 3.323363e-04
## [2156] 6.456820e-03 2.373831e-04 0.000000e+00 0.000000e+00 1.661682e-03
## [2161] 7.756554e-02 2.517530e-02 4.927869e-02 6.339234e-02 2.039631e-03
## [2166] 3.468862e-03 4.886183e-02 1.771651e-03 4.168590e-04 2.977564e-04
## [2171] 1.994968e-03 8.932692e-05 0.000000e+00 1.191026e-04 6.104006e-04
## [2176] 6.030879e-02 1.616345e-02 5.014145e-02 8.331805e-02 1.197293e-03
## [2181] 2.761497e-03 5.318296e-02 3.862234e-04 7.338245e-04 3.476011e-04
## [2186] 2.133884e-03 1.351782e-04 0.000000e+00 8.690027e-05 6.565798e-04
## [2191] 5.202717e-02 5.186252e-02 5.326199e-02 1.115867e-01 2.963573e-03
## [2196] 4.939288e-03 8.121013e-02 5.885985e-03 1.481786e-03 1.646429e-04
## [2201] 2.675448e-03 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
## [2206] 7.936374e-02 5.664030e-02 8.786400e-02 1.176570e-01 2.608988e-03
## [2211] 9.678505e-03 4.654099e-02 3.366437e-04 0.000000e+00 8.416092e-04
## [2216] 9.426023e-03 0.000000e+00 0.000000e+00 0.000000e+00 1.935701e-03
## [2221] 4.697200e-02 2.092468e-02 2.292164e-02 5.035815e-02 1.458650e-02
## [2226] 1.419579e-02 9.342305e-02 1.415238e-02 6.989364e-03 1.215542e-03
## [2231] 1.975255e-02 1.736488e-04 9.984806e-04 1.041893e-03 8.682440e-04
## [2236] 5.416478e-01 5.302503e-01 5.143989e-01 2.960190e-01 5.114324e-01
## [2241] 4.011679e-01 5.511148e-01 6.116001e-01 5.879770e-01 5.699077e-01
## [2246] 6.307963e-01 5.850985e-01 4.142786e-01 5.423342e-01 5.845656e-01
## [2251] 2.466405e-01 4.687474e-01 2.778176e-01 2.523225e-01 3.866625e-01
## [2256] 3.208959e-01 1.463034e-01 6.023522e-01 4.848586e-01 6.026325e-01
## [2261] 3.466451e-01 4.955733e-01 5.740210e-01 5.100161e-01 6.634094e-01
## [2266] 6.260574e-01 5.024889e-01 4.431644e-01 4.495186e-01 6.180158e-01
## [2271] 5.598283e-01 7.109864e-01 2.507532e-01 1.943574e-01 1.405985e-01
## [2276] 4.838036e-01 5.417191e-01 8.227821e-01 8.231488e-01 7.733392e-01
## [2281] 5.824601e-01 6.150202e-01 7.437008e-01 7.548036e-01 2.343026e-01
## [2286] 6.230306e-02 2.350913e-01 4.718890e-01 5.569126e-01 7.519612e-01
## [2291] 6.681839e-01 8.151271e-01 7.555046e-01 5.465133e-01 7.286164e-01
## [2296] 3.783965e-01 4.686934e-01 5.112275e-01 3.408583e-01 4.944075e-01
## [2301] 3.902806e-01 7.153057e-01 6.700190e-01 7.494587e-01 7.191606e-01
## [2306] 6.178213e-01 5.172851e-01 4.475631e-01 2.492114e-01 4.078463e-01
## [2311] 4.902829e-01 4.805770e-01 5.248625e-01 7.068749e-01 6.763554e-01
## [2316] 7.907971e-01 7.438862e-01 6.864811e-01 6.682738e-01 2.067274e-01
## [2321] 3.304607e-01 2.908438e-01 2.626501e-01 5.977994e-01 8.370517e-01
## [2326] 6.861218e-01 7.136619e-01 6.035219e-01 6.917169e-01 4.827784e-01
## [2331] 4.820823e-01 5.593350e-01 2.128705e-01 2.737879e-01 4.479538e-01
## [2336] 6.443964e-01 7.324740e-01 5.853716e-01 6.904398e-01 5.089220e-01
## [2341] 4.547673e-01 6.567244e-01 6.443640e-01 6.405407e-01 6.378967e-01
## [2346] 5.625531e-01 5.386474e-01 6.060455e-01 7.149639e-01 5.928994e-01
## [2351] 5.571850e-01 5.543635e-01 5.285568e-01 4.621426e-01 5.466610e-01
## [2356] 5.029523e-01 5.998643e-01 6.796403e-01 6.178746e-01 5.518006e-01
## [2361] 6.151537e-01 6.605457e-01 7.242441e-01 6.716507e-01 6.631007e-01
## [2366] 6.817506e-01 6.444678e-01 7.608996e-01 8.012140e-01 6.768496e-01
## [2371] 6.108833e-01 6.189205e-01 7.394982e-01 7.209765e-01 6.646319e-01
## [2376] 6.566576e-01 6.513591e-01 6.520824e-01 6.075108e-01 7.249177e-01
## [2381] 7.284463e-01 6.319407e-01 5.871065e-01 6.914261e-01
# you see all the observations with list()