% Read data y and x variables
% Read W matrix
% Determine N and T
% Create wx variables
for t=1:T
    t0=(t-1)*N;t1=t*N;
    wx(t0:t1,:)=W*x(t0:t1,:)
end
% Estimate GNS model by 
% Set info depending on the problem
info.model=3; % include spatial and time fixed effects
info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on
info.bc=1; % Bias correction
result = sac_panel_FE(y,[x wx],W,T,info); % change info if necessary
% vnames = ????
prt_sp(results,vnames,1); 
prt(result,vnames);  % print results
% Determine direct and indirect effects 
spat_model=1;
direct_indirect_effects_estimates(results,W,spat_model);