Mercurial > repos > siwaa > redelac_stics_t

comparison redelacSticsTool.R @ 20:475aa4451308 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/redelac commit b4077990f5df769fbae45689f5bb12dbc81b2e4e"
author siwaa
date Fri, 10 Jul 2026 11:43:17 +0000
parents 8848db582840
children 5dcba002e019
comparison
equal deleted inserted replaced
19:8848db582840 20:475aa4451308
224 224
225 for (j in 1:length(casesUpcSol)) { 225 for (j in 1:length(casesUpcSol)) {
226 226
227 successive_usms_tab_split[[j]]$Treatment <- NULL 227 successive_usms_tab_split[[j]]$Treatment <- NULL
228 successive_usms_tab_split[[j]]$Number_USM <- NULL 228 successive_usms_tab_split[[j]]$Number_USM <- NULL
229
230 }
231
232 foreach(
233 j = 1:length(casesUpcSol),
234 .options.future = list(packages =c("foreach", "tidyr", "arrow"))
235 ) %dordopar% {
236
237 splited <- strsplit(casesUpcSol[j], split = "-")[[1]] 229 splited <- strsplit(casesUpcSol[j], split = "-")[[1]]
238 230
239 upc <- splited[1] 231 upc <- splited[1]
240 232
241 sol <- splited[2] 233 sol <- splited[2]
234
235 simPCPlist <- list()
242 236
243 for (i in 1:nrow(successive_usms_tab_split[[j]])) { 237 for (i in 1:nrow(successive_usms_tab_split[[j]])) {
244 #i=1 238 #i=1
245 if (i == 1) { 239 if (i == 1) {
246 simPCP <- result[[j]][[i]] 240 simPCP <- result[[j]][[i]]
247 } else { 241 } else {
248 simPCP$sim_list <- c(simPCP$sim_list, result[[j]][[i]]$sim_list) 242 simPCP$sim_list <- c(simPCP$sim_list, result[[j]][[i]]$sim_list)
249 } 243 }
250 244 }
251 } 245
246 simPCPlist[[j]] <- simPCP
247 }
248
249 rm(casesUpcSol)
250
251 foreach(
252 j = 1:length(simPCPlist),
253 .options.future = list(packages =c("foreach", "tidyr", "arrow"))
254 ) %dordopar% {
252 255
253 #agregation des sorties dans un dataframe 256 #agregation des sorties dans un dataframe
254 sim_df <- 257 sim_df <-
255 rbindlist(simPCP$sim_list, use.names = TRUE, idcol = "usm") 258 rbindlist(simPCPlist[[j]]$sim_list, use.names = TRUE, idcol = "usm")
256 259
257 #ajout des colonnes explicatives + de l'horizon temporel 260 #ajout des colonnes explicatives + de l'horizon temporel
258 sim_df <- sim_df %>% 261 sim_df <- sim_df %>%
259 separate( 262 separate(
260 col = usm, 263 col = usm,