# A Tutorial on R

## Introduction

Welcome to the R tutorial. The tutorial is meant for people who already know something about programming in general.

### Aim of the tutorial

Starting from some fundamentals, learn R on a generic level. R has got so many different packages for many different approaches that it would be very difficult to cover all of them comprehensively. So I want to give you a nice start into learning R for whatever purpose you may need it for.

## Slicing different parts of objects

Run this code on your own computer to see different methods to slice objects.

dat <- data.frame( # Define a data frame with five rows and three columns.
A = 1:5,
B = c("c", "b", "d", "e", "c"),
Result = c(6, 45, 2, 4.5, 2)
)

dat		# The data frame dat.
dat\$A  		# The vector that forms column A.
dat[["A"]]	# The same as previous.
dat[[1]]	# The vector that forms the first column (the same as previous).
dat["A"] 	# The column with name A (this is a data frame with one column).
dat[1] 		# The first column  (with name A; the same as previous).
dat[c(1,3)]	# Data frame with the first and third columns.
dat[c("A", "Result")] 	The same as previous.
dat[2:4, 1:2] 			# Data frame with rows 2 to 4 and columns 1 to 2.
dat[2:4, "Result"] 		# The vector that is formed from rows 2 to 4 of column Result.
dat\$Result[2:4] 		# The same as previous.
dat[2:4, ]["Result"] 	# Data frame that first takes rows 2 to 4 and then column Result.
dat\$B		# Factor that forms column B. Note that factors are also vectors.
levels(dat\$B)	# Levels (i.e. possible values of the factor) of the previous factor.
as.numeric(dat\$B)	# The numeric position values of the previous factor.
levels(dat\$B)[dat\$B]	# The previous factor converted into a character vector.
c("b", "c", "d", "e")[c(2, 1, 3, 4, 2)]	# The same as previous

library(OpasnetUtils) # This package is needed to operate with ovariables.

odat <- Ovariable("odat", data = dat) # Ovariable that has dat as data.
odat@data 	# Data slot of odat.
odat 		# All slots of odat.
odat <- EvalOutput(odat) # Evaluate odat (i.e., calculate the output)
odat@output	# Output slot of odat.
result(odat)	# The result column of the output slot of odat.
summary(odat) # Summary of odat. If odat is probabilistic, summary includes mean and other statistics.
odat@marginal # Which columns in output are marginals?
colnames(odat@output) 	# Names of columns of the output of odat
colnames(odat@output)[odat@marginal]	# Names of marginal columns.
odat@output[odat@output\$B == "c" , ]	# Data frame from odat output of rows where column B has value c.
result(odat)[odat@output\$B == "c"]	# Vector from odat result of rows rows where column B has value c.

## Chapter 1 - the Basics of the R Language

• Basic syntax
• The creatures of R: vectors, lists, factors, data.frames, arrays
• Data types

R is an object oriented programming language.

The R console is a command line interpreter, typing out a command e.g. 2+2 returns 4. Mathematical operators (+,-,*,/,^) are used in their mathematical sense and the order of execution is the mathematical one. The operators can also be used in function form, e.g. '+'(1,1). More operators here. Brackets () can be used to control which expressions are evaluated first, e.g. (1+1)/2.

Variables and functions are saved in memory using either an arrow (<- or ->) or =, e.g. var1 <- 1. It is preferable to use the arrow operator.

Functions are used as function(parameter1, parameter2, parameter6 = value), e.g. mean(rnorm(10, sd = 2)). Checking syntax and other details for a function is easy by using ?function, e.g. ?mean. Lists of functions can be found using library(help="package"), where package is the container/library of functions of interest. E.g. library(help="base") or library(help="stats"). A semicolon (;) can be used to separate statements on the same line. E.g. a <- 2;a.

The basic structure most objects in R is a vector. Vectors can be atomic (contain only values), or recursive (vector of objects, basically). A vector in R is basically an ordered set of values or objects. Selecting (subscripting) the n:th element is done by vector[n] or vector[[n]]. The double bracket is used when subscripting from recursive vectors (lists and its subtypes) to extract the stored objects themselves rather than a recursive vector of the selected elements.

