[R] Trellis Plots: translating lattice xyplot() to ggplot()

[Rich Shepard]

On Fri, 10 Jul 2015, Hadley Wickham wrote:

> You shouldn't be explicitly opening a device in a knitr document.

Hadley,

   Didn't think so.

> I think maybe you should post a minimal document so we can figure out
> what's going wrong.

   Agreed. Attached are the raw data (carlin.csv) and a stripped down LyX
document with the knitr chunks.

   This is my first attempt to use knitr; I'm reading the knitr book and
that's where I got the impression that lattice graphics are not supported.

Rich
-------------- next part --------------
siteid,sampdate,Temp.h2o,Temp.air,Disc.cfs,Turb,SC,DO,pH,ANC,HCO3,CO3,Alk,Hard,TDS,TSS,N.tot,N.org,NH4,NO3,NO2,PO4,C,Ca,Mg,Na,K,Cl,SO4,F,Si,As,Ba,Be,Cd,Cr,Co,Cu,Fe,Pb,Mn,Mo,Ni,Ag,Sr,V,Zi,Al,Li,Se,CN,Hg
10321000,1965-10-01,,,98,,490,,8.4,,238,4,,150,,,,,,,,,,40,12,46,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1965-11-01,,,141,,507,,8.3,213,252,4,,170,,,,,,,,,,50,11,42,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1965-12-01,,,128,,564,,8.1,238,290,0,,190,359,,,,,,,,,56,13,47,6.4,16,37,0.5,40,,,,,,,,,,,,,,,,,,,,,
10321000,1966-01-01,,,122,,506,,8.2,212,259,0,,170,321,,,,,,,,,50,11,45,6.1,16,34,0.6,31,,,,,,,,,,,,,,,,,,,,,
10321000,1966-01-17,,,97,,551,,8.5,241,276,9,,200,,,,,,,,,,58,13,48,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-02-01,,,108,,506,,8.2,216,263,0,,180,,,,,,,,,,52,11,42,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-03-01,,,384,,516,,7.7,201,245,0,,160,,,,,,,,,,46,11,46,8,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-03-18,,,701,,628,,8.1,235,286,0,,190,,,,,,,,,,51,16,61,8.5,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-04-01,,,612,,470,,8,189,231,0,,160,301,,,,,,,,,43,12,40,6.3,14,37,0.5,34,,,,,,,,,,,,,,,,,,,,,
10321000,1966-04-11,,,541,,423,,8,172,210,0,,140,,,,,,,,,,40,10,34,5.4,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-04-18,,,355,,466,,8.1,193,235,0,,160,,,,,,,,,,43,12,39,5.9,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-05-01,,,167,,498,,8,208,253,0,,170,,,,,,,,,,46,13,42,6.8,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-05-06,,,260,,420,,8.1,181,221,0,,140,,,,,,,,,,41,10,33,5.4,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-07-01,,,7,,566,,8,226,275,0,,180,450,,,,,,,,,47,16,52,96,21,48,0.6,34,,,,,,,,,,,,,,,,,,,,,
10321000,1966-08-01,,,3.9,,550,,8.1,210,256,0,,170,,,,,,,,,,40,17,51,9.7,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-09-01,,,4.1,,536,,8.4,208,254,0,,170,,,,,,,,,,41,17,48,8.9,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-10-01,,,7.7,,546,,8.4,220,260,4,,190,340,,,,,,,,,48,16,46,8.5,19,43,0.5,27,,,,,,,,,,,,,,,,,,,,,
10321000,1966-11-01,,,16,,536,,8.4,224,263,5,,180,,,,,,,,,,49,15,45,7.6,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1966-12-01,,,31,,533,,8.6,222,248,11,,180,336,,,,,,,,,50,13,47,7.5,14,37,0.5,34,,,,,,,,,,,,,,,,,,,,,
10321000,1966-12-16,,,30,,615,,8.6,264,296,13,,210,,,,,,,,,,60,14,56,8.7,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-01-01,,,43,,517,,8.8,214,228,16,,170,333,,,,,,,,,50,12,48,7.3,19,40,0.5,28,,,,,,,,,,,,,,,,,,,,,
10321000,1967-02-01,,,72,,513,,8.7,217,234,15,,180,,,,,,,,,,51,12,46,6.9,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-02-12,,,113,,474,,8.3,189,230,3,,160,,,,,,,,,,46,11,41,6.8,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-03-01,,,232,,464,,8.5,174,212,7,,150,,,,,,,,,,44,10,41,6.6,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-03-25,,,422,,498,,8.5,179,218,6,,160,,,,,,,,,,48,10,43,6.5,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-04-01,,,333,,448,,8.5,164,200,5,,150,283,,,,,,,,,42,11,39,5.7,15,40,0.5,26,,,,,,,,,,,,,,,,,,,,,
10321000,1967-05-01,,,327,,451,,8.5,165,201,6,,150,,,,,,,,,,43,10,38,5.8,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-05-10,,,406,,371,,8.3,127,155,2,,110,,,,,,,,,,29,9.5,34,5.7,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-05-14,,,396,,435,,8.3,158,193,2,,160,,,,,,,,,,45,11,38,6.7,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-05-20,,,844,,365,,8,147,179,0,,150,,,,,,,,,,45,8.1,25,5.6,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-06-01,,,1560,,443,,8.4,189,230,4,,180,,,,,,,,,,55,10,35,6.7,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-07-01,,,1340,,386,,8.5,175,213,13,,150,255,,,,,,,,,47,8.6,27,4.8,7.4,14,0.5,28,,,,,,,,,,,,,,,,,,,,,
10321000,1967-07-11,,,651,,406,,8.3,187,228,2,,170,,,,,,,,,,55,9.1,27,5,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-07-21,,,254,,453,,8.2,207,252,0,,190,,,,,,,,,,58,10,35,6,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-08-01,,,48,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-08-01,,,48,,525,,8.2,223,272,0,,190,,,,,,,,,,53,13,47,7.4,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-09-01,,,18,,508,,8.1,210,256,0,,180,,,,,,,,,,46,15,45,7.9,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-10-01,,,20,,556,,8.4,221,255,7,,180,322,,,,,,,,,47,15,45,8.2,16,32,0.6,24,,,,,,,,,,,,,,,,,,,,,
10321000,1967-10-01,,,7.7,,546,,8.4,213,260,4,,190,340,,,,,0.02,,,,48,16,46,8.5,19,43,0.5,27,,,,,,,,,,,,,,,,,,,,,
10321000,1967-10-09,,,27,,612,,8.4,242,281,7,,190,,,,,,,,,,53,15,57,11,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-10-17,,,23,,586,,8.2,238,290,0,,190,,,,,,,,,,53,15,49,8.3,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-11-01,,,36,,582,,8.4,238,278,6,,190,,,,,,,,,,55,14,50,7.8,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-12-01,,,44,,565,,8.2,239,292,0,,170,,,,,,,,,,50,12,53,8.2,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-12-09,,,41,,626,,8.2,271,330,0,,210,,,,,,,,,,60,14,57,8.7,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1967-12-19,,,46,,553,,8.2,241,294,0,,180,,,,,,,,,,54,12,51,8.3,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-01-01,,,50,,519,,8.6,218,238,14,,170,334,,,,,,,,,45,13,50,7.9,16,41,0.6,28,,,,,,,,,,,,,,,,,,,,,
10321000,1968-01-15,,,48,,590,,8.3,251,298,4,,200,372,,,,,,,,,57,13,52,8.6,16,43,0.6,30,,,,,,,,,,,,,,,,,,,,,
10321000,1968-02-01,,,70,,406,,7.8,142,173,0,,110,,,,,,,,,,27,11,40,6.6,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-02-11,,,221,,480,,8.3,195,232,3,,160,,,,,,,,,,44,11,43,6.3,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-03-01,,,259,,503,,8.4,204,241,4,,160,,,,,,,,,,45,11,46,6.9,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-03-09,,,221,,471,,8.3,191,227,3,,150,,,,,,,,,,44,10,40,6.4,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-03-15,,,229,,479,,8,195,238,0,,160,,,,,,,,,,44,12,41,6.4,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-04-01,,,169,,415,,8.2,157,192,0,,130,271,,,,,,,,,36,10,38,5.8,18,35,0.5,27,,,,,,,,,,,,,,,,,,,,,
10321000,1968-05-01,,,160,,381,,8.4,163,193,3,,130,,,,,,,,,,38,8.6,33,5.5,,26,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-05-15,,,164,,380,,8,154,188,0,,120,235,,,,,,,,,33,8.7,36,5.8,12,26,0.4,19,,,,,,,,,,,,,,,,,,,,,
10321000,1968-06-01,,,498,,311,,7.8,138,168,0,,120,,,,,,,,,,37,5.8,20,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-06-07,,,914,,394,,8.2,174,212,0,,140,,,,,,,,,,42,8.5,31,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-06-15,,,985,,424,,7.7,184,224,0,,150,,,,,,,,,,45,9.2,33,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-07-01,,,140,,436,,8.2,195,238,0,,150,273,,,,,,,,,44,10,37,5.7,11,21,1.1,24,,,,,,,,,,,,,,,,,,,,,
10321000,1968-08-01,,,20,,520,,8.6,213,240,10,,170,,,,,,,,,,45,13,49,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-08-15,,,12,,566,,8.2,231,282,0,,180,,,,,,,,,,50,14,52,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-08-24,,,44,,473,,8.6,205,226,12,,160,,,,,,,,,,45,11,42,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-09-01,,,29,,512,,8.4,220,256,6,,170,,,,,,,,,,45,13,48,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1968-09-15,,,28,,536,,8.1,225,274,0,,170,,,,,,,,,,47,13,49,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-05-20,16,,126,,483,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-06-29,23.5,,149,,445,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-07-29,24.5,,14,,517,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-08-11,15.5,,14,,464,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-10-06,10.5,,17,,506,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-11-09,,,26,,496,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1977-12-19,0,,26,,551,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-01-24,0,,64,,527,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-02-27,6.5,,156,,541,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-03-24,10,,557,,407,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-04-20,12,,579,,389,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-07-07,17.5,,370,,384,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-08-22,19.5,,19,,550,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-09-21,16.5,,52,,498,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1978-11-21,4,,77,,490,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-01-04,0.5,,54,,619,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-02-05,0.5,,103,,495,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-03-02,3.5,,338,,536,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-04-26,13,,724,,359,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-04-26,13,,724,54,374,8.6,8.6,210,230,13,,140,273,90,0.6,0.46,0.03,,,,6.1,41,8.8,29,5.5,11,24,0.4,27,,,,,,,,,,,,,,,,,,,,,
10321000,1979-05-29,16.5,,2160,84,366,7.9,8.2,280,340,0,,130,234,53,1.1,1,0.01,,,,,39,8.3,28,7.1,9.7,23,0.4,28,6,100,,,,,,20,,20,,,,,,20,,,1,,0.10
10321000,1979-07-09,22,,204,13,418,8,8.6,215,240,11,,160,283,89,0.43,0.37,0.010,,,,5.1,46,11,33,5.6,9.3,21,0.4,27,,,,,,,,,,,,,,,,,,,,,
10321000,1979-07-09,22,,204,,377,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-07-30,26,,88,,448,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-07-30,26,,80,10,458,7.5,8.8,235,280,3,,160,299,83,0.66,0.62,0.010,,,,,44,11,35,7.1,11,24,0.4,25,8,90,,2,,3,,10,,1,,,,,,3,,,1,,0.10
10321000,1979-08-29,21.5,,29,,518,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1979-08-29,21.5,,19,13,509,8.5,8.4,185,220,3,,160,201,78,0.69,0.42,0.010,,,,5.4,46,11,48,9.6,21,39,0.6,26,,,,,,,,,,,,,,,,,,,,,
10321000,1979-10-26,12.5,,42,7,610,11,8.9,259,230,42,,190,395,,0.81,0.49,0.09,,,,,54,14,55,11,27,48,0.6,30,10,100,,,0,3,,10,0,,,,0,,,3,,,0,,0.2
10321000,1979-11-27,1,,74,10,521,12.2,8.6,204,220,14,,170,327,,0.59,0.39,0.01,,,,4,50,11,44,7.3,17,33,0.6,31,,,,,,,,,,,,,,,,,,,,,
10321000,1979-12-21,1,,96,6,514,12,8.5,,,,,170,,,1.2,0.92,0,,,,7.8,50,12,43,7.2,17,33,0.5,32,,,,,,,,,,,,,,,,,,,,,
10321000,1980-02-05,3,,227,21,522,11.2,7.9,205,250,0,,180,333,89,0.85,0.49,0.04,,,,,51,13,48,7.4,17,39,0.5,33,6,90,,1,0,3,,30,,,,0,0,,,3,,,0,,0.2
10321000,1980-02-28,8,,612,60,550,9.6,,197,240,0,,170,331,83,1.2,0.88,,,,,9.1,47,12,47,8.2,20,48,0.5,29,,,,,,,,,,,,,,,,,,,,,
10321000,1980-03-27,5,,299,19,505,,7.6,197,240,0,,160,307,,0.56,0.48,0.01,,,,13,46,12,42,6.4,15,37,0.4,30,,,,,,,,,,,,,,,,,,,,,
10321000,1980-04-24,12,,652,130,308,8.3,8.3,,,,,100,,81,1.8,1.6,0.08,,,,13,30,7.1,25,5.3,9.9,18,0.3,28,,,,,,,,,,,,,,,,,,,,,
10321000,1980-05-29,14.5,,3180,88,449,7.6,7.8,190,233,0,,140,281,,1.4,1.2,0.06,,,,,42,9.6,41,7.9,13,29,0.2,29,7,90,,1,0,3,10,20,0,10,,0,0,,,,,,0,,0.1
10321000,1980-06-23,20.5,,2070,23,362,6.7,,,,,,140,,,,,0.04,,,,8.9,42,9.2,24,5.2,7.9,11,0.3,32,,,,,,,,,,,,,,,,,,,,,
10321000,1980-07-29,24,,203,11,427,8.1,8.2,180,230,0,,160,257,,0.63,0.5,0.01,,,,,46,9.9,31,5.9,10,20,0.5,25,7,90,,,0,3,,40,,20,,0,0,,,10,,,0,,0
10321000,1980-08-25,19.5,,49,11,545,9.2,7.7,216,260,0,,170,342,,1.4,1.1,0,,,,10,50,12,47,8.2,20,40,0.8,26,,,,,,,,,,,,,,,,,,,,,
10321000,1980-09-30,17,,47,6,530,11,8.4,223,260,6,,170,336,,0.82,0.74,0,,,,4.4,47,12,50,8.8,23,36,0.6,25,,,,,,,,,,,,,,,,,,,,,
10321000,1980-11-21,5,,108,10,553,12.2,8.5,210,220,8,,180,346,,0.6,0.46,0.02,,,,,52,12,45,8.5,21,52,0.5,28,7,100,,1,0,3,,10,0,10,,0,0,,,3,,,0,,0
10321000,1981-01-28,3.5,,127,14,527,12,8.7,210,270,19,,170,355,,,,0.05,,,,,50,12,46,6.5,19,38,0.5,30,8,100,,1,0,3,0,30,,20,,0,0,,,10,,,0,,0
10321000,1981-03-25,9.5,,110,24,600,11.8,8.8,230,250,17,,200,364,,0.55,0.51,0.05,,,,4.9,56,14,53,6.8,20,47,0.5,28,,,,,,,,,,,,,,,,,,,,,
10321000,1981-05-20,13,,117,16,474,,8.6,200,230,7,,180,,,1,0.94,0.08,,,,,55,10,36,7.7,14,,0.5,27,9,,,1,0,3,,20,,20,,,0,,,,7,,0,,0
10321000,1981-07-21,25.5,,5.7,9,616,9.2,8.6,240,270,12,,200,382,,1,0.71,0.13,,,,,54,15,55,10,25,43,0.5,34,14,120,,1,,3,,10,,10,,,0,,,,,,0,,0.1
10321000,1981-09-17,21.5,,9.5,3,487,10.2,8.1,200,250,0,,170,292,,1.2,0.97,0.13,,,,,44,14,37,8,15,25,0.6,24,9,100,,1,0,3,,10,,20,,,0,,,10,,,1,,0
10321000,1981-12-04,3,,23,4,623,13,8.8,240,260,15,,210,406,,0.57,0.38,0.09,0.100,0.020,,,62,13,58,9.5,33,56,0.6,30,9,100,,1,10,3,1,30,1,20,,,1,,,3,,,1,,0.10
10321000,1982-01-27,1,,61,8,529,12.5,8.2,210,,,,180,342,,0.88,,,,,0.55,,52,11,43,7.7,23,42,0.5,32,,,,,,,,,,,,,,,,,,,,,
10321000,1982-03-12,7,,649,96,395,10.7,8.4,160,,,,130,248,87,0.70,0.48,0.13,,,0.31,,37,9.3,33,7.4,15,22,0.4,28,9,82,,1,10,,,,,,,,1,,,4,,,1,,0.10
10321000,1982-05-26,19,,1530,56,367,7,8.2,158,,,,130,235,70,1.3,1.1,0.12,0.100,0.020,0.31,,37,8.5,26,5.4,9.6,22,0.4,31,7,84,,1,10,1,,80,1,,,1,1,,,20,,,1,,0.10
10321000,1982-07-20,25,,349,26,418,7.5,8,190,,,,150,253,,1.7,1.5,0.14,,,0.31,,45,9.3,25,6,9,17,0.4,27,,,,,,,,,,,,,,,,,,,,,
10321000,1982-09-21,20.5,,56,6,510,,8.6,208,,,,160,318,,1.0,0.84,0.06,,,0.37,,49,10,46,8.5,20,32,0.5,27,9,100,,1,10,1,,,1,,,,1,,,10,,,1,,0.10
10321000,1982-11-24,,,281,25,486,12.7,8.6,213,,,,180,313,,1.1,0.87,0.13,,,0.09,,50,12,38,5.6,14,35,0.4,29,5,81,0.5,1,1,3,,,1,,10,1,1,300,6.0,10,10,30,1,,0.10
10321000,1983-01-25,2,,270,20,476,13.1,8.4,195,,,,160,294,,0.83,0.58,0.12,,,0.25,,47,11,37,5.5,16,30,0.4,30,,,,,,,,,,,,,,,,,,,,,
10321000,1983-03-06,5,,4800,200,308,,7.6,109,,,,79.9,175,79,2.1,1.6,0.33,,,0.521,,23,5.4,26,6.9,10,21,0.3,15,6,49,0.5,1,1,3,,40,,20,10,,1,150,6.0,,40,20,1,,0.1
10321000,1983-05-24,19.5,26,1500,110,359,7.2,8.2,153,,,,127,232,86,1.2,0.97,0.13,,,0.184,,36,9,26,5.4,12,20,0.3,30,6,73,0.5,1,1,3,,30,,,10,,1,220,6.0,40,40,20,1,,0.10
10321000,1983-07-20,19.5,22.5,499,26,369,7.7,8.4,176,,,,143,231,,1.0,0.77,0.13,,,0.245,,43,8.7,23,4.4,7.1,16,0.4,22,,,,,,,,,,,,,,,,,,,,,
10321000,1983-09-20,16,17,75,4,505,11,,212,,,,175,323,,0.9,0.80,0.100,,,0.368,,50,12,47,7.1,22,31,0.5,25,8,98,0.5,1,1,3,,3,,10,10,,1,390,6.0,,10,40,1,,0.10
10321000,1983-11-15,5.5,12.5,335,13,512,12.2,8.4,230,,,,193,333,88,1.1,0.87,0.13,,,0.123,,54,14,43,5.8,17,32,0.4,28,5,87,0.5,1,1,3,,,1,10,10,1,1,350,6.0,,10,40,1,,0.10
10321000,1984-01-24,0,3.5,369,32,532,10.7,7.7,222,,,,190,335,80,0.97,0.68,0.12,,,0.184,,53,14,40,5.6,20,37,0.4,31,,,,,,,,,,,,,,,,,,,,,
10321000,1984-03-14,3.5,8.5,2630,640,361,11,8.1,134,,,,111,215,80,2.6,2.4,0.13,,,0.245,,32,7.4,30,5.9,13,24,0.2,20,5,67,0.5,,1,3,,130,1,20,10,1,1,200,6.0,130,130,20,1,,0.10
10321000,1984-05-18,15,18.5,8130,180,352,8.4,8.1,140,,,,111,224,,2.1,1.9,0.07,,,0.153,,31,8,28,10,13,22,0.4,27,5,59,1,1,1,3,,,,10,10,1,1,190,6,3,10,20,1,,0.10
10321000,1984-07-25,24,34,1200,39,327,,,142,,,,126,196,,0.93,0.71,0.09,,,0.153,,38,7.5,19,3.7,6.9,14,0.3,21,,,,,,,,,,,,,,,,,,,,,
10321000,1984-09-27,13.5,20.5,160,8,519,10.6,8.5,205,,,,179,318,,0.40,0.27,0.03,,,0.245,,50,13,39,7.7,21,37,0.5,26,7,110,1,1,1,3,,,,10,10,,1,370,6.0,,10,30,1,,0.10
10321000,1984-11-30,2.5,9.5,335,15,600,14,8.4,239,,,,208,,,1.1,0.95,0.05,,,0.153,,55,17,48,6.9,,,0.4,30,5,110,0.5,1,,3,,,,20,10,,1,390,6.0,3,10,30,1,,0.10
10321000,1985-01-22,0.5,3.5,218,31,404,12.2,8.2,129,,,,144,294,,0.86,0.49,0.11,,,0.276,,41,10,29,4.3,19,29,0.3,21,,,,,,,,,,,,,,,,,,,,,
10321000,1985-03-27,3.5,4,1040,90,476,12.8,8.3,169,,,,166,305,88,1.5,1.2,0.07,,,0.215,,45,13,38,6.8,20,41,0.3,27,6,83,0.6,1.0,1,3.0,5,91,7,11,10.0,2,1,260,6.0,18,180,26,1,,0.1
10321000,1985-05-21,17.5,22.5,606,27,430,8.7,8.4,166,,,,151,268,90,0.60,0.42,0.08,,,0.031,,44,10,30,6.1,14,28,0.4,28,5,80,0.50,1.0,1,3.0,6,27,4,8,10.0,1,1,260,6.0,28,30,22,1,,0.1
10321000,1985-07-23,23.5,26,54,17,577,8.7,8.6,215,,,,176,358,98,0.70,0.57,0.03,,,0.276,,49,13,50,9.9,26,42,0.5,29,,,,,,,,,,,,,,,,,,,,,
10321000,1985-09-26,13.5,17,33,9,675,10.6,8.8,229,,,,188,394,,0.90,0.72,0.08,,,1.17,,52,14,66,13,40,57,0.6,27,12,140,0.50,1.0,1,3.0,1,4,1,11,10.0,1,1,590,6.0,14,10,67,1,,0.1
10321000,1985-11-19,0,-3.5,76,11,579,15.2,8.7,205,,,,190,361,,0.61,0.37,0.1,0.09,0.02,0.429,,56,13,50,7.1,26,51,0.6,33,8,110,0.50,1.0,1,3.0,3,9,1,56,10.0,2,1,440,6.0,10,10,37,1,,0.2
10321000,1986-01-22,1,5.5,247,50,451,12.6,8.3,165,,,,140,273,86,1.1,0.66,0.2,0.16,0.02,0.521,,39,11,36,7.1,16,38,0.4,26,,,,,,,,,,,,,,,,,,,,,
10321000,1986-03-27,12,21.5,822,47,399,9.8,8.2,146,,,,131,240,89,0.60,0.45,0.05,0.090,0.01,0.153,,36,10,29,5.3,16,29,0.4,28,7,68,0.50,1.0,1,3.0,4,47,1,8,10.0,1,1,230,6.0,7,50,18,1,,0.1
10321000,1986-05-22,12,11,942,40,384,8,8.4,149,,,,131,229,,1.0,0.85,0.01,0.100,0.010,0.184,,38,8.8,25,4.8,11,24,0.3,28,,,,,,,,,,,,,,,,,,,,,
10321000,1986-07-23,21,21,126,10,443,8.9,8.7,190,,,,159,276,,0.50,0.36,0.02,0.090,0.01,0.245,,47,10,33,6.3,14,30,0.4,23,,,,,,,,,,,,,,,,,,,,,
10321000,1986-09-23,15,22.5,29,1,612,12.3,8.8,231,,,,201,400,,,,0.02,0.100,0.010,0.797,,57,14,59,13,34,54,0.6,28,10,120,0.50,1.0,1,3.0,1,9,,12,10.0,1,1,510,6.0,9,10,55,1,,0.1
10321000,1986-11-25,4,6,152,10,475,11.6,8.3,198,,,,177,303,,0.30,0.19,0.010,0.100,0.010,0.061,,51,12,39,6,16,32,0.4,29,7,90,0.50,1.0,1,3.0,1,79,5,13,10.0,1,1,360,6.0,17,10,32,1,,0.1
10321000,1987-01-28,0.5,3.5,74,3,506,11.8,8.3,,235,4,,177,309,,0.60,0.48,0.010,0.100,0.010,0.092,,51,12,39,6.1,16,37,0.5,28,,,,,,,,,,,,,,,,,,,,,
10321000,1987-03-24,7.5,7.5,300,33,603,10.5,8.6,230,,,,209,392,79,0.56,0.44,0.01,0.06,0.010,0.123,,57,16,55,7.4,27,60,0.6,33,8,100,0.50,1.0,1,3.0,1,19,5,11,10.0,2,1,370,6.0,4,10,39,1,,0.5
10321000,1987-06-01,18.5,18,483,31,460,,8.5,193,,,,159,285,,1.3,1.1,0.02,0.100,0.010,0.123,,47,10,34,5.9,16,28,0.4,31,3,82,0.50,1,3,3.0,8,27,5,8,10.0,1,1,280,6.0,12,10,26,1,,0.1
10321000,1987-08-26,21,26.5,10,5,491,9.8,8.6,,,,,171,306,85,0.30,0.17,0.010,0.100,0.010,0.031,,45,14,41,8.6,21,30,0.8,24,9,100,0.50,1.0,2,3.0,2,4,5,12,10.0,1,1,490,6.0,7,10,45,1,,0.1
10321000,1987-11-03,14.5,16,8,3,521,,8.5,205,,,,173,307,98,0.60,0.47,0.02,0.100,0.010,0.031,,46,14,46,7.8,20,38,0.5,20,10,110,0.50,1.0,1,3.0,2,3.0,5,10,10.0,4,1,420,6.0,3.0,10,41,1,,0.1
10321000,1988-01-28,0.5,4,33,1,588,14,8.6,235,,,,201,369,,0.30,0.18,0.010,0.090,0.01,0.031,,59,13,48,8.3,22,46,0.5,31,,,,,,,,,,,,,,,,,,,,,
10321000,1988-03-31,8.5,11.5,234,33,409,,8.6,181,,,,162,322,98,0.60,0.47,0.03,0.100,0.010,0.153,,45,12,38,6.5,17,44,0.4,30,6,74,0.50,1.0,2,3.0,1,10,5,10,10.0,5,1,310,6.0,5,10,24,1,,0.1
10321000,1988-05-26,17.5,29.5,384,41,300,8.3,8.4,114,,,,109,195,98,0.40,0.28,0.02,0.100,0.010,0.061,,34,5.9,20,4.3,8.5,26,0.4,21,4,55,0.50,1.0,1,3.0,1,9,5,3,10.0,1,1,210,6.0,8,10,17,1,,0.1
10321000,1988-08-24,24.5,34.5,11,4,479,8.8,8.5,187,,,,171,298,85,1.0,0.90,0.010,0.100,0.010,0.031,,45,14,40,7.8,15,28,0.5,24,8,100,0.50,1.0,1,3.0,1,5,5,14,10.0,3,1,480,6.0,3.0,10,51,1,,0.1
10321000,1988-11-29,0.5,0.5,64,13,548,12.4,8.6,225,,,223,187,358,,1,0.36,0.04,0.62,0.01,0.276,,53,13,52,7.8,22,44,0.5,28,8,100,0.50,1.0,1,3.0,1,16,5,11,10.0,5,1,490,6.0,3,10,47,1,,0.1
10321000,1988-11-29,0.5,,,,545,12.4,8.6,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1988-11-29,0.5,,,,550,12.4,8.6,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1988-11-29,0.5,,,,545,12.4,8.6,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-01-26,0,-5.0,29,68,615,14.9,8.5,,,,270,217,383,,0.30,0.16,0.03,0.100,0.010,0.031,,64,14,51,8.9,17,36,0.3,30,,,,,,,,,,,,,,,,,,,,,
10321000,1989-04-08,14.5,21.5,905,50,395,8.4,8.3,136,,,136,128,240,80,0.80,0.65,0.04,0.100,0.010,0.031,,36,9.2,29,5.8,15,35,0.3,28,7,67,0.50,1.0,1,3.0,8,47,5,7,10.0,6,1,230,6.0,6,60,17,1,,0.1
10321000,1989-05-31,17,21,524,12,378,7.6,8.2,167,,,166,142,242,90,0.60,0.46,0.010,0.100,0.010,0.123,,42,9,27,5.2,8.9,21,0.3,28,5,69,0.50,1.0,1,3.0,4,13,1,6,10.0,2,1.0,250,6.0,5,20,20,1,,0.1
10321000,1989-05-31,17,,,,380,7.6,8.2,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-05-31,17,,,,378,7.6,8.1,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-05-31,17,,,,376,7.6,8.1,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-07-31,24.5,33.5,19,4,530,9.1,8.7,201,,,,163,320,72,0.60,0.50,0.010,0.100,0.010,0.031,,44,13,47,8.7,20,43,0.7,23,,,,,,,,,,,,,,,,,,,,,
10321000,1989-08-28,23,23.5,15,4,445,9.4,8.5,189,,,188,152,272,87,0.40,0.29,0.01,0.100,0.010,0.031,,41,12,36,8,14,27,0.5,21,8,98,0.50,1.0,1,3.0,1,6,1,8,10.0,1.0,1.0,440,6.0,4,10,43,1,,0.1
10321000,1989-08-28,22,,,,445,9.4,8.5,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-08-28,22,,,,445,9.4,8.5,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-08-28,22,,,,442,9.4,8.6,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,1989-11-29,0.5,2,49,12,502,13.2,8.3,208,,,240,175,350,68,0.50,0.39,0.010,0.100,0.010,0.061,,50,12,45,7.9,21,41,0.6,29,8,99,0.50,1.0,1,3.0,5,16,1,14,10.0,5,1.0,450,6.0,3.0,20,41,1,,0.1
10321000,1990-01-17,1,2.5,71,15,459,12.2,8.7,183,,,190,160,294,64,0.40,0.28,0.01,0.100,0.010,0.061,,46,11,38,5.5,17,36,0.5,27,,,,,,,,,,,,,,,,,,,,,
10321000,1990-03-30,11,17,307,26,420,9.9,8.5,175,,,168,149,272,99,0.90,0.77,0.010,0.100,0.010,0.092,,43,10,34,6.6,17,32,0.3,28,7,72,0.7,1.0,5.0,3.0,10.0,29,10,6,10,10.0,1,280,6.0,3,50,24,1,,0.1
10321000,1990-05-22,20,26,195,9,396,9.2,8.6,164,,,169,133,242,83,0.70,0.60,0.010,0.100,0.010,0.061,,40,7.9,29,5.9,13,24,0.1,21,5,62,0.50,1.0,1,3.0,2,11,1,2,10.0,1.0,1.0,260,6.0,8,10,22,1,,0.1
10321000,1990-06-28,20,28,238,10,378,8,8.7,174,,,173,139,235,93,0.50,0.37,0.03,0.100,0.010,0.031,,42,8.4,26,4.6,11,16,0.4,23,,,,,,,,,,,,,,,,,,,,,
10321000,1990-08-21,19.5,21,11,5,437,9.5,8.4,183,,,186,150,272,93,0.80,0.69,0.02,0.090,0.01,0.031,,40,12,34,7.3,14,29,1.3,22,7,84,0.50,1.0,1,3.0,2,3,1,4,10.0,1,1.0,430,6.0,3.0,20,42,1,,0.1
10321000,1990-11-01,7,4.5,14,3,576,13.6,8.6,209,,,210,179,302,61,0.60,0.49,0.010,0.100,0.010,0.061,,50,13,38,7.3,16,30,0.2,21,7,120,0.50,1.0,1,3.0,1,6,1,12,10.0,1.0,1.0,500,6.0,27,10,44,1,,0.1
10321000,1991-01-03,0,-5,15,3,623,13.2,7.8,264,,,265,219,397,43,0.4,0.25,0.03,0.1,0.010,0.031,,63,15,48,9.7,22,52,0.6,27,,,,,,,,,,,,,,,,,,,,,
10321000,1991-02-22,7,17,152,20,478,11.8,7.9,195,,,170,156,290,86,0.70,0.57,0.01,0.100,0.010,0.092,,46,9.9,39,6.3,14,46,0.4,25,6,79,0.50,1.0,1,3.0,3,25,1,11,10.0,1,1.0,330,6.0,5,40,30,1,,0.4
10321000,1991-04-29,11.5,8.5,119,14,488,10.4,8.7,199,,,193,159,295,99,0.56,0.48,0.03,0.05,0.01,0.031,,47,10,39,6.7,15,34,0.6,26,6,80,0.50,1.0,1,3.0,2,18,1,15,10.0,2,1.0,350,6.0,7,20,32,1,,0.1
10321000,1991-06-28,17,25.5,487,18,365,6.8,8.5,167,,,163,144,222,83,0.65,0.58,0.01,0.050,0.010,0.031,,45,7.6,22,3.8,8,16,0.1,22,,,,,,,,,,,,,,,,,,,,,
10321000,1991-08-28,17.5,21,16,3,455,11.2,8.5,199,,,190,160,279,97,0.45,0.40,0.010,0.050,0.010,0.031,,44,12,34,6.8,17,28,0.6,22,6,95,0.50,1.0,1,3.0,1,6,1,14,10.0,1,1.0,450,6.0,13,30,43,1,,0.1
10321000,1991-11-01,4,8,32,2,458,12.3,8.5,208,,,194,181,326,76,0.26,0.18,0.02,0.063,0.010,0.031,,51,13,40,7.7,21,48,0.6,27,,110,,,,3.0,,7,,12,10.0,1.0,1.0,430,6.0,,10,63,1,,
10321000,1991-12-28,0.5,4,68,5,506,12.2,8.2,230,,,159,189,308,65,0.25,0.19,0.010,0.050,0.010,0.031,,56,12,45,7.1,20,47,0.5,25,,,,,,,,,,,,,,,,,,,,0.01,
10321000,1992-02-26,7,16.5,215,49,412,11.2,8.4,184,,,185,154,274,97,0.55,0.48,0.010,0.050,0.010,0.061,,46,9.4,32,7.2,16,27,0.5,25,,71,,,,3.0,,64,,8,10.0,2,1.0,290,6.0,,120,26,1,,
10321000,1992-04-28,19.5,29,149,13,432,9.4,8.5,180,,,178,142,275,79,0.35,0.28,0.010,0.050,0.010,0.031,,42,9,39,6.6,18,31,0.5,22,,71,,,,3.0,,10,,5,10.0,1,1.0,330,6.0,,20,30,1,,
10321000,1992-08-27,13,14,11,4,451,8.8,8.4,191,,,190,153,292,95,0.45,0.38,0.02,0.050,0.010,0.031,,38,14,41,8.3,20,33,0.6,23,,110,,,,3.0,,8,,7,10.0,1.0,1.0,480,6.0,,10,49,1,0.01,
10321000,1992-10-27,11,11.5,16,2,474,9.4,8.4,208,,,208,170,299,72,0.35,0.29,0.010,0.050,0.010,0.031,,48,12,36,7.7,15,30,0.5,20,,120,,,,3.0,,6,,8,10.0,1.0,1.0,500,6.0,,20,45,1,0.01,
10321000,1993-02-26,1,-7,40,4,543,11.8,8.4,213,,,245,175,365,64,0.44,0.26,0.04,0.13,0.01,0.031,,50,12,50,7.5,22,49,0.4,26,,96,,,,3.0,,8,,12,10.0,1,1.0,500,6.0,,10,51,1,0.01,
10321000,1993-04-29,12,9.5,485,18,452,8.4,8.4,155,,,149,144,274,86,0.65,0.58,0.02,0.050,0.010,0.123,,41,10,35,6.7,20,46,0.4,25,,72,,,,3.0,,16,,6,10.0,1,1.0,270,6.0,,20,23,1,0.01,
10321000,1993-06-16,15,8.5,967,16,392,7.5,8.4,163,,,166,146,247,60,0.55,0.48,0.02,0.050,0.010,0.123,,44,8.8,28,4.9,11,22,0.4,24,5,63,1.00,1.00,1,1.00,1,,1.00,3,3,1,1.00,,,3,22,,,0.01,
10321000,1993-08-25,19.5,21,26,2,486,8.6,8.4,187,,,191,158,289,72,0.55,0.48,0.02,0.050,0.010,0.031,,45,11,38,7,17,37,0.5,18,,92,,,,3.0,,6,,16,10.0,1,1.0,410,6.0,,20,39,1,0.01,
10321000,1993-10-26,8,14,36,1,595,11.4,8.6,210,,,171,188,349,62,0.25,0.18,0.02,0.050,0.010,0.031,,52,14,58,8.7,26,63,0.7,23,,110,,,,3.0,,6,,7,10.0,1.0,1.0,530,6.0,,10,53,1,0.01,
10321000,1994-02-09,0,2,66,3,573,12.3,8.6,232,,,228,199,367,43,0.45,0.38,0.02,0.050,0.010,0.031,,58,13,50,7.9,22,52,0.6,26,,99,,,,3.0,,5,,22,10.0,1.0,1.0,480,6.0,,10,40,1,0.01,
10321000,1994-04-28,9.5,11,231,8,409,9.8,8.4,175,,,168,145,248,61,0.45,0.39,0.01,0.050,0.010,0.031,,43,9,31,5.8,13,26,0.4,19,,67,,,,3.0,,13,,8,10.0,1,1.0,290,6.0,,10,22,1,0.01,
10321000,1994-06-07,15,20.5,690,11,363,,8.4,171,,,160,139,218,62,0.55,0.48,0.02,0.050,0.010,0.092,,43,7.6,22,4.1,7.6,16,0.3,20,4,56,0.50,1.0,1.0,3.0,2,19,1,5,10.0,1.0,1.0,240,6.0,3.0,40,13,1,0.01,
10321000,1994-09-01,18.5,23,6.2,5,446,8.1,8.6,193,,,185,152,276,66,0.45,0.38,0.02,0.050,0.010,0.245,,41,12,36,8.1,16,24,0.5,24,,130,,,,3.0,,65,,160,10.0,1,1.0,410,6.0,,50,42,1,0.01,
10321000,1994-10-27,9,20,15,6,476,10.8,8.5,212,,,203,167,284,,,,0.02,0.050,0.010,0.031,,47,12,36,7.3,15,26,0.5,17,8,130,0.50,1.0,1.0,3.0,1.0,14,1,58,10.0,1,1.0,480,6.0,4,10,39,1,0.01,0.1
10321000,1994-12-21,0,6,18,7,630,12.6,8.5,245,,,240,200,383,54,0.25,0.20,0.02,0.050,0.010,0.061,,57,14,55,9.9,23,56,0.7,27,,,,,,,,,,,,,,,,,,,,0.01,
10321000,1995-03-02,6,8.5,245,18,434,10.2,8.3,170,,,165,151,269,95,0.55,0.50,0.02,0.050,0.010,0.092,,44,10,34,6.1,16,35,0.4,24,,78,,,,3.0,,25,,13,10.0,1.0,1.0,310,6.0,,20,27,1,0.01,
10321000,1995-04-27,11.5,17.5,453,17,411,9.3,8.4,160,,,142,147,243,84,0.25,0.20,0.02,0.050,0.010,0.061,,44,9,32,5,14,30,0.5,23,,72,,,,4,,60,,5,10.0,1,1.0,260,7,,90,24,1,0.01,
10321000,1995-06-28,21,26,2180,32,370,6.5,8.1,161,,,154,135,217,65,0.65,0.60,0.02,0.050,0.010,0.153,,41,8,22,4.3,9.3,15,0.3,25,,,,,,,,,,,,,,,,,,,,0.01,
10321000,1995-08-22,22,24,80,2,390,9,8.6,165,,,156,141,238,73,0.35,0.30,0.02,0.050,0.010,0.031,,42,8.7,30,5.5,14,26,0.4,17,,71,,,,3.0,,11,,10,10.0,1.0,1.0,320,6.0,,20,21,1,0.01,
10321000,1995-10-27,7.5,5.5,74,6,481,10.2,8.3,200,,,174,165,294,,,,0.02,0.050,0.010,0.031,,46,12,38,6.5,16,37,0.5,21,7,88,0.50,1.0,1.0,3.0,1.0,13,1,14,10.0,1.0,1.0,400,6.0,3.0,20,31,1,0.01,0.1
10321000,1996-04-19,8,16.5,1100,35,370,11.2,8.5,147,,,142,117,224,,,,0.02,0.050,0.010,0.123,,34,7.8,27,4.4,13,28,0.4,24,4,60,0.50,1.0,1.0,3.0,12,58,1,4,10.0,2,1.0,210,6.0,8,70,17,1,0.01,0.1
10321000,1996-09-09,16,20.5,24,1,448,9.2,8.5,198,,,186,158,261,,,,0.02,0.06,0.010,0.031,,45,11,30,6.5,14,24,0.5,18,6,88,0.50,1.0,1.0,3.0,4,11,1,8,10.0,1.0,1.0,430,6.0,4,20,33,1,0.01,0.1
10321000,1996-11-22,,,,12,,,,210,,,,175,328,,,,0.02,0.09,0.010,0.061,,50,12,43,7,20,46,0.6,22,7,90,0.50,1.0,1.0,3.0,8,24,1,6,10.0,1.0,1.0,410,6.0,3.0,40,30,1,0.01,0.1
10321000,1997-04-17,13,21,698,34,416,8.9,8.3,159,,,146,135,248,,,,0.02,0.050,0.010,0.123,,38.3,9.6,32.3,5.4,14.9,32.5,0.5,26,6,70.5,0.50,1.0,1.0,3.0,1.5,3.0,1,3.19,13.5,1.0,1.0,252,7.5,3.0,5,19.5,1,0.01,0.1
10321000,1997-08-25,22,32,140,2,313,9.5,8.7,134,,,131,101,182,,,,0.02,0.050,0.010,0.113,,30.7,5.97,20.8,4.3,9,15.5,0.3,16.6,6,47.4,0.50,1.0,1.0,3.0,1.0,4.3,1,8.53,10.0,1.0,1.0,225,6.0,8.1,5,18.1,1,0.01,0.1
10321000,1997-12-17,0.5,0,66,4,585,12.2,8.4,229,,,229,194,360,,,,0.02,0.050,0.010,0.055,,55.8,13.1,46.9,6.9,21,50.1,0.5,28.5,7,94.9,1.00,8.0,1.0,12.0,1.0,10.0,1,11.2,60.0,1.0,1.0,409,10.0,20.0,10,35.4,1,0.01,0.1
10321000,1998-04-07,7,6,864,44,413,,8.8,159,,,162,141,262,,,,0.02,0.067,0.010,0.159,,39,10.4,33,5.4,15.7,34.7,0.3,26.2,6,74.2,1.00,8.0,1.0,12.0,3.8,18.8,1,4.7,60.0,1.0,1.0,254,10.0,20.0,10,21,1,0.01,0.1
10321000,1998-08-27,18,24.5,62,2,532,9.7,8.5,212,,,200,170,314,,,,0.02,0.050,0.010,0.043,,47.8,12.2,44.2,7.7,19.2,40.3,0.6,21.6,10,97.4,1.00,8.0,1.0,12.0,1.9,10.0,1,9,60.0,1,1.0,401,10.0,20.0,10,37.2,1,0.01,0.1
10321000,1998-08-27,,,,2,,,,212,,,,170,321,,,,0.02,0.050,0.010,0.043,,47.8,12.2,44.2,7.7,19.2,40.3,0.6,21.6,10,97.4,1.00,8.0,1.0,12.0,1.9,10.0,1,9,60.0,1,1.0,401,10.0,20.0,10,37.2,1,0.01,0.1
10321000,2015-05-01,16.5,,80,,390,8.1,8.7,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
10321000,2015-06-26,21.6,,94,,422,3,8.2,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
-------------- next part --------------
#LyX 2.1 created this file. For more info see http://www.lyx.org/
\lyxformat 474
\begin_document
\begin_header
\textclass scrreprt
\begin_preamble
\date{}
\usepackage{textcomp,url,multicol}
%\setkomafont{sectioning}{\rmfamily}
\end_preamble
\options abstract=on
\use_default_options false
\begin_modules
natbibapa
knitr
\end_modules
\maintain_unincluded_children false
\language english
\language_package default
\inputencoding auto
\fontencoding global
\font_roman palatino
\font_sans default
\font_typewriter default
\font_math auto
\font_default_family default
\use_non_tex_fonts false
\font_sc false
\font_osf false
\font_sf_scale 100
\font_tt_scale 100
\graphics default
\default_output_format default
\output_sync 0
\bibtex_command bibtex
\index_command default
\paperfontsize default
\spacing single
\use_hyperref false
\papersize letterpaper
\use_geometry false
\use_package amsmath 1
\use_package amssymb 1
\use_package cancel 1
\use_package esint 0
\use_package mathdots 1
\use_package mathtools 1
\use_package mhchem 1
\use_package stackrel 1
\use_package stmaryrd 1
\use_package undertilde 1
\cite_engine natbib
\cite_engine_type authoryear
\biblio_style humannat
\use_bibtopic false
\use_indices false
\paperorientation portrait
\suppress_date false
\justification true
\use_refstyle 0
\index Index
\shortcut idx
\color #008000
\end_index
\secnumdepth 2
\tocdepth 2
\paragraph_separation indent
\paragraph_indentation default
\quotes_language english
\papercolumns 1
\papersides 2
\paperpagestyle default
\tracking_changes false
\output_changes false
\html_math_output 0
\html_css_as_file 0
\html_be_strict false
\end_header

