How to do 'effects coding' in Matlab

2020 3 2
1 个回答
回答已采纳
更新时间：2023 7 12
7 次查看（30 天）

0 个投票

Hi, right now I am writing a program for 'design of experiments' and I need to calculate the effects matrix for my factors and levels. The matrix should look like this with white = 1, grey = 0, black = − 1.
The only functions I found regarding effects coding in Matlab are 'fitlme' and 'fitglme', but I really don't understand how they work and how I get an effects matrix from this (like in the picture). For example how would that work for a 2^9 design plan?
Thanks in advance.

请先登录，再回答此问题。

 采纳的回答

Tobias Averbeck
Tobias Averbeck 2020-3-3

0 个投票

So, I found a solution, you have to use categorical input to get the -1,0,1 effects matrix, this is the code I used:
% Define the levels of the factors and size of design
level = [4 4 6 6];
reduced_design_size = 300;
% Construct the reduced design
[D,~] = rowexch(size(level,2),reduced_design_size,...
'linear','categorical',[1:size(level,2)],'levels',level);
% This is just to have a numerical response column for the
% formula of the lme
D(:,5) = randperm(size(D,1))';
D_tbl = mat2dataset(D);
% Convert the factors to 'categorical' to get a -1,0,1 effects matrix
D_tbl.D1 = nominal(D_tbl.D1);
D_tbl.D2 = nominal(D_tbl.D2);
D_tbl.D3 = nominal(D_tbl.D3);
D_tbl.D4 = nominal(D_tbl.D4);
% Calculate LME with just addition of factors as formula
% and response as the added column
lme = fitlme(D_tbl,'D5 ~ D1 + D2 + D3 + D4',...
'FitMethod','REML','DummyVarCoding','effects');
% Extract effects matrix
DM = designMatrix(lme,'Fixed');
% sort for better overview
DM = sortrows(DM,[1:size(DM,2)],'descend');

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Matthew
Matthew 2023-7-12
I had a general question about using effects coding in LMEs. Say I had 3 categorical variables I wanted to run effects coding on - running fitlme only returns the intercept and 2 of those 3 variables. How do I view the last one?

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Mixed Effects 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by