Jump to content


NLME myelosuppression model

  • Please log in to reply
4 replies to this topic

#1 Robert D. Faulkner

Robert D. Faulkner


  • Members
  • Pip
  • 2 posts

Posted 22 April 2014 - 03:57 PM

Working within NLME to implement a semi-mechanistic PK/PD model to describe the chemotherapy induced myelosuppression of various chemotherapeutics. Interested in using the PK/PD myelosuppression model described by Friberg et al. (J Clin Oncol, 20: 4713, 2002) which includes a proliferating compartment that is sensitive to the drug, 3 transit compartments reflecting cell maturation and a blood cell compartment. Want to know if it is possible to set up this model in NLME.


Thanks and kind regards


Bob [file name=Friberg_myelsupp_model_2002.pdf size=294601]http://www.pharsight.com/extranet/media/kunena/attachments/legacy/files/Friberg_myelsupp_model_2002.pdf[/file]

#2 Simon Davis

Simon Davis

    Advanced Member

  • Administrators
  • 1,111 posts

Posted 14 May 2014 - 09:03 AM

Bob, I'm sorry your question seems to have been overlooked - I am overdue to send out a digest.


  I don't see any immediate reason why this should not be possible in NLME, there are several examples with transit compartment models on the forum already that you should be able to find by typing 'transit' if that helps you get started.


e.g. http://www.pharsight...51&id=1911#1913


  Personally I would start with the graphical editor to describe the model on p4715 of the paper you attached and then edit it further in Textual mode. If you have some problems then post your project/code back here and we can try and take it from there.



#3 Samer Mouksassi

Samer Mouksassi

    Advanced Member

  • Members
  • PipPipPip
  • 90 posts

Posted 14 May 2014 - 03:38 PM

Dear Bob and Simon,

please find below working code I tested it in simulation mode with dummy one line input


id time CObs RATE ProlObs CircObs Aa1 ADDL II

1 0 100 13.88889 4.34 4.94 1000 7 168



it still needs work but it illustrate how e can code transit comparments with a feedback loop.

we can integrate the PK later on currently it has a biomarker linked to the ANC model from a recent publication: http://www.nature.co...psp201362a.html







Circ = BASE;


Transit1=KE*BASE/KP ;




deriv(BM = KIN*(1-EFF) - KOUT*BM )

deriv(Circ = -KE *Circ + KP*Transit3)

deriv(STEM = -KP *STEM + KP*STEM*(1-DRUG)*(BASE/Circ)**0.362)

deriv(Transit1= -KP *Transit1 + KP*STEM )

deriv(Transit2= -KP *Transit2 + KP*Transit1 )

deriv(Transit3= -KP *Transit3 + KP*Transit2 )


error(CEps2 = 0.001)

observe(CircObs = (Circ)+CEps2)




DRUG = 0.520*(-BMCHANGE)/(1.0099+(-BMCHANGE))

EFF = 1*(50/32.819)/(1.0099+(50/32.819))


stparm(KOUT = tvKout )

stparm(KIN = tvKout * BMBASE)

stparm(BMBASE =tvBM0 )

fixef(tvBM0 = c(,42554, ))

fixef(tvKout = c(,0.002719313, ))


FN = (BASE/Circ)**0.362




stparm(MTT = tvMTT * exp(nMTT))

stparm(kcirc = tvkcirc * exp(nkcirc))

stparm(BASE = tvCirc0* exp(nBASE) )



fixef(tvCirc0 = c(,4.94, ))

fixef(tvMTT = c(, 248, ))

fixef(tvkcirc = c(, 7, ))

ranef(diag( nMTT,nkcirc,nBASE) = c( 1,1,1))


#4 sk44



  • Members
  • Pip
  • 3 posts

Posted 28 January 2020 - 09:00 PM



Has anybody done Friberg  et al.2002 model in NLME Phoenix yet?





#5 bwendt@certara.com


    Advanced Member

  • Administrators
  • 188 posts

Posted 12 February 2020 - 02:54 PM

Hi SK,


not sure if this is the right one, but we have worked out a Friberg model as part of our PML School webinars. Please take a look here:




Let us know if there are specific questions.




0 user(s) are reading this topic

0 members, 0 guests, 0 anonymous users