Mercurial > repos > greg > ideas
comparison ideas.R @ 152:26c26cb32137 draft
Uploaded
author | greg |
---|---|
date | Fri, 12 Jan 2018 13:12:47 -0500 |
parents | 3762c27d820a |
children | b140097a984e |
comparison
equal
deleted
inserted
replaced
151:9d34f7e6d80c | 152:26c26cb32137 |
---|---|
324 } | 324 } |
325 cmd = paste(cmd, "-rseed", opt$rseed, sep=" "); | 325 cmd = paste(cmd, "-rseed", opt$rseed, sep=" "); |
326 cmd = paste(cmd, "-thread", opt$thread, sep=" "); | 326 cmd = paste(cmd, "-thread", opt$thread, sep=" "); |
327 | 327 |
328 if (is.null(opt$training_iterations)) { | 328 if (is.null(opt$training_iterations)) { |
329 cmd = paste(cmd, "-o", output_base_name, sep=" "); | 329 final_cmd = paste(cmd, "-o", output_base_name, sep=" "); |
330 cmd = add_output_redirect(cmd, opt$save_ideas_log, opt$output_log, default_log_name); | 330 final_cmd = add_output_redirect(final_cmd, opt$save_ideas_log, opt$output_log, default_log_name); |
331 run_cmd(cmd, opt$save_ideas_log, opt$output_log, default_log_name); | 331 run_cmd(final_cmd, opt$save_ideas_log, opt$output_log, default_log_name); |
332 } else { | 332 } else { |
333 output_para0 = paste(output_base_name, ".para0", sep=""); | |
334 output_profile0 = paste(output_base_name, ".profile0", sep=""); | |
333 for (i in 1:opt$training_iterations) { | 335 for (i in 1:opt$training_iterations) { |
334 cmd = paste(cmd, "-o", paste(output_base_name, ".tmp.", i, sep=""), sep=" "); | 336 final_cmd = paste(cmd, "-o", paste(output_base_name, ".tmp.", i, sep=""), sep=" "); |
335 cmd = add_output_redirect(cmd, opt$save_ideas_log, opt$output_log, default_log_name); | 337 final_cmd = add_output_redirect(final_cmd, opt$save_ideas_log, opt$output_log, default_log_name); |
336 run_cmd(cmd, opt$save_ideas_log, opt$output_log, default_log_name); | 338 run_cmd(final_cmd, opt$save_ideas_log, opt$output_log, default_log_name); |
337 } | 339 } |
338 tpara = combine_state(paste(output_base_name, ".tmp.", (1:opt$training_iterations), ".para", sep=""), mycut=0.5); | 340 tpara = combine_state(paste(output_base_name, ".tmp.", (1:opt$training_iterations), ".para", sep=""), mycut=0.5); |
339 para = tpara$para; | 341 para = tpara$para; |
340 write.table(tpara$profile, paste(output_base_name, ".profile0", sep=""), quote=F, row.names=F, col.names=F); | 342 write.table(tpara$profile, output_profile0, quote=F, row.names=F, col.names=F); |
341 para = apply(para, 1, function(x){paste(x, collapse=" ")}); | 343 para = apply(para, 1, function(x){paste(x, collapse=" ")}); |
342 para = c(readLines(paste(output_base_name, ".tmp.1.para", sep=""), n=1), para); | 344 para = c(readLines(paste(output_base_name, ".tmp.1.para", sep=""), n=1), para); |
343 output_para0 = paste(output_base_name, ".para0", sep=""); | |
344 writeLines(para, output_para0); | 345 writeLines(para, output_para0); |
345 } | 346 cmd = c(cmd, "-otherpara", output_para0, output_profile0); |
347 if (length(which(cmd == "-G")) == 0) { | |
348 cmd = c(cmd, "-G", length(para)-1); | |
349 } else { | |
350 tt = which(cmd == "-G"); | |
351 cmd[tt + 1] = length(para)-1; | |
352 } | |
353 tt = which(cmd == '-C'); | |
354 if(length(tt)>0) { | |
355 cmd = cmd[-c(tt, tt+1)]; | |
356 } | |
357 } |