categorical covariate with 5 levels

5 messages 5 people Latest: Jul 12, 2010

categorical covariate with 5 levels

From: Hong Li Date: July 09, 2010 technical
Dear all, I want to build a covariate model for the pk of a drug which is metabolized in the liver. One of the covariates is liver GVHD which has 5 levels ( 0,1,2,3,4). The only way I can think of is to group liver GVHD into two categories, 0 ( not GVHD) and 1 ( having GVHD). But there is interest to evaluate the effect of degree of GVHD on pk parameter, e.g. CL. Does anyone have any suggestions on building a categorical covariate model of more than two levels. Thanks for any input. Best regards, Hong Li, Ph.D. Postdoctoral Associate Department of Pharmaceutical Science University at Buffalo, SUNY 363 Hochstetter Hall Buffalo, NY 14260 Tel#: (716)645-4827

RE: categorical covariate with 5 levels

From: Ishidaysm Date: July 09, 2010 technical
Dear Hong, I usually use below script for covariates with more than two (n) levels. Please create n-1 indicator variables. ;GVHD1 = 1 for GVHD=1, 0 otherwise GVHD1=0 IF (GVHD.EQ.1) GVHD1=1 ;GVHD2 = 1 for GVHD=2, 0 otherwise GVHD2=0 IF (GVHD.EQ.2) GVHD2=1 ;GVHD3 = 1 for GVHD=3, 0 otherwise GVHD3=0 IF (GVHD.EQ.3) GVHD3=1 ;GVHD4 = 1 for GVHD=4, 0 otherwise GVHD4=0 IF (GVHD.EQ.4) GVHD4=1 CL = THETA(A) + THETA(B)*GVHD1 + THETA(C)*GVHD2 + THETA(D)*GVHD3+ THETA(E)*GVHD4 (Additive effect of GVHD on CL) Or CL = THETA(A) * (1+THETA(B)*GVHD1) * (1+THETA(C)*GVHD2) * (1+THETA(D)*GVHD3) * (1+THETA(E)*GVHD4) (Proportional shift in CL for GVHD) Regards, Yoshi _/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/ Yoshimasa Ishida  Chugai Pharmaceutical CO., LTD.  Clinical PK/PD Analysis Group  Clinical Research Planning Dept. _/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/
Quoted reply history
-----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Hong Li Sent: Friday, July 09, 2010 1:01 PM To: [email protected] Subject: [NMusers] categorical covariate with 5 levels Dear all, I want to build a covariate model for the pk of a drug which is metabolized in the liver. One of the covariates is liver GVHD which has 5 levels ( 0,1,2,3,4). The only way I can think of is to group liver GVHD into two categories, 0 ( not GVHD) and 1 ( having GVHD). But there is interest to evaluate the effect of degree of GVHD on pk parameter, e.g. CL. Does anyone have any suggestions on building a categorical covariate model of more than two levels. Thanks for any input. Best regards, Hong Li, Ph.D. Postdoctoral Associate Department of Pharmaceutical Science University at Buffalo, SUNY 363 Hochstetter Hall Buffalo, NY 14260 Tel#: (716)645-4827

RE: categorical covariate with 5 levels

From: Sam Liao Date: July 09, 2010 technical
I have a four levels covariate model. Hope this is useful to you. Sam Liao Pharmax Research ========================================================== RA1=0 RA2=0 RA3=0 IF(RACE.EQ.7) RA1=1 IF(RACE.EQ.8) RA2=1 IF(RACE.GT.8) RA3=1 TVCL= TVCL *(1 + THETA(33)*RA1+ THETA(34)*RA2+ THETA(35)*RA3)
Quoted reply history
-----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Hong Li Sent: Thursday, July 08, 2010 9:01 PM To: [email protected] Subject: [NMusers] categorical covariate with 5 levels Dear all, I want to build a covariate model for the pk of a drug which is metabolized in the liver. One of the covariates is liver GVHD which has 5 levels ( 0,1,2,3,4). The only way I can think of is to group liver GVHD into two categories, 0 ( not GVHD) and 1 ( having GVHD). But there is interest to evaluate the effect of degree of GVHD on pk parameter, e.g. CL. Does anyone have any suggestions on building a categorical covariate model of more than two levels. Thanks for any input. Best regards, Hong Li, Ph.D. Postdoctoral Associate Department of Pharmaceutical Science University at Buffalo, SUNY 363 Hochstetter Hall Buffalo, NY 14260 Tel#: (716)645-4827

RE: categorical covariate with 5 levels

