作者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