Rstudio

  Download and go through Week2Part2 video and follow the process I have to demo your dataset. Start from Week2_Part2.R script, change it according to your selected datset from Quiz2. Then submit to this link: 1. Your modified script. 2. A screenshot report showing the outputs from your rstudio after executing all activities. Week2_Part2.R script:    setwd(“C:/Users/ialsmadi/Desktop/University_of_Cumberlands/Lectures/Week2/RScripts”) getwd() # Import test data data<-read.csv("yearly_sales.csv") #A 5-number summary is a set of 5 descriptive statistics for summarizing a continuous univariate data set.  #It consists of the data set's: minimum, 1st quartile, median, 3rd quartile, maximum #Find the set, L, of data below the median. The 1st quartile is the median of L. #Find the set, U, of data above the median. The 3rd quartile is the median of U. print(summary(data)) anscombe<-read.csv("anscombe.csv") print(summary(anscombe)) sd(anscombe$X) var(anscombe$X) sd(anscombe$x1) var(anscombe$x1) sd(anscombe$x2) var(anscombe$x2) sd(anscombe$x3) var(anscombe$x3) sd(anscombe$x4) var(anscombe$x4) sd(anscombe$y1) var(anscombe$y1) sd(anscombe$y2) var(anscombe$y2) sd(anscombe$y3) var(anscombe$y3) ##-- now some "magic" to do the 4 regressions in a loop: ff <- y ~ x mods <- setNames(as.list(1:4), paste0("lm", 1:4)) for(i in 1:4) { ff[2:3] <- lapply(paste0(c("y","x"), i), as.name) ## or ff[[2]] <- as.name(paste0("y", i)) ## ff[[3]] <- as.name(paste0("x", i)) mods[[i]] <- lmi <- lm(ff, data = anscombe) print(anova(lmi)) } ## See how close they are (numerically!) sapply(mods, coef) lapply(mods, function(fm) coef(summary(fm))) ## Now, do what you should have done in the first place: PLOTS op <- par(mfrow = c(2, 2), mar = 0.1+c(4,4,1,1), oma = c(0, 0, 2, 0)) for(i in 1:4) { ff[2:3] <- lapply(paste0(c("y","x"), i), as.name) plot(ff, data = anscombe, col = "red", pch = 21, bg = "orange", cex = 1.2, xlim = c(3, 19), ylim = c(3, 13)) abline(mods[[i]], col = "blue") } mtext("Anscombe's 4 Regression data sets", outer = TRUE, cex = 1.5) par(op) plot(sort(data$num_of_orders)) hist(sort(data$num_of_orders)) plot(density(sort(data$num_of_orders))) plot(sort(data$gender)) hist(sort(data$sales_total)) plot(density(sort(data$sales_total))) library(lattice) densityplot(data$num_of_orders) # top plot # bottom plot as log10 is actually # easier to read, but this plot is in natural log densityplot(log(data$num_of_orders)) densityplot(data$sales_total) densityplot(log(data$sales_total)) hist(data$sales_total, breaks=100, main="Sales total",  xlab="sales", col="gray") # draw a line for the media abline(v = median(data$sales_total), col = "magenta", lwd = 4) # use rug() function to see the actual datapoints rug(data$sales_total) #Boxplots can be created for individual variables or for variables by group.  #The format is boxplot(x, data=), where x is a formula and data= denotes the data frame providing  #the data. boxplot(data$sales_total,data=data, main="Dis by Sales",  xlab="Sales", ylab="Total") # Boxplot of MPG by Car Cylinders, using one of R built-in datasets  boxplot(mpg~cyl,data=mtcars, main="Car Milage Data",  xlab="Number of Cylinders", ylab="Miles Per Gallon") #in our boxplot above, we might want to draw a horizontal line at 12 where the national standard is. abline(h = 12) boxplot(data$sales_total,data=data, main="Total sales Bplot",  xlab="Sales", ylab="Total") # Dot chart of a single numeric vector dotchart(mtcars$mpg, labels = row.names(mtcars), cex = 0.6, xlab = "mpg") #install.packages("ROCR") #library(ROCR) # Simple Scatterplot attach(mtcars) plot(wt, mpg, main="Scatterplot Example",  xlab="Car Weight ", ylab="Miles Per Gallon ", pch=19) #The R function abline() can be used to add vertical, horizontal or regression lines to a graph plot(data$sales_total, data$gender) # Add fit lines abline(lm(data$sales_total~ data$num_of_orders), col="red") # regression line (y~x)  lines(lowess(data$sales_total, data$num_of_orders), col="blue") # lowess line (x,y) # Basic Scatterplot Matrix pairs(data) pairs(data[0:2]) # Scatterplot Matrices from the car Package install.packages("car") library(car) install.packages("ggplot2") library(ggplot2) quit()