Jump to content


Photo

Two-Stage Design and EMA’s fixed effects


  • Please log in to reply
2 replies to this topic

#1 Helmut Schütz

Helmut Schütz

    Advanced Member

  • Members
  • PipPipPip
  • 316 posts
  • LocationVienna, Austria

Posted 03 November 2012 - 05:51 PM

Hi,

a goody from the BE-wizard. I used to setup Potvin’s model of pooled data (stage 1&2) as:
Fixed: Sequence+Stage+Period(Stage)+Treatment
Random: Subject(Sequence*Stage)

Results agree with reported ones (Example 2) of the paper:
94.12% CI: 88.45–116.38%

Unfortunately EMA wants all effects fixed (which works in the conventioanl setup of stage 1) but in the pooled analysis I got:
“ERROR 11050: Parsing error: The containing term must be in model.”
The core output stops after the partial test; no LSMeans are calculated.

I could run the model in LME – but only if I didn’t ask for the LSMeans (otherwise I got the same error like from the BE wizard).
Results validated against R and SAS (http://forum.bebac.a...try.php?id=9462).

Do you think it will be possible to catch the error (I guess LME passes the results to BE)? Most users are not familiar with LME and are lost when they try to use the BE wizard.

Example project attached.
[file name=Potvin2.phxproj size=2567353]http://www.pharsight.com/extranet/media/kunena/attachments/legacy/files/Potvin2.phxproj[/file]

P.S.: Of course LSMeans are not necessary for the BE analysis. Some regulators will be surprised if the are missing anyhow.

Attached Files


Edited by Simon Davis, 08 January 2019 - 02:39 PM.

 Best regards,
Helmut

https://forum.bebac.at/


#2 Linda Hughes

Linda Hughes

    Advanced Member

  • Members
  • PipPipPip
  • 73 posts

Posted 05 November 2012 - 09:51 PM

Hi Helmut,

 

I'm just getting a chance to look at this. I changed the last BE model in your project to:

 

Sequence+Stage+Period(Stage)+Treatment+Sequence*Stage+Subject(Sequence*Stage)

 

so that the "containing term" is in the model, and the error message is no longer given.

 

The CI_User_Lower and CI_User_Upper in the BE output are now:

 

88.445196, 116.37701

 

The Ratio_%REf is:

 

101.45436

 

These are the same as the point estimate and back-transformed CI's in the "pooled fixed backtransform" object.

 

Does that resolve your problem?

 

Linda



#3 Helmut Schütz

Helmut Schütz

    Advanced Member

  • Members
  • PipPipPip
  • 316 posts
  • LocationVienna, Austria

Posted 06 November 2012 - 12:28 AM

Hi Linda!

Sequence+Stage+Period(Stage)+Treatment+Sequence*Stage+Subject(Sequence*Stage)
so that the "containing term" is in the model, and the error message is no longer given.

I see. You have to introduce "Sequence*Stage" in order to use the nested subject-term.

The CI_User_Lower and CI_User_Upper in the BE output are now […]
These are the same as the point estimate and back-transformed CI's in the "pooled fixed backtransform" object.

Fine. Agrees with SAS and R. ;-)

Does that resolve your problem?

Yes and no. I wouldn't have guessed this setup. The paper reads: "… model including sequence, stage, period(stage), treatment, subject(sequence × stage) …"

The SAS-code is: treatment period(stage) sequence stage subject(stage*sequence)
The R code is: 0+treatment+period%in%stage+sequence+stage+subjectj%in%(stage*sequence)

Maybe you post this code at the support site as well. Two-stage designs are acceptable in Canada, the US, and the EU. Russia is expected to follow soon.

THX for the solution; I'll post it at my forum.

 Best regards,
Helmut

https://forum.bebac.at/





0 user(s) are reading this topic

0 members, 0 guests, 0 anonymous users