There are 3 basic data types: text (class = character), numbers (class = numeric) and logical. Data types in R consist of these basic data types and their more elaborate derivatives. E.g. a factor is a character vector (a vector that consists of textual elements) stored as a numeric vector, where each number represents a unique element of the character vector. The unique elements are stored as a levels attribute of the factor object. Integers are a special case of the numeric class, they are handled as numeric except in storage; normally all numbers are numeric, a special case is the 'a:b' operator which produces a vector of integers from a to b. An atomic vector can only contain values of a single data type.

 a <- 1:6 a b <- c("Q", "W", "E", "R", "T", "Y") b d <- data.frame(Row = a, Letter = b) d b <- c(b, c("W", "W", "T")) b class(b) b <- as.factor(b) b class(b) b <- as.numeric(b) b class(b) d[ , "Letter"] d\$Letter d["Letter"] d[d\$Row == 2, ] d[d\$Row %in% 2:4, ] print(d) library(xtable) print(xtable(d), type = 'html') library(OpasnetUtils) oprint(d)

Attributes are used to simulate more complex data structures. An atomic vector can be given a dimensions attribute dim (which is a numeric vector containing the lengths of the dimensions), to turn it into an array (a matrix is an array with length(dim) = 2). Dimensions can be given names in the dimnames attribute which is a list of named character vectors. Because arrays and matrices are atomic vectors by nature they can only contain values of one data type.

 a <- 1:27;dim(a) <- c(3,3,3);dimnames(a) <- list(dim1=1:3, dim2=1:3, dim3=1:3);a b <- array(1:27, dim = c(3,3,3), dimnames = list(dim1=1:3, dim2=1:3, dim3=1:3));b

Recursive vectors; lists and data.frames (which are lists with elements of equal length atomic vectors, character vectors are by default converted to factors) can have values of different data types since they consist of different objects. The data.frame is perhaps the most common object type in R. It resembles the basic rectangular table format.

 a <- 1:9; b <- factor(rep(c("a","b","c"),3)); d <- factor(rep(c("q","w","e"), each = 3)) df <- data.frame(a,b,d);df df2 <- data.frame(a = 18:10, b = factor(rep(c("a","b","c"),3)), d = factor(rep(c("q","w","e"), each = 3)));df2 l <- list(df, df2, k = 1:2);l l[[1]] l[[3]] l[3] l[[1]][5,1] l[[1]][5,] l[[1]][,1]

In general arrays are produced and used when data is summarized (summed or averaged over some marginals), for computational purposes the data.frame is superior.

Classes in R refer to either the data type of a simple atomic vector, or the object type of a more complex object. Functions in R may have different methods for handling inputs of different classes and this may sometimes confuse newcomers; e.g. some functions take a factor input as only a numeric vector instead of a character vector. Many functions try to coerce their input into the format they can operate with by using functions like as.character and as.numeric.

 class(1);class(1:4);class("a");class(TRUE);class(as.factor("a","a","b")) 1:4

Many R functions are vectorized, meaning that a function can take one or more vectors as input to produce a vector as output. I.e. 1:5 + 10:6 produces 11 11 11 11 11. If the vectors are of different length, the shorter one is usually recycled to match the length of the longest vector. E.g. 1:2 + rep(4,5) produces 5 6 5 6 5 with a warning calling attention to the arguments different lengths. Vectorized operations are incredibly fast and should be used whenever possible in place of for, while or other loops.

## Chapter 2 - Getting your data into R

 library(OpasnetUtils) e <- opbase.data("Op_en5103") oprint(e) e <- tidy(e) oprint(e)

To upload data, you can use Table2Base tables on a wiki page:

A Tutorial on R: Difference between revisions(kton CO2e)
ObsRowSectorYearAmountPer.person
11Consumers' use of electricity 2010126.961.31
22Electric heating 201032.650.34
33District heating 2010321.183.31
44Separate heating 201041.100.42
55Traffic 2010164.471.70
66Agriculture 201024.360.25
77Waste management 201020.550.21
88Total 2010731.277.54

