SAS Proc Mixed coding, replicate design [Regulatives / Guidelines]
Dear Matz,
the deeper the thoughts, the numerously the questions .
SUB(SEQ) is the random subject term in the syntax of PROC GLM in case of a 2x2 cross over.
It is included in the PROC MIXED code but by different syntax.
If you need comparable code you can go with:
But with that you loose the disparate variance terms (between and within) for the formulations under study.
The FA0(2) structure is chosen to ensure a positive definite covariance matrix, what ever this means to us as run-of-the-mill statisticians (in german: Feld-Wald-und-Wiesen-) .
Roughly spoken it ensures variance terms in the model that are >0, a very reasonable assumption for such terms I think .
Type=CSH (heterogenous compound symmetry) is mentioned in the FDA guidance also.
Both parameterizations should be theoretically the same.
Even type UN (unstructured) may be used, but suffer from sometimes negative estimates of variance components and is therefore not recommended in the FDA guidance.
This is too wide a subject and goes beyond the themes discussed here. Keywords are 'likelihood ratio test' and 'information criteria'. Ask the great mahatma Goooooogle.
Seems the personal view of the author(s?) of the FDA guidance.
The book
B Jones and MG Kenward
Design and Analysis of Cross-over Trials
Chapman & Hall/CRC, Boca Raton (2nd ed. 2003)
uses DDFM=KR (abbreviations, abbreviations, abbreviations ...) most time. Not so astonishing in looking at the second author .
So if you need deeper insight, I recommend this book and once again (see above) the papers by
Patterson and Jones
GSK BDS Technical Report 2002 – 01 (part I)
http://biometrics.com/wp-content/uploads/2009/06/tr2002-01.pdf
GSK BDS Technical Report 2002 – 05 (part II)
http://biometrics.com/wp-content/uploads/2002/12/tr2002-051.pdf
--
Edit: Updated URLs. [Helmut]
❝ One doubt ...
the deeper the thoughts, the numerously the questions .
❝ How can I include a SUB(SEQ) term in the above SAS code?.
SUB(SEQ) is the random subject term in the syntax of PROC GLM in case of a 2x2 cross over.
It is included in the PROC MIXED code but by different syntax.
If you need comparable code you can go with:
PROC MIXED data=yours;
class treatment period sequence subject;
model lnPKParam=treatment period sequence / DDFM=KR;
random subject(sequence) subject*treatment;
run;
But with that you loose the disparate variance terms (between and within) for the formulations under study.
❝ Why the covariance structure FA0(2) is used in this code?.Why not the other covariance structures?
The FA0(2) structure is chosen to ensure a positive definite covariance matrix, what ever this means to us as run-of-the-mill statisticians (in german: Feld-Wald-und-Wiesen-) .
Roughly spoken it ensures variance terms in the model that are >0, a very reasonable assumption for such terms I think .
Type=CSH (heterogenous compound symmetry) is mentioned in the FDA guidance also.
Both parameterizations should be theoretically the same.
Even type UN (unstructured) may be used, but suffer from sometimes negative estimates of variance components and is therefore not recommended in the FDA guidance.
❝ Is there any method to choose the covariance structure?
This is too wide a subject and goes beyond the themes discussed here. Keywords are 'likelihood ratio test' and 'information criteria'. Ask the great mahatma Goooooogle.
❝ Why DDFM=satterth ? Why not DDFM=kenwardroger?.
Seems the personal view of the author(s?) of the FDA guidance.
The book
B Jones and MG Kenward
Design and Analysis of Cross-over Trials
Chapman & Hall/CRC, Boca Raton (2nd ed. 2003)
uses DDFM=KR (abbreviations, abbreviations, abbreviations ...) most time. Not so astonishing in looking at the second author .
So if you need deeper insight, I recommend this book and once again (see above) the papers by
Patterson and Jones
GSK BDS Technical Report 2002 – 01 (part I)
http://biometrics.com/wp-content/uploads/2009/06/tr2002-01.pdf
GSK BDS Technical Report 2002 – 05 (part II)
http://biometrics.com/wp-content/uploads/2002/12/tr2002-051.pdf
--
Edit: Updated URLs. [Helmut]
—
Regards,
Detlew
Regards,
Detlew
Complete thread:
- random effect, fix effect in the guidance Appendix E LS 2008-02-05 18:53 [Regulatives / Guidelines]
- random effect, fix effect in the guidance Appendix E d_labes 2008-02-06 13:34
- random effect, fix effect in the guidance Appendix E LS 2008-02-06 17:58
- random effect, fix effect in the guidance Appendix E d_labes 2008-02-07 10:14
- random effect, fix effect in the guidance Appendix E SKR 2008-03-06 10:54
- Satterthwaite degrees of freedom Helmut 2008-03-06 13:57
- Satterthwaite degrees of freedom d_labes 2008-03-10 09:59
- Satterthwaite DF; SAS Helmut 2008-03-10 12:24
- Satterthwaite DF; SAS d_labes 2008-03-10 14:50
- Satterthwaite DF; SAS mathews 2008-03-11 07:12
- SAS Proc Mixed coding, replicate designd_labes 2008-03-11 09:43
- SAS Proc Mixed coding, replicate design Sathya 2008-09-17 07:21
- Replicate design evaluation d_labes 2008-09-17 10:49
- Replicate design evaluation Sathya 2008-09-18 06:43
- Replicate design evaluation d_labes 2008-09-17 10:49
- SAS Proc Mixed coding, replicate design Sathya 2008-09-17 07:21
- SAS Proc Mixed coding, replicate designd_labes 2008-03-11 09:43
- Satterthwaite DF; SAS mathews 2008-03-11 07:12
- Satterthwaite DF; SAS d_labes 2008-03-10 14:50
- Satterthwaite DF; SAS Helmut 2008-03-10 12:24
- Satterthwaite degrees of freedom d_labes 2008-03-10 09:59
- Satterthwaite degrees of freedom Helmut 2008-03-06 13:57
- random effect, fix effect in the guidance Appendix E SKR 2008-03-06 10:54
- random effect, fix effect in the guidance Appendix E d_labes 2008-02-07 10:14
- random effect, fix effect in the guidance Appendix E LS 2008-02-06 17:58
- random effect, fix effect in the guidance Appendix E d_labes 2008-02-06 13:34