Thanks Peter,
I still wonder why it thinks it's unbalanced...
The se's of the contrasts are different than the se's of the means,
which is the point of se=T in model.tables (type "means") I would have
thought. No big deal though, the following code makes a nice table with
the se's of the means to use in further CI's be it t's or
Bonferroni:
kidney.mse <- kidney.anova[length(row.names(kidney.anova)),3]
means=se=means.a=se.a=means.b=se.b=rep(0,6)
kidney.means <- data.frame(means,se,means.a,se.a,means.b,se.b,a,b)
for (i in 1:6){
kidney.means[i,1] = mean(kidney$dhosp.trans[kidney$a.b==i])
if (i < 4){
if (i < 3){
kidney.means[i,3] = mean(kidney$dhosp.trans[kidney$duration==i])
kidney.means[i,4] =
kidney.mse/length(kidney$duration[kidney$duration==i])
}
kidney.means[i,5] = mean(kidney$dhosp.trans[kidney$wgain==i])
kidney.means[i,6] = kidney.mse/length(kidney$wgain[kidney$wgain==i])
kidney.means[i,7] = 1
kidney.means[i,8] = i
}
if (i > 3){
kidney.means[i,7] = 2
kidney.means[i,8] = i-3
}
kidney.means[i,2] = kidney.mse/length(kidney$dhosp.trans[kidney$a.b==i])
}
> kidney.means
means se means.a se.a means.b se.b a b
1 0.4434824 0.01012537 0.7867166 0.003375125 0.4208548 0.005062687 1 1
2 0.8099675 0.01012537 0.6151953 0.003375125 0.6954658 0.005062687 1 2
3 1.1066998 0.01012537 0.0000000 0.000000000 0.9865472 0.005062687 1 3
4 0.3982271 0.01012537 0.0000000 0.000000000 0.0000000 0.000000000 2 1
5 0.5809641 0.01012537 0.0000000 0.000000000 0.0000000 0.000000000 2 2
6 0.8663947 0.01012537 0.0000000 0.000000000 0.0000000 0.000000000 2 3
Thanks again
Peter Alspach wrote:
>Damián
>
>I asked a similar question a few months ago (3 August 2004):
>
>
>
>>temp.aov <- aov(S~rep+trt1*trt2*trt3, data=dummy.data)
>>model.tables(temp.aov, type='mean', se=T)
>>
>>Returns the means, but states "Design is unbalanced - use
se.contrasts
>>for se's" which is a little surprising since the design is
balanced.
>>
>>
>
>To which Prof Ripley replied: If you used the default treatment contrasts,
it is not. Try Helmert
>contrasts with aov().
>
>If I recall correctly, following Prof Ripley's suggestion led aov() to
accept the design was balanced, but model.tables() still did not (but that could
have been my error). However, se.contrast() worked.
>
>Cheers ........
>
>Peter Alspach
>
>
>
>
>
>>>>Damián Cirelli <damian.cirelli@maine.edu> 02/12/04
09:50:13 >>>
>>>>
>>>>
>Hi all,
>I'm new to R and have the following problem:
>I have a 2 factor design (a has 2 levels, b has 3 levels). I have an
>object kidney.aov which is an aov(y ~ a*b), and when I ask for
>model.tables(kidney.avo, se=T) I get the following message along with
>the table of effects:
>
>Design is unbalanced - use se.contrast() for se's
>
>but the design is NOT unbalanced... each fator level combination has the
>same n
>
>I' d appreciate any help.
>Thanks.
>
>______________________________________________
>R-help@stat.math.ethz.ch mailing list
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
>
>
>______________________________________________________
>
>The contents of this e-mail are privileged and/or confidential to the
>named recipient and are not to be used by any other person and/or
>organisation. If you have received this e-mail in error, please notify
>the sender and delete all material pertaining to this e-mail.
>______________________________________________________
>
>
>
[[alternative HTML version deleted]]