![[logo] US EPA](../gif/logo_epaseal.gif){kind=logo}
Example Analysis With EPA Region 10, EMAP West Wadeable Stream Data
The material for this example was prepared and presented at a workshop titled: Aquatic Monitoring Design and Analysis Workshop with an emphasis on EMAP Western Pilot, October 18-19, 2005. Instructors: Tony Olsen & John Van Sickle. The data files, analytical procedures, and analysis results are for illustrative purposes and should not be interpreted or cited as final results for the EMAP Western Pilot. The workshop included several sessions and the following are included:
- Population Extent Estimation
- Estimating CDFs, densities, percentiles, means
- Comparing two CDFs
Population Extent Estimation
| Purpose | Extent estimation for Region 10 and States. Data files include final design information and only have one record for each site from the probability sample. Originial design weights and adjusted weights are included in the data file. Extent Staus estimates are based on all evaluated sites, i.e. sites assigned to target, non-target status. Used Extent estimates are based on all sites with status assigned during field visits. |
| Input file | status_wgt_all7.csv (125 KB) |
| R Script file | extent_estimates.R (4 KB) |
| Original weights adjusted to sample size of 852
Sum adjwts= 243,622 km |
Extent estimates for site evaluated as to Target (n=646) and NonTarget (n=206)
for Region 10 and States:
|
| Results R Console [extent_R.doc (MS WORD 40 KB)] | Extent estimates based on Site Status: T (451), NT (172) and LD (n=114), NT
(n=172), PB (n=43), TS (n=294) Total (n=623)
|
All Sample Sites with Determination of Indicator Category (Good, Fair, Poor)
| Purpose | Condition class estimation for Region 10 and States useful for assessment document. Estimates based on Site Status: T (451), NT (172) and LD (n=114), NT (n=172), PB (n=43), TS (n=294) Total (n=623) |
| Input File | condclass_wgt3.csv (47 KB) |
| R Script | assessment-class.R (3 KB) |
| Sum adjwts= 243,622 km | Estimates based on Indicator Categories (Good, Fair, Poor) for biological and
chemical indicators.
|
| Results R Console [assessment_R.doc (MS WORD, 36KB)] | condition_assessments.csv (9 KB) |
Estimating CDFs
| Purpose | Estimate the cdf's for Chemical Indicators pH and Nitrate for Region and by State |
| Input file | chem_wgt3.csv (41 KB) |
| R Script file | cdf_chem.R (4 KB) , plus
R scripts: cont.ash.r, ash1.wgt.r, cont.plot.r, plot.cont.r, plotcdf.r These R scripts must be in the current directory selected for R. |
| Results R Console [chem_cdf_R.doc (MS WORD, 5 pp, 53 KB)] |
|
Comparing two CDFs
| Purpose | Compare distribution and cdf's for two subpopulations. Chemical Indicators, pH and NO3, for level 3 ecoregions Mountain and Xeric within Region 10 |
| Input file | chem_wgt3.csv (41 KB) |
| R Script file | chem_cdf_ compare.R (7
KB), plus R scripts: cont.ash.r,
ash1.wgt.r, cont.plot.r, plot.cont.r, plotcdf.r These R scripts must be in the current directory selected for R. |
| Results R Console [chem_cdf_comp_R.doc (MS WORD, 8pp, 73KB) |
|