@@ -96,18 +96,15 @@ colnames(observed)[1] <- "ID"
9696# select available samples
9797observed_samples <- colnames(observed )[2 : ncol(observed )]
9898print(paste( " Observed samples:" , paste(observed_samples ,collapse = " ," )))
99- SAMPLES <- observed_samples
10099
101100# Read expected abundance table
102- if ( file.exists(expabundFILE ) ) {
103- exp = read.table( expabundFILE , header = TRUE , sep = " \t " , stringsAsFactors = FALSE , check.names = FALSE , strip.white = TRUE )
104- colnames(exp )[1 ] <- " ID"
105- # extract samples to analyse
106- exp_samples <- colnames(exp )[2 : ncol(exp )]
107- print(paste( " Expected samples:" , paste( exp_samples ,collapse = " ," )))
108- SAMPLES <- exp_samples [exp_samples %in% observed_samples ]
109- print(paste( " Investigate samples:" , paste( SAMPLES ,collapse = " ," )))
110- }
101+ exp = read.table( expabundFILE , header = TRUE , sep = " \t " , stringsAsFactors = FALSE , check.names = FALSE , strip.white = TRUE )
102+ colnames(exp )[1 ] <- " ID"
103+ # extract samples to analyse
104+ exp_samples <- colnames(exp )[2 : ncol(exp )]
105+ print(paste( " Expected samples:" , paste( exp_samples ,collapse = " ," )))
106+ SAMPLES <- exp_samples [exp_samples %in% observed_samples ]
107+ print(paste( " Investigate samples:" , paste( SAMPLES ,collapse = " ," )))
111108
112109# check if there are any samples to analyse
113110if (length(SAMPLES ) == 0 ) {
0 commit comments