[R] nesting multiple for loops

From: Ben Hunter <bjameshunter_at_gmail.com>
Date: Thu, 14 Apr 2011 16:34:22 -0700

Hi everyone.

I am quite frustrated that this doesn't work, as all the functions within work fine by themselves. I'd also like any pointers to how to avoid 'for' loops in my code. I understand it's less than desirable, but I'm still quite new and use them a lot.

I have a few wide datasets (90 to 120) with long column names, each name contains a number of different 'markers'. Each could be considered a factor variable within the column name. Their are two categories of factors, we'll call them f1 and f2.

The data frame names look something like this:  'ace_van' , 'boy_van', 'car_xes' , 'ace_xes', 'dog_wall' , 'car_zounds'

f1 <- c('ace', 'boy', 'car', 'dog')
f2 <- c('van', 'wall', 'xes', 'zounds') # actual vectors are length 6 and 7, so I don't want to individually sum the 42 combinations.

> var.table <- function(data, vec1, vec2)

  table <- as.data.frame(matrix(nrow = length(vec1), ncol = length(vec2)), row.names = vec1)
  names(table) <- vec2
  for (i in 1:length(vec1))
    for (j in 1:length(vec2))
      indices <- intersect(grep(vec1[j], names(data), value = TRUE), grep(cats[i], names(data), value = TRUE))

      table[i,j] <- sum(data[ ,indices])     }

> var.table(mydf, f1, f2)


Error in FUN(X[[1L]], ...) :
only defined on a data frame with all numeric variables

Every entry in mydf is an integer with no missing values.

Thanks a ton.


        [[alternative HTML version deleted]]

R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 14 Apr 2011 - 23:37:09 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 15 Apr 2011 - 10:30:31 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive