grand average and UnfoldMakie

[EmmanuelleKris]

Hello,
It is not clear to me how we can compute a grand average of the models across subjects (a model per subject) in order for the data to be in a format compatible with UnfoldMakie.
Does Unfold have a specific function for this?

[ReneSkukies]

Hey!
If you have the coeftable() or effects() of your models you can combine all the tidy DataFrames that you get to a big data frame with an additional :subject collumn. Then you can use the split-apply-combine functionality of DataFrames to calculate the grand average like this:

using DataFrames
using Statistics

mean_df = combine(groupby(tidy_df, [:time, :coefname]), :estimate => mean => :estimate)

There is also a bit of a longer tutorial in our UnfoldBIDS documentation

The tidy dataframe you get from this should be compatible with UnfoldMakie, but you might either have to specify your yhat column or rename it in your DataFrame. But @behinger or @vladdez could help with that.

[behinger]

if you prefer R-style tidier syntax (havent tested, but should work):

using Tidier
@chain tidy_df begin
  @group_by(time,coefname)
  @summarize(estimate = mean(estimate))
end

Estimate should be fine as a column name. @vladdez @ReneSkukies let's think together where such Tutorials / Functionality should "live" in, tutorial, extra functions, extra package UnfoldGroup, UnfoldStats? Include in UnfoldMakie directly? And link at helpful locations...

[ReneSkukies]

Ah, I meant that combine(groupby(tidy_df, [:time, :coefname]), :estimate => mean) will result in a column named :estimate_mean.

You can either change the outputcolumn by changing the command to
combine(groupby(..., :estimate => mean => :estimate),

or provide your new collumn to
plot_erp(... ; mapping = (y = :estimate_mean))

[EmmanuelleKris]

Thank you very much, I hadn't thought to look at Unfold Bids...

[vladdez]

@vladdez @ReneSkukies let's think together where such Tutorials / Functionality should "live" in, tutorial, extra functions, extra package UnfoldGroup, UnfoldStats? Include in UnfoldMakie directly? And link at helpful locations...

Can imagine an example in plot_erp page.

[ReneSkukies]

@vladdez @ReneSkukies let's think together where such Tutorials / Functionality should "live" in, tutorial, extra functions, extra package UnfoldGroup, UnfoldStats? Include in UnfoldMakie directly? And link at helpful locations...

Can imagine an example in plot_erp page.

For functionality, I think UnfoldStats would be a good place. We could then also provide more fancy versions than taking the mean/ median.