\begin_body

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

load('.RData')
\begin_inset Argument 1
status open

\begin_layout Plain Layout
echo=FALSE
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_layout Standard
The first step in analyzing water chemistry data for CWA compliance is reading
 it into the analytical software.
\end_layout

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

carlin <- read.csv("./carlin.csv", header = TRUE, sep = ",", stringsAsFactors
 = F)
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Next, check that the data are what you expect to see and convert dates from
 factors.
 
\end_layout

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

carlin$sampdate <- as.Date(carlin$sampdate)
\end_layout

\begin_layout Plain Layout

str(carlin)
\end_layout

\end_inset


\end_layout

\begin_layout Standard
The above table shows the structure of the data in a R data.frame format.
 
\end_layout

\begin_layout Standard
Another perspective on the data is the summary of each of the data columns.
\end_layout

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

summary(carlin)
\end_layout

\end_inset

 
\end_layout

\begin_layout Standard
Plots of these distributions are the next step because the graphic conveys
 all the written information, plus much more insight into the data characteristi
ecs, more easily for decision-makers.
\end_layout

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

carlin.1 <- subset(carlin, select = siteid:CO3)
\end_layout

\begin_layout Plain Layout

carlin.2 <- subset(carlin, select = c(siteid, sampdate, Alk:Ca))
\end_layout

