作者piq (享)
看板Statistics
標題[程式] R中的if條件
時間Tue Mar 19 16:35:07 2013
[軟體程式類別]:
R
[程式問題]:
if條件問題
[軟體熟悉度]:
低(1~3個月)
[問題敘述]:
v是一個矩陣
a是一個我們任意給的數(已知)
今天我希望v中的某行(假設第一行) 中的各元素做是否等於a的判斷運算
我希望能在不用for下 完成我想要的運算
但是一直想不出
[程式範例]:
以下為大略範例
>test<-function(v,a){
+if(v[,1]==a){執行A運算}
+else{執行B運算}}
>
Warning message:
In if (a == v[, 1]) { :
條件的長度 > 1，因此只能用其第一元素
===
想請問高手能有辦法步用for的狀況下
比對v[,1]中的元素是否等於a做邏輯運算嗎?
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 219.85.83.97
1F：推 ching0629:for (i in 1:length(v[,1])) { 03/19 16:39
2F：→ ching0629:if (v[i,1]==a) {A運算} 03/19 16:39
3F：→ ching0629:else {B運算} 03/19 16:40
4F：→ ching0629:} 03/19 16:40
5F：→ ching0629:沒看到不用for,我還以為指定用for 03/19 16:41
6F：→ ching0629:要做A運算的DATA=v[v[,1]==a,] 03/19 16:43
7F：→ ching0629:要做B運算的DATA=v[v[,1]!=a,] 03/19 16:43
8F：→ ching0629:像上面那樣分兩個DATA處理吧 03/19 16:44
9F：→ piq:看不太懂樓上DATA的想法 DATA是指邏輯算式? 03/19 16:54
10F：推 Wush978:從問題的敘述中我不確定 all(a == v[,1]) 是不是你要的 03/19 17:34
11F：推 Yukirin:直接下 matrix[,1] = a 會返回一個邏輯向量 03/19 17:38
12F：→ Yukirin:還是要全都等於a才做第一個算式? 03/19 17:39
13F：→ gsuper:apply(v,2,function(vector){vector==a}) 03/19 19:47
14F：→ gsuper:apply(v==a,2,which) 03/19 19:47
15F：→ piq:應該說是立即判斷 如果==a 就執行A !=a 就執行B 03/19 20:12
16F：→ piq:就是v[i,1]==a 執行A 反之B 是否有可以不用for的方法? 03/19 20:31
17F：→ gsuper:寫是可以寫 用雙apply可以取代雙回圈 但重點是你的分析 03/19 22:31
18F：→ gsuper:會花多少資源 用apply不會比回圈快 03/19 22:31
19F：→ gsuper:apply(v,2,function(x){apply(as.matrix(x),1,fx)} 03/19 22:32
20F：→ gsuper:把你的分析包在 fx 這個 function 裡面就可以 03/19 22:32
21F：→ piq:我現在改用一樓的方法 但是跑非常久 其實自己也知道最簡單的 03/20 08:55
22F：→ piq:方法是for 但是v[,1]的長度非常大 所以不是很想要用這招 03/20 08:56
23F：推 HeroNoah:這是你想要的 ? 03/20 11:03
24F：→ HeroNoah:x <- array( 03/20 11:03
25F：→ HeroNoah: sample(1 : 5, 30, replace = TRUE), 03/20 11:03
26F：→ HeroNoah: dim = c(10, 3) 03/20 11:03
27F：→ HeroNoah:) 03/20 11:04
28F：→ HeroNoah:AAA <- function(x) paste(x, "A運算") 03/20 11:04
29F：→ HeroNoah:BBB <- function(x) paste(x, "B計畫") 03/20 11:04
30F：→ HeroNoah:y <- ifelse(x[, 1] > 3, AAA(x), BBB(x)) 03/20 11:04
31F：→ HeroNoah:y 03/20 11:04
32F：→ HeroNoah:要等於a, 就把 > 3 改成 == a 這樣就是了 03/20 11:09
33F：推 ching0629:為什麼不分成兩個檔案跑呢? 我覺得你描述的清楚一點會 03/20 17:12
34F：→ ching0629:對你更有幫助,我覺得我提供的第二個方法也很理想 03/20 17:12
35F：→ ching0629:我的意思是說,你先把原始檔分成要做A運算的檔案及要做B 03/20 17:13
36F：→ ching0629:運算的檔案,然後這兩個檔案分別在每列做運算就好了 03/20 17:13