% Este programa destina-se a resolver planos fatoriais 2^3
% As respostas tem que entrar pela seguinte ordem, rxyz, ryz...
alfa=0.05;
rxyz=9; ryz=10.5; rxy=12; ry=14.5; rxz=13; rz=16.5; rx=18; r1=22;
% alfa=input('Qual o alfa pretendido ? = ')
% rxyz=input('Qual a resposta rxyz ? = ')
% ryz=input('Qual a resposta ryz ? = ')
% rxy=input('Qual a resposta rxy ? = ')
% ry=input('Qual a resposta ry ? = ')
% rxz=input('Qual a resposta rxz ? = ')
% rz=input('Qual a resposta rz ? = ')
% rx=input('Qual a resposta rx ? = ')
% r1=input('Qual a resposta r1 ? = ')
respostas=[rxyz;ryz;rxy;ry;rxz;rz;rx;r1];
t=[1 -1 1 -1 1 -1 1 -1; 1 1 1 1 -1 -1 -1 -1; 1 1 -1 -1 1 1 -1 -1; 1 -1 1 -1 -1 1 -1 1; 1 -1 -1 1 1 -1 -1 1; 1 1 -1 -1 -1 -1 1 1; 1 -1 -1 1 -1 1 1 -1];
Efeitos=['Ex  ';'Ez  ';'Ey  '; 'Exy '; 'Exz ';'Eyz '; 'Exyz'];
Ef=t*respostas;
SQ=((Ef.^2)./8);
GL=[1;1;1;1;1;1;1];
QM=SQ./GL;
w=min(QM);
RV=QM/w;
f=finv(1-alfa,1,1)
ftab=[f; f; f; f; f; f; f];
signf=[];
for i=1:length(RV)
if RV(i)>=ftab(i)
signf=[signf;'***'];
else
signf=[signf;'   '];
end
end
format short G;MatEf=[Ef SQ GL QM RV ftab]
format
Final=[Efeitos signf]
display ('*** - efeitos significativos')