[R] I know this is wrong, but why?

Peter Dalgaard BSA p.dalgaard at biostat.ku.dk
Fri Aug 2 23:41:08 CEST 2002


Stuart Luppescu <s-luppescu at uchicago.edu> writes:

> My mind has become corrupted by having to use SAS too much. I wanted to
> calculate the difference of elements of two vectors if the signs are the
> same. I tried it the corrupted way first:
> 
> if (effects1[,3] < 0 & effects0[,3] < 0) {
>   imprv1 <- effects0[,3] - effects1[,3]
> } else  if  (effects1[,3] > 0 & effects0[,3] > 0) {
>   imprv1 <- effects1[,3] - effects0[,3]
> } else {
>   imprv1 <- NA
> 
> Then I realized it was giving me wrong results, (for one thing, I was
> getting no NAs) and I should be doing it ``The R Way'' like this:
> 
> imprv1 <- rep(NA, 267)
> imprv1[effects1[,3] < 0 & effects0[,3] < 0] <-  effects0[,3] - \
> effects1[,3]
> imprv1[effects1[,3] > 0 & effects0[,3] > 0] <-  effects1[,3] - \
> effects0[,3]
> 
> But what I don't understand is why the first method doesn't work, or at
> least why it didn't give me any warning or error messages. Can someone
> enlighten me?

The main thing is that R generally deals in entire vectors and that
if-statements will only make one decision, never do one thing for some
elements and something else for the rest. In  if (effects1[,3] < 0 &
effects0[,3] < 0)  the condition is a vector. but only the first
element is used, as in 

> X <- c(T,T,F)
> if (X) X else !X
[1]  TRUE  TRUE FALSE
> X <- c(F,T,F)
> if (X) X else !X
[1]  TRUE FALSE  TRUE

for by-element decisions consider ifelse(), or, as you did, indexing.

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list