About

R is a language and environment for statistical computing and graphics. R provides a wide variety of statistical (linear and nonlinear modelling, classical statistical tests, time-series analysis, classification, clustering, .) and graphical techniques, and is highly extensible.

This notebook is a tutorial on how to use R.

Setup

Remember to always set your working directory to the source file location. Go to ‘Session’, scroll down to ‘Set Working Directory’, and click ‘To Source File Location’. Read carefully the below and follow the instructions to complete the tasks and answer any questions. Submit your work to RPubs as detailed in previous notes.

Note

For your assignment you may be using different data sets than what is included here. Always read carefully the instructions on Sakai. For clarity, tasks/questions to be completed/answered are highlighted in red color and numbered according to their particular placement in the task section. Quite often you will need to add your own code chunk.

Execute all code chunks, preview, publish, and submit link on Sakai.

Basics Operations

First we will begin with a few basic operations.

Variable assignment

We assign values to variables using the assignment operator ‘=’. Another form of assignment, more general, is the ‘<-’ operator. A variable allows you to store values or an object (e.g. a function).

x = 128
y = 16
z <- 5
vars = c(2,4,8,16,32) # Creates a vector list using the generic combine function 'c'

x # display value of variable x

[1] 128

z # displays value of variable z

[1] 5

vars[1] #This calls the first value in the vector vars

[1] 2

vars[2] #This calls the second value in the vector vars

[1] 4

vars[1:3] #This calls the first through third values in the vector vars

[1] 2 4 8

vars #This calls the vector list

[1]  2  4  8 16 32

Common Arithmetic Operations

Below shows some simple arithmetic operations.

12*6

[1] 72

128/16

[1] 8

9^2

[1] 81

Basic Data Types

R works with numerous data types. Some of the most basic types are: numeric,integers, logical (Boolean-TRUE/FALSE) and characters (string-"TEXT").

#Type: Character                   
#Example:"TRUE",'23.4'

v = "TRUE"                       
class(v)

[1] "character"

#Type: Numeric                
#Example: 12.3,5

v = 23.5                  
class(v)

[1] "numeric"

              
#Type: Logical    
#Example: TRUE,FALSE

v = TRUE
class(v)

[1] "logical"

#Type: Factor (nominal, categorical)
#Example: m f m f m

v = as.factor(c("m", "f", "m"))
class(v)

[1] "factor"

Functions

R Functions are invoked by its name, followed by the parenthesis, and zero or more arguments.

# The following applies the function 'c' (seen earlier) to combine three numeric values into a vector 
c(1,2,3)

[1] 1 2 3

# Example of function mean() to calcule the mean of three values
mean(c(5,6,7))

[1] 6

# Square root of a number
sqrt(99)

[1] 9.949874

Importing Data and Variable Assignment

# Here we are reading a file of type csv (comma seperated values) typical of many Excel files
il_income = read.csv(file = "data/il_income.csv")
top_il_income = read.csv(file = "data/top_il_income.csv")

Arithmetic Operations with Data

We can extract values from the dataset to perform calculations by referencing the proper elements of a list

DuPage = top_il_income$per_capita_income[1]
Lake = top_il_income$per_capita_income[2]
DuPage-Lake

[1] 472

DuPage+Lake

[1] 77390

(DuPage+Lake)/2

[1] 38695

McHenry = top_il_income$per_capita_income[3]
      Sangamon = top_il_income$per_capita_income[10] 
      McHenry-Sangamon

[1] 2524

      McHenry+Sangamon

[1] 63712

      (McHenry+Sangamon)/2

[1] 31856

Basic Statistics

mean(il_income$per_capita_income)
median(il_income$per_capita_income)
quantile(il_income$per_capita_income)
# Summary 
summary(il_income)

mean(top_il_income$per_capita_income)

[1] 32918.5

median(top_il_income$per_capita_income)

[1] 31430

quantile(top_il_income$per_capita_income)

      0%      25%      50%      75%     100% 
30594.00 30743.75 31430.00 33103.25 38931.00

# Summary 
summary(top_il_income)

      rank           county  per_capita_income
 Min.   : 2.00   DuPage :1   Min.   :30594    
 1st Qu.: 4.25   Kane   :1   1st Qu.:30744    
 Median :12.00   Kendall:1   Median :31430    
 Mean   :27.10   Lake   :1   Mean   :32918    
 3rd Qu.:41.00   McHenry:1   3rd Qu.:33103    
 Max.   :90.00   McLean :1   Max.   :38931    
                 (Other):4                    
   population         region   
 Min.   :  7032   Min.   :2.0  
 1st Qu.: 36920   1st Qu.:2.0  
 Median :194782   Median :3.0  
 Mean   :334866   Mean   :3.2  
 3rd Qu.:648159   3rd Qu.:4.0  
 Max.   :933736   Max.   :5.0

