Ben,
Thanks for the suggestions.

> It's a little dodgy to use year (with only 3 levels) as a random
> effect -- you might try with a fixed effect and see whether you get
> different results.

The reason I would like to use year as a random effect is because there are
large differences between years (this is a post-fire study on vegetation).
Does this justify setting year as a random effect?

> I don't understand why you're using aovlmer.fnc with                     >
which=c("TreatC","TreatUC") -- it looks to me as though this
> function is intended to be used to test hypotheses involving         >
multiple parameters simultaneously equalling zero.  Testing both > levels of
the treatment in this way doesn't really make sense (to   > me).
Furthermore, I think you'd need "Treat-1" rather than "Treat" > as your
fixed effect to do this.

Thanks for pointing out that I don't really need to use the aovlmer.fnc
function. After reading through the languageR document more carefully in
light of your comments it's clear that you're right, I really only need the
pvals.fnc function.

> If you just want to test for differences between the two treatment
> levels, what's wrong with the results of 'mcmc0' below?

I tried the code that you included in your response and it came up with the
same error that I've been getting all along:

Error in pvals.fnc(lmer.seeded, nsim = 10000) :
  trying to get slot "coefs" from an object (class "summaryDefault")
that is not an S4 object

Were you able to get results using this code? Any ideas on why I'm getting
this error when I try the pvals.fnc function?

-- 
Daniel Zvirzdin
Plant and Wildlife Sciences Department
Brigham Young University
(801) 422-2211 (w)
(801) 472-6006 (c)

	[[alternative HTML version deleted]]


