Motor evoked potentials (MEPs) caused by transcranial magnetic stimulation (TMS) provide a possibility of noninvasively mapping cortical muscle representations for clinical and research purposes. The interpretation of such results is complicated by the high variability in MEPs and the lack of a standard optimal mapping protocol. Comparing protocols requires the determination of the accuracy of estimated representation parameters (such as the area), which is problematic without ground truth data. We addressed this problem and obtained two main results: (1) the development of a bootstrapping-based approach for estimating the within-session variability and bias of representation parameters and (2) estimations of the area and amplitude-weighted area accuracies for motor representations using this approach. The method consists in the simulation of TMS mapping results by subsampling MEPs from a single map with a large number of stimuli. We studied the extensor digitorum communis (EDC) and flexor digitorum superficialis (FDS) muscle maps of 15 healthy subjects processed using Voronoi diagrams. We calculated the (decreasing) dependency of the errors in the area and weighted area on the number of stimuli. This result can be used to choose a number of stimuli sufficient for studying the effects of a given size (e.g., the protocol with 150 stimuli leads to relative errors of 7% for the area and 11% for the weighted area in 90% of the maps). The approach is applicable to other parameters (e.g., the center of gravity) and other map processing methods, such as spline interpolation.
Keywords: Cortical muscle representation mapping; Motor cortex mapping; Motor evoked potential; Statistical bootstrapping; TMS; Voronoi diagram.