# # hivnet.q # # ------------------------------------------------------------ # # PURPOSE: analysis of IC data # # AUTHOR: P. Heagerty # # DATE: 00/05/09 revised 08March2002 # # ------------------------------------------------------------ # 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("