# # cfkids-read.q # # ------------------------------------------------------------ # # PURPOSE: build the data # # AUTHOR: P. Heagerty # # DATE: 00/07/07 # # ------------------------------------------------------------ # data <- read.table( "NewCFkids.data", header=F ) # cfkids <- data.frame( id = data[,1], fev1 = data[,2], age = data[,3], female = as.integer( data[,4]==2 ), pseudoA = as.integer( data[,5]==3 ), f508 = 3-data[,6], pancreat = as.integer( data[,7]==2 ) ) # ##### Add the age-at-entry # source("cluster.min.q") # age0 <- cluster.min( cfkids$id, cfkids$age ) cfkids$age0 <- age0 # ##### Define age-since-entry # ageL <- cfkids$age - age0 cfkids$ageL <- ageL rm( age0, ageL ) # ##### Write data for SAS input # write( file="NewCFkids-SAS.data", t(cfkids), ncol=ncol(cfkids) ) # # end-of-file...