Rでは、NAではなく加算の値をどのように得ることができますか

Question

次の例では、NAではなく加算の値を求めています。どうやってやるの？

set.seed(7354) 
a <- rbinom(1,1,.5) 
x <- ifelse(a==0, rexp(1, 1/50), NA) 
y <- ifelse(a==1, rexp(1, 1/100), NA) 
b <- ifelse(a==1, rbinom(1,1,.5), NA) 
z1 <- ifelse(b==1, rexp(1, 1/100), NA) 
z2 <- ifelse(b==0, rexp(1, 1/190), NA) 
s <- (1-a)*x + a *(y + b * z1 + (1-b) * z2) 

希望する出力はs = 968.8501です。

Answer 1

この

s = c((1-a)*x, a*y, a*b*z1, a*(1-b)*z2) 
s = sum(s, na.rm = T) 

これは、あなたは、NASを持つようになるんが、そうではありません項の合計には影響しません。

あなたはシミュレーションやブートストラップループの一部としてこれをやっている場合は

Answer 2

set.seed(7354) 

a <- rbinom(1,1,.5) 
x <- (a==0) * rexp(1, 1/50) 
y <- (a==1) * rexp(1, 1/100) 
b <- (a==1) * rbinom(1,1,.5) 
z1 <- (b==1) * rexp(1, 1/100) 
z2 <- (b==0) * rexp(1, 1/190) 

(1-a)*x + a *(y + b * z1 + (1-b) * z2) 
# [1] 968.8501 

、これはずっとifelse Sの束よりも速くなります。