*** Summary statistics for the S-ESTREND data. *** Constituent: SO4 Station: 19 112 observations from 12/08/1982 to 07/06/1998 106 less than values and 0 zero values ( 94.6 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 10 10 10 10 10 10 16 Station: 20 127 observations from 11/16/1979 to 07/06/1998 97 less than values and 0 zero values ( 76.4 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 1.2 10 10 10 10 11.4 22 Station: 21 151 observations from 06/23/1967 to 07/06/1998 70 less than values and 0 zero values ( 46.4 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 3 6 10 10 11 20 39 Station: 22 127 observations from 09/27/1967 to 07/06/1998 74 less than values and 0 zero values ( 58.3 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 4 10 10 10 11 16 24 Station: 23 121 observations from 11/16/1979 to 07/06/1998 2 less than values and 4 zero values ( 5 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 0 22 33 44 80 80 94 Station: 27 107 observations from 11/16/1979 to 07/06/1998 13 less than values and 0 zero values ( 12.1 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 6 10 14 18 21 26 40.2 Station: 32 58 observations from 11/16/1979 to 06/05/1991 51 less than values and 0 zero values ( 87.9 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 2.6 10 10 10 10 10 12 Station: 38 165 observations from 06/23/1967 to 01/16/2002 98 less than values and 1 zero values ( 60 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 0 5.24 6.4 10 10 10 19 Station: 46 108 observations from 12/08/1982 to 07/06/1998 105 less than values and 0 zero values ( 97.2 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 10 10 10 10 10 10 23 Station: 47 117 observations from 11/16/1979 to 07/06/1998 80 less than values and 0 zero values ( 68.4 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 2.2 10 10 10 10 14.4 27 Station: 53 58 observations from 11/16/1979 to 05/26/1998 44 less than values and 1 zero values ( 77.6 %) Statistics of the raw data 0% 10% 25% 50% 75% 90% 100% 0 8.21 10 10 10 10 15 Summary report saved in obed.all.SO4.summary. *** S-ESTREND Processing status. *** SO4 19 data OK 20 data OK 21 data OK 22 data OK 23 data OK 27 data OK 32 data OK 38 data OK 46 data OK 47 data OK 53 data OK *** S-ESTREND Seasonal Analysis *** Constituent: SO4 Station: 19 Record length -- 15.6 years, samples 12/08/1982 to 07/06/1998 Period of seasonal analysis: 1983 to 1998 Beginning-ending record number of years = 6 Middle record number of years = 10 All possible within-season comparisons can be made for the 2-season definition. Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 100 83.33 100 100 MI 100 90.00 100 80 Station: 20 Record length -- 18.6 years, samples 11/16/1979 to 07/06/1998 Period of seasonal analysis: 1980 to 1998 Beginning-ending record number of years = 8 Middle record number of years = 11 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 75 75 MI 100 100 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 75 62.50 75 75.00 MI 100 90.91 100 81.82 Station: 21 Record length -- 31 years, samples 06/23/1967 to 07/06/1998 Period of seasonal analysis: 1967 to 1998 Beginning-ending record number of years = 14 Middle record number of years = 18 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 64.29 71.43 MI 66.67 61.11 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 64.29 50.00 64.29 64.29 MI 55.56 61.11 61.11 55.56 Station: 22 Record length -- 30.8 years, samples 09/27/1967 to 07/06/1998 Period of seasonal analysis: 1968 to 1998 Beginning-ending record number of years = 12 Middle record number of years = 19 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 66.67 66.67 MI 52.63 52.63 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 66.67 50.00 58.33 58.33 MI 52.63 47.37 52.63 42.11 Station: 23 Record length -- 18.6 years, samples 11/16/1979 to 07/06/1998 Period of seasonal analysis: 1980 to 1998 Beginning-ending record number of years = 8 Middle record number of years = 11 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 75 75 MI 100 100 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 75 62.50 62.5 75 MI 100 90.91 100.0 100 Station: 27 Record length -- 18.6 years, samples 11/16/1979 to 07/06/1998 Period of seasonal analysis: 1980 to 1998 Beginning-ending record number of years = 8 Middle record number of years = 11 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 75.00 75.00 MI 81.82 90.91 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 75.00 62.50 75.00 75.00 MI 72.73 72.73 81.82 72.73 Station: 32 Record length -- 11.6 years, samples 11/16/1979 to 06/05/1991 Period of seasonal analysis: 1980 to 1991 Beginning-ending record number of years = 6 Middle record number of years = 6 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 50 66.67 MI 100 100.00 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 50 16.67 66.67 16.67 MI 100 100.00 100.00 83.33 Station: 38 Record length -- 34.6 years, samples 06/23/1967 to 01/16/2002 Period of seasonal analysis: 1967 to 2001 Beginning-ending record number of years = 14 Middle record number of years = 21 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 42.86 50.0 MI 61.90 61.9 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 35.71 42.86 50.0 42.86 MI 61.90 57.14 61.9 57.14 Station: 46 Record length -- 15.6 years, samples 12/08/1982 to 07/06/1998 Period of seasonal analysis: 1983 to 1998 Beginning-ending record number of years = 6 Middle record number of years = 10 All possible within-season comparisons can be made for the 2-season definition. Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 100 83.33 100 100 MI 100 80.00 100 70 Station: 47 Record length -- 18.6 years, samples 11/16/1979 to 07/06/1998 Period of seasonal analysis: 1980 to 1998 Beginning-ending record number of years = 8 Middle record number of years = 11 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 75 75 MI 100 100 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 75 62.50 75 75.00 MI 100 90.91 100 81.82 Station: 53 Record length -- 18.5 years, samples 11/16/1979 to 05/26/1998 Period of seasonal analysis: 1980 to 1998 Beginning-ending record number of years = 8 Middle record number of years = 11 Percent of possible within-season comparisons for the 2-season definition. 1 2 BE 25.00 50.00 MI 63.64 63.64 Percent of possible within-season comparisons for the 4-season definition. 1 2 3 4 BE 25.00 25.00 37.50 37.50 MI 63.64 45.45 63.64 36.36 Percent possible comparisons within a season adequate to define trend: 50 Percent of seasons satisfying above criterion to qualify for seasonal definition: 80 *** S-ESTREND Seasonal selections *** Constituent: SO4 Station BE MI select 1 19 4 4 4 2 20 4 4 4 3 21 4 4 4 4 22 4 2 4 5 23 4 4 4 6 27 4 4 4 7 32 2 4 2 8 38 1 4 1 9 46 4 4 4 10 47 4 4 4 11 53 1 2 1 *** S-ESTREND Flow-adjustment models *** Constituent: SO4 Station: 19 Percentage of censored values ( 94.6428571428571 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 20 Percentage of censored values ( 76.3779527559055 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 21 Percentage of censored values ( 46.3576158940397 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 22 Percentage of censored values ( 58.2677165354331 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 23 No flow adjustment Station: 27 Percentage of censored values ( 12.1495327102804 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 32 Percentage of censored values ( 87.9310344827586 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 38 Percentage of censored values ( 60 ) exceeds 5 1 censoring levels were found [1] 10 Zero values (nondetected values with unrecorded reporting limit) were found No flow adjustment Station: 46 Percentage of censored values ( 97.2222222222222 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 47 Percentage of censored values ( 68.3760683760684 ) exceeds 5 1 censoring levels were found [1] 10 No flow adjustment Station: 53 Percentage of censored values ( 77.5862068965517 ) exceeds 5 1 censoring levels were found [1] 10 Zero values (nondetected values with unrecorded reporting limit) were found No flow adjustment *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 19 $tobit: [1] "Fewer than 10 detections in record" Station: 20 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 0.870 0.3886 2.24 0.0251100 Year 0.090 0.0293 3.07 0.0021426 Log(scale) -0.655 0.1513 -4.33 0.0000148 Scale= 0.519 Gaussian distribution Loglik(model)= -36.8 Loglik(intercept only)= -48.5 Chisq= 23.33 on 1 degrees of freedom, p= 1.4e-006 Number of Newton-Raphson Iterations: 5 n= 118 Correlation of Coefficients: (Intercept) Year Year -0.970 Log(scale) -0.656 0.558 Station: 21 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 1.532 0.1607 9.54 1.49e-021 Year 0.081 0.0155 5.23 1.68e-007 Log(scale) -0.550 0.1130 -4.87 1.11e-006 Scale= 0.577 Gaussian distribution Loglik(model)= -71.5 Loglik(intercept only)= -87 Chisq= 31 on 1 degrees of freedom, p= 2.6e-008 Number of Newton-Raphson Iterations: 3 n= 114 Correlation of Coefficients: (Intercept) Year Year -0.883 Log(scale) -0.529 0.351 Station: 22 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 1.7579 0.127 13.87 8.99e-044 Year 0.0474 0.011 4.29 1.77e-005 Log(scale) -0.9385 0.136 -6.88 6.00e-012 Scale= 0.391 Gaussian distribution Loglik(model)= -48.3 Loglik(intercept only)= -59.7 Chisq= 22.83 on 1 degrees of freedom, p= 1.8e-006 Number of Newton-Raphson Iterations: 3 n= 109 Correlation of Coefficients: (Intercept) Year Year -0.881 Log(scale) -0.631 0.417 Station: 23 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 3.75626 0.1062 35.373 4.48e-274 Year 0.00302 0.0123 0.246 8.05e-001 Log(scale) -0.48815 0.0694 -7.038 1.95e-012 Scale= 0.614 Gaussian distribution Loglik(model)= -109.7 Loglik(intercept only)= -109.7 Chisq= 0.06 on 1 degrees of freedom, p= 0.81 Number of Newton-Raphson Iterations: 2 n= 114 Correlation of Coefficients: (Intercept) Year Year -0.840 Log(scale) -0.010 -0.002 Station: 27 $tobit: [1] "Too sparse observations in ending fifth, latest: 12/03/1997" Station: 32 $tobit: [1] "Fewer than 10 detections in record" Station: 38 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 1.6895 0.1340 12.61 1.92e-036 Year 0.0135 0.0112 1.20 2.30e-001 Log(scale) -0.9702 0.1161 -8.35 6.62e-017 Scale= 0.379 Gaussian distribution Loglik(model)= -33.1 Loglik(intercept only)= -33.9 Chisq= 1.59 on 1 degrees of freedom, p= 0.21 Number of Newton-Raphson Iterations: 4 n= 125 Correlation of Coefficients: (Intercept) Year Year -0.880 Log(scale) -0.547 0.316 Station: 46 $tobit: [1] "Fewer than 10 detections in record" Station: 47 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 1.3523 0.218 6.19 5.87e-010 Year 0.0738 0.018 4.09 4.27e-005 Log(scale) -0.7107 0.146 -4.86 1.17e-006 Scale= 0.491 Gaussian distribution Loglik(model)= -45.4 Loglik(intercept only)= -58.5 Chisq= 26.21 on 1 degrees of freedom, p= 3.1e-007 Number of Newton-Raphson Iterations: 4 n= 107 Correlation of Coefficients: (Intercept) Year Year -0.926 Log(scale) -0.663 0.502 Station: 53 $tobit: [1] "Fewer than 10 detections in record" *** S-ESTREND Trend Summary *** Constituent: SO4 STAID Type Begyear Numyears Numseasons N Trend Trend.pct P.value 1 19 NA NA NA NA NA NA NA NA 2 20 Tobit 1983 15 NA 118 NA 9.4227105 2.142637e-003 3 21 Tobit 1983 15 NA 114 NA 8.4409330 1.681465e-007 4 22 Tobit 1983 15 NA 109 NA 4.8558391 1.770989e-005 5 23 Tobit 1983 15 NA 114 NA 0.3025809 8.054184e-001 6 27 NA NA NA NA NA NA NA NA 7 32 NA NA NA NA NA NA NA NA 8 38 Tobit 1983 15 NA 125 NA 1.3551404 2.302720e-001 9 46 NA NA NA NA NA NA NA NA 10 47 Tobit 1983 15 NA 107 NA 7.6618771 4.272036e-005 11 53 NA NA NA NA NA NA NA NA Sig 1 NA 2 up 3 up 4 up 5 none 6 NA 7 NA 8 none 9 NA 10 up 11 NA *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 27 $tobit: [1] "Too sparse observations in ending fifth, latest: 12/18/1996" *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 27 $tobit: [1] "Too sparse observations in ending fifth, latest: 12/12/1995" *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 27 $tobit: [1] "Too sparse observations in ending fifth, latest: 12/06/1993" *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 27 Call: survReg(formula = formula.temp, data = data, na.action = na.exclude, dist = "gaussian", model = T) Value Std. Error z p (Intercept) 2.8172 0.0923 30.51 1.71e-204 Year -0.0222 0.0187 -1.19 2.36e-001 Log(scale) -1.0441 0.1050 -9.94 2.75e-023 Scale= 0.352 Gaussian distribution Loglik(model)= -27.6 Loglik(intercept only)= -28.3 Chisq= 1.38 on 1 degrees of freedom, p= 0.24 Number of Newton-Raphson Iterations: 2 n= 57 Correlation of Coefficients: (Intercept) Year Year -0.858 Log(scale) -0.047 0.008 *** S-ESTREND Trend Summary *** Station: 27 Constituent Type Begyear Numyears Numseasons N Trend Trend.pct P.value 1 SO4 Tobit 1982 10 NA 57 NA -2.194902 0.2358117 Sig 1 none *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 32 $tobit: [1] "Fewer than 10 detections in record" *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 46 $tobit: [1] "Fewer than 10 detections in record" *** S-ESTREND Trend Test results *** Constituent: SO4 Station: 53 $tobit: [1] "Fewer than 10 detections in record"