Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

assignment 2 #245

Open
wants to merge 1 commit into
base: master
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
91 changes: 68 additions & 23 deletions Assignment 2-2020.Rmd
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
---
title: "Assignment 2"
author: "Charles Lang"
author: "Shuying Xiong"
date: "September 24, 2020"
output: html_document
---
Expand Down Expand Up @@ -69,7 +69,8 @@ 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))
D3 <- D1 %>% group_by(year) %>%
summarise(mean_key = mean(key.points))

plot(D3$year, D3$mean_key, type = "l", lty = "dashed")

Expand All @@ -92,25 +93,38 @@ pairs(D5)

```{r}
#rnorm(100, 75, 15) creates a random sample with a mean of 75 and standard deviation of 20

#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,20)
hist (score, breaks=30)
S1 <- data.frame(score)
library(dplyr)
S1 <- filter(S1,score<=100)
hist(S1$score)
S2 <- data.frame(rep(100,11))
names(S2) <- "score"
S3 <- bind_rows(S1,S2)
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)
```

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.

label <- letters[1:10]
S3$breaks <- cut(S3$score, breaks=10,labels=label)
```

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.
Expand All @@ -123,56 +137,65 @@ library(RColorBrewer)
#Make RColorBrewer palette available to R and assign to your bins

#Use named palette in histogram

display.brewer.all()
S3$colors <- brewer.pal(10,"Set3")
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

interest.col <- brewer.pal(4,"Dark2")
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$col1 <- ifelse(S3$interest == "music","red","green")
```


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}

plot(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}

plot(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.

```{r}
library(dplyr)
library(tidyr)
```
### Instructions

1. Insert a new code block
2. Create a data frame from the `swirl-data.csv` file called `DF1`

1. Insert a new code block
2. Create a data frame from the `swirl-data.csv` file called `DF1`

```{r}
DF1 <- read.csv("swirl-data.csv", header=TRUE)
```
The variables are:

`course_name` - the name of the R course the student attempted
Expand All @@ -185,19 +208,41 @@ The variables are:
`hash` - anonymyzed student ID

3. Create a new data frame that only includes the variables `hash`, `lesson_name` and `attempt` called `DF2`

```{r}
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(attempt=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

6. Convert `DF3` to this format

```{r}
DF3 <- na.omit(DF3)
DF3b <- spread(DF3, lesson_name,attempt)
```
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 <- ifelse(DF4$correct==TRUE,1,ifelse(is.na(DF4$correct)==TRUE,NA,0))
```

9. Create a new data frame called `DF5` that provides a mean score for each student on each course

```{r}
DF5 <- DF4 %>% group_by (hash,lesson_name) %>% summarise (mean.correct = mean(correct, na.rm=TRUE))
```
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$correct <- ifelse(DF6$correct == TRUE, 1,0)
DF6$datetime <- as.POSIXlt (DF6$datetime,origin="1970-01-01 00:00:00 UTC")
DF6$datetime <- strftime (DF6$datetime, format="%b:%e")
DF7 <- DF6 %>% group_by(datetime) %>% summarise (av.correct=mean(correct,na.rm=TRUE))
DF7$datetime <- strptime(DF7$datetime,format="%b:%e")
plot(DF7$datetime,DF7$av.correct)
```
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.
Loading