From 3587293928088fa1c99602e8daab21a0485bb05b Mon Sep 17 00:00:00 2001 From: Logan Date: Sun, 24 Sep 2017 20:15:05 -0400 Subject: [PATCH 1/5] Finished R Code for parts 1 and 2 --- Tutorial5Challenge.R | 32 ++++++++++++++++++++++++++++++++ 1 file changed, 32 insertions(+) create mode 100644 Tutorial5Challenge.R diff --git a/Tutorial5Challenge.R b/Tutorial5Challenge.R new file mode 100644 index 0000000..804af88 --- /dev/null +++ b/Tutorial5Challenge.R @@ -0,0 +1,32 @@ +#Reading in the Data (make sure your working directory is correct) +wages=read.csv(file="wages.csv", header = TRUE, stringsAsFactors = FALSE) +attach(wages) + +#Part 1 +data1 <- wages[order(gender,yearsExperience),] +data1 <- (data1[,1:2]) +UNIQUE <- unique(data1) +write.table(UNIQUE, file = "Tutorial5Output.dat", sep = " ") + +#Part 2 +data2 <- wages[order(wage),] +min.wage <- data2[1,] +max.wage <- data2[nrow(wages),] +bottom = nrow(wages) -9 +top = nrow(wages) +top.10 <- data2[bottom:top,] + +counter = 0 +for (i in 1:10){ + if (top.10[i,]$gender == "female") { + counter = counter + 1 + } +} +print("The person with the lowest wage is") +min.wage +print("The person with the highest wage is") +max.wage + +print("The number of females in the top 10 earners is") +counter + From ae56fe353d3884d784736ffe7335afb920b6fb76 Mon Sep 17 00:00:00 2001 From: Logan Date: Wed, 27 Sep 2017 10:47:51 -0400 Subject: [PATCH 2/5] The In-Class Practice Problems --- ExtraPractice.R | 12 ++++++++++++ 1 file changed, 12 insertions(+) create mode 100644 ExtraPractice.R diff --git a/ExtraPractice.R b/ExtraPractice.R new file mode 100644 index 0000000..5668b6d --- /dev/null +++ b/ExtraPractice.R @@ -0,0 +1,12 @@ +#Monday Activity: Wage for all males and all females +fwage = 0 +mwage = 0 + +for (i in 1:nrow(wages)){ + if (wages[i,1] == "female"){ + fwage = fwage + wages[i,]$wage + } + else if (wages[i,1] == "male"){ + mwage = mwage + wages[i,]$wage + } +} \ No newline at end of file From 4f3f935a445e0ef8e6f2b2a2e3b7146931e1007e Mon Sep 17 00:00:00 2001 From: Logan Date: Wed, 27 Sep 2017 20:53:15 -0400 Subject: [PATCH 3/5] Finished the DNA Assignment --- ExtraPractice.R | 67 +++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 67 insertions(+) create mode 100644 ExtraPractice.R diff --git a/ExtraPractice.R b/ExtraPractice.R new file mode 100644 index 0000000..62da2d3 --- /dev/null +++ b/ExtraPractice.R @@ -0,0 +1,67 @@ +#Monday Activity: Wage for all males and all females +fwage = 0 +mwage = 0 + +for (i in 1:nrow(wages)){ + if (wages[i,1] == "female"){ + fwage = fwage + wages[i,]$wage + } + else if (wages[i,1] == "male"){ + mwage = mwage + wages[i,]$wage + } +} + +#DNA Script + +#Initializing Data vectors to empty +sequenceID <- NULL +sequenceLENGTH <- NULL +percentCG <- NULL +meltTEMP <- NULL + +#Reading the File +fileName <- "Lecture11.fasta" +readFile<- file(fileName,open="r") #Makng fileName readable +linn <-readLines(readFile) #Setting linn to the lines of conn + + +#Looping through contents of file +j = 1 +k = 1 +for (i in 1:length(linn)){ + linn_split <- strsplit(linn[i], "")[[1]] #Splits line into series of characters to be accessed + + #Conditional to see if the ith line is a sequence ID + if(linn_split[1] == ">"){ + sequenceID[j] = substring(linn[i],5) + j = j + 1 + + #Calculations for the base pairs + } else { + + #Getting length of sequence + sequenceLENGTH[k] = nchar(linn[i], type = "chars") + + #Calculating percentCG + CG = 0 + for (base in linn_split) { + if (base == "G" || base == "C") { + CG = CG + 1 + } + } + percentCG[k] = CG/sequenceLENGTH[k]*100 + + #Calculating melting temp + if (sequenceLENGTH[k] > 14){ + meltTEMP[k] = -9999 + } else { + meltTEMP[k] = 4*CG + 2*(sequenceLENGTH[k] - CG) + } + + k = k + 1 + } +} +close(readFile) + +#Putting everything into a table called DNA.dat +DNA.dat <- data.frame(sequenceID, sequenceLENGTH, percentCG, meltTEMP) \ No newline at end of file From 765b525599040f6c71f3c9cc5a57048d34511340 Mon Sep 17 00:00:00 2001 From: Logan Date: Wed, 27 Sep 2017 20:54:51 -0400 Subject: [PATCH 4/5] Fixed Typos --- ExtraPractice.R | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/ExtraPractice.R b/ExtraPractice.R index 0db1946..9644078 100644 --- a/ExtraPractice.R +++ b/ExtraPractice.R @@ -9,7 +9,7 @@ for (i in 1:nrow(wages)){ else if (wages[i,1] == "male"){ mwage = mwage + wages[i,]$wage } -<<<<<<< HEAD + } #DNA Script @@ -66,6 +66,3 @@ close(readFile) #Putting everything into a table called DNA.dat DNA.dat <- data.frame(sequenceID, sequenceLENGTH, percentCG, meltTEMP) -======= -} ->>>>>>> 2a45a565f982abf8216f12c364dafc420e9f4307 From 0989cd51eb87def974e3740513b942ee737b5826 Mon Sep 17 00:00:00 2001 From: Logan Date: Thu, 28 Sep 2017 13:51:00 -0400 Subject: [PATCH 5/5] Finished --- Tutorial5Challenge.R | 26 ++++++++++++++++++++++++-- 1 file changed, 24 insertions(+), 2 deletions(-) diff --git a/Tutorial5Challenge.R b/Tutorial5Challenge.R index 804af88..8702066 100644 --- a/Tutorial5Challenge.R +++ b/Tutorial5Challenge.R @@ -23,10 +23,32 @@ for (i in 1:10){ } } print("The person with the lowest wage is") -min.wage +min.wage[,c(1,2,4)] print("The person with the highest wage is") -max.wage +max.wage[,c(1,2,4)] print("The number of females in the top 10 earners is") counter +#Part 3 +highSchoolWages <- NULL +collegeWages <- NULL + +j = 1 +k = 1 +for (i in 1:nrow(wages)){ + if (wages[i,3] == "12"){ + highSchoolWages[j] = wages [i,4] + j = j+1 + }else if (wages[i,3] == "16"){ + collegeWages[k] = wages[i,4] + k = k+1 + } +} + +print("The minimum wage of a high school graduate is:") +min(highSchoolWages) +print("The minimum wage of a college graduate is:") +min(collegeWages) +print("The difference in these wages is") +min(collegeWages) - min(highSchoolWages) \ No newline at end of file