diff --git a/MATLAB/beefblup.m b/MATLAB/beefblup.m
index 896cb07..08ff208 100644
--- a/MATLAB/beefblup.m
+++ b/MATLAB/beefblup.m
@@ -90,4 +90,18 @@ end
 disp(' ')
 disp('If no animal matching this description exists, the results may appear')
 disp('outlandish, but are still as correct as the accuracy suggests')
-disp(' ')
\ No newline at end of file
+disp(' ')
+
+% Form the fixed effect matrix
+X = zeros(numanimals, sum(numgroups)-length(numgroups)+1);
+X(:,1) = ones(1, numanimals);
+
+% Iterate through each group
+for i = 1:length(normal)
+    % Find the traits that are present in this trait
+    traits = uniquecell(data(:,i+5));
+    
+    % Remove the "normal" version from the analysis
+    normalindex = find(strcmp(traits, normal{i}));
+    traits(normalindex)  = [];
+end
\ No newline at end of file