Permalink
Switch branches/tags
Nothing to show
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
38 lines (29 sloc) 1013 Bytes
### combine CP/CM output
### M. Ragonnet
### 5th Feb 2014
## CM is cluster matcher
## CP is cluster picker
## the files output from each program should be in the same folder
setwd("")
path <- getwd()
print(path)
## read the files
#read the cm file
txtFolder <- dir(path)
inds <- grep("clustInfo", txtFolder)
clustInfoFile <- txtFolder[inds]
cm <- read.csv(clustInfoFile, header=TRUE)
#read the CP log file
cp_output <- readLines(list.files(pattern="*clusterPicks_log.txt"))
cp_output2 <- cp_output[20:(length(cp_output)-1)]
writeLines(cp_output2, "CP_output.txt")
cp <- read.table("CP_output.txt", sep="\t", header=TRUE)
colnames(cp)[1] <- "Clust_ID"
# merge the files into one superfile and write to folder
cluster_table <- merge(cp,cm,by=c("Clust_ID"), all=TRUE)
for (line in 1:length(cluster_table[,1])){
if(cluster_table$NumberOfTips [line] != cluster_table$Num_Seqs[line]){
print("There was a problem combining the files!")
}
}
write.csv (cluster_table, file ="cpcm.csv", row.names = FALSE)