Bogus PROC MIXED/FDA [Software]
Grosser Meister!
Thats the question!
Just to cite myself: "I have the strong believe ...". See above.
Without really understanding the AlGore rhythms for REML as a "statistical amateur" (stolen from Helmut
) the overall variance matrix, which is part of the likelihood, which is maximised or equivalently -log likelihood is minimised ...
is
where the intra-individual variance parameter go into R matrix and the inter-individual go into G. There is no possibility within Proc MIXED, and I suppose also in other real mixed model software, to parameterize a interdependence between both. And such an interdependence we have.
Just a short outline (according to McNally's paper, citation above):
From the values of Test -> total variability of T = var(between-T)+var(within-T)
From the subject mean values of Ref -> total variability of R = var(between-R)+var(within-R)/2
From R-R -> within variability of R = var(within-R)
From T-R -> var(diff) = var(SxF) + var(within-T) + 0.5*var(within-R)
red: not separable between and within terms because of
var(SxF)=var(between-T)+var(between-R)-2*cov(between-T,between-R)
Thats the basis of my strong believe.
❝ But how ?
Thats the question!
Just to cite myself: "I have the strong believe ...". See above.
Without really understanding the AlGore rhythms for REML as a "statistical amateur" (stolen from Helmut


V= Z*G*Zt + R
where the intra-individual variance parameter go into R matrix and the inter-individual go into G. There is no possibility within Proc MIXED, and I suppose also in other real mixed model software, to parameterize a interdependence between both. And such an interdependence we have.
Just a short outline (according to McNally's paper, citation above):
From the values of Test -> total variability of T = var(between-T)+var(within-T)
From the subject mean values of Ref -> total variability of R = var(between-R)+var(within-R)/2
From R-R -> within variability of R = var(within-R)
From T-R -> var(diff) = var(SxF) + var(within-T) + 0.5*var(within-R)
red: not separable between and within terms because of
var(SxF)=var(between-T)+var(between-R)-2*cov(between-T,between-R)
Thats the basis of my strong believe.
—
Regards,
Detlew
Regards,
Detlew
Complete thread:
- Convergence problem in PROC MIXED dixit 2010-03-21 10:19
- Convergence problem in PROC MIXED ElMaestro 2010-03-21 11:35
- Convergence problem in PROC MIXED dixit 2010-03-22 06:26
- Convergence problem in PROC MIXED ElMaestro 2010-03-22 08:21
- Convergence problem in PROC MIXED dixit 2010-03-23 19:15
- Convergence problem in PROC MIXED ElMaestro 2010-03-23 20:29
- Convergence problem in PROC MIXED dixit 2010-03-24 15:09
- Convergence problem in PROC MIXED ElMaestro 2010-03-25 07:44
- Bogus PROC MIXED/FDA d_labes 2010-03-25 09:15
- Bogus PROC MIXED/FDA ElMaestro 2010-03-25 10:42
- Bogus PROC MIXED/FDAd_labes 2010-03-26 09:45
- Bogus PROC MIXED/FDA ElMaestro 2010-03-25 10:42
- Convergence problem in PROC MIXED dixit 2010-03-24 15:09
- Convergence problem in PROC MIXED shri 2010-04-05 08:01
- Convergence problem in PROC MIXED ElMaestro 2010-03-23 20:29
- Convergence problem in PROC MIXED dixit 2010-03-23 19:15
- partial replicate? d_labes 2010-03-22 09:50
- Convergence problem in PROC MIXED ElMaestro 2010-03-22 08:21
- Convergence problem in PROC MIXED dixit 2010-03-22 06:26
- Convergence problem in PROC MIXED ElMaestro 2010-03-21 11:35