[R] adding variable into dataframe by indice

From: Adrian Katschke <adrian_at_atstatconsulting.com>
Date: Thu 09 Feb 2006 - 04:40:43 EST


  R-Helpers,    

  I am trying to insert a value into a dataframe. This value is a proportion calculated by counting the number of those individuals with that value and then inserting the proportion at the end of the dataframe to only those individuals with the given value. The problem I am running into is that the proportions are not being attached to only those individuals with the specified value for that proportion.    

  Below is an example of the code that I am using. The data is made up for the dataframe. Should give you an idea, but the original has 'NA' in many rows. The original data is what is reported in the output below.    

    #Read in Data
  age.int <- data.frame(IND_ID = seq(1, 140, 10), rs1042364 = sample( c("(1,1)","(1,2)","(2,2)"),14,replace = T),   first_drink = sample(5:17,14,replace = T))    

    asubs112 <- subset(age.int, rs1042364 != "(2,2)")    

    ages112 <- sort(unique(na.omit(asubs112$first_drink)))    

  for ( i in ages112) {
    indce <- which(na.omit(asubs112$first_drink == i))     prop <- length(indce)/nrow(asubs112)     asubs112[indce,4] <- prop
    asubs112[indce,]
  }    

  Below is the output that I get from the script above. Notice the proportion for the first NA but not any of the others. Not sure what I am doing wrong, any suggestions are a big help.    

  TIA,
  Adrian    

   asubs112[1:50,]

      IND_ID rs1042364 first_drink age_int          V5
4   10008007     (1,2)          NA      16 0.003891051
6   10013012     (1,2)          13      14 0.116731518
7   10015006     (1,2)          12      17 0.105058366
8   10015007     (1,1)          12      16 0.105058366
10  10021009     (1,2)          NA      15          NA
14  10039036     (1,2)          NA      15          NA
15  10039037     (1,2)          NA      13          NA
17  10045005     (1,2)          13      17 0.116731518
18  10045014     (1,2)          13      14 0.116731518
21  10055022     (1,2)          NA      15          NA






	[[alternative HTML version deleted]]

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Feb 09 06:08:54 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:26 EST