##### 2) Repeat here the above basic statistics code chunk using instead the data from the file top_il_income

Vectors

Defining a Vector

A sequence of data elements of the same basic type is defined as a vector.

# vector of numeric values
c(2, 3, 5, 8)

# vector of logical values.
c(TRUE, FALSE, TRUE)

# vector of character strings.
c("A", "B", "B-", "C", "D")

Lists

Defining a List

Lists, as opposed to vectors, can hold components of different types.

scores = c(80, 75, 55)  # vector of numeric values                   
grades = c("B", "C", "D-")  # vector of character strings.          

office_hours = c(TRUE, FALSE, FALSE) # vector of logical values.
student = list(scores,grades,office_hours) # list of vectors
student

List Slicing

We can retrieve components of the list with the single square bracket [] operator.

student[1]     
student[2]
student[3]

# first two components of the list
student[1:2]

Member Reference

Using the double square bracket [[]] operator we can reference a member of the list directly. Using one bracket [] would still reference the list but will not allow you to extract a particular member of the list.

student[[1]] # Components of the Scores Vector

First element of the Scores vector

student[[1]][1]

First three elements of the Scores vector

student[[1]][1:3]

student[[1]][2]

##### 3) Repeat here the above code chunk to extract instead the second element of the grades vector

Named List Members

It’s possible to assign names to list members and reference them by names instead of by numeric indexes.

student = list(myscores = scores, mygrades = grades , myoffice_hours = office_hours) 

student
student$myscores
student$mygrades
student$myoffice_hours

Matrices

All columns in a matrix must have the same data type and the same length.

Create a numeric matrix of 5 rows and 4 columns made of sequential numbers 1:20

x_mat = matrix(1:20, nrow=5, ncol=4)
x_mat

Retrieve the 4th column of matrix

x_mat[,4]

Retrieve the 3rd row of matrix

x_mat[3,]

Retrieve rows 2,3,4 of columns 1,2,3

x_mat[2:4,1:3]

x_mat[3,3]

[1] 13

##### 4) Repeat here the above code chunk to extract instead the third row and third column of the matrix

Data Frames

A data frame is more general than a matrix, in that different columns can have different data types (numeric, character, logic, factor). It is a powerful way to work with mixed data structures.

Defining a Data Frame

When we need to store data in table form, we use data frames, which are created by combining lists of vectors of equal length. The variables of a data set are the columns and the observations are the rows.

The str() function helps us to display the internal structure of any R data structure or object to make sure that it’s correct.

str(il_income)

Creating a Data Frame

Snapshot of the solar system.

name = c("Earth", "Mars", "Jupiter")
type = c("Terrestrial","Terrestrial", "Gas giant")
diameter = c(1, 0.532, 11.209)
rotation = c(1, 1.03, 0.41)
rings = c(FALSE, FALSE, TRUE)

Now, by combining the vectors of equal size, we can create a data frame object.

planets_df = data.frame(name,type,diameter,rotation,rings)
planets_df

Suggested Exercises & Resources

Exercises

Datacamp - Learn Data Science from your browser: https://www.datacamp.com/courses/free-introduction-to-r
R-tutor - An R intro to stats that explains basic R concepts: http://www.r-tutor.com/r-introduction

Data Sources

Data samples used in this worksheet were downloaded from the U.S. Census Bureau American FactFinder site.

“SELECTED ECONOMIC CHARACTERISTICS 2006-2010 American Community Survey 5-Year Estimates” - U.S. Census Bureau. Retrieved 2016-09-09: https://factfinder.census.gov/faces/nav/jsf/pages/index.xhtml

