BookmarkSubscribeRSS Feed
deleted_user
Not applicable
Based on the limited information I have provided, Mr. D (a dedicated contributor at this forum) has kindly provided me a base code for my ‘drug test’ scenario. I trust his base code, while I am not certain of the modifications I have made to it. Please take a look and see if the code below looks alright. Thanks in advance.

Scenario:
A specific drug X is administered to patients at various dosages (DOSE), and the change in patients’ heart rates is monitored hourly (PATIENT, HEARTRATE, and TIME). Men and women are tested separately (GENDER). Also, prior to assigning patients to a certain dosage, three weight categories (BLOCK) are created.


proc mixed data=mydata;
class dose patient time gender block;
model heartrate= dose|time|gender;
random block;
repeated time /subject=patient(dose*gender*block) type=un;
run;


Any suggestions/comments are greatly appreciated. Especially, I am not sure how to handle the block effect in the proc mixed repeated procedure.
4 REPLIES 4
SteveDenham
Jade | Level 19
The statements you have will give marginal means averaging over blocks. This seems logical to me, as the realization of blocks in this particular study is a sample from all possible blocks.

I think your measurements are at equal intervals, and if there are more than three or four, you may wish to use a time-dependent covariate structure, rather than the unstructured matrix, simply because the number of parameters to be estimated is smaller. You could check the AIC for this.

Finally, you may want BLUPs for heartrate for the three blocks. These differe from the marginal means (lsmeans), but can be obtained from estimate statements. See the documentation for examples.

Good luck.

Steve Denham
deleted_user
Not applicable
Thanks a lot Steve!
deleted_user
Not applicable
Hi again. I am trying to modify the above ‘drug test’ scenario to accommodate “unequal” measurement intervals. Using the SAS book as a reference, I made the following modifications:

- Base line heartrate measurement at T=0 is made into a separate entry as BASEHEARTRATE
- New entry TIME1 is created which is identical to TIME to be used in the REPEATED statement
- Chanes indicated in the upper letters

proc mixed data=mydata ORDER=DATA;
class dose patient time gender block;
model heartrate= dose|time|gender BASEHEARTRATE;
random block;
repeated time /subject=patient(dose*gender*block) type= SP(POW)(TIME1);
run;

This seems to run o.k., but is this right? Also, I don’t understand:
- why you need ORDER=DATA?
- why you need BASEHEARTRATE for “unequal” intervals but not for the “equal” intervals?
(==> Without BASEHEARTRATE in the MODEL statement, does it not account for the pre-treatment differences (i.e., differences in the individual heartrate at T=0?))
deleted_user
Not applicable
Correction: - Changes are indicated in uppercase letters.

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 4 replies
  • 914 views
  • 0 likes
  • 2 in conversation