ifelse无法正常运行中的R

Question

我试图使用命令ifelse如下

temp <- function(p){ 
ifelse(p < 0.5 *(1 + 0.5), (1 + 0.5) * qnorm(p/(1 +0.5)), (1 - 0.5) * qnorm((p - 0.5)/(1 - 0.5))) 
} 

功能与标量值效果很好，例如

> temp(0.1) 
[1] -2.251629 
> temp(0.9) 
[1] 0.4208106 

创建R中的矢量功能但不适用于矢量：

> temp(c(0.1,0.9)) 
[1] -2.2516289 0.4208106 
Warning message: 
In qnorm((p - 0.5)/(1 - 0.5)) : NaNs produced 

奇怪的是它返回正确的答案，b表示警告。

我在做什么错？看起来ifelse正在评估向量p的所有条目中的两个函数，这应该通过此命令来避免。

Answer 1

ifelse基本上是这样：

p<- c(.1,.9) 
a<-(1 + 0.5) * qnorm(p/(1 +0.5)) 
b<- (1 - 0.5) * qnorm((p - 0.5)/(1 - 0.5)) 
c<-NULL 
c[which(p < 0.5 *(1 + 0.5))] <-a[which(p < 0.5 *(1 + 0.5))] 
c[which(!(p < 0.5 *(1 + 0.5)))] <-b[which(!(p < 0.5 *(1 + 0.5)))] 

也就是说，它创造了“是”一个载体和“不”的向量。它创建的'否'矢量会引发警告。

本文档中的例子暗示了这一点。

x <- c(6:-4) 
sqrt(x) #- gives warning 
sqrt(ifelse(x >= 0, x, NA)) # no warning 

## Note: the following also gives the warning ! 
ifelse(x >= 0, sqrt(x), NA)