Module 8

 #Module 8 Assignment

##Question 1

##Data input for stress level vectors


HighStress <- c(10,9,8,9,10,8)


ModerateStress <- c(8,10,6,7,8,8)


LowStress <- c(4,6,6,4,2,2)


#generate data frame to help incorpirate all three variables


stress <- data.frame(cbind(HighStress, ModerateStress, LowStress))


stressStack <- stack(stress)


#Next, I will run an ANOVA test on the dataframe created from the inital data for each stress level


oneway.test(values~ind,data = stressStack)



## Remember this assumes a non-equal variance (note to self: See personal notepad as to why)


#Next WE will perform an ANOVA test o the data derived from the initial Stress level vectors


aovResult <- aov(values~ind,data = stressStack)


summary(aovResult)


##Key for results in Console

#df = degrees of freedom

#SumSq = Sum of the Squares from the variable

#Mean sq = sum sq / df

#F-Value = mean sq var / mean sq residual

#p-value associated with the f-value denoted by Pr(>F)


# Because the p-value is less than a 0.05 

##The results are signifigant and the null hypothesis is rejected







#Question 2

##First I need to set up to zelazo dataset


library(ISwR)

data("zelazo")

##First line of buisiness involves creating a dataframe and employing cbind for the columns

##We will then conduct an independant t-test with the dataframe and then an ANOVA test


##The data set zelazo is a list of 4 vectors 

##We need to convert the list into dataframecolumns


zdf <- as.data.frame(do.call(cbind, zelazo))


#stack zdf


zdfstsack <- stack(zdf)


#Independant T-Test on active and passive


t.test(zelazo$active, zelazo$passive)


#Defining the test results to a variable to examine the summary results


A.S.Stat <- t.test(zelazo$active, zelazo$passive)


summary(A.S.Stat)


##Next, an ANOVA test with the datafram zdf


aov(values~ind,data = zdfstack)


##This dataset ran warnings and would not compile or give any entries

##For some reason my computer is having issues internally with this package 


Comments

Post a Comment

Popular posts from this blog

R Package: pfStat

Image
 Final Project: pfStat, a personal finance mini-suite in R GitHub URL: https://github.com/tahousex/pfStat What is pfStat? This package is entitled pfStat, and is a small package with a few functions dedicated to helping a user calculate different personal finance metrics. The functions include the ability to derive Returns on Interest, Tax Rates based on Income, Taxable Income and Home Equity, and the purpose behind choosing these kinds of metrics was to emphasize the importance of knowing where one puts their wallet. It goes without being entirely sanctimonious that the dollar does not go as far as it once did, and knowing where one stands amongst their own finances and parameters in the financial world will assist them in making more calculated decisions.  The Functions: The taxrate() income receives input from the user income and state tax percentage represented as a decimal percentage (0-1)   and gives the user output of a dollar amount of the taxes owed based on the input The taxa
Read more

Module 7 Assignment

Module 7 Assignment  R Code: # Module 7 Assignment #1.1 #First define the data set required, I/O values x <- c(16, 17, 13, 18, 12, 14, 19, 11, 11, 10) y <- c(63, 81, 56, 91, 47, 57, 76, 72, 62, 48) #Then we must create the relationship (correlation) cor(x,y) #This value 0.728 that is in the console tells us a few things #This relationship is correlated and is increasing linearly #I may choose to plot these findings using plot() plot(x,y) #See the lower righthand window for a composite plot for this data, as the regression trednline would positively increase through the points #1.2 #This settion of code will get the value for the coefficients of the linear regression for this data Regression_model=lm(y~x) Regression_model$coefficients #2.1 #This problem will go over the aquisition of similar values with a different data set discharge <- c(3.600,1.800,3.333,2.283,4.533,2.883) waiting <- c(79,54,74,62,85,55) plot(discharge,waiting) #2.1 Regression_model_1=lm(waiting~discharge)
Read more

The Tampa Feasibility Report featuring R based Visualizations

Image
  Visualizations for Professional Writing's Final Deliverable These visualizations were taken from FDOT data for our local district on the subject of transit usage, and crash data that was munged into relative ratios that had a regression performed on these variables within.  Pie Chart The above figure is a pie chart using data collected that estimates the higher end limits of hourly numbers of passengers that could utilize each method of public transportation during peak hours Transit Regression This figure here shows the monthly totals for crashes, fatalities and serious injuries from accidents recorded within FDOT data and the regression line method was used against the 'free_y' inputs which allowed there to be a 95 percent confidence regression on the other two variables from the one in question showing certain trends within and between variables. Bus versus Pubic Transport with 95 percent regression line fitted The same style of regression was sued here but this time i
Read more