<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Cher Alain,
<div class=""><br class="">
</div>
<div class="">Thanks for reporting this! The SQUID jump indeed is quite prominent. In the absence of the ground truth, I quickly computed a volumetric image of kurtosis, using data from which the jump was cut out (prior to filtering). Interestingly, the spatial
image does not look altogether too different (apart from the scale) between the two situations (left panels in the attached screenshots always refer to the ‘original’ results, right panels to the results obtained with the jump cut out). These spatial images
were obtained with covariance matrices that obviously were quite different. I think that the beamformer to a large extent protects against the disruptive effect of the SQUID jump (at least spatially) because this jump does not localize well in brain space.
On the other hand, it had an effect on the reconstructed kurtosis value throughout the volume.</div>
<div class=""><br class="">
</div>
<div class="">However, I agree that the data that is subjected to any downstream analysis should be as clean as possible. I will make some small cosmetic changes to the tutorial to do justice to this.</div>
<div class=""><br class="">
</div>
<div class="">With respect to any ground truth data: </div>
<div class=""><br class="">
</div>
<div class="">In the context of the 2022 BIOMAG epilepsy data analysis competition I adjusted the kurtosis pipeline to be applied to the competition data. During the conference we discussed the possibility to make the data (+ prize winning analysis pipelines,
among which the kurtosis pipeline was one, whoohoo!) publicly available in some way. Despite the initial enthusiasm (of the competitors + organizers(: Jean-Michel Badier and Christian Bénar)) this ended up on everybody’s backburner unfortunately. Yet, I am
pretty sure that the competition data have been labeled. </div>
<div class=""><br class="">
</div>
<div class="">Best wishes,</div>
<div class="">Jan-Mathijs</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><img apple-inline="yes" id="8207AD58-CD16-4521-893F-F8898C7C595F" src="cid:47D17A13-8703-4B3C-86CA-D31893331724" class=""></div>
<div class=""><img apple-inline="yes" id="B73F389A-CE4E-43BF-98BA-085BF23F5C03" src="cid:AE32304A-8001-4A10-AA48-876F890F58D5" class=""></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 24 Mar 2024, at 19:00, Alain de Cheveigne via fieldtrip <<a href="mailto:fieldtrip@science.ru.nl" class="">fieldtrip@science.ru.nl</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">
<div class="BodyFragment"><font size="2" class=""><span style="font-size:11pt;" class="">
<div class="PlainText">Hi all, <br class="">
<br class="">
I've been going through the 'Virtual channel analysis of epilepsy MEG data' tutorial (<a href="https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.fieldtriptoolbox.org%2Ftutorial%2Fepilepsy%2F&data=05%7C02%7Cfieldtrip%40science.ru.nl%7C2f712284a8a341bec31c08dc4d843a59%7C084578d9400d4a5aa7c7e76ca47af400%7C1%7C0%7C638470477509485173%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=OlOcYeraMNfryIIj7tCBAVzcz3zeudR0bZX%2BJtC0YDo%3D&reserved=0" originalsrc="https://www.fieldtriptoolbox.org/tutorial/epilepsy/" shash="FnXaJEevi2MtZ0Vnh2id1OWGa2co/jn8II84C0c4Ks3pu6MVAbZUZ2PJVrRUv7zjhxacz0zNv2VFW8m/ulQZibjdtTAvWGFx3aEbycJjWlUvGND+95haBK41BHAHVcR8xicoGXacK/SHJRpUzRMFe3a29grK3AmGnkhvr21mc9E=" originalsrc="https://www.fieldtriptoolbox.org/tutorial/epilepsy/" shash="Y0tIPTlCvDRIMkE0pebPGkzJiU366BW/G3qsJVo3PomBeQbtHiJcyv14z1knEZtYP7OZMQgnFirzivYwQnpd+gBtQmxq0tHGitHoTPhSqHz/C40TrfHNKqs1MbSRmMT6RNzmH8duXZJTSdPJ3/4eBMrK7NkI4N+nMORk4By4+CI=" class="">https://www.fieldtriptoolbox.org/tutorial/epilepsy/</a>),
and I noticed that the CTF data for case3 has what looks like a flux jump, see below. The code applies a high-pass filter that turns this into a spike which might drive the locations of high kurtosis in the beamformer images.
<br class="">
</div>
</span></font></div>
<div class="BodyFragment"><font size="2" class=""><span style="font-size:11pt;" class="">
<div class="PlainText"><br class="">
Other examples seem OK, but it's hard to judge the outcome without ground truth. <br class="">
<br class="">
Is anyone aware of MEG/MRI data with labels that could be plugged into the tutorial? The ideal would be data labeled for the type of activity that a surgeon would really like to localize, e.g. interictal or pre-ictal onset.<br class="">
<br class="">
Alain<br class="">
<br class="">
</div>
</span></font></div>
<div class="BodyFragment"><font size="2" class=""><span style="font-size:11pt;" class="">
<div class="PlainText">_______________________________________________<br class="">
fieldtrip mailing list<br class="">
<a href="https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fmailman.science.ru.nl%2Fmailman%2Flistinfo%2Ffieldtrip&data=05%7C02%7Cfieldtrip%40science.ru.nl%7C2f712284a8a341bec31c08dc4d843a59%7C084578d9400d4a5aa7c7e76ca47af400%7C1%7C0%7C638470477509485173%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=jM9jd5Kta%2B3%2BOWZsG2Yj6g3AEgRb2Fh77W%2BZixVFhWc%3D&reserved=0" originalsrc="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" shash="Woi3TB6J7RwhEsJvwtRUmBPiLSzTksD6KmXIDw2/Oi+4dslNNIufLtOpo7nOvb/mqyU84SUUbCWZeLmp7kp7bYk/zxrTogLLaSBx/a4ooF2aLCRpTfpiFnezjSyuwPwKV7fnzAqBjPFZfef1eb24bSw4905bEC8YX56J4345BWQ=" originalsrc="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" shash="rIjt8kMBfPAMA5Q/m2TqgQhi/jfLr3/gZ60NlVGNcFKJ/xnXNunniBBN0H+V+I0cX78wihrcCIVpGJpjPU62mTch6MLA92j16J6MUSMNlsLAtv2RiJhMhfrjMsJ6HuWzUvd1A4WJp7ZSjBW/wrMLySmGviyZDuTCGTnhpSEf/cI=" class="">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br class="">
<a href="https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdoi.org%2F10.1371%2Fjournal.pcbi.1002202&data=05%7C02%7Cfieldtrip%40science.ru.nl%7C2f712284a8a341bec31c08dc4d843a59%7C084578d9400d4a5aa7c7e76ca47af400%7C1%7C0%7C638470477509641411%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C0%7C%7C%7C&sdata=FkUj3VulsFNiOg6R9XJ3dU62ilAlZLoxZCW5%2BzZvrNM%3D&reserved=0" originalsrc="https://doi.org/10.1371/journal.pcbi.1002202" shash="ttQZyYyQSkSDj8HDobFc4qyIZOUfOvTiOvh0vcCa63ncpFphrciCWJxgalDInno1kPYtSVwG+TIzgOk7TnRaU9WwX7xM7ffzF8oYdwoIE/v3hRE3JHB8nRDeVtnlFRWzB5MQmPQoI6jYwomUDigP5krz7dKfbIlLG0mISpBkCV8=" originalsrc="https://doi.org/10.1371/journal.pcbi.1002202" shash="FVXKHu9Yg/VDObJASVxNAb8r9x8Srqqopd9/c5+AjXem9ZkTloYwchWhoxXdNkdmizso9KzPq7ob2WNi0z93tIj1c388GCSNT0w93dzEyWXzlaJm3KG6TlAqp/ucBY0wVXIZCpY6gxT2ZHVRbfoMhi5OcSLQs1bcYAQ3uIena6k=" class="">https://doi.org/10.1371/journal.pcbi.1002202</a><br class="">
</div>
</span></font></div>
</div>
<span id="cid:d6983c96-0c5e-473c-9ce4-e00ef02ec867@EURPRD10.PROD.OUTLOOK.COM"><PastedGraphic-2.pdf></span></div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>