Teaching: 45
Exercises: 0
Questions
Objectives
To make the best of the R language, you’ll need a strong understanding of the basic data types and data structures and how to operate on them.
Data structures are very important to understand because these are the objects you will manipulate on a day-to-day basis in R. Dealing with object conversions is one of the most common sources of frustration for beginners.
Everything in R is an object.
R has 6 basic data types. (In addition to the five listed below, there is also raw which will not be discussed in this workshop.)
Elements of these data types may be combined to form data structures, such as atomic vectors. When we call a vector atomic, we mean that the vector only holds data of a single data type. Below are examples of atomic character vectors, numeric vectors, integer vectors, etc.
"a"
, "swc"
2
, 15.5
2L
(the L
tells R to store this as an integer)TRUE
, FALSE
1+4i
(complex numbers with real and imaginary parts)R provides many functions to examine features of vectors and other objects, for example
class()
- what kind of object is it (high-level)?typeof()
- what is the object’s data type (low-level)?length()
- how long is it? What about two dimensional objects?attributes()
- does it have any metadata?# Example
"dataset"
x <-typeof(x)
[1] "character"
attributes(x)
NULL
1:10
y <- y
[1] 1 2 3 4 5 6 7 8 9 10
typeof(y)
[1] "integer"
length(y)
[1] 10
as.numeric(y)
z <- z
[1] 1 2 3 4 5 6 7 8 9 10
typeof(z)
[1] "double"
R has many data structures. These include
A vector is the most common and basic data structure in R and is pretty much the workhorse of R. Technically, vectors can be one of two types:
although the term “vector” most commonly refers to the atomic types not to lists.
A vector is a collection of elements that are most commonly of mode character
, logical
, integer
or numeric
.
You can create an empty vector with vector()
. (By default the mode is logical
. You can be more explicit as shown in the examples below.) It is more common to use direct constructors such as character()
, numeric()
, etc.
vector() # an empty 'logical' (the default) vector
logical(0)
vector("character", length = 5) # a vector of mode 'character' with 5 elements
[1] "" "" "" "" ""
character(5) # the same thing, but using the constructor directly
[1] "" "" "" "" ""
numeric(5) # a numeric vector with 5 elements
[1] 0 0 0 0 0
logical(5) # a logical vector with 5 elements
[1] FALSE FALSE FALSE FALSE FALSE
You can also create vectors by directly specifying their content. R will then guess the appropriate mode of storage for the vector. For instance:
c(1, 2, 3) x <-
will create a vector x
of mode numeric
. These are the most common kind, and are treated as double precision real numbers. If you wanted to explicitly create integers, you need to add an L
to each element (or coerce to the integer type using as.integer()
).
c(1L, 2L, 3L) x1 <-
Using TRUE
and FALSE
will create a vector of mode logical
:
c(TRUE, TRUE, FALSE, FALSE) y <-
While using quoted text will create a vector of mode character
:
c("Sarah", "Tracy", "Jon") z <-
The functions typeof()
, length()
, class()
and str()
provide useful information about your vectors and R objects in general.
typeof(z)
[1] "character"
length(z)
[1] 3
class(z)
[1] "character"
str(z)
chr [1:3] "Sarah" "Tracy" "Jon"
The function c()
(for combine) can also be used to add elements to a vector.
c(z, "Annette")
z <- z
[1] "Sarah" "Tracy" "Jon" "Annette"
c("Greg", z)
z <- z
[1] "Greg" "Sarah" "Tracy" "Jon" "Annette"
You can create vectors as a sequence of numbers.
1:10
series <-seq(10)
[1] 1 2 3 4 5 6 7 8 9 10
seq(from = 1, to = 10, by = 0.1)
[1] 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4
[16] 2.5 2.6 2.7 2.8 2.9 3.0 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9
[31] 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9 5.0 5.1 5.2 5.3 5.4
[46] 5.5 5.6 5.7 5.8 5.9 6.0 6.1 6.2 6.3 6.4 6.5 6.6 6.7 6.8 6.9
[61] 7.0 7.1 7.2 7.3 7.4 7.5 7.6 7.7 7.8 7.9 8.0 8.1 8.2 8.3 8.4
[76] 8.5 8.6 8.7 8.8 8.9 9.0 9.1 9.2 9.3 9.4 9.5 9.6 9.7 9.8 9.9
[91] 10.0
R supports missing data in vectors. They are represented as NA
(Not Available) and can be used for all the vector types covered in this lesson:
c(0.5, NA, 0.7)
x <- c(TRUE, FALSE, NA)
x <- c("a", NA, "c", "d", "e")
x <- c(1+5i, 2-3i, NA) x <-
The function is.na()
indicates the elements of the vectors that represent missing data, and the function anyNA()
returns TRUE
if the vector contains any missing values:
c("a", NA, "c", "d", NA)
x <- c("a", "b", "c", "d", "e")
y <-is.na(x)
[1] FALSE TRUE FALSE FALSE TRUE
is.na(y)
[1] FALSE FALSE FALSE FALSE FALSE
anyNA(x)
[1] TRUE
anyNA(y)
[1] FALSE
Inf
is infinity. You can have either positive or negative infinity.
1/0
[1] Inf
NaN
means Not a Number. It’s an undefined value.
0/0
[1] NaN
R will create a resulting vector with a mode that can most easily accommodate all the elements it contains. This conversion between modes of storage is called “coercion”. When R converts the mode of storage based on its content, it is referred to as “implicit coercion”. For instance, can you guess what the following do (without running them first)?
c(1.7, "a")
xx <- c(TRUE, 2)
xx <- c("a", TRUE) xx <-
You can also control how vectors are coerced explicitly using the as.<class_name>()
functions:
as.numeric("1")
[1] 1
as.character(1:2)
[1] "1" "2"
Do you see a property that’s common to all these vectors above?
All vectors are one-dimensional and each element is of the same type.
Objects can have attributes. Attributes are part of the object. These include:
You can also glean other attribute-like information such as length (works on vectors and lists) or number of characters (for character strings).
length(1:10)
[1] 10
nchar("Software Carpentry")
[1] 18
In R matrices are an extension of the numeric or character vectors. They are not a separate type of object but simply an atomic vector with dimensions; the number of rows and columns. As with atomic vectors, the elements of a matrix must be of the same data type.
matrix(nrow = 2, ncol = 2)
m <- m
[,1] [,2]
[1,] NA NA
[2,] NA NA
dim(m)
[1] 2 2
You can check that matrices are vectors with a class attribute of matrix
by using class()
and typeof()
.
matrix(c(1:3))
m <-class(m)
[1] "matrix" "array"
typeof(m)
[1] "integer"
While class()
shows that m is a matrix, typeof()
shows that fundamentally the matrix is an integer vector.
Consider the following matrix:
matrix(
FOURS <-c(4, 4, 4, 4),
nrow = 2,
ncol = 2)
Given that typeof(FOURS[1])
returns "double"
, what would you expect typeof(FOURS)
to return? How do you know this is the case even without running this code?
Hint Can matrices be composed of elements of different data types?
We know that typeof(FOURS)
will also return "double"
since matrices are made of elements of the same data type. Note that you could do something like as.character(FOURS)
if you needed the elements of FOURS
as characters.
Matrices in R are filled column-wise.
matrix(1:6, nrow = 2, ncol = 3) m <-
Other ways to construct a matrix
1:10
m <-dim(m) <- c(2, 5)
This takes a vector and transforms it into a matrix with 2 rows and 5 columns.
Another way is to bind columns or rows using rbind()
and cbind()
(“row bind” and “column bind”, respectively).
1:3
x <- 10:12
y <-cbind(x, y)
x y
[1,] 1 10
[2,] 2 11
[3,] 3 12
rbind(x, y)
[,1] [,2] [,3]
x 1 2 3
y 10 11 12
You can also use the byrow
argument to specify how the matrix is filled. From R’s own documentation:
matrix(c(1, 2, 3, 11, 12, 13),
mdat <-nrow = 2,
ncol = 3,
byrow = TRUE)
mdat
[,1] [,2] [,3]
[1,] 1 2 3
[2,] 11 12 13
Elements of a matrix can be referenced by specifying the index along each dimension (e.g. “row” and “column”) in single square brackets.
2, 3] mdat[
[1] 13
In R lists act as containers. Unlike atomic vectors, the contents of a list are not restricted to a single mode and can encompass any mixture of data types. Lists are sometimes called generic vectors, because the elements of a list can by of any type of R object, even lists containing further lists. This property makes them fundamentally different from atomic vectors.
A list is a special type of vector. Each element can be a different type.
Create lists using list()
or coerce other objects using as.list()
. An empty list of the required length can be created using vector()
list(1, "a", TRUE, 1+4i)
x <- x
[[1]]
[1] 1
[[2]]
[1] "a"
[[3]]
[1] TRUE
[[4]]
[1] 1+4i
vector("list", length = 5) # empty list
x <-length(x)
[1] 5
The content of elements of a list can be retrieved by using double square brackets.
1]] x[[
NULL
Vectors can be coerced to lists as follows:
1:10
x <- as.list(x)
x <-length(x)
[1] 10
x[1]
?x[[1]]
?{r examine-lists-1} class(x[1])
{r examine-lists-2} class(x[[1]])
Elements of a list can be named (i.e. lists can have the names
attribute)
list(a = "Karthik Ram", b = 1:10, data = head(mtcars))
xlist <- xlist
$a
[1] "Karthik Ram"
$b
[1] 1 2 3 4 5 6 7 8 9 10
$data
mpg cyl disp hp drat wt qsec vs am gear carb
Mazda RX4 21.0 6 160 110 3.90 2.620 16.46 0 1 4 4
Mazda RX4 Wag 21.0 6 160 110 3.90 2.875 17.02 0 1 4 4
Datsun 710 22.8 4 108 93 3.85 2.320 18.61 1 1 4 1
Hornet 4 Drive 21.4 6 258 110 3.08 3.215 19.44 1 0 3 1
Hornet Sportabout 18.7 8 360 175 3.15 3.440 17.02 0 0 3 2
Valiant 18.1 6 225 105 2.76 3.460 20.22 1 0 3 1
names(xlist)
[1] "a" "b" "data"
{r examine-named-lists-1} length(xlist)
{r examine-named-lists-2} str(xlist)
Lists can be extremely useful inside functions. Because the functions in R are able to return only a single object, you can “staple” together lots of different kinds of results into a single object that a function can return.
A list does not print to the console like a vector. Instead, each element of the list starts on a new line.
Elements are indexed by double brackets. Single brackets will still return a(nother) list. If the elements of a list are named, they can be referenced by the $
notation (i.e. xlist$data
).
A data frame is a very important data type in R. It’s pretty much the de facto data structure for most tabular data and what we use for statistics.
A data frame is a special type of list where every element of the list has same length (i.e. data frame is a “rectangular” list).
Data frames can have additional attributes such as rownames()
, which can be useful for annotating data, like subject_id
or sample_id
. But most of the time they are not used.
Some additional information on data frames:
read.csv()
and read.table()
, i.e. when importing the data into R.data.frame()
function.nrow(dat)
and ncol(dat)
, respectively.To create data frames by hand:
data.frame(id = letters[1:10], x = 1:10, y = 11:20)
dat <- dat
id x y
1 a 1 11
2 b 2 12
3 c 3 13
4 d 4 14
5 e 5 15
6 f 6 16
7 g 7 17
8 h 8 18
9 i 9 19
10 j 10 20
See that it is actually a special list:
is.list(dat)
[1] TRUE
class(dat)
[1] "data.frame"
Because data frames are rectangular, elements of data frame can be referenced by specifying the row and the column index in single square brackets (similar to matrix).
1, 3] dat[
[1] 11
As data frames are also lists, it is possible to refer to columns (which are elements of such list) using the list notation, i.e. either double square brackets or a $
.
"y"]] dat[[
[1] 11 12 13 14 15 16 17 18 19 20
$y dat
[1] 11 12 13 14 15 16 17 18 19 20
The following table summarizes the one-dimensional and two-dimensional data structures in R in relation to diversity of data types they can contain.
Dimensions | Homogenous | Heterogeneous |
---|---|---|
1-D | atomic vector | list |
2-D | matrix | data frame |
Knowing that data frames are lists, can columns be of different type?
What type of structure do you expect to see when you explore the structure of the iris
data frame? Hint: Use str()
.
The Sepal.Length, Sepal.Width, Petal.Length and Petal.Width columns are all numeric types, while Species is a Factor. Lists can have elements of different types. Since a Data Frame is just a special type of list, it can have columns of differing type (although, remember that type must be consistent within each column!).
str(iris)
'data.frame': 150 obs. of 5 variables:
$ Sepal.Length: num 5.1 4.9 4.7 4.6 5 5.4 4.6 5 4.4 4.9 ...
$ Sepal.Width : num 3.5 3 3.2 3.1 3.6 3.9 3.4 3.4 2.9 3.1 ...
$ Petal.Length: num 1.4 1.4 1.3 1.5 1.4 1.7 1.4 1.5 1.4 1.5 ...
$ Petal.Width : num 0.2 0.2 0.2 0.2 0.2 0.4 0.3 0.2 0.2 0.1 ...
$ Species : Factor w/ 3 levels "setosa","versicolor",..: 1 1 1 1 1 1 1 1 1 1 ...
{% include links.md %}