\begin_layout Plain Layout

carlin.3 <- subset(carlin, select = c(siteid, sampdate, Mg:Cr))
\end_layout

\begin_layout Plain Layout

carlin.4 <- subset(carlin, select = c(siteid, sampdate, Co:Hg))
\end_layout

\end_inset


\end_layout

\begin_layout Standard
To plot each chemical constituent's concentration as a function of collection
 date the data format needs to be reshaped from wide to long:
\end_layout

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

carlin.1.melt <- melt(carlin.1, na.rm = F, id.vars = c('siteid', 'sampdate'))
\end_layout

\begin_layout Plain Layout

carlin.2.melt <- melt(carlin.2, na.rm = F, id.vars = c('siteid', 'sampdate'))
\end_layout

\begin_layout Plain Layout

carlin.3.melt <- melt(carlin.3, na.rm = F, id.vars = c('siteid', 'sampdate'))
\end_layout

\begin_layout Plain Layout

carlin.4.melt <- melt(carlin.4, na.rm = F, id.vars = c('siteid', 'sampdate'))
\end_layout

\end_inset


\end_layout

\begin_layout Standard
The X-Y plots use only measured data; missing data are not included.
\end_layout

\begin_layout Standard
\begin_inset Flex Chunk
status open

\begin_layout Plain Layout

pdf('carlin-1-descriptive.pdf')
\end_layout

\begin_layout Plain Layout

print(xyplot(value ~ sampdate | variable, data=carlin.1.melt, rm.na = T))
\end_layout

\end_inset


\end_layout

\end_body
\end_document