% PURPOSE: An example of using sem() on a small data set (eire.dat)
%                                    
% Assumes that eire.mat is in the workspace

% y = Blood_Group_A data;
% x = pale indicator variable
% W = (sparse) Cliff-Ord weight matrix

n = length(y);
X = [ones(n,1),x]; %add intercept to data

%Construct "names" vector

names = strvcat('Blood_Group','Constant','Pale');
          
% Do exact sem regression

info.lflag = 0; % Full computation of log determinant

res = sem_ts(y,X,W,info);

%res = sem(y,X,W,info);

prt(res,names); %print results to screen