# Reformat Example -- reformat.r source code.

# Read data from input file.
orig = read.table("reading-paired.txt", header=TRUE)

# Print original dataset.
cat("Orig Dataset\n")
print(orig)

# Create Reformat dataset.
reformat = data.frame(Scores=c(orig$Standard, orig$New), 
                      Method=c(rep("std", length(orig$Standard)), 
                               rep("new", length(orig$New))))

# Print Reformat dataset. 
cat("Reformat Dataset\n")  
print(reformat)

# Redirect graphics to pdf file.
pdf("sidebysidebp.pdf")

# Create side-by-side boxplots.
boxplot(Scores ~ Method, data=reformat)

# Turn off graphics redirection.
dev.off( )