# # hivnet-DataSetup.q # # ------------------------------------------------------------ # # PURPOSE: analysis of IC data # # AUTHOR: P. Heagerty # # DATE: 00/05/09 # # ------------------------------------------------------------ # data <- read.table( "HivnetWide.dat", header=F ) # vps.data <- data.frame( id = data[,1], risk.group = factor( data[,2], levels=1:4, labels=c("MSM","MaleIDU","WAHR","WAHR+IDU") ), education = factor( data[,3], levels=1:6, labels=c(" 0 ) # ################################################## # # ##### demographics / risk # vps.stacked$risk.group <- factor( rep(data[,2],rep(3,nsubjects)), levels=1:4, labels=c("MSM","MaleIDU","WAHR","WAHR+IDU") ) vps.stacked$education <- factor( rep(data[,3],rep(3,nsubjects)), levels=1:6, labels=c("