Mercurial > repos > siwaa > carboseq_s

diff carboseqUnsplit.R @ 8:5fdc9d0c223b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://forgemia.inra.fr/carboseq/record-projet-carboseq/-/tree/43e7f4ebfe4e552bf684936df0a050cb6def5aab/toos/CarboSeqSimulator commit 43e7f4ebfe4e552bf684936df0a050cb6def5aab-dirty"
author siwaa
date Fri, 30 Aug 2024 13:29:48 +0000
parents 3e7978af2ba2
children 2af4b1e46d90
line wrap: on
line diff
--- a/carboseqUnsplit.R	Thu Aug 29 14:58:19 2024 +0000
+++ b/carboseqUnsplit.R	Fri Aug 30 13:29:48 2024 +0000
@@ -1,13 +1,26 @@
 
-resU <- list(SOC=data.frame(), forcings=data.frame())
+resU <- data.frame()
+
+listFiles <- list.files('./inputs')
+
+# restoring the order first
+
+orderVec <- c()
 
 for (i in list.files('./inputs')) {
     load(file.path('.', 'inputs',i))
-    resU$SOC <- rbind(resU$SOC, res$SOC);
-    resU$forcings <- rbind(resU$forcings, res$forcings)
+    orderVec <- c(orderVec, res$ID[1]);
+    rm(res)
+    }
+
+listFilesReOrder <- listFiles[order(orderVec)]
+
+for (i in listFilesReOrder) {
+    load(file.path('.', 'inputs',i))
+    resU<- rbind(resU, res);
     rm(res)
     }
 
 res <-resU
 
-save(res, file = 'csopraRes.Rdata')
\ No newline at end of file
+save(res, file = 'csopraRes.Rdata')	
\ No newline at end of file