LS0tCnRpdGxlOiAiQlNBRDM0MyBGYWxsIDIwMTggTGFiIFdvcmtzaGVldCAwMSIKYXV0aG9yOiAiU2t5bGVyIFN1bW1lcnNvbiIKZGF0ZTogIlNlcHRlbWViZXIgMTJ0aCAyMDE4IgpvdXRwdXQ6CiAgaHRtbF9ub3RlYm9vazogZGVmYXVsdAogIHBkZl9kb2N1bWVudDogZGVmYXVsdApzdWJ0aXRsZTogSW50cm9kdWN0aW9uIHRvIFIgKGJzYWQtbGFiMDEpCi0tLQoKIyMjIEFib3V0CgpSIGlzIGEgbGFuZ3VhZ2UgYW5kIGVudmlyb25tZW50IGZvciBzdGF0aXN0aWNhbCBjb21wdXRpbmcgYW5kIGdyYXBoaWNzLiBSIHByb3ZpZGVzIGEgd2lkZSB2YXJpZXR5IG9mIHN0YXRpc3RpY2FsIChsaW5lYXIgYW5kIG5vbmxpbmVhciBtb2RlbGxpbmcsIGNsYXNzaWNhbCBzdGF0aXN0aWNhbCB0ZXN0cywgdGltZS1zZXJpZXMgYW5hbHlzaXMsIGNsYXNzaWZpY2F0aW9uLCBjbHVzdGVyaW5nLCAuKSBhbmQgZ3JhcGhpY2FsIHRlY2huaXF1ZXMsIGFuZCBpcyBoaWdobHkgZXh0ZW5zaWJsZS4gCgpUaGlzIG5vdGVib29rIGlzIGEgdHV0b3JpYWwgb24gaG93IHRvIHVzZSBSLgoKCiMjIyBTZXR1cAoKUmVtZW1iZXIgdG8gYWx3YXlzIHNldCB5b3VyIHdvcmtpbmcgZGlyZWN0b3J5IHRvIHRoZSBzb3VyY2UgZmlsZSBsb2NhdGlvbi4gR28gdG8gJ1Nlc3Npb24nLCBzY3JvbGwgZG93biB0byAnU2V0IFdvcmtpbmcgRGlyZWN0b3J5JywgYW5kIGNsaWNrICdUbyBTb3VyY2UgRmlsZSBMb2NhdGlvbicuIFJlYWQgY2FyZWZ1bGx5IHRoZSBiZWxvdyBhbmQgZm9sbG93IHRoZSBpbnN0cnVjdGlvbnMgdG8gY29tcGxldGUgdGhlIHRhc2tzIGFuZCBhbnN3ZXIgYW55IHF1ZXN0aW9ucy4gIFN1Ym1pdCB5b3VyIHdvcmsgdG8gUlB1YnMgYXMgZGV0YWlsZWQgaW4gcHJldmlvdXMgbm90ZXMuIAoKIyMjIE5vdGUKCkZvciB5b3VyIGFzc2lnbm1lbnQgeW91IG1heSBiZSB1c2luZyBkaWZmZXJlbnQgZGF0YSBzZXRzIHRoYW4gd2hhdCBpcyBpbmNsdWRlZCBoZXJlLiBBbHdheXMgcmVhZCBjYXJlZnVsbHkgdGhlIGluc3RydWN0aW9ucyBvbiBTYWthaS4gIEZvciBjbGFyaXR5LCB0YXNrcy9xdWVzdGlvbnMgdG8gYmUgY29tcGxldGVkL2Fuc3dlcmVkIGFyZSBoaWdobGlnaHRlZCBpbiByZWQgY29sb3IgYW5kIG51bWJlcmVkIGFjY29yZGluZyB0byB0aGVpciBwYXJ0aWN1bGFyIHBsYWNlbWVudCBpbiB0aGUgdGFzayBzZWN0aW9uLiAgUXVpdGUgb2Z0ZW4geW91IHdpbGwgbmVlZCB0byBhZGQgeW91ciBvd24gY29kZSBjaHVuay4KCkV4ZWN1dGUgYWxsIGNvZGUgY2h1bmtzLCBwcmV2aWV3LCBwdWJsaXNoLCBhbmQgc3VibWl0IGxpbmsgb24gU2FrYWkuCgotLS0tLS0tLS0tCgohW10oaW1ncy9pbWcwMS5wbmcpCgoKIyBCYXNpY3MgT3BlcmF0aW9ucyAKRmlyc3Qgd2Ugd2lsbCBiZWdpbiB3aXRoIGEgZmV3IGJhc2ljIG9wZXJhdGlvbnMuIAoKIyMgVmFyaWFibGUgYXNzaWdubWVudApXZSBhc3NpZ24gdmFsdWVzIHRvIHZhcmlhYmxlcyB1c2luZyB0aGUgYXNzaWdubWVudCBvcGVyYXRvciAnPScuIEFub3RoZXIgZm9ybSBvZiBhc3NpZ25tZW50LCBtb3JlIGdlbmVyYWwsIGlzIHRoZSAnPC0nIG9wZXJhdG9yLiAKQSB2YXJpYWJsZSBhbGxvd3MgeW91IHRvIHN0b3JlIHZhbHVlcyBvciBhbiBvYmplY3QgKGUuZy4gYSBmdW5jdGlvbikuCgpgYGB7cn0KeCA9IDEyOAp5ID0gMTYKeiA8LSA1CnZhcnMgPSBjKDIsNCw4LDE2LDMyKSAjIENyZWF0ZXMgYSB2ZWN0b3IgbGlzdCB1c2luZyB0aGUgZ2VuZXJpYyBjb21iaW5lIGZ1bmN0aW9uICdjJyAKYGBgCgpgYGB7cn0KeCAjIGRpc3BsYXkgdmFsdWUgb2YgdmFyaWFibGUgeAp6ICMgZGlzcGxheXMgdmFsdWUgb2YgdmFyaWFibGUgegpgYGAKCmBgYHtyfQp2YXJzWzFdICNUaGlzIGNhbGxzIHRoZSBmaXJzdCB2YWx1ZSBpbiB0aGUgdmVjdG9yIHZhcnMKYGBgCgpgYGB7cn0KdmFyc1syXSAjVGhpcyBjYWxscyB0aGUgc2Vjb25kIHZhbHVlIGluIHRoZSB2ZWN0b3IgdmFycwpgYGAKCmBgYHtyfQp2YXJzWzE6M10gI1RoaXMgY2FsbHMgdGhlIGZpcnN0IHRocm91Z2ggdGhpcmQgdmFsdWVzIGluIHRoZSB2ZWN0b3IgdmFycwpgYGAKCmBgYHtyfQp2YXJzICNUaGlzIGNhbGxzIHRoZSB2ZWN0b3IgbGlzdApgYGAKCiMjIENvbW1vbiBBcml0aG1ldGljIE9wZXJhdGlvbnMKQmVsb3cgc2hvd3Mgc29tZSBzaW1wbGUgYXJpdGhtZXRpYyBvcGVyYXRpb25zLgpgYGB7cixldmFsPVRSVUV9CjEyKjYKMTI4LzE2CjleMgpgYGAKCiMjIEJhc2ljIERhdGEgVHlwZXMKUiB3b3JrcyB3aXRoIG51bWVyb3VzIGRhdGEgdHlwZXMuIFNvbWUgb2YgdGhlIG1vc3QgYmFzaWMgdHlwZXMgYXJlOiBudW1lcmljLGludGVnZXJzLCBsb2dpY2FsIChCb29sZWFuLWBUUlVFL0ZBTFNFYCkgYW5kIGNoYXJhY3RlcnMgKHN0cmluZy1gIlRFWFQiYCkuCgpgYGB7cixldmFsPVRSVUV9CiNUeXBlOiBDaGFyYWN0ZXIgICAgICAgICAgICAgICAgICAgCiNFeGFtcGxlOiJUUlVFIiwnMjMuNCcKCnYgPSAiVFJVRSIgICAgICAgICAgICAgICAgICAgICAgIApjbGFzcyh2KSAgICAgICAgICAgICAgICAgICAgICAgICAgIAoKI1R5cGU6IE51bWVyaWMgICAgICAgICAgICAgICAgCiNFeGFtcGxlOiAxMi4zLDUKCnYgPSAyMy41ICAgICAgICAgICAgICAgICAgCmNsYXNzKHYpICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgIAojVHlwZTogTG9naWNhbCAgICAKI0V4YW1wbGU6IFRSVUUsRkFMU0UKCnYgPSBUUlVFCmNsYXNzKHYpCgojVHlwZTogRmFjdG9yIChub21pbmFsLCBjYXRlZ29yaWNhbCkKI0V4YW1wbGU6IG0gZiBtIGYgbQoKdiA9IGFzLmZhY3RvcihjKCJtIiwgImYiLCAibSIpKQpjbGFzcyh2KQpgYGAKCiMjICBGdW5jdGlvbnMKClIgRnVuY3Rpb25zIGFyZSBpbnZva2VkIGJ5IGl0cyBuYW1lLCBmb2xsb3dlZCBieSB0aGUgcGFyZW50aGVzaXMsIGFuZCB6ZXJvIG9yIG1vcmUgYXJndW1lbnRzLiAKYGBge3J9CiMgVGhlIGZvbGxvd2luZyBhcHBsaWVzIHRoZSBmdW5jdGlvbiAnYycgKHNlZW4gZWFybGllcikgdG8gY29tYmluZSB0aHJlZSBudW1lcmljIHZhbHVlcyBpbnRvIGEgdmVjdG9yIApjKDEsMiwzKQoKIyBFeGFtcGxlIG9mIGZ1bmN0aW9uIG1lYW4oKSB0byBjYWxjdWxlIHRoZSBtZWFuIG9mIHRocmVlIHZhbHVlcwptZWFuKGMoNSw2LDcpKQoKIyBTcXVhcmUgcm9vdCBvZiBhIG51bWJlcgpzcXJ0KDk5KQpgYGAKCiMjIEltcG9ydGluZyBEYXRhIGFuZCBWYXJpYWJsZSBBc3NpZ25tZW50CgpgYGB7cixldmFsPVRSVUV9CiMgSGVyZSB3ZSBhcmUgcmVhZGluZyBhIGZpbGUgb2YgdHlwZSBjc3YgKGNvbW1hIHNlcGVyYXRlZCB2YWx1ZXMpIHR5cGljYWwgb2YgbWFueSBFeGNlbCBmaWxlcwppbF9pbmNvbWUgPSByZWFkLmNzdihmaWxlID0gImRhdGEvaWxfaW5jb21lLmNzdiIpCnRvcF9pbF9pbmNvbWUgPSByZWFkLmNzdihmaWxlID0gImRhdGEvdG9wX2lsX2luY29tZS5jc3YiKQpgYGAKCgojIyBBcml0aG1ldGljIE9wZXJhdGlvbnMgd2l0aCBEYXRhCgpXZSBjYW4gZXh0cmFjdCB2YWx1ZXMgZnJvbSB0aGUgZGF0YXNldCB0byBwZXJmb3JtIGNhbGN1bGF0aW9ucyBieSByZWZlcmVuY2luZyB0aGUgcHJvcGVyIGVsZW1lbnRzIG9mIGEgbGlzdApgYGB7cixldmFsPVRSVUV9CkR1UGFnZSA9IHRvcF9pbF9pbmNvbWUkcGVyX2NhcGl0YV9pbmNvbWVbMV0KTGFrZSA9IHRvcF9pbF9pbmNvbWUkcGVyX2NhcGl0YV9pbmNvbWVbMl0KRHVQYWdlLUxha2UKRHVQYWdlK0xha2UKKER1UGFnZStMYWtlKS8yCmBgYAoKYGBge3IsZXZhbD1UUlVFfQpNY0hlbnJ5ID0gdG9wX2lsX2luY29tZSRwZXJfY2FwaXRhX2luY29tZVszXQogICAgICBTYW5nYW1vbiA9IHRvcF9pbF9pbmNvbWUkcGVyX2NhcGl0YV9pbmNvbWVbMTBdIAogICAgICBNY0hlbnJ5LVNhbmdhbW9uCiAgICAgIE1jSGVucnkrU2FuZ2Ftb24gCiAgICAgIChNY0hlbnJ5K1NhbmdhbW9uKS8yCmBgYAoKCgojIyBCYXNpYyBTdGF0aXN0aWNzCgpgYGB7cixldmFsPVRSVUV9Cm1lYW4oaWxfaW5jb21lJHBlcl9jYXBpdGFfaW5jb21lKQptZWRpYW4oaWxfaW5jb21lJHBlcl9jYXBpdGFfaW5jb21lKQpxdWFudGlsZShpbF9pbmNvbWUkcGVyX2NhcGl0YV9pbmNvbWUpCiMgU3VtbWFyeSAKc3VtbWFyeShpbF9pbmNvbWUpCmBgYAoKYGBge3IsZXZhbD1UUlVFfQptZWFuKHRvcF9pbF9pbmNvbWUkcGVyX2NhcGl0YV9pbmNvbWUpCm1lZGlhbih0b3BfaWxfaW5jb21lJHBlcl9jYXBpdGFfaW5jb21lKQpxdWFudGlsZSh0b3BfaWxfaW5jb21lJHBlcl9jYXBpdGFfaW5jb21lKQojIFN1bW1hcnkgCnN1bW1hcnkodG9wX2lsX2luY29tZSkKYGBgCgo8c3BhbiBzdHlsZT0iY29sb3I6cmVkIj4KIyMjIyMgMikgUmVwZWF0IGhlcmUgdGhlIGFib3ZlIGJhc2ljIHN0YXRpc3RpY3MgY29kZSBjaHVuayB1c2luZyBpbnN0ZWFkIHRoZSBkYXRhIGZyb20gdGhlIGZpbGUgdG9wX2lsX2luY29tZQo8L3NwYW4+CgojIFZlY3RvcnMKCiMjIERlZmluaW5nIGEgVmVjdG9yCgpBIHNlcXVlbmNlIG9mIGRhdGEgZWxlbWVudHMgb2YgdGhlIHNhbWUgYmFzaWMgdHlwZSBpcyBkZWZpbmVkIGFzIGEgdmVjdG9yLgpgYGB7cixldmFsPVRSVUV9CiMgdmVjdG9yIG9mIG51bWVyaWMgdmFsdWVzCmMoMiwgMywgNSwgOCkKCiMgdmVjdG9yIG9mIGxvZ2ljYWwgdmFsdWVzLgpjKFRSVUUsIEZBTFNFLCBUUlVFKQoKIyB2ZWN0b3Igb2YgY2hhcmFjdGVyIHN0cmluZ3MuCmMoIkEiLCAiQiIsICJCLSIsICJDIiwgIkQiKQpgYGAKCiMgTGlzdHMKCiMjIERlZmluaW5nIGEgTGlzdAoKTGlzdHMsIGFzIG9wcG9zZWQgdG8gdmVjdG9ycywgY2FuIGhvbGQgY29tcG9uZW50cyBvZiBkaWZmZXJlbnQgdHlwZXMuCgpgYGB7cixldmFsPVRSVUV9CnNjb3JlcyA9IGMoODAsIDc1LCA1NSkgICMgdmVjdG9yIG9mIG51bWVyaWMgdmFsdWVzICAgICAgICAgICAgICAgICAgIApncmFkZXMgPSBjKCJCIiwgIkMiLCAiRC0iKSAgIyB2ZWN0b3Igb2YgY2hhcmFjdGVyIHN0cmluZ3MuICAgICAgICAgIAoKb2ZmaWNlX2hvdXJzID0gYyhUUlVFLCBGQUxTRSwgRkFMU0UpICMgdmVjdG9yIG9mIGxvZ2ljYWwgdmFsdWVzLgpzdHVkZW50ID0gbGlzdChzY29yZXMsZ3JhZGVzLG9mZmljZV9ob3VycykgIyBsaXN0IG9mIHZlY3RvcnMKc3R1ZGVudApgYGAKCiMjIExpc3QgU2xpY2luZyAKCldlIGNhbiByZXRyaWV2ZSBjb21wb25lbnRzIG9mIHRoZSBsaXN0IHdpdGggdGhlIHNpbmdsZSBzcXVhcmUgYnJhY2tldCBgW11gIG9wZXJhdG9yLgoKCmBgYHtyLGV2YWw9VFJVRX0Kc3R1ZGVudFsxXSAgICAgCnN0dWRlbnRbMl0Kc3R1ZGVudFszXQoKIyBmaXJzdCB0d28gY29tcG9uZW50cyBvZiB0aGUgbGlzdApzdHVkZW50WzE6Ml0KYGBgCgojIyBNZW1iZXIgUmVmZXJlbmNlCgpVc2luZyB0aGUgZG91YmxlIHNxdWFyZSBicmFja2V0IGBbW11dYCBvcGVyYXRvciB3ZSBjYW4gcmVmZXJlbmNlIGEgbWVtYmVyIG9mIHRoZSBsaXN0IGRpcmVjdGx5LiBVc2luZyBvbmUgYnJhY2tldCBbXSB3b3VsZCBzdGlsbCByZWZlcmVuY2UgdGhlIGxpc3QgYnV0IHdpbGwgbm90IGFsbG93IHlvdSB0byBleHRyYWN0IGEgcGFydGljdWxhciBtZW1iZXIgb2YgdGhlIGxpc3QuCgpgYGB7cixldmFsPVRSVUV9CnN0dWRlbnRbWzFdXSAjIENvbXBvbmVudHMgb2YgdGhlIFNjb3JlcyBWZWN0b3IKYGBgCgoqRmlyc3QgZWxlbWVudCBvZiB0aGUgU2NvcmVzIHZlY3RvcioKCmBgYHtyLGV2YWw9VFJVRX0Kc3R1ZGVudFtbMV1dWzFdCmBgYAoKCipGaXJzdCB0aHJlZSBlbGVtZW50cyBvZiB0aGUgU2NvcmVzIHZlY3RvcioKCmBgYHtyLGV2YWw9VFJVRX0Kc3R1ZGVudFtbMV1dWzE6M10KYGBgCmBgYHtyLGV2YWw9VFJVRX0Kc3R1ZGVudFtbMV1dWzJdCmBgYAoKPHNwYW4gc3R5bGU9ImNvbG9yOnJlZCI+CiMjIyMjIDMpIFJlcGVhdCBoZXJlIHRoZSBhYm92ZSBjb2RlIGNodW5rIHRvIGV4dHJhY3QgaW5zdGVhZCB0aGUgc2Vjb25kIGVsZW1lbnQgb2YgdGhlIGdyYWRlcyB2ZWN0b3IKPC9zcGFuPgoKIyMgTmFtZWQgTGlzdCBNZW1iZXJzCgpJdCdzIHBvc3NpYmxlIHRvIGFzc2lnbiBuYW1lcyB0byBsaXN0IG1lbWJlcnMgYW5kIHJlZmVyZW5jZSB0aGVtIGJ5IG5hbWVzIGluc3RlYWQgb2YgYnkgbnVtZXJpYyBpbmRleGVzLgoKYGBge3IsZXZhbD1UUlVFfQpzdHVkZW50ID0gbGlzdChteXNjb3JlcyA9IHNjb3JlcywgbXlncmFkZXMgPSBncmFkZXMgLCBteW9mZmljZV9ob3VycyA9IG9mZmljZV9ob3VycykgCgpzdHVkZW50CnN0dWRlbnQkbXlzY29yZXMKc3R1ZGVudCRteWdyYWRlcwpzdHVkZW50JG15b2ZmaWNlX2hvdXJzCmBgYAoKIyBNYXRyaWNlcwoKQWxsIGNvbHVtbnMgaW4gYSBtYXRyaXggbXVzdCBoYXZlIHRoZSBzYW1lIGRhdGEgdHlwZSBhbmQgdGhlIHNhbWUgbGVuZ3RoLgoKKkNyZWF0ZSBhIG51bWVyaWMgbWF0cml4IG9mIDUgcm93cyBhbmQgNCBjb2x1bW5zIG1hZGUgb2Ygc2VxdWVudGlhbCBudW1iZXJzIDE6MjAqCmBgYHtyfQp4X21hdCA9IG1hdHJpeCgxOjIwLCBucm93PTUsIG5jb2w9NCkKeF9tYXQKYGBgCgoqUmV0cmlldmUgdGhlIDR0aCBjb2x1bW4gb2YgbWF0cml4KgpgYGB7cn0KeF9tYXRbLDRdCmBgYAoKCipSZXRyaWV2ZSB0aGUgM3JkIHJvdyBvZiBtYXRyaXgqCmBgYHtyfQp4X21hdFszLF0KYGBgCgoqUmV0cmlldmUgcm93cyAyLDMsNCBvZiBjb2x1bW5zIDEsMiwzKgpgYGB7cn0KeF9tYXRbMjo0LDE6M10KYGBgCmBgYHtyfQp4X21hdFszLDNdCmBgYAoKPHNwYW4gc3R5bGU9ImNvbG9yOnJlZCI+CiMjIyMjIDQpIFJlcGVhdCBoZXJlIHRoZSBhYm92ZSBjb2RlIGNodW5rIHRvIGV4dHJhY3QgaW5zdGVhZCB0aGUgdGhpcmQgcm93IGFuZCB0aGlyZCBjb2x1bW4gb2YgdGhlIG1hdHJpeAo8L3NwYW4+CgoKIyBEYXRhIEZyYW1lcwoKQSBkYXRhIGZyYW1lIGlzIG1vcmUgZ2VuZXJhbCB0aGFuIGEgbWF0cml4LCBpbiB0aGF0IGRpZmZlcmVudCBjb2x1bW5zIGNhbiBoYXZlIGRpZmZlcmVudCBkYXRhIHR5cGVzIChudW1lcmljLCBjaGFyYWN0ZXIsIGxvZ2ljLCBmYWN0b3IpLiBJdCBpcyBhIHBvd2VyZnVsIHdheSB0byB3b3JrIHdpdGggbWl4ZWQgZGF0YSBzdHJ1Y3R1cmVzLgoKIyMgRGVmaW5pbmcgYSBEYXRhIEZyYW1lCgpXaGVuIHdlIG5lZWQgdG8gc3RvcmUgZGF0YSBpbiB0YWJsZSBmb3JtLCB3ZSB1c2UgZGF0YSBmcmFtZXMsIHdoaWNoIGFyZSBjcmVhdGVkIGJ5IGNvbWJpbmluZyBsaXN0cyBvZiB2ZWN0b3JzIG9mIGVxdWFsIGxlbmd0aC4gVGhlIHZhcmlhYmxlcyBvZiBhIGRhdGEgc2V0IGFyZSB0aGUgY29sdW1ucyBhbmQgdGhlIG9ic2VydmF0aW9ucyBhcmUgdGhlIHJvd3MuIAoKClRoZSBgc3RyKClgIGZ1bmN0aW9uIGhlbHBzIHVzIHRvIGRpc3BsYXkgdGhlIGludGVybmFsIHN0cnVjdHVyZSBvZiBhbnkgUiBkYXRhIHN0cnVjdHVyZSBvciBvYmplY3QgdG8gbWFrZSBzdXJlIHRoYXQgaXQncyBjb3JyZWN0LgoKCmBgYHtyLGV2YWw9VFJVRX0Kc3RyKGlsX2luY29tZSkKYGBgCgojIyBDcmVhdGluZyBhIERhdGEgRnJhbWUKClNuYXBzaG90IG9mIHRoZSBzb2xhciBzeXN0ZW0uCgpgYGB7ciwgZXZhbD1UUlVFfQpuYW1lID0gYygiRWFydGgiLCAiTWFycyIsICJKdXBpdGVyIikKdHlwZSA9IGMoIlRlcnJlc3RyaWFsIiwiVGVycmVzdHJpYWwiLCAiR2FzIGdpYW50IikKZGlhbWV0ZXIgPSBjKDEsIDAuNTMyLCAxMS4yMDkpCnJvdGF0aW9uID0gYygxLCAxLjAzLCAwLjQxKQpyaW5ncyA9IGMoRkFMU0UsIEZBTFNFLCBUUlVFKQpgYGAKCk5vdywgYnkgY29tYmluaW5nIHRoZSB2ZWN0b3JzIG9mIGVxdWFsIHNpemUsIHdlIGNhbiBjcmVhdGUgYSBkYXRhIGZyYW1lIG9iamVjdC4KCmBgYHtyLCBldmFsPVRSVUV9CnBsYW5ldHNfZGYgPSBkYXRhLmZyYW1lKG5hbWUsdHlwZSxkaWFtZXRlcixyb3RhdGlvbixyaW5ncykKcGxhbmV0c19kZgoKYGBgCgojIFN1Z2dlc3RlZCBFeGVyY2lzZXMgJiBSZXNvdXJjZXMKCiMjIEV4ZXJjaXNlcwoKKiBEYXRhY2FtcCAtIExlYXJuIERhdGEgU2NpZW5jZSBmcm9tIHlvdXIgYnJvd3NlcjogCmh0dHBzOi8vd3d3LmRhdGFjYW1wLmNvbS9jb3Vyc2VzL2ZyZWUtaW50cm9kdWN0aW9uLXRvLXIKCiogUi10dXRvciAtIEFuIFIgaW50cm8gdG8gc3RhdHMgdGhhdCBleHBsYWlucyBiYXNpYyBSIGNvbmNlcHRzOgpodHRwOi8vd3d3LnItdHV0b3IuY29tL3ItaW50cm9kdWN0aW9uCgojIyBEYXRhIFNvdXJjZXMKRGF0YSBzYW1wbGVzIHVzZWQgaW4gdGhpcyB3b3Jrc2hlZXQgd2VyZSBkb3dubG9hZGVkIGZyb20gdGhlIFUuUy4gQ2Vuc3VzIEJ1cmVhdSBBbWVyaWNhbiBGYWN0RmluZGVyIHNpdGUuCgoqICJTRUxFQ1RFRCBFQ09OT01JQyBDSEFSQUNURVJJU1RJQ1MgMjAwNi0yMDEwIEFtZXJpY2FuIENvbW11bml0eSBTdXJ2ZXkgNS1ZZWFyIEVzdGltYXRlcyIgLSBVLlMuIENlbnN1cyBCdXJlYXUuIFJldHJpZXZlZCAyMDE2LTA5LTA5OgpodHRwczovL2ZhY3RmaW5kZXIuY2Vuc3VzLmdvdi9mYWNlcy9uYXYvanNmL3BhZ2VzL2luZGV4LnhodG1sCg==

BSAD343 Fall 2018 Lab Worksheet 01

Introduction to R (bsad-lab01)

Skyler Summerson

Septemeber 12th 2018