[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
Re: st: SP/RM[1,2] EMS
David Airey <firstname.lastname@example.org> wonders about what he thinks
is a difference between the example I show at
and a similar problem he is working on.
The crux of the question is:
> Ken's command is:
> . anova score noise / subject|noise period noise*period /
> period*subject|noise dial noise*dial / dial*subject|noise
> period*dial noise*period*dial , repeated(period dial)
> My data are a little differently named but basically the same:
> Ken has noise, I have "S"
> Ken has dial, I have "I"
> Ken has period, I have "P"
> Ken has subject, I have "A"
> For EMSs I get:
> A|S = piA + E
> S = piA + apiS + E
> P = asiP + iPA|S + E
> I = aspI + pIA|S + E
> SP = aiSP + iPA|S + E
> SI = apSI + pIA|S + E
> PI = asPI + E
> PA|S = iPA|S + E
> IA|S = pIA|S + E
> SPI = PIA|S + E
> PIA|S = PIA|S + E
> E = E
> So my F ratios to test are:
> MS(S) / MS(A|S)
> MS(P) / MS(PA|S) <-- different --> Ken does MS(P) / MS(E)
> MS(I) / MS(IA|S) <-- different --> Ken does MS(I) / MS(E)
> MS(SP) / MS(PA|S)
> MS(SI) / MS(IA|S)
> MS(PI) / MS(E)
> MS(SPI) / MS(E)
> My question is Ken tests each within subject factor over the error.
> What is different between testing over the error as Ken has done, vs
> testing over the terms I have determined above? It seems these terms
> are like the error as well.
David has misread the example in the FAQ. He says that it uses
MS(E) for the test of P and I (his letters, or period and dial
the corresponding names for the example in the FAQ).
Look carefully at the -anova- command and at the output shown in
term tested using
If I am following the mapping from his names to the FAQ example
names correctly, it appears that we are in agreement.
By the way, what I present in the FAQ as Stata output of -anova-
agrees with what Winer, Brown, & Michels show in Table 7.16 where
they show the ANOVA table for the example dataset.
Ken Higbee email@example.com
* For searches and help try: