библиотека ( здесь ) # получить лист данных из URL- адреса Scripps <- c ( "http://scrippsco2.ucsd.edu/assets/data/atmospheric/stations/flask_co2/monthly/monthly_flask_co2_nzd.csv" ) файл <- c ( "co2_nzd .csv " ) download.file ( url , file ) rm ( url , file ) url <- c ( "http://scrippsco2.ucsd.edu/assets/data/atmospheric/stations/in_situ_co2/monthly/monthly_in_situ_co2_mlo.csv" ) файл <- c ( "co2_mlo.csv" ) download.file ( url , файл )# читать данные Baring Head, пропуская первые 57 строк bhd <- read.csv ( "co2_nzd.csv" , skip = 57 , header = FALSE , sep = "," , dec = "." , stringsAsFactors = FALSE , strip .white = TRUE , na.strings = -99.99 ) # подмножество только даты и измерений co2 (10-й столбец, который включает заполненные данные, поэтому нет NAs bhd <- bhd [ 7 : 504 , c ( -1 , -2 , -3 , -5 , -6 , -7 , -8 , -9 )] names ( bhd ) <- c ( "Дата" , "CO2" ) str ( bhd ) 'data.frame' : 498 набл. Из 2 переменных : $ Дата : число 1978 1978 1978 1978 1978 ... $ CO2 : число 332 332 332 332 334 ...# читать данные Мауна-Лоа, пропуская первые 291 строку, поэтому первая строка имеет размер 1977 07 mlo <- read.csv ( "co2_mlo.csv" , skip = 291 , header = FALSE , sep = "," , dec = "." , na.strings = -99.99 ) # подмножество только даты и co2 mlo <- mlo [ 1 : 498 , c ( -1 , -2 , -3 , -6 , -7 , -8 , -9 , -10 )] names ( mlo ) <- c ( "Дата" , "CO2" ) str ( mlo ) 'data.frame' : 498 набл. из 2 переменных : $ Дата : число 1978 1978 1978 1978 1978 ... $ CO2 : число 335 333 332 331 332 ...# создать график в формате SVG svg ( filename = "Baringhead_co2_720-540.svg" , width = 8 , height = 6 , pointsize = 14 , onefile = FALSE , family = "sans" , bg = "white" ) par ( mar = c ( 3.1 , 3.1 , 1 , 1 ) +0.1 ) plot ( mlo «[[Date]]» , mlo «[[CO2]]]» , ylim = c ( 325 , 415 ), tck = 0,01 , оси = F , ann = F , las = 1 , pch = 20 , cex = 0,75 , type = "o" , col = "darkgray" , lwd = 1 ) ось ( сторона = 1 , tck = 0,01 , at = NULL , label = NULL , tick = T , lwd = 0 , lwd.tick = 1 ) ось ( сторона = 2 , tck = 0,01 , at = NULL , label = NULL , tick = T , lwd = 0 , lwd.tick = 1 , las = 1 ) grid () box () lines ( bhd "[[Date]]" , bhd "[[CO2]]" , col = 2 , lwd = 2 , lty = 1 ) mtext ( side = 2 , cex = 1 , line = -1,5 , выражение ( paste ( «Частей углекислого газа на миллион» ))) mtext ( side = 3 , cex = 1.5 , line = -2 , expression ( paste ( «Atmospheric C» , O [ 2 ], «Baring Head 1977 в 2018" году ))) Mtext ( сторона = 1 , линия = -2,8 , CEX = 1 , выражение ( вставка ( "Data: Скрипс C" , O [ 2 ], "Программа " ))) mtext ( side = 1 , cex = 0.7 , line = -1.3 , " http://scrippsco2.ucsd.edu/assets/data/atmospheric/stations/in_situ_co2/monthly/monthly_in_situ_co2_mlo.csv\nhttp:// scrippsco2.ucsd.edu/assets/data/atmospheric/stations/flask_co2/monthly/monthly_flask_co2_nzd.csv " ) legend ( 1980 , 400 , bty = " n " , c ( " Mauna Loa Hawaii " , " Baring Head New Zealand " ) , lwd = c ( 1 , 2 ), pch = c ( 20 , NA ), lty = 1 , col = c ( "темно-серый" , 2 )) dev.off ()