BIOL 339 CUNY Borough of Manhattan Community College Data Manipulation Worksheet

BIOL 339 CUNY Borough of Manhattan Community College Data Manipulation Worksheet
Sample Answer for BIOL 339 CUNY Borough of Manhattan Community College Data Manipulation Worksheet Included After Question
Description

has two scientific goals: 1) to learn how to summarize data and 2) to help you gain a better sense of the terrestrial weather variables that impact the behavior and physiology of animals. As we pursue these goals, we will add to our knowledge of RStudio, specifically improving how to make figures that look more professional (sources -Chapter 3 of R for Data Science and ggplot book, https://ggplot2-book.org/index.html) and how to manipulate data frames so we can make the visualizations that are needed (material in chapter 5 of R for Data Science). This homework requires one figure with a legend (caption) in Part 1 and 2 figures with legends for Part 2.  Uploading your R code with the document containing the figures is often helpful. Below are more details about the assignment.
BIOL 339 CUNY Borough of Manhattan Community College Data Manipulation Worksheet
Part 1. This first part is to give you experience with the dplyr options to filter rows, arrange rows, select columns, create a new column with mutate and summarize subsets of the data frame.  The book chapter 5 and the code I have uploaded with the homework give examples to try.  The requirement here is to use the group_by() statement to group data by day, week or month, calculate the mean of one variable and one at least one other statistical variable such as mode, median, max, min, standard deviation and standard error then plot the data.  If you want to use standard deviation or standard error look at examples using geom_error to make the figures. In addition to including a legend, compare you compare the figure you have make to a standard boxplot and comment about which you prefer and why.
Stevenson 2023 HW04 Graph interpretations & Data manipulation Comparative Animal Physiology Biology 337-339 Due Monday Feb 13 Midnight – upload the assignment to Blackboard HW04 has two scientific goals: 1) to learn how to summarize data and 2) to help you gain a better sense of the terrestrial weather variables that impact the behavior and physiology of animals. As we pursue these goals, we will add to our knowledge of RStudio, specifically improving how to make figures that look more professional (sources -Chapter 3 of R for Data Science and ggplot book, https://ggplot2book.org/index.html) and how to manipulate data frames so we can make the visualizations that are needed (material in chapter 5 of R for Data Science). This homework requires one figure with a legend (caption) in Part 1 and 2 figures with legends for Part 2. Uploading your R code with the document containing the figures is often helpful. Below are more details about the assignment. Part 1. This first part is to give you experience with the dplyr options to filter rows, arrange rows, select columns, create a new column with mutate and summarize subsets of the data frame. The book chapter 5 and the code I have uploaded with the homework give examples to try. The requirement here is to use the group_by() statement to group data by day, week or month, calculate the mean of one variable and one at least one other statistical variable such as mode, median, max, min, standard deviation and standard error then plot the data. If you want to use standard deviation or standard error look at examples using geom_error to make the figures. In addition to including a legend, compare you compare the figure you have make to a standard boxplot and comment about which you prefer and why. Part 2. The relationships among terrestrial weather variables. Here you need to make two or more professional looking figures (plot + figure legend) that enlighten the relationship among variables. Weather is driven by a series of processes that are now so well understood and modeled that meteorologists can make accurate short- and long-term forecasts. The goal here is much simpler – to make sense of patterns we see in the data from the Plum Island LTER site. We can use our everyday experiences as our guide. To orient yourself, remember there are three prominent time scales of weather change. (Ignoring other drivers for now.) the daily, annual and high and low cycles of weather cycle changes (days to weeks). 1) The daily cycle of solar input driven by from the rotation of the earth gives us the daily cycle of solar energy input 2) The tilt in the earth and its 3651/4 day long trip around the sun are largely responsible for giving us the annual cycle. The earth’s tilt and path around the solar determine the number of daylight hours and the path of the sun through the sky. It reaches its highest point above the horizon at solar noon each day and its highest point in the year at the summer solstice in June (19-21th each year). These variations are well studied and change in a very predictable way. See simulations ( https://ccnmtl.github.io/astro-simulations/sun-motion-simulator/ , http://andrewmarsh.com/apps/staging/sunpath3d.html. 3) The uneven heating of the earth by the sun causes thermal gradients that drive the movement of air masses in the atmosphere (fronts, cyclones, the jet stream, etc.) and water currents in the ocean together causing our weather. These forces bring complex exposures to cold and hot and wet and dry air masses that circulate across the globe. In Massachusetts these changes usually occur over several days or weeks. Despite the complexity patterns of weather caused by the three drivers listed above, experience and common sense teaches there are many known relationships. Here are some you may have thought about, heard or observed before. Winds usually blow from the west. Winds from the west and south bring warmer air and winds from the north bring cooler air. When we get a nor’easter storm strong winds off the ocean bring moisture as rain or snow. Precipitation is associated with clouds and thus low sunshine or radiation. Precipitation is associated with high relative humidity. Relative humidity goes up at night and can lead to the formation of dew if it gets to 100%. The rise of air temperature lags behind the solar radiation. The maximum daily air temperature is usually several hours later than the maximum radiation and the maximum average daily temperatures occur several months after the maximum daily solar input. High pressures are associated with clear skies and sunny weather. Falling pressures are an indication of storms coming. Solar radiation data (Pyranometer readings) correspond closely with Photosynthetically Active Radiation (PAR) readings. Based on the ideas of relationships stated in the previous paragraph or ideas of your own, create a minimum of two figures (more are fine) to test the relationships. Are the relationships true in the data set we are using? Try to make the figures as much like figures in published papers as possible. Successful figures communicate the results clearly and in an aesthetically pleasing manner. Based on ggplot’s grammar of graphics this means making the data visible and using the theme’s functions to communicate your information clearly. From among the research papers you have read for the course, find and study two or three figures that were published in recent years. Study the graphic elements of the figures. (Note the background style, lack of titles, the numerical range of axes, whether or not they use log axes, the units of the axis labels, the location and relative sizes of the tick and axis labels, the number of panels in the figure, etc.). Use these published figures as guides to what you create. For the overall look use themes such as theme_classic(), theme_bw() and theme_prism() (you need the ggprism library installed and loaded for theme_prism() to work). Making labels with exponents, Greek letter and subscripts can be done. I put an example in the code. http://abacus.bates.edu/~ganderso/biology/resources/writing/HTWtablefigs.html All axis parameters such as the number of tick marks, the number, size and font of the tick mark labels and the labels can be specified if the default do not work well. One important suggestion is to use xlab(0 and ylab() to make your x and y labels be to set the background. If you cannot get elements of the figure the way you want, please comment on your homework what you could not do. One element that we have touch on in the first R homework were axis labels. You can specific axis labels using ylab() and xlab() or labs(y= “y axis name”, x = “x axis name”). You can change the font size with a statement such as + theme(axis.text.x = element_text(size = 20)) or with more terms You also need to write a clear figure legend (called a caption sometimes) for each figure with three objectives in mind 1) Describes that relationship between the variables in the plot. Describe in words what the viewer should be able to see if they had not seen the plot. 2) Make sure to include units and numerical values to allow the reader to compare what is in your words with what is in the plot 3) Provide an interpretation of the data for the reader. What is the take home message? Are the patterns in the plot standard or unexpected? Plum Island Ecosystems LTER DATABASE DATASET ID: Metacat Package ID Year Released to Public Distribution URL for file DATASET TITLE: MON-PR-Met15Min2019.01 knb-lter-pie.542.1 2020 http://ecosystems.mbl.edu/PIE/data-archive/MON/MON-PR-Met15Min2019.html PIE LTER year 2019, meteorological data, 15 minute intervals, from the PIE LTER Marshview Farm weather station located in Newbury, MA INVESTIGATOR INFORMATION : Investigator 1 First Name Anne Last Name Giblin Address line 1 Ecosystems Center Address line 2 MBL Address line 3 7 MBL St City Woods Hole State MA Zip Code 02543 Country USA OTHERS: Investigator 2 Hap Garritt DATA FILE INFORMATION: Data File URL http://ecosystems.mbl.edu/PIE/data-archive/MON/data/MON-PR-Met15min2019.csv Data File Name MON-PR-Met15min2019 Beginning Date 01-Jan-2019 End Date 31-Dec-2019 Number of Data Records 35040 Other Files to Reference Availability Status Type 1 Quality Control Information Maintenance Description Final data set for Year 2019 Log of Changes: Version 01: January 29, 2020, data and metadata updates to comply with importation to DEIMS7 and LTER Data Portal. Used MarcrosExportEML_HTML (working)pie_excel2007_Sep2019.xlsm 9/3/19 5:19 PM for QA/QC to EML 2.1.0 MON-PR-MFM Met Station RESEARCH LOCATION: Geographic Description PIE meteorological station at MBL, Marshview Farm, Newbury, MA. Location Bounding Box West Bounding Coordinate East Bounding Coordinate North Bounding Coordinate South Bounding Coordinate OR if single point location Latitude 42.756971 Longitude -70.891366 Elevation TAXONOMIC COVERAGE: Taxonomic Protocols Organisms studied Site 2 Example: Spartina; Spartina patens; Carex aquatilis var. aquatilis; Carex KEYWORD INFORMATION KEYWORDS: PIE LTER, primary production, Massachusetts, Newbury, climate, meteorology, temperature, precipitation, rain, humidity, solar, radiation, P KeywordThesaurus ABSTRACT: Year 2019 meteorological measurements at MBL Marshview Farm of air temperature, humidity, precipitation, solar radiation, pho active radiation (PAR), wind speed and direction and barometric pressure. Sensors conduct measurements every 5 secs and meas reported as averages or totals for 15 minute intervals. 15 minute averages are reported for air temperature, humidity, solar radiation, PAR, wind speed and direction and barometric 15 minute totals are reported for precipitation. METHODS: The data file is a summary of daily weather data from the PIE LTER weather station located at MBL Marshview Farm in Newbury, The weather station is solar powered using a Solarex MSX-20 20 watt solar module including a 12V charger and regulator and 12 AH battery. The system was supplied by Campbell Scientific Inc., The data logger is supplied by Campbell Scientific CR10X-1M data logger with 1 megabyte memory storage module and PC208W Windows logging software. The enclosure and tower is a Campbell Scientific ENC 12/14 (12” X 14” weather proof enclosure), CM10 (10’ tripod, grounding k kit. Data telemetry/Communication at Remote Site/Weather station, Governor’s Academy – RF 310M, RF Modem for use w/Maxon radios – RF 310, Maxon radio for 148-174 MHZ – 158.34 MHZ – Yagi antenna – Antennex Y1503 3DB 150-174 MHz at Base station, Rowley field station – RF 310B Base station for use w/Maxon radio – FG1563, Omni directional antenna 155-162 MHZ, 3 DB gain Sensors used: Temperature and Humidity Vaisala HMP45C temperature and humidity sensor with solar radiation shield. Specifications: +/- 1% accuracy over Precipitation gauge Texas Electronic TE525WS-L, 8” rain gage with CS705 precipitation adapter for snow fall. Specifications: +/- 1% accuracy up to Pyranometer Solar radiation, 400-1100 nm, Licor LI200X-L pyranometer, mounted at 3 meter height Specifications: +/- 3% typical accuracy Solar radiation, 400-1100 nm, Licor LI200X-L pyranometer, mounted at 3 meter height Specifications: +/- 3% typical accuracy PAR Photosynthetically active radiation, 400-700 nm, Licor LI190SB quantum sensor, mounted at 3 meter height Specifications: +/- 5% calibration, National Institute of Standards Technology Wind speed and direction RM Young 05103 Wind monitor, mounted at 3 meter height Barometric pressure Vaisala PTB110 Specifications: +/- 1 mB @ -20oC to +45o C. April 16, 2019 Precipitation snow adapter removed May 22, 2019 Temp/RH HMP45C radiation shield cleaned of debris lodged between panels, PTB101B barometer sensor repalced with December 6, 2019 Precipitation snow adapter put on for Winter season Sampling and/or Lab Protocols Protocol Title URL of online Protocol OR Protocol Document VARIABLE DESCRIPTIONS: Variable Name Variable Description Units Date Time date (YYYY-MM-DD) time of day (hh:mm), (eastern standard time EST), 24 hour format Julian julian day number Precip total precipitation, wet only, rain and snow (mm) millimeter Pyranometer solar radiation (kW m-2), 15 minute average kilowattPerMeterSquared -2 -1 PAR photosynthetically active radiation (µM m sec ), 15 minute average micromolePerMeterSquaredPerSecond Temp air temperature (Celsius), 15 minute average celsius RH relative humidity (percent), 15 minute average percent Wind wind speed (m sec-1), 15 minute average meterPerSecond WindDir wind direction (degrees), 15 minute average degree BAR barometer (mbar), 15 minute average millibar Comments Comments, notes about data Investigator 3 Site 3 ns; Carex aquatilis var. aquatilis; Carex atlantica ssp. atlantica ecipitation, solar radiation, photosynthetically urements every 5 secs and measurements are and direction and barometric pressure. Marshview Farm in Newbury, MA. charger and regulator and 12 AH lead acid storage module and PC208W Windows data M10 (10’ tripod, grounding k it and CM10 guy wire % accuracy over – 40 to 60 °C temperature range tions: +/- 1% accuracy up to 1 in./hr. ter height 1B barometer sensor repalced with PTB110 datetime Missing DateTim Value e Format Code YYYYMM-DD datetime hh:mm Measurement Scale Code Number Informati Type on interval integer NA ratio real NA ratio real NA ratio real NA interval real NA ratio real NA ratio real NA interval real NA ratio real NA nominal NA Missing Value Code Explanation Do Not Modify. This is the lists for the drop MeasurementScale Number Type Unit Name NA = not available NA = not available NA = not available NA = not available NA = not available NA = not available NA = not available NA = not available NA = not available NA = not available datetime integer ampere interval natural amperePerMeter nominal real amperePerSquareMeter ordinal whole angstrom ratio atmosphere bar calorie celsius centigram centimeter centimeterPerYear centimetersPerSecond centimolesChargePerKi centisecond coulomb cubicCentimetersPerCu cubicFeetPerSecond cubicMeter cubicMeterPerKilogram cubicMetersPerSecond cubicMicrometersPerGr decibar decigram decimeter decisecond degree dimensionless disintegrationsPerMinute farad gram gramsPer0.04SquareMe gramsPerCentimeterSqu gramsPerCubicCentime gramsPerGram gramsPerHectarePerDa gramsPerLiter gramsPerLiterPerDay gramsPerMeterSquared gramsPerMeterSquared gramsPerMilliliter gramsPerNumber gramsPerSquareMeter gramsPerYear gray hectare hectoPascal henry hertz hour joule joulesPerCentimeterSqu joulesPerCentimeterSqu kelvin kilogram kilogramPerCubicMeter kilogramsPerHectare kilogramsPerHectarePe kilogramsPerMeterSqua kilogramsPerMeterSqua kilogramsPerSecond kilogramsPerSquareMet kilohertz kiloliter kilometer kilometersPerHour kilopascal kilosecond kilovolt kilowatt kilowattPerMeterSquare liter litersPerHectare litersPerSecond litersPerSquareMeter lumen lux megagram megahertz megameter megapascal megasecond megavolt megawatt meter metersPerDay metersPerGram metersPerSecond metersPerSecondSquar metersSquaredPerDay metersSquaredPerSeco microCuriePerMicroMole microEinsteinsPerSquar microEinsteinsPerSquar microequivalentsPerLite microgram microgramsPerCubicCe microgramsPerGram microgramsPerLiter microgramsPerMilliliter microliter micrometer microMolesCarbonPerM microMolesPerKilogram microMolesPerLiter microMolesPerSquareM microMolesPerSquareM micron microsecond microsiemensPerCentim millibar milligram milliGramPerKilogram milliGramPerSegment milligramsPerCubicMete milligramsPerLiter milliGramsPerMilliLiter milligramsPerMillimeter milligramsPerSquareMe milligramsPerSquareMe millihertz milliliter milliliterPerLiter millimeter millimetersPerNumber millimetersPerSecond millimetersPerSegment millimolesPerGram millimolesPerMole millimolesPerSquareMe millisecond millivolt milliwatt minute molality molarity mole molePerCubicMeter molesPerGram molesPerKilogram molesPerKilogramPerSe molesPerMeterSqaureP nanogram nanometer nanomolesPerGramPer nanosecond newton nominalDay nominalYear number numberPerGram numberPerKilometerSqu numberPerMeterCubed numberPerMeterSquare numberPerMilliliter numberPerSquareCenti ohm ohmMeter partsPerMillion partsPerThousand pascal percent picoMolesPerLiter picoMolesPerLiterPerHo PSU radian second serialDateNumberYear0 siemen siemensPerMeter squareCentimeters squareCentimetersPer0 squareCentimetersPerG squareKilometers squareMeter squareMeterPerKilogram squareMeterPerNumber squareMeterPerSquareM squareMillimeters volt watt s the lists for the drop-downs. ubicCentimetersPerCubicCentimeters ramsPerCentimeterSquaredPerSecond ramsPerMeterSquaredPerYear oulesPerCentimeterSquaredPerDay oulesPerCentimeterSquaredPerHour logramsPerMeterSquaredPerSecond logramsPerMeterSquaredPerYear microEinsteinsPerSquareMeterPerSecond microMolesCarbonPerMicroMolePhoton microMolesPerSquareMeterPerMinute microMolesPerSquareMeterPerSecond milligramsPerSquareMeterPerDay millimolesPerSquareMeterPerHour umberPerSquareCentimeterPerHour quareCentimetersPer0.04SquareMeter Date Time Julian Precip Pyranometer PAR Temp RH Wind WindDir BAR 2019-01-01 00:15 1 1.016 0.000 0.174 2.52 102 2.06 141.3 1012 2019-01-01 00:30 1 1.016 0.000 0.806 3.12 102 2.52 146.2 1012 2019-01-01 00:45 1 1.270 0.000 0.523 3.63 102 3.20 142.7 1012 2019-01-01 01:00 1 1.016 0.000 0.834 3.99 101 2.49 139.4 1011 2019-01-01 01:15 1 1.524 0.000 0.273 4.19 101 2.53 142.2 1010 2019-01-01 01:30 1 1.016 0.000 0.042 4.34 101 1.90 140.3 1010 2019-01-01 01:45 1 0.762 0.000 0.000 4.64 101 2.76 144.1 1009 2019-01-01 02:00 1 0.508 0.000 0.000 4.86 101 2.14 137.8 1008 2019-01-01 02:15 1 0.762 0.000 0.000 5.05 101 2.22 139.7 1007 2019-01-01 02:30 1 0.508 0.000 0.000 5.29 101 2.80 145.1 1005 2019-01-01 02:45 1 0.254 0.000 0.000 5.46 102 2.79 134.9 1004 2019-01-01 03:00 1 0.000 0.000 0.000 5.63 102 2.23 149.4 1003 2019-01-01 03:15 1 0.254 0.000 0.000 5.84 102 2.70 141.6 1003 2019-01-01 03:30 1 0.000 0.000 0.221 6.02 102 2.55 146.1 1002 2019-01-01 03:45 1 0.000 0.000 0.584 6.17 102 1.88 148.4 1001 2019-01-01 04:00 1 0.254 0.000 0.000 6.35 102 2.16 160.1 1001 2019-01-01 04:15 1 0.000 0.000 0.000 6.49 102 2.72 150.2 1000 2019-01-01 04:30 1 0.254 0.000 0.000 6.53 102 1.40 164.2 999 2019-01-01 04:45 1 0.254 0.000 0.000 6.55 102 0.88 204.5 999 2019-01-01 05:00 1 0.254 0.000 0.000 6.53 102 0.65 183.8 998 2019-01-01 05:15 1 0.000 0.000 0.000 6.48 102 0.69 206.7 997 2019-01-01 05:30 1 0.000 0.000 0.000 6.44 102 0.64 241.5 997 2019-01-01 05:45 1 0.254 0.000 0.000 6.34 102 0.52 232.7 996 2019-01-01 06:00 1 0.000 0.000 0.000 6.29 103 0.92 241.5 996 2019-01-01 06:15 1 0.000 0.000 0.042 6.27 103 1.06 230.6 996 2019-01-01 06:30 1 0.000 0.000 0.085 6.41 103 1.16 257.5 996 2019-01-01 06:45 1 0.000 0.000 0.019 6.44 103 1.18 257.6 996 2019-01-01 07:00 1 0.000 0.000 0.118 6.34 103 0.78 241.5 996 2019-01-01 07:15 1 0.000 0.001 1.728 5.97 103 1.44 233.6 996 2019-01-01 07:30 1 0.000 0.004 8.980 5.58 103 1.21 245.5 996 2019-01-01 07:45 1 0.000 0.011 25.620 5.16 103 1.49 230.3 996 2019-01-01 08:00 1 0.000 0.013 35.890 4.93 103 1.64 255.3 996 2019-01-01 08:15 1 0.000 0.017 56.350 4.76 103 1.31 247.3 996 2019-01-01 08:30 1 0.000 0.208 357.400 5.11 103 1.82 248.6 996 2019-01-01 08:45 1 0.000 0.060 114.800 5.71 103 1.52 245.5 996 2019-01-01 09:00 1 0.000 0.030 72.500 6.05 103 0.81 241.7 996 2019-01-01 09:15 1 0.000 0.100 219.800 6.75 103 1.17 253.0 996 2019-01-01 09:30 1 0.000 0.201 400.500 7.71 103 1.02 260.5 997 2019-01-01 09:45 1 0.000 0.191 384.600 8.86 103 1.90 234.2 997 2019-01-01 10:00 1 0.254 0.340 623.700 9.71 98 2.92 252.9 997 2019-01-01 10:15 1 0.508 0.362 672.700 10.1 88 3.27 245.2 997 2019-01-01 10:30 1 0.000 0.376 707.000 10.9 80 3.01 251.2 997 2019-01-01 10:45 1 0.254 0.390 737.000 11.7 74 3.55 250.2 997 2019-01-01 11:00 1 0.254 0.401 767.000 12.4 69 3.06 256.3 998 2019-01-01 11:15 1 0.254 0.412 805.000 12.7 67 2.97 264.3 998 2019-01-01 11:30 1 0.254 0.427 863.000 12.7 64 3.34 276.1 998 2019-01-01 11:45 1 0.254 0.415 846.000 12.6 63 3.71 285.8 999 2019-01-01 12:00 1 0.000 0.315 659.100 12.5 62 3.90 283.4 999 2019-01-01 12:15 1 0.508 0.315 663.000 12.1 59 4.68 280.1 1000 2019-01-01 12:30 1 0.000 0.400 821.000 12.1 58 3.67 287.2 1000 2019-01-01 12:45 1 0.000 0.453 909.000 12.2 58 4.10 297.7 1001 2019-01-01 13:00 2019-01-01 13:15 2019-01-01 13:30 2019-01-01 13:45 2019-01-01 14:00 2019-01-01 14:15 2019-01-01 14:30 2019-01-01 14:45 2019-01-01 15:00 2019-01-01 15:15 2019-01-01 15:30 2019-01-01 15:45 2019-01-01 16:00 2019-01-01 16:15 2019-01-01 16:30 2019-01-01 16:45 2019-01-01 17:00 2019-01-01 17:15 2019-01-01 17:30 2019-01-01 17:45 2019-01-01 18:00 2019-01-01 18:15 2019-01-01 18:30 2019-01-01 18:45 2019-01-01 19:00 2019-01-01 19:15 2019-01-01 19:30 2019-01-01 19:45 2019-01-01 20:00 2019-01-01 20:15 2019-01-01 20:30 2019-01-01 20:45 2019-01-01 21:00 2019-01-01 21:15 2019-01-01 21:30 2019-01-01 21:45 2019-01-01 22:00 2019-01-01 22:15 2019-01-01 22:30 2019-01-01 22:45 2019-01-01 23:00 2019-01-01 23:15 2019-01-01 23:30 2019-01-01 23:45 2019-01-01 00:00 2019-01-02 00:15 2019-01-02 00:30 2019-01-02 00:45 2019-01-02 01:00 2019-01-02 01:15 2019-01-02 01:30 2019-01-02 01:45 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 0.000 0.000 0.254 0.254 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.432 0.348 0.138 0.231 0.138 0.134 0.111 0.071 0.057 0.146 0.099 0.036 0.015 0.009 0.003 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 856.000 699.300 302.400 489.500 311.900 299.100 262.100 186.600 156.600 279.600 192.400 90.800 51.410 27.380 7.680 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 12.1 11.6 10.6 10 9.71 9.31 8.55 7.95 7.45 7.15 6.93 6.27 5.79 5.5 5.11 4.74 4.43 4.21 4.03 3.86 3.68 3.59 3.49 3.47 3.25 2.98 2.64 2.26 1.98 1.7 1.53 1.27 1.13 1.1 1.02 0.86 0.69 0.43 0.23 -0.08 -0.27 -0.45 -0.66 -0.95 -1.17 -1.29 -1.61 -1.64 -1.71 -1.92 -2.08 -2.25 58 58 60 60 59 60 58 58 57 58 58 58 59 58 58 59 59 60 62 63 64 64 64 63 63 64 64 62 62 63 64 65 65 64 61 60 60 61 61 61 61 61 62 63 63 63 65 62 61 61 61 62 4.26 5.11 5.92 5.31 4.56 3.59 4.35 4.74 5.90 4.65 4.29 4.77 4.70 4.97 4.68 4.24 4.21 4.18 4.05 3.44 3.45 3.62 3.03 4.11 3.68 2.62 3.32 4.15 3.46 2.62 2.27 1.70 2.81 3.09 3.55 2.87 2.14 2.40 2.29 1.88 1.97 2.01 2.13 1.82 1.68 1.33 1.23 2.55 2.13 1.97 1.86 2.04 297.6 1001 310.1 1002 305.9 1002 305.5 1003 299.4 1003 304.9 1004 302.1 1004 302.5 1005 292.2 1005 296.0 1006 309.5 1006 306.0 1007 309.6 1007 303.0 1008 299.5 1008 305.6 1009 305.3 1009 307.1 1010 303.3 1010 304.2 1011 302.9 1011 300.5 1011 301.4 1012 301.4 1012 296.3 1012 303.7 1013 303.3 1014 305.8 1014 307.7 1014 298.9 1014 310.1 1015 298.2 1015 301.6 1015 312.2 1015 306.2 1015 294.6 1016 297.8 1016 292.3 1016 311.1 1017 293.2 1017 287.2 1017 289.5 1017 300.7 1017 290.9 1017 297.0 1018 307.9 1018 297.5 1018 303.2 1018 307.6 1019 309.1 1019 297.1 1019 306.2 1020 2019-01-02 02:00 2019-01-02 02:15 2019-01-02 02:30 2019-01-02 02:45 2019-01-02 03:00 2019-01-02 03:15 2019-01-02 03:30 2019-01-02 03:45 2019-01-02 04:00 2019-01-02 04:15 2019-01-02 04:30 2019-01-02 04:45 2019-01-02 05:00 2019-01-02 05:15 2019-01-02 05:30 2019-01-02 05:45 2019-01-02 06:00 2019-01-02 06:15 2019-01-02 06:30 2019-01-02 06:45 2019-01-02 07:00 2019-01-02 07:15 2019-01-02 07:30 2019-01-02 07:45 2019-01-02 08:00 2019-01-02 08:15 2019-01-02 08:30 2019-01-02 08:45 2019-01-02 09:00 2019-01-02 09:15 2019-01-02 09:30 2019-01-02 09:45 2019-01-02 10:00 2019-01-02 10:15 2019-01-02 10:30 2019-01-02 10:45 2019-01-02 11:00 2019-01-02 11:15 2019-01-02 11:30 2019-01-02 11:45 2019-01-02 12:00 2019-01-02 12:15 2019-01-02 12:30 2019-01-02 12:45 2019-01-02 13:00 2019-01-02 13:15 2019-01-02 13:30 2019-01-02 13:45 2019-01-02 14:00 2019-01-02 14:15 2019-01-02 14:30 2019-01-02 14:45 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.004 0.013 0.022 0.044 0.051 0.092 0.240 0.261 0.293 0.322 0.348 0.372 0.395 0.407 0.421 0.429 0.444 0.439 0.435 0.433 0.419 0.378 0.387 0.325 0.375 0.286 0.261 0.192 0.156 0.131 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 3.232 10.980 27.770 49.440 94.000 117.400 192.100 415.100 467.100 532.800 593.100 648.500 698.200 744.000 771.000 800.000 819.000 848.000 843.000 836.000 833.000 809.000 736.000 749.000 637.500 723.000 567.500 514.600 383.400 316.400 265.400 -2.38 -2.52 -2.72 -2.84 -3 -3.13 -3.17 -3.24 -3.33 -3.49 -3.62 -3.69 -3.77 -3.9 -4.04 -4.14 -4.2 -4.21 -4.19 -4.19 -4.23 -4.24 -4.31 -4.31 -4.22 -4.12 -4.01 -3.98 -3.31 -2.62 -2.19 -1.86 -1.53 -1.16 -0.98 -0.53 -0.17 0.02 0.1 0.17 0.48 0.62 0.65 1.14 1.34 1.14 1.1 1.18 1.18 1.32 1.08 0.83 62 62 62 61 61 61 60 61 61 61 60 61 61 60 61 61 62 61 62 61 61 61 61 60 58 57 56 55 51 48 44 42 39 36 34 32 31 29 30 28 28 27 27 28 27 27 26 24 24 25 25 26 2.21 2.18 2.09 2.43 1.66 1.56 2.36 1.76 1.86 2.24 2.74 2.12 1.81 1.84 1.80 1.62 1.40 1.74 1.28 1.96 1.99 1.75 1.55 1.56 1.37 1.24 0.70 0.79 1.35 1.24 1.45 1.76 1.30 1.85 1.19 1.44 1.33 1.56 2.02 1.89 1.47 1.81 1.65 1.14 1.22 1.68 1.35 1.05 0.95 1.13 1.04 0.70 310.8 1020 307.4 1020 303.9 1020 305.9 1020 290.0 1020 296.5 1021 299.6 1021 301.7 1021 300.7 1021 297.2 1021 313.1 1021 305.2 1022 308.4 1022 311.3 1022 305.6 1022 313.7 1023 301.2 1023 309.1 1024 301.1 1024 312.7 1025 305.6 1025 308.6 1025 313.7 1025 313.1 1025 297.7 1025 309.1 1025 272.6 1026 287.3 1026 298.4 1026 295.5 1026 289.1 1026 270.0 1026 247.7 1026 153.9 1026 288.4 1026 267.5 1026 271.9 1025 274.2 1025 241.6 1025 284.7 1025 268.2 1025 257.3 1025 277.7 1025 263.7 1025 266.3 1025 284.0 1025 271.7 1025 216.7 1024 265.9 1024 248.9 1024 233.5 1024 246.5 1024 2019-01-02 15:00 2019-01-02 15:15 2019-01-02 15:30 2019-01-02 15:45 2019-01-02 16:00 2019-01-02 16:15 2019-01-02 16:30 2019-01-02 16:45 2019-01-02 17:00 2019-01-02 17:15 2019-01-02 17:30 2019-01-02 17:45 2019-01-02 18:00 2019-01-02 18:15 2019-01-02 18:30 2019-01-02 18:45 2019-01-02 19:00 2019-01-02 19:15 2019-01-02 19:30 2019-01-02 19:45 2019-01-02 20:00 2019-01-02 20:15 2019-01-02 20:30 2019-01-02 20:45 2019-01-02 21:00 2019-01-02 21:15 2019-01-02 21:30 2019-01-02 21:45 2019-01-02 22:00 2019-01-02 22:15 2019-01-02 22:30 2019-01-02 22:45 2019-01-02 23:00 2019-01-02 23:15 2019-01-02 23:30 2019-01-02 23:45 2019-01-02 00:00 2019-01-03 00:15 2019-01-03 00:30 2019-01-03 00:45 2019-01-03 01:00 2019-01-03 01:15 2019-01-03 01:30 2019-01-03 01:45 2019-01-03 02:00 2019-01-03 02:15 2019-01-03 02:30 2019-01-03 02:45 2019-01-03 03:00 2019-01-03 03:15 2019-01-03 03:30 2019-01-03 03:45 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.145 0.106 0.065 0.045 0.024 0.008 0.002 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 277.500 207.000 131.700 89.600 49.720 19.300 6.041 0.660 0.000 0.005 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.83 0.72 0.48 0.18 -0.39 -1.23 -2.02 -2.78 -3.66 -4.26 -4.42 -4.84 -5.26 -5.4 -5.33 -5.38 -5.84 -6.08 -6.23 -5.98 -6.03 -6.07 -5.77 -5.33 -4.93 -4.58 -4.65 -4.83 -5 -5.21 -5.19 -5.17 -4.99 -4.67 -4.34 -4.2 -4.08 -3.93 -3.79 -3.62 -3.35 -2.95 -2.08 -1.04 -0.49 -0.08 0.49 0.83 0.98 1.25 1.48 1.52 26 25 26 28 32 35 40 44 49 52 59 63 63 63 66 67 69 71 72 75 76 76 78 78 79 80 80 81 82 82 84 84 84 85 85 86 87 88 88 88 89 89 85 82 81 82 82 81 81 81 83 84 0.89 0.76 0.47 0.22 0.19 0.13 0.04 0.00 0.05 0.03 0.35 0.08 0.04 0.02 0.24 0.00 0.02 0.00 0.13 0.14 0.00 0.00 0.00 0.04 0.11 0.07 0.00 0.02 0.00 0.00 0.02 0.00 0.08 0.18 0.01 0.06 0.00 0.02 0.00 0.02 0.16 0.51 1.00 1.21 1.02 1.26 0.98 0.56 0.43 1.08 1.00 1.03 248.1 1024 200.5 1024 255.9 1024 249.3 1024 193.8 1024 182.1 1024 194.9 1024 194.9 1023 171.0 1023 148.0 1024 93.2 1024 123.7 1023 109.8 1023 297.4 1023 214.5 1023 12.9 1024 61.0 1024 97.2 1023 82.6 1023 115.0 1023 148.7 1023 148.8 1023 108.0 1023 240.0 1022 225.7 1022 100.8 1022 115.9 1021 120.5 1021 122.8 1021 131.0 1021 236.8 1020 233.0 1020 277.5 1020 280.2 1020 275.6 1020 266.9 1020 218.5 1020 157.1 1019 170.9 1019 256.1 1019 149.3 1018 118.4 1018 128.8 1018 131.0 1017 118.3 1017 133.8 1016 154.0 1016 195.2 1016 197.0 1016 175.3 1015 156.8 1015 151.2 1014 2019-01-03 04:00 2019-01-03 04:15 2019-01-03 04:30 2019-01-03 04:45 2019-01-03 05:00 2019-01-03 05:15 2019-01-03 05:30 2019-01-03 05:45 2019-01-03 06:00 2019-01-03 06:15 2019-01-03 06:30 2019-01-03 06:45 2019-01-03 07:00 2019-01-03 07:15 2019-01-03 07:30 2019-01-03 07:45 2019-01-03 08:00 2019-01-03 08:15 2019-01-03 08:30 2019-01-03 08:45 2019-01-03 09:00 2019-01-03 09:15 2019-01-03 09:30 2019-01-03 09:45 2019-01-03 10:00 2019-01-03 10:15 2019-01-03 10:30 2019-01-03 10:45 2019-01-03 11:00 2019-01-03 11:15 2019-01-03 11:30 2019-01-03 11:45 2019-01-03 12:00 2019-01-03 12:15 2019-01-03 12:30 2019-01-03 12:45 2019-01-03 13:00 2019-01-03 13:15 2019-01-03 13:30 2019-01-03 13:45 2019-01-03 14:00 2019-01-03 14:15 2019-01-03 14:30 2019-01-03 14:45 2019-01-03 15:00 2019-01-03 15:15 2019-01-03 15:30 2019-01-03 15:45 2019-01-03 16:0