From: Yaming Hang Date: July 09, 2010 technical
Hi Hong, Other approaches I can think of are: 1.to plot the EBE of CL vs. different GVHD categories to determine whether you need a different parameter for each level or indeed you can (and how to) combine several groups together; 2.since your GVHD variable is an ordinal categorical variable instead of nominal groups like male and female, you may be able to introduce some functional form (i.e. treating the covariate as continuous) Assigning each category with a different parameter is a very intuitive approach, but the underlying assumption is that each group is different and no clear pattern exist for the relationship between CL and GVHD. This might be over-parameterized. Hope this helps, Yaming
Quoted reply history
-----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Sam Liao Sent: Friday, July 09, 2010 1:01 AM To: [email protected]; [email protected] Subject: RE: [NMusers] categorical covariate with 5 levels I have a four levels covariate model. Hope this is useful to you. Sam Liao Pharmax Research ========================================================== RA1=0 RA2=0 RA3=0 IF(RACE.EQ.7) RA1=1 IF(RACE.EQ.8) RA2=1 IF(RACE.GT.8) RA3=1 TVCL= TVCL *(1 + THETA(33)*RA1+ THETA(34)*RA2+ THETA(35)*RA3) -----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Hong Li Sent: Thursday, July 08, 2010 9:01 PM To: [email protected] Subject: [NMusers] categorical covariate with 5 levels Dear all, I want to build a covariate model for the pk of a drug which is metabolized in the liver. One of the covariates is liver GVHD which has 5 levels ( 0,1,2,3,4). The only way I can think of is to group liver GVHD into two categories, 0 ( not GVHD) and 1 ( having GVHD). But there is interest to evaluate the effect of degree of GVHD on pk parameter, e.g. CL. Does anyone have any suggestions on building a categorical covariate model of more than two levels. Thanks for any input. Best regards, Hong Li, Ph.D. Postdoctoral Associate Department of Pharmaceutical Science University at Buffalo, SUNY 363 Hochstetter Hall Buffalo, NY 14260 Tel#: (716)645-4827 Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.

RE: categorical covariate with 5 levels

From: Jeff Wald Date: July 12, 2010 technical
Another consideration is that GVHD is likely to have a rank-ordered effect on clearance (though I agree with the earlier suggestion that 5 levels may be over-parameterized). You could try making a cumulative sum of the effects on clearance extending on the samples from Sam and Yoshi. This is essentially a piecewise-linear, monotonic function. This is not going so far as having a continuous function, but it is intuitive and would probably simplify identifying which parameters you really need. Once you plot out this empirical function, a parsimonious, continuous relationship may present itself to you. Jeff
Quoted reply history
-----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Hang, Yaming Sent: Friday, July 09, 2010 10:00 AM To: Sam Liao; [email protected]; [email protected] Subject: RE: [NMusers] categorical covariate with 5 levels Hi Hong, Other approaches I can think of are: 1.to plot the EBE of CL vs. different GVHD categories to determine whether you need a different parameter for each level or indeed you can (and how to) combine several groups together; 2.since your GVHD variable is an ordinal categorical variable instead of nominal groups like male and female, you may be able to introduce some functional form (i.e. treating the covariate as continuous) Assigning each category with a different parameter is a very intuitive approach, but the underlying assumption is that each group is different and no clear pattern exist for the relationship between CL and GVHD. This might be over-parameterized. Hope this helps, Yaming -----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Sam Liao Sent: Friday, July 09, 2010 1:01 AM To: [email protected]; [email protected] Subject: RE: [NMusers] categorical covariate with 5 levels I have a four levels covariate model. Hope this is useful to you. Sam Liao Pharmax Research ========================================================== RA1=0 RA2=0 RA3=0 IF(RACE.EQ.7) RA1=1 IF(RACE.EQ.8) RA2=1 IF(RACE.GT.8) RA3=1 TVCL= TVCL *(1 + THETA(33)*RA1+ THETA(34)*RA2+ THETA(35)*RA3) -----Original Message----- From: [email protected] [mailto:[email protected]] On Behalf Of Hong Li Sent: Thursday, July 08, 2010 9:01 PM To: [email protected] Subject: [NMusers] categorical covariate with 5 levels Dear all, I want to build a covariate model for the pk of a drug which is metabolized in the liver. One of the covariates is liver GVHD which has 5 levels ( 0,1,2,3,4). The only way I can think of is to group liver GVHD into two categories, 0 ( not GVHD) and 1 ( having GVHD). But there is interest to evaluate the effect of degree of GVHD on pk parameter, e.g. CL. Does anyone have any suggestions on building a categorical covariate model of more than two levels. Thanks for any input. Best regards, Hong Li, Ph.D. Postdoctoral Associate Department of Pharmaceutical Science University at Buffalo, SUNY 363 Hochstetter Hall Buffalo, NY 14260 Tel#: (716)645-4827 Notice: This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station, New Jersey, USA 08889), and/or its affiliates Direct contact information for affiliates is available at http://www.merck.com/contact/contacts.html) that may be confidential, proprietary copyrighted and/or legally privileged. It is intended solely for the use of the individual or entity named on this message. If you are not the intended recipient, and have received this message in error, please notify us immediately by reply e-mail and then delete it from your system.