Data <- na.omit(read.table(file="C:/Documents and Settings/steeleb/My Documents/Math 543/GlobalTB.txt",header=T)) name <- names(Data) head(Data,10) n <- dim(Data)[1] p <- dim(Data)[2]-2 Country <- Data$Country Region <- Data$Region regions <- levels(Region) Year <- 1982:2005 rates <- matrix(unlist(Data[,3:dim(Data)[2]]),nrow=n) head(rates) rates <- as.vector(t(rates)) head(rates) # check the ordering Year.2 <- rep(Year,times=n) Country.2 <- rep(Country,each=p) Region.2 <- rep(Region,each=p) data <- groupedData(rates~Year.2|Country.2,data=data.frame(rates,Year.2,Country.2),outer~Region.2) head(data) unique(getGroups(data))