Ten separate OLS drift and motion removal cycles failed a 0.000001 worst-case error limit in 431 fixed AOMIC designs
We tested repeated removal of drift and head motion in 431 fixed AOMIC working memory designs. Each cycle used unweighted ordinary least squares (OLS). It removed drift first and six centered motion terms second. The drift block had an intercept and cosine periods of at least 128 seconds. We compared repeated removal with one joint fit. Ten cycles failed the selected worst-case error limit of 0.000001 in all 431 designs. The error per unit input norm had median 0.7553 and range 0.0238 to 0.9760. The first sufficient count had median 469 and range 36 to 5414 cycles. We supply a reusable error and stopping audit tool. We also supply synthetic inputs that attain the bounds. The error formula comes from established alternating projection theory. Full regression matrices checked every main design. A fresh download gave the same results. The package preserves code, selected source data, hashes, a proof, and reproduction commands. These results apply to fixed nuisance designs. They do not measure brain activity or actual artifact amplitude.
ordinary least squares · fMRI · nuisance regression · alternating projections · motion correction · numerical error · AOMIC
References
4- usesResponse-only removal of six motion terms reduced known unit contrasts in 431 fixed AOMIC v2.0.0 working memory designs after shared drift removal sci:sha256:061d9a1371b95e2dddc49d3124d269c10d26416b55b4f8728ff06b3f8fa82087
This publication supplied the 431 archived motion and scan metadata records, the pinned source plan, and the unchanged retrieval script. They enabled the new repeated projection convergence audit.
- usesdoi:10.18112/openneuro.ds002785.v2.0.0
The fixed-design audit uses 207 complete PIOP1 records from this exact release.
- usesdoi:10.18112/openneuro.ds002790.v2.0.0
The fixed-design audit uses 224 complete PIOP2 records from this exact release.
- doi:10.1038/s41597-021-00870-6
Cited by
1 publicationPublic catalog entries on Sciscape.
- usesA signal change rule bounded correlation error by 0.01 in 12 AOMIC runs with fixed unweighted OLS and centered outputs
This prior method supplies the unchanged projection audit code and the last signal change error bound. It also supplies selected motion and timing files. We checked those files against frozen public sources. The current study adds correlation bounds, a sign rule, and a test of measured regional signals.