[R-sig-ME] unreplicated repeated measures
Christine Griffiths
Christine.Griffiths at bristol.ac.uk
Tue Jun 23 19:29:42 CEST 2009
Dear R users,
Sorry if this question is not applicable to this site. I am having problems
analysing unreplicated repeated measures. I calculated food web properties
for three treatments over time (10 months). It is unreplicated in that I
only have one observation per month per treatment. My problem is that I am
interested in how a food web property varies between Treatments and over
time.
Originally I had tried using lmer:
m4<-lmer(lncon~Treatment*month+(1|month),data=dataset)
but this provides the following error for which I have not found an
explanation to on the R site.
Error in mer_finalize(ans) : Calculated PWRSS for a LMM is negative
I suspect it is because I had treated month as a factor and consequently
due to the lack of replication and trying to interact these two categorical
variables it fails. I can overcome this problem by treating month as a
continuous variable from which I calculated confidence intervals using MCMC
method. However I am not sure how these are being calculated and if this is
accurate, given I have no mean per Treatment as such. Is it acceptable to
use month as a continuous variable?
Alternatively, I tried using repeated measures ANOVA, aov, to model the
data. I am cautious to use this method as it indicates significant
differences which are not apparent from the plotted raw data. Given that my
data lacks replication, I am wary of this method.
I have investigated time series analysis, but I am reluctant to venture
down this route.
Any reassurance or advice as to the best technique given my data would be
greatly appreciated.
Many thanks,
Christine
----------------------
Christine
More information about the R-sig-mixed-models
mailing list