Proc Mixed FA0(1) for partial replicate [RSABE / ABEL]
Hi Wolfgang,
If I recall correctly, and mind you I could be completely wrong, you can omit the consideration of correlated random effects and still get the right estimates of the s2w,r, s2b,r and compound sum for test.
Effectively, in my naïve notation where I would write q=rsb,r sb,t it would correspond to setting r=0; then the restricted max likelihood of the optimised V matrix changes, but the determination of s2w,r and s2b,r (and s2w,t+s2b,t, where applicable) do not.
This might still be in line with the random effects model, if I get Chow & Liu right. After all, the model does not rely on correlated random effects but justs allows to calculate correlations if we find it appropriate.
Not a statistician and not a SAS user so don't know any syntax that fits anything.
❝ With this notation, the remedy would formally be equivalent to setting q=1 - but this is more difficult to justify, and also I would not know how to specify it in SAS.
❝ So I guess we have to stick to FA0(1)
and the corresponding notation.
If I recall correctly, and mind you I could be completely wrong, you can omit the consideration of correlated random effects and still get the right estimates of the s2w,r, s2b,r and compound sum for test.
Effectively, in my naïve notation where I would write q=rsb,r sb,t it would correspond to setting r=0; then the restricted max likelihood of the optimised V matrix changes, but the determination of s2w,r and s2b,r (and s2w,t+s2b,t, where applicable) do not.
This might still be in line with the random effects model, if I get Chow & Liu right. After all, the model does not rely on correlated random effects but justs allows to calculate correlations if we find it appropriate.
Not a statistician and not a SAS user so don't know any syntax that fits anything.
—
Pass or fail!
ElMaestro
Pass or fail!
ElMaestro
Complete thread:
- RSAB model overspecified Wolfgang 2013-08-05 14:37 [RSABE / ABEL]
- RSAB model overspecified jag009 2013-08-05 15:35
- Proc Mixed FA0(1) for partial replicate d_labes 2013-08-05 15:59
- Proc Mixed FA0(1) for partial replicate Wolfgang 2013-08-05 16:24
- Proc Mixed FA0(1) for partial replicate ElMaestro 2013-08-05 17:19
- Proc Mixed FA0(1) for partial replicate Wolfgang 2013-08-06 08:35
- Proc Mixed FA0(1) for partial replicateElMaestro 2013-08-06 10:34
- Proc Mixed FA0(1) for partial replicate Wolfgang 2013-08-06 08:35
- Proc Mixed FA0(1) for partial replicate d_labes 2013-08-06 10:30
- FA0(1) for partial replicate in Phoenix Helmut 2013-08-07 15:13
- FA0(1) for partial replicate in SAS d_labes 2013-08-08 09:42
- FA0(1) for partial replicate in Phoenix lhughes 2013-08-30 15:17
- FA0(1) in PHX: works (& more questions) Helmut 2013-08-30 17:07
- FA0(1) for partial replicate in Phoenix Helmut 2013-08-07 15:13
- Proc Mixed FA0(1) for partial replicate ElMaestro 2013-08-05 17:19
- Proc Mixed FA0(1) for partial replicate Wolfgang 2013-08-05 16:24