
miammt = [...
1  0.1334  0.2194;  ...
2  0.1214  0.1929;  ...
3  0.1276  0.1833;  ...
4  0.1152  0.1879;  ...
5  0.1047  0.2529; ...
6  0.1312  0.2527; ...
7  0.1222  0.2595; ...
8  0.2569  0.2726; ...
9  0.2101  0.2234; ...
10 0.1852  0.2216; ...
11 0.1798  0.1905; ...
12 0.1049  0.1444; ...
13 0.2649  0.2841; ...
14 0.2383  0.2731];

miatreat = miammt(1:7,2);
miacont  = miammt(1:7,3);
mmttreat = miammt(8:14,2);
mmtcont  = miammt(8:14,3);

miadif = miatreat - miacont;
mmtdif = mmttreat - mmtcont;

tmia = mean(miadif)/(std(miadif)/7)
tmmt = mean(mmtdif)/(std(mmtdif)/7)

pvalmia = tcdf( tmia, 6)
pvalmmt = tcdf( tmmt, 6)


t = mean(miadif - mmtdif)/sqrt( var(miadif)/7 + var(mmtdif)/7 )


