In ICTatRTI/PersonAlytics: Analytics for N-of-1 and small N intensive longitudinal designs, idiographic clinical trials (ICT), and interrupted time series

The Ovary data included in nlme can be used to illustrate analysis steps from an ideographic clinical trial (ICT).

The Ovary Data

see ?Ovary.

Modifications Needed

# change data to standard data.frame, and create a phase variable
Ovary <- as.data.frame(nlme::Ovary)
Ovary$Mare <- factor(Ovary$Mare, ordered = FALSE)
Ovary$Phase <- as.numeric(Ovary$Time > .5)
Ovary$Time <- sin(2*pi*Ovary$Time) # the model won't converge without this!
#par(mfrow=c(2,2))
#hist(Ovary$Time, xlim=c(-1,1.5))
#plot(Ovary$Time, Ovary$Time)
#hist(Ovary$Time, xlim=c(-1,1.5))
#par(mfrow=c(1,1))
# figure 5.10 page 240
#ggplot2::ggplot(Ovary, aes(x = Time, y = follicles)) + geom_line() + geom_point() +
#  facet_grid( ~ Mare)
## rescaled time
#ggplot2::ggplot(Ovary, aes(x = Time, y = follicles)) + geom_line() + geom_point() +
#  facet_grid( ~ Mare)
frm <- as.formula('follicles ~ Phase + Time + Phase*Time')

This vignette will be finished following the code in 'R:\PaCCT\Users\stueller\Pinheiro & Bates\r\Ovary_as_ICT.r'.

ICTatRTI/PersonAlytics documentation built on Dec. 13, 2024, 11:06 p.m.