P55.R

In fredfeng/MorpheusData: Data sets for Morpheus

# making table data sets
library(dplyr)
library(tidyr)
library(MorpheusData)

#############benchmark 55
dat <- read.table(text=
"
   Geotype   Strategy Year.1 Year.2 
A     Demand      1      5     
A Strategy_1      2      6    
A Strategy_2      3      7   
B     Demand      8      8  
B Strategy_1      9      9 
B Strategy_2     10     10
", header=T)



write.csv(dat, "data-raw/p55_input1.csv", row.names=FALSE)

df_out = dat %>%
    gather(key, value, - Geotype, - Strategy) %>%
    filter(Strategy!="Demand") %>% group_by(Geotype, key) %>%
      summarize(sum = sum(value)) %>% spread(key, sum)

write.csv(df_out, "data-raw/p55_output1.csv", row.names=FALSE)

p55_output1 <- read.csv("data-raw/p55_output1.csv", check.names = FALSE)
fctr.cols <- sapply(p55_output1, is.factor)
int.cols <- sapply(p55_output1, is.integer)

p55_output1[, fctr.cols] <- sapply(p55_output1[, fctr.cols], as.character)
p55_output1[, int.cols] <- sapply(p55_output1[, int.cols], as.numeric)
save(p55_output1, file = "data/p55_output1.rdata")

p55_input1 <- read.csv("data-raw/p55_input1.csv", check.names = FALSE)
fctr.cols <- sapply(p55_input1, is.factor)
int.cols <- sapply(p55_input1, is.integer)

p55_input1[, fctr.cols] <- sapply(p55_input1[, fctr.cols], as.character)
p55_input1[, int.cols] <- sapply(p55_input1[, int.cols], as.numeric)
save(p55_input1, file = "data/p55_input1.rdata")