#
# Import under R data from the Nucl folder at 
# http://pbil.univ-lyon1.fr/members/lobry/repro/cabios96/Nucl/
#
path <- "http://pbil.univ-lyon1.fr/members/lobry/repro/cabios96/Nucl"
#
# Bacillus subtilis data:
#
Bslea <- read.table(paste(path, "Bslea.frc", sep = "/"))
rownames(Bslea) <- readLines(paste(path, "Bslea.lst", sep = "/"))
Bslag <- read.table(paste(path, "Bslag.frc", sep = "/"))
rownames(Bslag) <- readLines(paste(path, "Bslag.lst", sep = "/"))
Bs <- rbind(Bslea, Bslag)
Bsfac <- factor(rep(c("lea", "lag"), c(nrow(Bslea), nrow(Bslag))))
names(Bs) <- read.table(paste(path, "codons", sep = "/"))$V1
#
# Escherichia coli data:
#
Eclea <- read.table(paste(path, "Eclea.frc", sep = "/"))
rownames(Eclea) <- readLines(paste(path, "Eclea.lst", sep = "/"))
Eclag <- read.table(paste(path, "Eclag.frc", sep = "/"))
rownames(Eclag) <- readLines(paste(path, "Eclag.lst", sep = "/"))
Ec <- rbind(Eclea, Eclag)
Ecfac <- factor(rep(c("lea", "lag"), c(nrow(Eclea), nrow(Eclag))))
names(Ec) <- read.table(paste(path, "codons", sep = "/"))$V1
#
# Haemophilus influenzae data:
#
Hilea <- read.table(paste(path, "Hilea.frc", sep = "/"))
rownames(Hilea) <- readLines(paste(path, "Hilea.lst", sep = "/"))
Hilag <- read.table(paste(path, "Hilag.frc", sep = "/"))
rownames(Hilag) <- readLines(paste(path, "Hilag.lst", sep = "/"))
Hi <- rbind(Hilea, Hilag)
Hifac <- factor(rep(c("lea", "lag"), c(nrow(Hilea), nrow(Hilag))))
names(Hi) <- read.table(paste(path, "codons", sep = "/"))$V1
#
# Mycoplasma genitalium data (variant genetic code):
#
Mglea <- read.table(paste(path, "Mglea.frc", sep = "/"))
rownames(Mglea) <- readLines(paste(path, "Mglea.lst", sep = "/"))
Mglag <- read.table(paste(path, "Mglag.frc", sep = "/"))
rownames(Mglag) <- readLines(paste(path, "Mglag.lst", sep = "/"))
Mg <- rbind(Mglea, Mglag)
Mgfac <- factor(rep(c("lea", "lag"), c(nrow(Mglea), nrow(Mglag))))
names(Mg) <- read.table(paste(path, "codonsmg", sep = "/"))$V1