[R] 1st script

[memilanuk]

So... wrote my first script, rather than just using the interactive 
console.  I think I got everything working more or less the way I want, 
but I'm sure there's a ton of room for improvement.  Specifically in the 
way of automation - but thats where I kind of ran out of steam.  Any 
suggestions would be much appreciated.

case_data.r

# Import CSV file into a data frame.
case_weights <- read.csv(file = "case_weights.csv")

# For each row, take the number in the Weight column and replicate it
# as many times as there are in each count column.
LC09 <- rep(case_weights$Weight, case_weights$LC09)
LC10 <- rep(case_weights$Weight, case_weights$LC10)
LP14b1 <- rep(case_weights$Weight, case_weights$LP14b1)
LP14b2 <- rep(case_weights$Weight, case_weights$LP14b2)

# Determine the longest vector, to help with the next step.
max.len <- max(length(LC09), length(LC10), length(LP14b1),
		length(LP14b2))

# Pad each vector with NA so they are all the same length and will
# go in a data frame.
LC09 <- c(LC09, rep(NA, max.len - length(LC09)))
LC10 <- c(LC10, rep(NA, max.len - length(LC10)))
LP14b1 <- c(LP14b1, rep(NA, max.len - length(LP14b1)))
LP14b2 <- c(LP14b2, rep(NA, max.len - length(LP14b2)))

# Stick everything back into one data frame.
case_dat <- data.frame(LC09)
case_dat$LC10 <- LC10
case_dat$LP14b1 <- LP14b1
case_dat$LP14b2 <- LP14b2

# Stuff said data frame back into a CSV for use elsewhere (plot.ly).
write.csv(case_dat, file = "expanded_case_weights.csv")

# Boxplot it
boxplot(case_dat, varwidth = TRUE, notch = TRUE, horizontal = TRUE,
         main = "Case Weights", xlab = "Weight (grains)",
         ylab = "Batch", las = 1, names = c("LC09", "LC10", "LP14b1",
	"LP14b2"))



-- 
Shiny!  Let's be bad guys.

Reach me @ memilanuk (at) gmail dot com