0.8.1: added significant_clusters.txt and stack_names.txt

Author: slamballais

NEWS.md

@@ -4,7 +4,11 @@
 * If the estimated smoothness is below 1, we now increase it to 1 to avoid problems down the line.
 * Fixed a bug in qdecr_clusters where it assumes that there is always at least 1 cluster significant.
 * "w-g.pct" files can now be used as a measure by specifying "qdecr_w_g.pct" (underscore instead of hyphen).
-* Added `fwhm` argument to `qdecr_fastlm`.
+* Added `fwhm` argument to `qdecr_fastlm`, which was missing before.
+
+## New features
+* Added `cwp_thr` argument to `qdecr_fastlm` and `qdecr` to set the further cluster-wise p-value adjustment (default is 0.025 due to having 2 hemispheres, thus 0.05 / 2).
+* Automatically output two extra files: "significant_clusters.txt" (the output of `summary(vw, annot = TRUE)`) and "stack_names.txt" (the output of `stacks(vw)` and the corresponding stack numbers)
 
 # QDECR 0.8.0: Momo
 

R/qdecr_fastlm.R

@@ -112,6 +112,12 @@ vw <- qdecr(id = id,
 
 vw$describe$call <- rbind(vw$describe$call, c("call", "qdecr_fastlm call", paste(trimws(deparse(match.call())), collapse = "")))
 
+stacks_df <- data.frame(stack_number = seq_along(stacks(vw)), stack_name = stacks(vw))
+write.table(stacks_df, file.path(vw$paths$dir_out, "stack_names.txt"), quote = FALSE, row.names = FALSE, sep = "\t")
+
+summary_df <- summary(vw, annot = TRUE)
+write.table(summary_df, file.path(vw$paths$dir_out, "significant_clusters.txt"), quote = FALSE, row.names = FALSE, sep = "\t")
+
 if (save) qdecr_save(vw, save_data = save_data)
 vw
 }