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)