Or you can go to upload functionality to upload an Excel or CSV file:

• Importing and exporting data from/to files
• Working with data: subscription, merge, apply, reshape, conversion between data.frame and array
• OpasnetUtils?

Example data: File:ArkS280.csv. Be careful when converting your Excel sheets to .csv, Microsoft Office is an idiotic piece of software and it sometimes writes empty cells in places where you have once edited something. To fix this open the exported .csv in the Open Office equivalent of Excel and overwrite the previous file.

Example:

test <- read.table("M:/R koulutus/arkS280.csv", sep = ";", dec = ",", header = TRUE)
test

The output is in the data.frame format. We can select row(s) and column(s) by using subscription. The data.frame is the most flexible format when it comes to data exploration and subscription. Since data.frames are essentially lists, we can use list1[[x]] to select a column x, where x can be a numerical vector (it can be longer than one, experiment with nested lists if you're interested). If the list is named (data.frames always are) we can select its elements by using list1[[x]], where x is a character vector of length one, or list1\$col1, where col1 is simply the name of the list element (column in a data.frame). After selecting an object from a list we can subscript from it again, e.g. list1[[1]][1] returns the first element of the first object stored in list1, list1\$col1[1] is similar; list1[[1]][[1]]... can be used for nested lists. A data.frame is special in that it is also subscribeable as a two dimensional array: df1[x,y] returns the x:th value of the row y, both can be vectors of any length and of any basic data type (either numeric, character or logical). Either of the x and y can be left blank so that a vector is returned; if a vertical slice is extracted, the result is an atomic vector if only one column was selected (same as selecting the object from the list); if a horizontal slice is extracted, the result is a data.frame. Arrays with more dimensions can be subscribed from in similar fashion e.g. arr1[x,y,z,...], where x, y, z and so on can be vectors of any basic data type.

To better utilize subscription you should learn about logical operators. <, >, <=, >=, ==, and !=, are pretty self-explanatory. More advanced ones include the and (&) and or (|). %in% is also pretty useful. To obtain a inversion use (statement)==FALSE or !(statement). These operators are vectorized, however the expression on the right side must be of length one (comparing all values of a vector on the left side to a single value/expression is allowed, but element by element comparison is not allowed). grep can be used to find regular expressions. All logical operators return a logical vector.

Example:

test <- read.table("M:/R koulutus/arkS280.csv", sep = ";", dec = ",", header = TRUE)
test[test\$Suklaa0>=5,]
test[test\$Ryhmä=="ip"]

R has got very powerful data manipulation facilities. Data in data.frames usually consist of a few factors and numeric vectors. The factors are usually indices to the data, e.g. in a population data there could be indices for Age, Year and Place. Unique cells in the data would be identified by a unique combination of these factors. This format can be used very similarly to an array. We could sum or take a mean over the levels of specified factors using tapply (data.frame variant of apply). xtabs creates a contingency table from a data.frame with cross-classifying factors, this is similar to an array with some extras. table is similar to xtabs but simpler and only does counts. reshape is a function that transforms data.frames with a single numeric vector into a data.frame with multiple numeric vectors specified by one or more indexing factors and vice versa. merge is a function that merges two data.frames by finding matching (indexing) vectors, any extra vectors in either data are carried over to the resulting data.frame. To simplify variable selection, components of data.frames can be attached to the general namespace, i.e. attach(data) would enable calling the component vectors of data directly: e.g. vec1 instead of data[,"vec1"].

There are many packages for database connections. One for use with the Opasnet Base is the OpasnetUtils, which uses the RODBC package for the actual connection.

A comprehensive guide to importing and exporting data in R can be found on http://cran.r-project.org/doc/manuals/R-data.html.

Useful data packages

• fmi: direct access to Finnish Meteorological Institute database: all monitoring station data, modelled weather. Note! You need a personal key from FMI to use it.

Opening sources from Opasnet and own computer

 library(OpasnetUtils) # Read a csv file from your own computer d <- read.csv("C:/Users/Neukkari/AppData/Local/Temp/ArkS280.csv", sep = ";", dec = ",") # Write a data.frame to a csv file on your own computer write.table(rob, "C:/Users/Neukkari/Documents/koe.csv", sep = ",", dec = ".") # Download a csv file from Opasnet wiki d <- opasnet.csv("/f/fc/ArkS280.csv", wiki = "opasnet_en") # Download a table from Opasnet Base e <- opbase.data("Op_en5103") rob <- opbase.data("Op_fi5339", subset = "Robottiautot ja matkojenyhdistely") sarc <- opbase.data("op_en2721", subset = "Koe") # Useful package to access data from Finnish Meteorological Institute FMI library(fmi)

## Chapter 3: Graphs

Examples of making graphs with ggplot2 package.

 library(OpasnetUtils) library(ggplot2) # graphical package library(reshape2) # package for melt function # Create a data.frame with 3 columns of data a <- data.frame( A = 1:100, B = rnorm(100), C = rnorm(100, 2, 1), D = rnorm(100, -1, 3) ) # Melt the data columns into one and create an explanatory column for source. a <- melt(a, id.vars = "A", variable.name = "Source") # Draw a line graph. You can also use geom_point() for point graphs and a lot of other alternatives. ggplot(a, aes(x = A, y = value, colour = Source)) + geom_line()+ theme_gray(base_size = 24)+ labs(x = "Number", y = "Result", title = "Great graph" )

## Chapter 4

Getting help

 a <- 1:10 print(a) print ?print

Plotting fancy plots

 n <- 1000 population <- data.frame( Sex = rep(c("Male", "Female"), each = n), Height = c( rnorm(n, 178, 18), rnorm(n, 168, 15) ) ) head(population) plot(population\$Sex, population\$Height) tapply(population\$Height, population\$Sex, mean) tapply(population\$Height, population\$Sex, sd)

Fancy map plots:

 What is the value for Finland?: library(xtable) data <- data.frame( Country = c("AT", "BE", "BG", "CH", "CY", "CZ", "DE", "DK", "EE", "ES", "FI", "FR", "GR", "HU", "IE", "IT", "LT", "LU", "LV", "MT", "NL", "NO", "PL", "PT", "RO", "SE", "SI", "SK", "UK"), Result = c(42, 78, 33, 57, 82, 66, 40, 65, 93, 50, 37, 74, 93, 26, 27, 15, 83, 36, 34, 89, 45, 96, 23, 39, 40, 22, 58, 20, 10) ) data[data\$Country == "FI", "Result"] <- fi cat("Results by country (using country codes).\n") print(xtable(data), type = 'html') eumap(data\$Result)

 library(rgdal) library(maptools) library(RColorBrewer) library(classInt) library(OpasnetUtils) shp<-readOGR('PG:host=localhost user=postgres dbname=spatial_db','watson_wkt') ###################### eumap plots values of EU29 countries on map. Parameters: ######### data: must be a vector with length 29 (country values in this order: ## "AT", "BE", "BG", "CH", "CY", "CZ", "DE", "DK", "EE", "ES", "FI", "FR", "GR", "HU", ## "IE", "IT", "LT", "LU", "LV", "MT", "NL", "NO", "PL", "PT", "RO", "SE", "SI", "SK", "UK" ######### nclr: number of colours to be used. Default: 24 eumap <- function(data = 1:29, nclr = 24) { shp@data\$value_inhalation <- data plotvar <- shp@data\$value_inhalation nclr <- 24 plotclr <- brewer.pal(nclr,"BuPu") class <- classIntervals(plotvar,nclr,style="quantile") colcode <- findColours(class,plotclr) epsg4326String <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs") proj4string(shp)<- ("+init=epsg:3035") shp2 <- spTransform(shp,epsg4326String) out<-sapply( slot(shp2,"polygons"), function(x){ kmlPolygon( x, name=as(shp2,"data.frame")[slot(x,"ID"),"country_code"], col=colcode[[((as.numeric(slot(x,"ID"))+1))]], lwd=1, border='black', description=paste("Value:",as(shp2,"data.frame")[slot(x,"ID"),"value_inhalation"]) ) } ) data<-paste( paste( kmlPolygon( kmlname="This will be layer name", kmldescription="More info about layer here" )\$header, collapse="\n" ), paste(unlist(out["style",]), collapse="\n"), paste(unlist(out["content",]), collapse="\n"), paste(kmlPolygon()\$footer, collapse="\n"), sep='' ) google.show_kml_data_on_maps(data) }

• The magic
• Packages
• Fancy plots

R is completely modular, i.e. all functions in R come in packages (libraries). The basic installation of R comes with some 10 packages, which define most of R's basic functionalities. Installing new packages is easy through a top bar menu in the R GUI. Alternatively if you know what you're doing you can use the install.package function directly. Only the some basic packages are loaded into memory during R startup, though those settings can be altered. Specific packages can be loaded using library(packagename).

A useful package for plotting gorgeous graphs is ggplot2. Information on it can be found here. Help about graphics:

• Modeling
• Probability distributions
• VOI analysis on R ==
• BRUGS (Open BUGS on R)
• RJAGS (Just another Gibbs sampler)
• ff (on disk objects)
• ...

## Cool tricks with R

### Multiplying data and adding depth to a time dimension (e.g. minutes) using string selections and regular expressions (completely vectorized)

Time format "1.1.2011 00:00:00", the length of this string varies.

 # First we'll create some random data data <- data.frame(Time=paste("1.1.2011 ", c(paste(0, 0:9, sep = ""), as.character(10:23)), ":00:00", sep = ""), Conc = rnorm(24,10,2)) data # Then multiply the number of rows by 60 data <- data[rep(1:nrow(data), each = 60),] #select from data 60 of each row, overwrite data rownames(data) <- 1:nrow(data) #make rownames sensible data # Change minutes to "00" to "59", currently repeated "00" temp <- as.character(data[,"Time"]) svec <- regexpr(":00:", temp) # gives the positions of ":00:" in the strings, output is a vector so the strings may be of different length e.g. "11.1.2011 00:00:00" vs "1.1.2011 00:00:00" substring(temp,svec+1,svec+2) <- c(paste(0, 0:9, sep = ""), as.character(10:59)) # left side selects parts of the strings (vectorized) based on the svec -- # and the right side substitutes a vector of strings to replace the selection data[,"Time"] <- temp data

### Applying a function in a data.frame with multiple rows with values

Example data from Chapter 2.

test <- read.table("M:/R koulutus/arkS280.csv", sep = ";", dec = ",", header = TRUE)

# The basic case where we want to take mean of a set of observations, indexed by something
tapply(test[,"Paino2"], test[,c("Annos", "Ryhmä")], mean) # mean of Paino2 by Annos and Ryhmä

# Tricky version where we take means of several sets of observations indexed by some indices
testf <- function(X, INDEX2, FUN2) tapply(X, INDEX2, FUN2) # define a custom function that does the above to a set of observations X
lapply(test[,7:15], testf, INDEX2 = test[,c("Annos", "Ryhmä")], FUN2 = mean) # apply our test function to all selected columns

# Could also be done with a loop
output <- list() # define a variable so it can be used inside the loop
for (i in 7:15) { # loop for i so that it takes the values given in a vector, 7:15 in this case
output[[length(output)+1]] <- tapply(test[,i], test[,c("Annos", "Ryhmä")], mean) # apply the above to column defined by i and put into the output list
}
names(output) <- colnames(test)[7:15] # give names to the list objects from the column names of the original table
output # this is identical to the output from the tricky version. Both methods actually use loops so the main difference here is the syntax.

## HELP

• How to merge two veeery big datasets together by R? Pauliina is going to workout it in here:

Temperature and population in Europe. Codes are very heavy and free to rewrite.

• If I want to attach new column with rownumbers, how can I do it (What is the code for that)?
• Would it not suffice to use the rownames attribute? e.g. rownames(test) <- 1:nrow(test)
• Yeah, I'll managed to do that with this code:
rownames(data) = c(1:nrow(data))#names rows by numbers of the row
rownames(data) = rownames(data, prefix = "ID")#name "ID" for the new column