dt1=data.frame(ID=c(1,1,2,2),
OCC=c(1,2,1,2),
## ID level
eta1=c(1,1,3,3),
## occasion level
eta2=c(1,3,1,5),
## not used
eta3=0
)
## model level
findCovs(dt1)
## ID level
findCovs(dt1,"ID")
## acual ID level
findVars(findCovs(dt1,"ID"))
## occasion level
findCovs(findVars(dt1,"ID"),c("ID","OCC"))
## Based on a "real data example"
if (FALSE) {
dat <- NMscanData(system.file("examples/nonmem/xgxr001.lst", package = "NMdata"))
findCovs(dat,by="ID")
### Without an ID column we get non-varying columns
findCovs(dat)
}
Run the code above in your browser using DataLab