Hello R-Help users,
 
I have a data frame named fd, a sample of which looks like this:
 




cbi_A

cbi_B

cbi_B1

cbi_B2

cbi_C

cbi_D

cbi_E

cbi_F


2.183451

1.047546

NA

NA

NA

NA

0.428528

NA


0.795837

0.510152

0.510152

NA

NA

NA

NA

NA


0.795837

1.149577

0.843485

1.122334

NA

NA

NA

NA


1.885522

2.203959

NA

NA

3.020202

3.040506

0.428528

0.224467


2.877257

1.353637

NA

NA

3.020202

NA

0.836649

NA


1.479441

2.816141

NA

NA

3.020202

NA

3.040506

3.040506


1.836547

1.659729

1.203959

2.401184

3.020202

NA

NA

NA


2.069177

1.870625

1.768595

2.989593

NA

NA

NA

NA


2.046985

1.203959

1.203959

NA

NA

NA

NA

NA


1.469238

NA

1.278849

NA

NA

NA

NA

NA









I want to create a new column (fd$cbi_tot) that is an average of some of the 
other columns (chosen based on the conditions as indicated in the code below).
 
ifelse(is.na(fd$cbi_C) & is.na(fd$cbi_D) & is.na(fd$cbi_E) & is.na(fd$cbi_F) & 
is.na(fd$cbi_B1) & !is.na(fd$cbi_B2), 
cbi_totlist <- cbind(fd$cbi_A, fd$cbi_B2),
ifelse(is.na(fd$cbi_C) & is.na(fd$cbi_D) & is.na(fd$cbi_E) & is.na(fd$cbi_F) & 
is.na(fd$cbi_B2) & !is.na(fd$cbi_B1),
cbi_totlist <- cbind(fd$cbi_A, fd$cbi_B1),
ifelse(is.na(fd$cbi_C) & is.na(fd$cbi_D) & is.na(fd$cbi_E) & is.na(fd$cbi_F) & 
!is.na(fd$cbi_B1) & !is.na(fd$cbi_B2), 
cbi_totlist <- cbind(fd$cbi_A, fd$cbi_B1, fd$cbi_B2),
cbi_totlist <- cbind(fd$cbi_A, fd$cbi_B, fd$cbi_C, fd$cbi_D, fd$cbi_E, 
fd$cbi_F))))
 
fd$cbi_tot <- apply(cbi_totlist, 1, mean, na.rm = TRUE)
 
I do not get an error and the new column (fd$cbi_tot) is created, but when I 
check the results, the first condition works properly, as does the last (the 
catch-all), but the middle two conditions seem to default to the catch-all. For 
example, in my data, look at row 3, where fd$cbi_C, fd$cbi_D, fd$cbi_E and 
fd$cbi_F  are all NA, and fd$cbi_B1 and fd$cbi_B2 have values. In this case, R 
uses the cbi_totlist that includes fd$cbi_B, whereas I would like it to use the 
list that includes fd$cbi_B1 and fd$cbi_B2 (and I thought I indicated that in 
the conditional statements).   Anybody have ideas about what I am doing wrong 
here? 
 
Thanks much for any suggestions!
 
                                          
        [[alternative HTML version deleted]]

______________________________________________
R-help@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.

Reply via email to