用在R如果功能

我正在運行中的R短if功能，但我得到以下警告消息：用在R如果功能

In if ((runif(50, 0, 1) < 0.69)) { : the condition has length > 1 and only the first element will be used`.

我的R初步的瞭解使我相信，runif生成向量但if產量一個單一的價值，所以我認爲這是問題。在這裏有沒有簡單的替代？

另外我想要最終產品是兩個參數的矩陣組合，但我不確定是否將兩個方案中的rnorm函數放置爲50是正確的。

Test <- 
if((runif(50, 0, 1)<0.69)) { 
rnorm(50, 25, 4) 
} else { 
rnorm(50, 28, 4.3) 
}

來源

2012-08-01 YesSure

你想在最後得到什麼？一個50件物品的矢量？ – Pop 2012-08-01 10:59:29

if()預計whetever在括號內以評估長度爲1的邏輯向量（TRUE或FALSE）。如果矢量長度超過1，則警告說只使用矢量的第一個元素。

如果ifelse()它作爲一種替代方案：

ifelse(test, TRUE_CASE, FALSE_CASE)

爲了您的例子中，我們有：

set.seed(1) 
ifelse(runif(50, 0, 1) < 0.69, rnorm(50, 25, 4), rnorm(50, 28, 4.3))

或

set.seed(1) 
test <- runif(50, 0, 1) 
Tcase <- rnorm(50, 25, 4) 
Fcase <- rnorm(50, 28, 4.3) 
ifelse(test, Tcase, Fcase)

無論是OK的ifelse()生成兩個之前做TRUE_CASE和FALSE_CASE的載體比較。

示例輸出：

R> ifelse(test < 0.69, Tcase, Fcase) 
[1] 24.77549 24.37682 19.11699 28.31967 26.67177 25.55472 27.41873 26.55069 
[9] 24.78478 19.49176 23.34002 23.42284 24.76275 29.40010 29.14852 24.34191 
[17] 33.19383 32.98973 27.22665 34.82338 30.40149 26.45833 28.07413 24.55062 
[25] 28.52443 26.59242 22.55189 26.36448 28.67952 30.73209 32.92160 23.53111 
[33] 20.82346 27.27888 35.23336 34.60647 26.01493 27.75896 25.20201 22.02691 
[41] 26.31093 17.78017 26.79981 25.61301 33.69045 25.82438 22.16021 27.44291 
[49] 27.22791 27.56918

另一種方法是做，這是隻生成所需的值

want <- test < 0.69 
res <- test ## copy vector of correct length 
res[want] <- rnorm(S <- sum(want), 25, 4) 
res[!want] <- rnorm(length(res) - S, 28, 4.3) 

R> res 
[1] 27.85067 24.70574 24.84946 22.04696 22.27336 36.03795 29.82793 23.70292 
[9] 25.24064 22.64442 27.12598 18.92642 26.22623 18.85420 26.97382 23.79610 
[17] 32.55148 31.81162 22.88688 25.33725 37.48624 22.39162 24.77241 17.34256 
[25] 29.70633 18.34011 23.14588 20.53632 26.90338 21.99672 33.34867 25.06958 
[33] 19.85480 18.43758 21.87467 26.80075 27.37908 24.92576 28.89241 23.72773 
[41] 37.92431 21.28255 28.45495 19.05016 20.69923 29.96509 29.00012 22.51493 
[49] 27.66824 26.56380

這顯然給了不同的數字與前面的，但僅僅是因爲我們產生更少隨機數字。對於大多數問題，我懷疑這會比ifelse()更有效率。

來源

2012-08-01 10:59:40

非常好，謝謝。我剛剛偶然發現了ifelse – YesSure 2012-08-01 11:12:21

@ Gavin只是爲了完整說明ifelse正在做什麼：使用runif在0和1之間生成的50個隨機數中的第一個，如果該值小於0.69，則從正態分佈（50,25,4），如果它更大，則從（50,28,4.3）生成一個值。然後最終的輸出就是將這些值合併成一個50個值的列表/矩陣。我的摘要是否對所執行的程序正確？ – YesSure 2012-08-01 11:28:43

是的，除了最後一位;結果是50個值的數字向量。我已經添加了一個輸出示例。 – 2012-08-01 11:42:15

需要any或all，其將邏輯值的向量成一個邏輯（當然，實際上是一個長度爲1的向量）。

來源

2012-08-01 10:59:35

用在R如果功能

回答

相關問題