Thanks to visit codestin.com
Credit goes to github.com

Skip to content
Open
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
36 changes: 36 additions & 0 deletions week2/programming1.R
Original file line number Diff line number Diff line change
@@ -1,3 +1,39 @@
source("http://d396qusza40orc.cloudfront.net/rprog%2Fscripts%2Fsubmitscript1.R")
submit()

> setwd("C:/Users/Dad/Documents/R/JHop")
> ## develope file reading
> file <- read.csv("./specdata/023.csv")
> ## create data frame
> df <- data.frame(Date=as.Date(character()),
+ sulfate=numeric(),
+ nitrate=numeric(),
+ ID=integer())
> ## variable name for directory - combine three dirrerent
> directory <- "specdata"
> id <- 23
> file <- read.csv(paste("./", directory, "/", sprintf("%03d", id), ".csv", sep=""))
> df<-rbind(df, file)
> id<-1
> file1<- read.csv(paste("./", directory, "/", sprintf("%03d", id), ".csv", sep=""))
> df<-rbind(df, file1)
> id<-2
> file2<- read.csv(paste("./", directory, "/", sprintf("%03d", id), ".csv", sep=""))
> df<-rbind(df, file2)
> ## - mean
>
> mean(file$nitrate, na.rm=TRUE) #calculate mean
[1] 1.280833
>
> #develop column name variable to calculate mean
> pollutant <- "nitrate"
> avector <- as.vector(file[pollutant]) #extract column to vector
> avector2 <-avector[!is.na(avector)] #remove and N/A from vector
> mean(avector2)
[1] 1.280833
>
> #combining statemente to reduce coding
> pollutant <- "nitrate"
> avector <- as.vector(file[pollutant]) #extract column to vector
> mean(avector[!is.na(avector)]) #remove any N/A and calculate mean
[1] 1.280833