diff --git a/module3/exo3/exercice_en.Rmd b/module3/exo3/exercice_en.Rmd index 13b258ddd0da29bc3bf08c64b6a1db742f6d5409..a85afaf78e17c3d04946c8adbfd6dd0f512b1a63 100644 --- a/module3/exo3/exercice_en.Rmd +++ b/module3/exo3/exercice_en.Rmd @@ -1,33 +1,51 @@ --- -title: "Your title" -author: "Your name" -date: "Today's date" +title: "CO2 Keeling" +author: "David Elser" +date: "05-03-2021" output: html_document --- +## Data preprocessing fetch the Data file and check if file already exists +The data is available from the Web site of the Scripps Institute "https://scrippsco2.ucsd.edu/data/atmospheric_co2/primary_mlo_co2_record.html" -```{r setup, include=FALSE} -knitr::opts_chunk$set(echo = TRUE) -``` +Soure script available at : https://github.com/DemirhanOzelTrojan/Keeling_Predict/blob/master/main.R +install.packages('forecast', dependencies = TRUE) +install.packages("imputeTS") + +```{r} -## Some explanations +install.packages("imputeTS") +library(forecast) +library(imputeTS) +``` -This is an R Markdown document that you can easily export to HTML, PDF, and MS Word formats. For more information on R Markdown, see . -When you click on the button **Knit**, the document will be compiled in order to re-execute the R code and to include the results into the final document. As we have shown in the video, R code is inserted as follows: +```{r} +data_url = "https://scrippsco2.ucsd.edu/assets/data/atmospheric/stations/in_situ_co2/monthly/monthly_in_situ_co2_mlo.csv" +data_file = "CO2.csv" -```{r cars} -summary(cars) +if (!file.exists(data_file)) { + download.file(data_url, data_file, method="auto") +} ``` -It is also straightforward to include figures. For example: +Skip the frist 5 lines of the file -```{r pressure, echo=FALSE} -plot(pressure) +```{r} +data = read.csv(data_file, header=TRUE, skip=5) +``` + +Let's have a look at what we got: +```{r} +head(data) +tail(data) ``` -Note the parameter `echo = FALSE` that indicates that the code will not appear in the final version of the document. We recommend not to use this parameter in the context of this MOOC, because we want your data analyses to be perfectly transparent and reproducible. -Since the results are not stored in Rmd files, you should generate an HTML or PDF version of your exercises and commit them. Otherwise reading and checking your analysis will be difficult for anyone else but you. +```{r} +data$X.ppm.[data$X.ppm. == -99.99] <- NA +data$X.ppm. <- na.interpolation(data$X.ppm.) +head(data) +tail(data) +``` -Now it's your turn! You can delete all this information and replace it by your computational document.