forked from core-methods-in-edm/assignment2
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Assignment 2-2020.Rmd
263 lines (185 loc) · 8.11 KB
/
Assignment 2-2020.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
---
title: "Assignment 2"
author: "JINGSHU ZHANG"
date: "OCTOBER 3, 2020"
output: html_document
---
#Part I
## Data Wrangling
In the hackathon a project was proposed to collect data from student video watching, a sample of this data is available in the file video-data.csv.
stid = student id
year = year student watched video
participation = whether or not the student opened the video
watch.time = how long the student watched the video for
confusion.points = how many times a student rewatched a section of a video
key,points = how many times a student skipped or increased the speed of a video
```{r}
#Install the 'tidyverse' package or if that does not work, install the 'dplyr' and 'tidyr' packages.
#Load the package(s) you just installed
library(tidyverse)
library(tidyr)
library(dplyr)
D1 <- read.csv("video-data.csv", header = TRUE)
#Create a data frame that only contains the years 2018
D2 <- filter(D1, year == 2018)
```
## Histograms
```{r}
#Generate a histogram of the watch time for the year 2018
hist(D2$watch.time)
#Change the number of breaks to 100, do you get the same impression?
hist(D2$watch.time, breaks = 100)
#Cut the y-axis off at 10
hist(D2$watch.time, breaks = 100, ylim = c(0,10))
#Restore the y-axis and change the breaks so that they are 0-5, 5-20, 20-25, 25-35
hist(D2$watch.time, breaks = c(0,5,20,25,35))
```
## Plots
```{r}
#Plot the number of confusion points against the watch time
plot(D1$confusion.points, D1$watch.time)
#Create two variables x & y
x <- c(1,3,2,7,6,4,4)
y <- c(2,4,2,3,2,4,3)
#Create a table from x & y
table1 <- table(x,y)
#Display the table as a Barplot
barplot(table1)
#Create a data frame of the average total key points for each year and plot the two against each other as a lines
D3 <- D1 %>% group_by(year) %>% summarise(mean_key = mean(key.points))
plot(D3$year, D3$mean_key, type = "l", lty = "dashed")
#Create a boxplot of total enrollment for three students
D4 <- filter(D1, stid == 4|stid == 20| stid == 22)
#The drop levels command will remove all the schools from the variable with no data
D4 <- droplevels(D4)
boxplot(D4$watch.time~D4$stid, xlab = "Student", ylab = "Watch Time")
```
## Pairs
```{r}
#Use matrix notation to select columns 2, 5, 6, and 7
D5 <- D1[,c(2,5,6,7)]
#Draw a matrix of plots for every combination of variables
pairs(D5)
```
## Part II
1. Create a simulated data set containing 100 students, each with a score from 1-100 representing performance in an educational game. The scores should tend to cluster around 75. Also, each student should be given a classification that reflects one of four interest groups: sport, music, nature, literature.
```{r}
#rnorm(100, 75, 15) creates a random sample with a mean of 75 and standard deviation of 15
#filter() can be used to set a maximum and minimum value
#round() rounds numbers to whole number values
#sample() draws a random samples from the groups vector according to a uniform distribution
score<-rnorm(100,75,15)
hist(score,breaks = 30)
S1<-data.frame(score)
library(dplyr)
S1<-filter(S1,score<=100)
hist(S1$score)
S2<-data.frame(rep(100,5))
names(S2)<-"score"
S3 <- bind_rows(S1,S2)
S3$score<-round(S3$score,0)
```
2. Using base R commands, draw a histogram of the scores. Change the breaks in your histogram until you think they best represent your data.
```{r}
hist(S3$score, breaks = 10)
```
3. Create a new variable that groups the scores according to the breaks in your histogram.
```{r}
#cut() divides the range of scores into intervals and codes the values in scores according to which interval they fall. We use a vector called `letters` as the labels, `letters` is a vector made up of the letters of the alphabet.
lable<-letters[1:10]
S3$letterbreaks<-cut(S3$score,breaks = 10,labels = lable)
```
4. Now using the colorbrewer package (RColorBrewer; http://colorbrewer2.org/#type=sequential&scheme=BuGn&n=3) design a pallette and assign it to the groups in your data on the histogram.
```{r}
library(RColorBrewer)
#Let's look at the available palettes in RColorBrewer
display.brewer.all()
#The top section of palettes are sequential, the middle section are qualitative, and the lower section are diverging.
#Make RColorBrewer palette available to R and assign to your bins
S3$colors <- brewer.pal(10,"Set3")
#Use named palette in histogram
hist(S3$score, col = S3$colors)
```
5. Create a boxplot that visualizes the scores for each interest group and color each interest group a different color.
```{r}
#Make a vector of the colors from RColorBrewer
S3$score<-ifelse(S3$score == 1,100,S3$score)
S3$score <- round(S3$score,0)
interest <- c("sport","music","nature","literature")
S3$interest <- sample(interest,100,replace = TRUE)
S3$stid <- seq(1,100,1)
interest.col <- brewer.pal(4,"Set3")
boxplot(score ~ interest, S3, col = interest.col)
```
6. Now simulate a new variable that describes the number of logins that students made to the educational game. They should vary from 1-25.
```{r}
S3$login <- sample(1:25,100,replace = TRUE)
```
7. Plot the relationships between logins and scores. Give the plot a title and color the dots according to interest group.
```{r}
plot(S3$login, S3$score, col=S3$colors, main = "students login VS. scores")
S3$coll<-ifelse(S3$interest=="music","blue","red")
```
8. R contains several inbuilt data sets, one of these in called AirPassengers. Plot a line graph of the the airline passengers over time using this data set.
```{r}
Airpassengers<- data.frame(AirPassengers)
plot(Airpassengers,main = "Airline Passengers over time",ylab = "airpassengers")
```
9. Using another inbuilt data set, iris, plot the relationships between all of the variables in the data set. Which of these relationships is it appropraiet to run a correlation on?
```{r}
pairs(iris)
```
# Part III - Analyzing Swirl
## Data
In this repository you will find data describing Swirl activity from the class so far this semester. Please connect RStudio to this repository.
### Instructions
1. Insert a new code block
```{r}
```
2. Create a data frame from the `swirl-data.csv` file called `DF1`
The variables are:
`course_name` - the name of the R course the student attempted
`lesson_name` - the lesson name
`question_number` - the question number attempted
`correct` - whether the question was answered correctly
`attempt` - how many times the student attempted the question
`skipped` - whether the student skipped the question
`datetime` - the date and time the student attempted the question
`hash` - anonymyzed student ID
```{r}
DF1<- read.csv("swirl-data.csv", header = TRUE)
```
3. Create a new data frame that only includes the variables `hash`, `lesson_name` and `attempt` called `DF2`
```{r}
library(dplyr)
DF2<-select(DF1,hash,lesson_name, attempt)
```
4. Use the `group_by` function to create a data frame that sums all the attempts for each `hash` by each `lesson_name` called `DF3`
```{r}
DF3<-DF2 %>% group_by(hash, lesson_name) %>% summarise(attempts = sum(attempt))
```
5. On a scrap piece of paper draw what you think `DF3` would look like if all the lesson names were column names
```{r}
table<-table(DF3$lesson_name,col=DF3$lesson_name)
```
6. Convert `DF3` to this format
```{r}
```
7. Create a new data frame from `DF1` called `DF4` that only includes the variables `hash`, `lesson_name` and `correct`
```{r}
DF4<-select(DF1,hash,lesson_name,correct)
```
8. Convert the `correct` variable so that `TRUE` is coded as the **number** `1` and `FALSE` is coded as `0`
```{r}
DF4$correct <- factor(x=DF4$correct,levels=c('FALSE','TRUE'),labels = c(0,1),exclude = NA)
```
9. Create a new data frame called `DF5` that provides a mean score for each student on each course
```{r}
DF5<- DF1 %>% group_by(course_name,hash) %>% summarise(mean_score = mean(score))
```
10. **Extra credit** Convert the `datetime` variable into month-day-year format and create a new data frame (`DF6`) that shows the average correct for each day
```{r}
DF6<-select(DF1,correct,datetime)
DF6$datetime<-as.Date("2000-01-01")
```
Finally use the knitr function to generate an html document from your work. Commit, Push and Pull Request your work back to the main branch of the repository. Make sure you include both the .Rmd file and the .html file.