owl/decoyScoring/score_corr_plots.R



DECOYSETS<-c("4state_reduced", "fisa", "fisa_casp3", "hg_structal", "ig_structal", "ig_structal_hires", "lattice_ssfit")


for (set in DECOYSETS) {
        files<-list.files(path = ".", pattern = paste(set,"_.*.scores",sep=""))
        
        cols=0  
        if (length(files)%%2==0) {
                cols=length(files)/2
        } else {
                cols=(length(files)+1)/2
        }
        par(mfrow = c(cols, cols))
        
        for (file in files) {
                data<-read.table(file)
                plot(data$V2,data$V3,xlab="score",ylab="rmsd",main=file)
        }

}
Revision:	1006
Committed:	Wed Mar 31 13:00:35 2010 UTC (10 years, 10 months ago) by hstehr
File size:	467 byte(s)
Log Message:	refactoring: created packages runners.blast, runners.gromacs, structure.graphs, moved actionTools to util, structure.decoyScoring to owl.decoyScoring, tinker to runners.tinker
Line	User	Rev	File contents
1	duarte	957
2
3			DECOYSETS<-c("4state_reduced", "fisa", "fisa_casp3", "hg_structal", "ig_structal", "ig_structal_hires", "lattice_ssfit")
4
5
6			for (set in DECOYSETS) {
7			files<-list.files(path = ".", pattern = paste(set,"_.*.scores",sep=""))
8
9			cols=0
10			if (length(files)%%2==0) {
11			cols=length(files)/2
12			} else {
13			cols=(length(files)+1)/2
14			}
15			par(mfrow = c(cols, cols))
16
17			for (file in files) {
18			data<-read.table(file)
19			plot(data$V2,data$V3,xlab="score",ylab="rmsd",main=file)
20			}
21
22			}