[FieldTrip] Plantar pressure data

[Schoffelen, J.M. (Jan Mathijs)]

Dear Lucija,

The fact that you get a p-value of 1 throughout (in combination with the cfg that you have used for ft_timelockstatistics) indicates that either:
-there are no data points that survive the a priori threshold for clustering, or
-the observed clusters in the data (after thresholding the T-statistic at the nominal level) yield a cluster-based test statistic that is smaller than any of the test statistics obtained for the null distribution.

With respect to the error you get when computing the difference of the average across subjects: you don’t specify it, but I suspect that the number of subjects with and without ulcer is different. But on the off change that this is not the case, I recommend to compute the average across subjects first (per condition), and then use ft_math to do the subtraction. Of course, you’d need to check (and ensure) that the data matrices are compatible, i.e. with respect to the channels and (number of) time points.

Good luck and keep up the good work,
Jan-Mathijs


On 15 Jul 2023, at 15:18, Lucija via fieldtrip <fieldtrip at science.ru.nl<mailto:fieldtrip at science.ru.nl>> wrote:

Dear community,

I'm Lucija and currently I'm analysing the plantar pressure data collected with Pedar-X system. I'm using SPM-method to analyse the data. The data files are ascii. and fgt. files.  The data size is 90x40x number of subjects. For my research, I'm comparing two groups of different sizes with independent sample t-test.  I tried using function ft_timelockstatistics to perform statistical analysis.
This is the code I used:
cfg = [];
cfg.channel     = 'all';
cfg.neighbours = neighbours;
cfg.latency     = [1 norm_frames];
cfg.avgovertime = 'no';
cfg.parameter   = 'avg';
cfg.method      = 'montecarlo';
cfg.statistic   = 'ft_statfun_indepsamplesT';
cfg.alpha       = 0.05;
cfg.correctm    = 'cluster';
cfg.correcttail = 'prob';
cfg.numrandomization = 500;
cfg.labels = neighbours;

cfg.design(1,1:(Nsub_ulcer+Nsub_noulcer))  = [ones(1,Nsub_ulcer) 2*ones(1,Nsub_noulcer)];
cfg.ivar                = 1;

stat = ft_timelockstatistics(cfg, DataC1{:}, DataC2{:});

After running the statistics, I am getting p-values equal to 1 for many sensors. I understand there might be no significant differences because of the variation between the groups, but the fact that the probabilities are almost all 1 seems incorrect to me. Does anyone see the error in the code that might have caused this problem?

I have also tried plotting the data by calculating the average difference between the groups. I have used this code for plotting:
cfg = []
cfg.parameter = 'avg'
cfg.operation = 'subtract'

for i = 1:Nsub_ulcer
data{i} = ft_math(cfg, DataC1{i})
 end
for k = 1: Nsub_noulcer
 data{k} = ft_math(cfg, DataC2{k})
 end

cfg.gridscale = 90

cfg.projection  = 'orthographic'
  cfg.viewpoint = 'left'
cfg.style     = 'both';
cfg.layout    = 'ast_insole_xw';
cfg.highlight = 'on';
cfg.highlightchannel = find(stat.mask);
cfg.comment   = 'no';
cfg.colorbar           = 'yes'
figure; ft_topoplotER(cfg, data{:})

I am getting an error with the cfg.operation where arrays have incompatible sizes for this operation. Does anyone know the solution to this error?

Any help would be appreciated!

Kind regards,
Lucija Ninčević
Master student at Vrije Universiteit Amsterdam


_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
https://urldefense.com/v3/__https://doi.org/10.1371/journal.pcbi.1002202__;!!HJOPV4FYYWzcc1jazlU!6hZzyflII3grBizkj-kgyIhpW3NK5XjEfx7qLa_Mo9NuWdOIpJVMYLZy7Ap_DTvMlbLrm6ys5IGX7l8zvGKvLNYp504eZhfweOowLQ$

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20230807/6fc320af/attachment.htm>