Due in large part to my nonexistent statistics vocabulary, I'm not sure
what the name is for the type of analysis I'm trying to do.  I have,
however, started to cobble something together (below) that is undoubtedly
reinventing the wheel (i.e., someone has probably already written something
to do this type of analysis, especially within R).  I've compiled 1-year's
worth of data (which is a subset of a much larger dataset, but extends from
10/1/97 to 9/30/98), that includes daily mean and standard deviations for
flow (Q, Qsd) and specific conductance (EC, ECsd).  My goal is to calculate
total load for the year and include a 95% CI on that value.  The load is
calculated by multiplying each day's Q and EC and summing over the year.
 Assuming normally-distributed errors, the calculated load for this data
set is 565,240 tons/yr.  The last line of code below is where I have left
off.  For each day in the full 25 year data set, I want to generate ~500
realizations of each day's load (Q*SC), sum over the daily realization
values to arrive at 500 realizations of total annual load.  So, my question
is, what is the package that already does this? or, alternatively, could
someone help me complete the code below based on what I just described?
-Eric


a <- read.csv(textConnection("date,Q,QAcc,Qsd,SC,SCAcc,SCsd
10/1/1997,654,0.15,50.05,1020.00,0.15,78.06
10/2/1997,592,0.15,45.31,1070.00,0.15,81.89
10/3/1997,560,0.15,42.86,1100.00,0.15,84.18
10/4/1997,541,0.15,41.40,1110.00,0.15,84.95
10/5/1997,545,0.15,41.71,1100.00,0.15,84.18
10/6/1997,528,0.15,40.41,1110.00,0.15,84.95
10/7/1997,452,0.15,34.59,1120.00,0.15,85.71
10/8/1997,458,0.15,35.05,NA,0.15,NA
10/9/1997,524,0.15,40.10,NA,0.15,NA
10/10/1997,580,0.15,44.39,NA,0.15,NA
10/11/1997,600,0.15,45.92,NA,0.15,NA
10/12/1997,716,0.15,54.80,NA,0.15,NA
10/13/1997,731,0.15,55.94,NA,0.15,NA
10/14/1997,550,0.15,42.09,NA,0.15,NA
10/15/1997,400,0.15,30.61,NA,0.15,NA
10/16/1997,413,0.15,31.61,NA,0.15,NA
10/17/1997,403,0.15,30.84,1050.00,0.15,80.36
10/18/1997,330,0.15,25.26,1100.00,0.15,84.18
10/19/1997,222,0.15,16.99,1210.00,0.15,92.60
10/20/1997,198,0.15,15.15,NA,0.15,NA
10/21/1997,197,0.15,15.08,NA,0.15,NA
10/22/1997,204,0.15,15.61,1400.00,0.15,107.14
10/23/1997,204,0.15,15.61,1410.00,0.15,107.91
10/24/1997,229,0.15,17.53,1390.00,0.15,106.38
10/25/1997,162,0.15,12.40,1220.00,0.15,93.37
10/26/1997,313,0.15,23.95,1230.00,0.15,94.13
10/27/1997,784,0.15,60.00,987.00,0.15,75.54
10/28/1997,691,0.15,52.88,1010.00,0.15,77.30
10/29/1997,939,0.15,71.86,974.00,0.15,74.54
10/30/1997,953,0.15,72.93,966.00,0.15,73.93
10/31/1997,1020,0.15,78.06,NA,0.15,NA
11/1/1997,1090,0.15,83.42,1000.00,0.15,76.53
11/2/1997,1070,0.15,81.89,991.00,0.15,75.84
11/3/1997,963,0.15,73.70,958.00,0.15,73.32
11/4/1997,896,0.15,68.57,984.00,0.15,75.31
11/5/1997,823,0.15,62.98,1020.00,0.15,78.06
11/6/1997,749,0.15,57.32,1070.00,0.15,81.89
11/7/1997,698,0.15,53.42,1100.00,0.15,84.18
11/8/1997,704,0.15,53.88,1140.00,0.15,87.24
11/9/1997,752,0.15,57.55,1100.00,0.15,84.18
11/10/1997,793,0.15,60.69,1080.00,0.15,82.65
11/11/1997,821,0.15,62.83,1070.00,0.15,81.89
11/12/1997,822,0.15,62.91,1080.00,0.15,82.65
11/13/1997,931,0.15,71.25,1040.00,0.15,79.59
11/14/1997,923,0.15,70.64,1020.00,0.15,78.06
11/15/1997,932,0.15,71.33,1010.00,0.15,77.30
11/16/1997,924,0.15,70.71,987.00,0.15,75.54
11/17/1997,928,0.15,71.02,984.00,0.15,75.31
11/18/1997,947,0.15,72.47,982.00,0.15,75.15
11/19/1997,786,0.15,60.15,1020.00,0.15,78.06
11/20/1997,729,0.15,55.79,1050.00,0.15,80.36
11/21/1997,771,0.15,59.01,1040.00,0.15,79.59
11/22/1997,745,0.15,57.02,1060.00,0.15,81.12
11/23/1997,463,0.15,35.43,1160.00,0.15,88.78
11/24/1997,397,0.15,30.38,1220.00,0.15,93.37
11/25/1997,848,0.15,64.90,1030.00,0.15,78.83
11/26/1997,878,0.15,67.19,971.00,0.15,74.31
11/27/1997,857,0.15,65.59,969.00,0.15,74.16
11/28/1997,939,0.15,71.86,965.00,0.15,73.85
11/29/1997,1130,0.15,86.48,958.00,0.15,73.32
11/30/1997,1040,0.15,79.59,1080.00,0.15,82.65
12/1/1997,956,0.15,73.16,1060.00,0.15,81.12
12/2/1997,1020,0.15,78.06,1040.00,0.15,79.59
12/3/1997,1150,0.15,88.01,944.00,0.15,72.24
12/4/1997,1120,0.15,85.71,933.00,0.15,71.40
12/5/1997,1040,0.15,79.59,928.00,0.15,71.02
12/6/1997,816,0.15,62.45,1020.00,0.15,78.06
12/7/1997,797,0.15,60.99,1060.00,0.15,81.12
12/8/1997,835,0.15,63.90,1070.00,0.15,81.89
12/9/1997,980,0.15,75.00,1020.00,0.15,78.06
12/10/1997,971,0.15,74.31,1020.00,0.15,78.06
12/11/1997,937,0.15,71.71,1030.00,0.15,78.83
12/12/1997,880,0.15,67.35,1010.00,0.15,77.30
12/13/1997,895,0.15,68.49,1050.00,0.15,80.36
12/14/1997,914,0.15,69.95,1030.00,0.15,78.83
12/15/1997,966,0.15,73.93,1020.00,0.15,78.06
12/16/1997,913,0.15,69.87,1010.00,0.15,77.30
12/17/1997,730,0.15,55.87,1160.00,0.15,88.78
12/18/1997,610,0.15,46.68,1290.00,0.15,98.72
12/19/1997,591,0.15,45.23,1350.00,0.15,103.32
12/20/1997,589,0.15,45.08,1360.00,0.15,104.08
12/21/1997,560,0.15,42.86,1370.00,0.15,104.85
12/22/1997,527,0.15,40.33,1380.00,0.15,105.61
12/23/1997,532,0.15,40.71,1390.00,0.15,106.38
12/24/1997,535,0.15,40.94,1400.00,0.15,107.14
12/25/1997,558,0.15,42.70,1390.00,0.15,106.38
12/26/1997,535,0.15,40.94,1380.00,0.15,105.61
12/27/1997,504,0.15,38.57,1410.00,0.15,107.91
12/28/1997,476,0.15,36.43,1400.00,0.15,107.14
12/29/1997,484,0.15,37.04,1410.00,0.15,107.91
12/30/1997,519,0.15,39.72,1400.00,0.15,107.14
12/31/1997,539,0.15,41.25,1410.00,0.15,107.91
1/1/1998,552,0.15,42.24,1400.00,0.15,107.14
1/2/1998,570,0.15,43.62,1400.00,0.15,107.14
1/3/1998,597,0.15,45.69,1420.00,0.15,108.67
1/4/1998,630,0.15,48.21,1400.00,0.15,107.14
1/5/1998,633,0.15,48.44,1390.00,0.15,106.38
1/6/1998,632,0.15,48.37,1390.00,0.15,106.38
1/7/1998,601,0.15,45.99,1390.00,0.15,106.38
1/8/1998,553,0.15,42.32,1430.00,0.15,109.44
1/9/1998,531,0.15,40.64,1440.00,0.15,110.20
1/10/1998,538,0.15,41.17,1430.00,0.15,109.44
1/11/1998,581,0.15,44.46,1400.00,0.15,107.14
1/12/1998,597,0.15,45.69,1400.00,0.15,107.14
1/13/1998,602,0.15,46.07,1390.00,0.15,106.38
1/14/1998,596,0.15,45.61,1390.00,0.15,106.38
1/15/1998,596,0.15,45.61,1380.00,0.15,105.61
1/16/1998,632,0.15,48.37,1390.00,0.15,106.38
1/17/1998,597,0.15,45.69,1390.00,0.15,106.38
1/18/1998,609,0.15,46.61,1400.00,0.15,107.14
1/19/1998,593,0.15,45.38,1400.00,0.15,107.14
1/20/1998,547,0.15,41.86,1410.00,0.15,107.91
1/21/1998,547,0.15,41.86,1400.00,0.15,107.14
1/22/1998,478,0.15,36.58,1400.00,0.15,107.14
1/23/1998,497,0.15,38.04,1400.00,0.15,107.14
1/24/1998,503,0.15,38.49,1390.00,0.15,106.38
1/25/1998,515,0.15,39.41,1370.00,0.15,104.85
1/26/1998,508,0.15,38.88,1370.00,0.15,104.85
1/27/1998,520,0.15,39.80,1350.00,0.15,103.32
1/28/1998,530,0.15,40.56,1340.00,0.15,102.55
1/29/1998,490,0.15,37.50,1320.00,0.15,101.02
1/30/1998,567,0.15,43.39,1320.00,0.15,101.02
1/31/1998,558,0.15,42.70,1340.00,0.15,102.55
2/1/1998,505,0.15,38.65,1340.00,0.15,102.55
2/2/1998,523,0.15,40.03,1370.00,0.15,104.85
2/3/1998,558,0.15,42.70,1380.00,0.15,105.61
2/4/1998,536,0.15,41.02,1400.00,0.15,107.14
2/5/1998,521,0.15,39.87,1400.00,0.15,107.14
2/6/1998,507,0.15,38.80,1400.00,0.15,107.14
2/7/1998,488,0.15,37.35,1390.00,0.15,106.38
2/8/1998,489,0.15,37.42,1390.00,0.15,106.38
2/9/1998,498,0.15,38.11,1390.00,0.15,106.38
2/10/1998,492,0.15,37.65,1390.00,0.15,106.38
2/11/1998,499,0.15,38.19,1380.00,0.15,105.61
2/12/1998,494,0.15,37.81,1380.00,0.15,105.61
2/13/1998,488,0.15,37.35,1390.00,0.15,106.38
2/14/1998,563,0.15,43.09,1380.00,0.15,105.61
2/15/1998,525,0.15,40.18,1300.00,0.15,99.49
2/16/1998,470,0.15,35.97,1360.00,0.15,104.08
2/17/1998,470,0.15,35.97,1380.00,0.15,105.61
2/18/1998,525,0.15,40.18,1380.00,0.15,105.61
2/19/1998,489,0.15,37.42,1370.00,0.15,104.85
2/20/1998,525,0.15,40.18,1370.00,0.15,104.85
2/21/1998,501,0.15,38.34,1380.00,0.15,105.61
2/22/1998,491,0.15,37.58,1390.00,0.15,106.38
2/23/1998,476,0.15,36.43,1400.00,0.15,107.14
2/24/1998,324,0.15,24.80,1460.00,0.15,111.73
2/25/1998,303,0.15,23.19,1520.00,0.15,116.33
2/26/1998,324,0.15,24.80,1500.00,0.15,114.80
2/27/1998,380,0.15,29.08,1350.00,0.15,103.32
2/28/1998,436,0.15,33.37,1240.00,0.15,94.90
3/1/1998,440,0.15,33.67,1230.00,0.15,94.13
3/2/1998,433,0.15,33.14,1220.00,0.15,93.37
3/3/1998,460,0.15,35.20,1220.00,0.15,93.37
3/4/1998,654,0.15,50.05,1180.00,0.15,90.31
3/5/1998,600,0.15,45.92,1190.00,0.15,91.07
3/6/1998,579,0.15,44.31,1250.00,0.15,95.66
3/7/1998,544,0.15,41.63,1280.00,0.15,97.96
3/8/1998,536,0.15,41.02,1280.00,0.15,97.96
3/9/1998,536,0.15,41.02,1270.00,0.15,97.19
3/10/1998,547,0.15,41.86,1270.00,0.15,97.19
3/11/1998,638,0.15,48.83,1210.00,0.15,92.60
3/12/1998,680,0.15,52.04,1130.00,0.15,86.48
3/13/1998,704,0.15,53.88,1100.00,0.15,84.18
3/14/1998,663,0.15,50.74,1100.00,0.15,84.18
3/15/1998,684,0.15,52.35,1090.00,0.15,83.42
3/16/1998,712,0.15,54.49,1080.00,0.15,82.65
3/17/1998,804,0.15,61.53,1050.00,0.15,80.36
3/18/1998,843,0.15,64.52,1040.00,0.15,79.59
3/19/1998,982,0.15,75.15,968.00,0.15,74.08
3/20/1998,1130,0.15,86.48,973.00,0.15,74.46
3/21/1998,1340,0.15,102.55,973.00,0.15,74.46
3/22/1998,1300,0.15,99.49,959.00,0.15,73.39
3/23/1998,1140,0.15,87.24,1010.00,0.15,77.30
3/24/1998,1060,0.15,81.12,1010.00,0.15,77.30
3/25/1998,791,0.15,60.54,1090.00,0.15,83.42
3/26/1998,934,0.15,71.48,1060.00,0.15,81.12
3/27/1998,1070,0.15,81.89,963.00,0.15,73.70
3/28/1998,2140,0.15,163.78,NA,0.15,NA
3/29/1998,1940,0.15,148.47,765.00,0.15,58.55
3/30/1998,1810,0.15,138.52,766.00,0.15,58.62
3/31/1998,1570,0.15,120.15,788.00,0.15,60.31
4/1/1998,1350,0.15,103.32,866.00,0.15,66.28
4/2/1998,1250,0.15,95.66,NA,0.15,NA
4/3/1998,1240,0.15,94.90,NA,0.15,NA
4/4/1998,1260,0.15,96.43,NA,0.15,NA
4/5/1998,1180,0.15,90.31,NA,0.15,NA
4/6/1998,1050,0.15,80.36,1000.00,0.15,76.53
4/7/1998,818,0.15,62.60,1040.00,0.15,79.59
4/8/1998,865,0.15,66.20,NA,0.15,NA
4/9/1998,1010,0.15,77.30,NA,0.15,NA
4/10/1998,1050,0.15,80.36,NA,0.15,NA
4/11/1998,2010,0.15,153.83,NA,0.15,NA
4/12/1998,2720,0.15,208.16,NA,0.15,NA
4/13/1998,2680,0.15,205.10,NA,0.15,NA
4/14/1998,2610,0.15,199.74,704.00,0.15,53.88
4/15/1998,1290,0.15,98.72,768.00,0.15,58.78
4/16/1998,825,0.15,63.14,962.00,0.15,73.62
4/17/1998,1120,0.15,85.71,937.00,0.15,71.71
4/18/1998,984,0.15,75.31,963.00,0.15,73.70
4/19/1998,1040,0.15,79.59,969.00,0.15,74.16
4/20/1998,1110,0.15,84.95,947.00,0.15,72.47
4/21/1998,1100,0.15,84.18,949.00,0.15,72.63
4/22/1998,1040,0.15,79.59,949.00,0.15,72.63
4/23/1998,964,0.15,73.78,950.00,0.15,72.70
4/24/1998,1010,0.15,77.30,933.00,0.15,71.40
4/25/1998,1210,0.15,92.60,859.00,0.15,65.74
4/26/1998,1850,0.15,141.58,780.00,0.15,59.69
4/27/1998,2740,0.15,209.69,749.00,0.15,57.32
4/28/1998,2780,0.15,212.76,780.00,0.15,59.69
4/29/1998,2850,0.15,218.11,761.00,0.15,58.24
4/30/1998,2650,0.15,202.81,729.00,0.15,55.79
5/1/1998,1920,0.15,146.94,736.00,0.15,56.33
5/2/1998,1290,0.15,98.72,825.00,0.15,63.14
5/3/1998,1220,0.15,93.37,832.00,0.15,63.67
5/4/1998,1330,0.15,101.79,832.00,0.15,63.67
5/5/1998,1470,0.15,112.50,818.00,0.15,62.60
5/6/1998,1820,0.15,139.29,824.00,0.15,63.06
5/7/1998,2200,0.15,168.37,795.00,0.15,60.84
5/8/1998,2250,0.15,172.19,796.00,0.15,60.92
5/9/1998,2360,0.15,180.61,830.00,0.15,63.52
5/10/1998,2380,0.15,182.14,815.00,0.15,62.37
5/11/1998,2090,0.15,159.95,831.00,0.15,63.60
5/12/1998,1980,0.15,151.53,863.00,0.15,66.05
5/13/1998,1800,0.15,137.76,874.00,0.15,66.89
5/14/1998,1400,0.15,107.14,849.00,0.15,64.97
5/15/1998,1400,0.15,107.14,807.00,0.15,61.76
5/16/1998,1380,0.15,105.61,801.00,0.15,61.30
5/17/1998,1410,0.15,107.91,782.00,0.15,59.85
5/18/1998,1280,0.15,97.96,789.00,0.15,60.38
5/19/1998,1250,0.15,95.66,786.00,0.15,60.15
5/20/1998,1390,0.15,106.38,761.00,0.15,58.24
5/21/1998,1490,0.15,114.03,764.00,0.15,58.47
5/22/1998,1570,0.15,120.15,745.00,0.15,57.02
5/23/1998,1680,0.15,128.57,721.00,0.15,55.18
5/24/1998,1610,0.15,123.21,727.00,0.15,55.64
5/25/1998,1640,0.15,125.51,753.00,0.15,57.63
5/26/1998,1380,0.15,105.61,802.00,0.15,61.38
5/27/1998,1180,0.15,90.31,795.00,0.15,60.84
5/28/1998,1180,0.15,90.31,792.00,0.15,60.61
5/29/1998,1020,0.15,78.06,811.00,0.15,62.07
5/30/1998,923,0.15,70.64,840.00,0.15,64.29
5/31/1998,1000,0.15,76.53,814.00,0.15,62.30
6/1/1998,1170,0.15,89.54,786.00,0.15,60.15
6/2/1998,1370,0.15,104.85,765.00,0.15,58.55
6/3/1998,1660,0.15,127.04,735.00,0.15,56.25
6/4/1998,1790,0.15,136.99,727.00,0.15,55.64
6/5/1998,2120,0.15,162.24,713.00,0.15,54.57
6/6/1998,2880,0.15,220.41,700.00,0.15,53.57
6/7/1998,3170,0.15,242.60,693.00,0.15,53.04
6/8/1998,2550,0.15,195.15,712.00,0.15,54.49
6/9/1998,1250,0.15,95.66,773.00,0.15,59.16
6/10/1998,595,0.15,45.54,964.00,0.15,73.78
6/11/1998,559,0.15,42.78,1040.00,0.15,79.59
6/12/1998,612,0.15,46.84,1060.00,0.15,81.12
6/13/1998,563,0.15,43.09,1090.00,0.15,83.42
6/14/1998,616,0.15,47.14,998.00,0.15,76.38
6/15/1998,679,0.15,51.96,942.00,0.15,72.09
6/16/1998,1040,0.15,79.59,868.00,0.15,66.43
6/17/1998,1020,0.15,78.06,847.00,0.15,64.82
6/18/1998,889,0.15,68.04,860.00,0.15,65.82
6/19/1998,792,0.15,60.61,886.00,0.15,67.81
6/20/1998,744,0.15,56.94,913.00,0.15,69.87
6/21/1998,680,0.15,52.04,927.00,0.15,70.94
6/22/1998,664,0.15,50.82,930.00,0.15,71.17
6/23/1998,703,0.15,53.80,929.00,0.15,71.10
6/24/1998,699,0.15,53.49,879.00,0.15,67.27
6/25/1998,682,0.15,52.19,867.00,0.15,66.35
6/26/1998,856,0.15,65.51,816.00,0.15,62.45
6/27/1998,993,0.15,75.99,777.00,0.15,59.46
6/28/1998,1130,0.15,86.48,709.00,0.15,54.26
6/29/1998,1290,0.15,98.72,671.00,0.15,51.35
6/30/1998,1290,0.15,98.72,653.00,0.15,49.97
7/1/1998,1370,0.15,104.85,651.00,0.15,49.82
7/2/1998,1460,0.15,111.73,731.00,0.15,55.94
7/3/1998,1540,0.15,117.86,687.00,0.15,52.58
7/4/1998,1520,0.15,116.33,632.00,0.15,48.37
7/5/1998,1450,0.15,110.97,647.00,0.15,49.52
7/6/1998,1070,0.15,81.89,673.00,0.15,51.51
7/7/1998,1160,0.15,88.78,653.00,0.15,49.97
7/8/1998,1360,0.15,104.08,642.00,0.15,49.13
7/9/1998,1190,0.15,91.07,660.00,0.15,50.51
7/10/1998,1140,0.15,87.24,681.00,0.15,52.12
7/11/1998,1380,0.15,105.61,663.00,0.15,50.74
7/12/1998,1550,0.15,118.62,638.00,0.15,48.83
7/13/1998,1020,0.15,78.06,640.00,0.15,48.98
7/14/1998,912,0.15,69.80,655.00,0.15,50.13
7/15/1998,995,0.15,76.15,636.00,0.15,48.67
7/16/1998,955,0.15,73.09,641.00,0.15,49.06
7/17/1998,924,0.15,70.71,635.00,0.15,48.60
7/18/1998,961,0.15,73.55,663.00,0.15,50.74
7/19/1998,824,0.15,63.06,695.00,0.15,53.19
7/20/1998,656,0.15,50.20,768.00,0.15,58.78
7/21/1998,596,0.15,45.61,746.00,0.15,57.09
7/22/1998,535,0.15,40.94,777.00,0.15,59.46
7/23/1998,579,0.15,44.31,763.00,0.15,58.39
7/24/1998,938,0.15,71.79,719.00,0.15,55.03
7/25/1998,993,0.15,75.99,704.00,0.15,53.88
7/26/1998,1220,0.15,93.37,711.00,0.15,54.41
7/27/1998,1360,0.15,104.08,681.00,0.15,52.12
7/28/1998,1200,0.15,91.84,697.00,0.15,53.34
7/29/1998,1340,0.15,102.55,630.00,0.15,48.21
7/30/1998,3200,0.15,244.90,602.00,0.15,46.07
7/31/1998,2850,0.15,218.11,628.00,0.15,48.06
8/1/1998,2570,0.15,196.68,684.00,0.15,52.35
8/2/1998,1880,0.15,143.88,710.00,0.15,54.34
8/3/1998,1900,0.15,145.41,684.00,0.15,52.35
8/4/1998,2450,0.15,187.50,756.00,0.15,57.86
8/5/1998,2760,0.15,211.22,717.00,0.15,54.87
8/6/1998,1320,0.15,101.02,837.00,0.15,64.06
8/7/1998,1030,0.15,78.83,934.00,0.15,71.48
8/8/1998,1080,0.15,82.65,883.00,0.15,67.58
8/9/1998,744,0.15,56.94,1010.00,0.15,77.30
8/10/1998,1250,0.15,95.66,910.00,0.15,69.64
8/11/1998,1360,0.15,104.08,927.00,0.15,70.94
8/12/1998,1670,0.15,127.81,889.00,0.15,68.04
8/13/1998,1210,0.15,92.60,1010.00,0.15,77.30
8/14/1998,668,0.15,51.12,1010.00,0.15,77.30
8/15/1998,1180,0.15,90.31,911.00,0.15,69.72
8/16/1998,1450,0.15,110.97,871.00,0.15,66.66
8/17/1998,1350,0.15,103.32,843.00,0.15,64.52
8/18/1998,1260,0.15,96.43,853.00,0.15,65.28
8/19/1998,936,0.15,71.63,897.00,0.15,68.65
8/20/1998,672,0.15,51.43,1000.00,0.15,76.53
8/21/1998,714,0.15,54.64,982.00,0.15,75.15
8/22/1998,771,0.15,59.01,974.00,0.15,74.54
8/23/1998,598,0.15,45.77,1000.00,0.15,76.53
8/24/1998,488,0.15,37.35,NA,0.15,NA
8/25/1998,528,0.15,40.41,NA,0.15,NA
8/26/1998,812,0.15,62.14,936.00,0.15,71.63
8/27/1998,720,0.15,55.10,993.00,0.15,75.99
8/28/1998,653,0.15,49.97,944.00,0.15,72.24
8/29/1998,796,0.15,60.92,892.00,0.15,68.27
8/30/1998,689,0.15,52.73,901.00,0.15,68.95
8/31/1998,511,0.15,39.11,1020.00,0.15,78.06
9/1/1998,410,0.15,31.38,1150.00,0.15,88.01
9/2/1998,374,0.15,28.62,1180.00,0.15,90.31
9/3/1998,444,0.15,33.98,1170.00,0.15,89.54
9/4/1998,473,0.15,36.20,1150.00,0.15,88.01
9/5/1998,727,0.15,55.64,951.00,0.15,72.78
9/6/1998,1130,0.15,86.48,735.00,0.15,56.25
9/7/1998,1090,0.15,83.42,737.00,0.15,56.40
9/8/1998,1020,0.15,78.06,773.00,0.15,59.16
9/9/1998,916,0.15,70.10,793.00,0.15,60.69
9/10/1998,859,0.15,65.74,820.00,0.15,62.76
9/11/1998,595,0.15,45.54,900.00,0.15,68.88
9/12/1998,458,0.15,35.05,1060.00,0.15,81.12
9/13/1998,429,0.15,32.83,1160.00,0.15,88.78
9/14/1998,418,0.15,31.99,1110.00,0.15,84.95
9/15/1998,493,0.15,37.73,1150.00,0.15,88.01
9/16/1998,438,0.15,33.52,1170.00,0.15,89.54
9/17/1998,406,0.15,31.07,1150.00,0.15,88.01
9/18/1998,397,0.15,30.38,1140.00,0.15,87.24
9/19/1998,429,0.15,32.83,1110.00,0.15,84.95
9/20/1998,450,0.15,34.44,1080.00,0.15,82.65
9/21/1998,449,0.15,34.36,1060.00,0.15,81.12
9/22/1998,450,0.15,34.44,1090.00,0.15,83.42
9/23/1998,478,0.15,36.58,1060.00,0.15,81.12
9/24/1998,531,0.15,40.64,1030.00,0.15,78.83
9/25/1998,573,0.15,43.85,999.00,0.15,76.45
9/26/1998,533,0.15,40.79,1040.00,0.15,79.59
9/27/1998,564,0.15,43.16,1020.00,0.15,78.06
9/28/1998,434,0.15,33.21,1130.00,0.15,86.48
9/29/1998,399,0.15,30.54,1180.00,0.15,90.31
9/30/1998,359,0.15,27.47,1230.00,0.15,94.13"),header=T)
closeAllConnections()
a[,1] <- transform(a$date,as.Date(a[,1],"%m/%d/%Y"))

#convert cfs to cfd
Cat.Dam.ConvertQ <- function(Q) (Q * 86400)
#convert SC to TDS using locally regressed data (and do unit conversion)
Cat.Dam.ConvertSC <- function(SC) ((0.8652 * SC - 145.43) *
(28.3168/1000000/907.185))

# estimated annul load, omits missing days
sum(Cat.Dam.ConvertQ(a$Q) * Cat.Dam.Convert(a$SC), na.rm=TRUE)
# 565239.3

# 10 realizations for the first day
Cat.Dam.ConvertQ(rnorm(10,a$Q[1],a$Qsd[1])) *
Cat.Dam.ConvertSC(rnorm(10,a$SC[1],a$SCsd[1]))

	[[alternative HTML version deleted]]


