Hi David,

Thanks. Yes this partitions the error variance, but not the random effect
variance. That part is working fine (and in fact matches what I get from
running both subsets). But still curious about the random effects being
different.

Kris


On Tue, Jul 30, 2013 at 8:55 PM, David Duffy <David.Duffy@qimr.edu.au>wrote:

> On Wed, 31 Jul 2013, Kristofor Voss wrote:
>
>  So when I run the model I first ran (i.e. the one you said is correct),
>>
>
> I don't use lme very often, but I believe that if you use
>
> weights = varIdent(form=~1|factor))
>
> the output should include something like:
>
> Variance function:
> Structure: Different standard deviations per stratum
> Formula: ̃ 1 | Factor
> Parameter estimates:
> Restored Non-restored
> 1 0.42537




-- 
Kristofor A. Voss
Ph.D. Candidate
Bernhardt Lab
Nicholas School of the Environment
Duke University
P.O. Box 90328
Durham, NC 27708
619-249-3614
kristofor.voss@duke.edu
Office: FFSC 3304

	[[alternative HTML version deleted]]

