Displaying 1 result from an estimated 1 matches for "sdres".
Did you mean:
adres
2010 Apr 29
1
How to estimate the residual SD for each sample separately in mixed-effects model?
...the SDs using b<-VarCorr(c.lme). There:
- the 'attr(b$Assay, "stddev")' is the assay-to-assay SD component for
each sample (SDaa)
- the 'attr(b$Run, "stddev")' is the run-to-run component (SDrr)
- the 'attr(b, "sc")' i.e. the residual (SDres), would be the
within-run component, but it's a single number for all the samples.
* The problem:
- how to estimate the 'within-run' component (SDres) for each sample
separately, as the two other components?
* Solutions tried:
- subtracting SDaa and SDrr from total SD - sometime...