如果运算不复杂，又找不到好的vectorize的方式 就可以请Rcpp出场了 以下搭配inline套件服用 ※ 引述《piq (享)》之铭言： : 感谢大家的建议 : 其实一开始最简单的方法是用for : 但是因为真实在跑分析时 : v[,1]的长度大约是17万6000个 : 也就是说v是一个176000x20up的矩阵 假设是v是numeric : 如果这时候用for(i 1:length(v[,1]))是可以跑... : 但是会跑非常久~我并不想要这样 : 而A,B两运算是要把v[,1]重新编程一个号码 只是规则不同 : 有想过用ifelse但似乎是失败了 : ifelse(v[,1]==a,A,B) : 分析上 A运算是：v[i,1]-a : B运算是：12-a+v[i,1] library(Rcpp) library(inline) v <- matrix(rnorm(176000*20), 176000, 20) v[,1] <- sample(1:10, 176000, TRUE) f_R <- function(v, a) { for(i in 1:nrow(v)) { if (v[i,1] == a) { v[i,1] <- v[i,1] - a } else { v[i,1] <- 12 - a + v[i,1] } } return(v) } f_Rcpp <- cxxfunction(sig=c(Rv = "numeric", Ra = "numeric"), plugin="Rcpp", body=' NumericMatrix v(Rv); double a(as<double>(Ra)); for(int i = 0;i < v.nrow();i++) { if (v(i,0) == a) { v(i,0) = v(i,0) - a; } else { v(i,0) = 12 - a + v(i,0); } } ') system.time({v_R <- f_R(v, 1)}) # 1.228秒 system.time(f_Rcpp(v, 1)) # Rcpp版本会直接改v本身, 0.001秒 all.equal(v_R, v) # 一样 : 所以我当初是写成: ifelse(v[,1]==a,v[i,1]-a,12-a+v[i,1]) : 但是不能跑.. : 莫非是ifelse里面不能用逻辑算式? : 还有是recodes指令... : 我的R只接跟我说没有这个指令 囧 : 所以我在想是否还有其他的方法~"~ : : [软体程式类别]: : : R : : [程式问题]: : : if条件问题 : : [软体熟悉度]: : : 低(1~3个月) : : [问题叙述]: : : v是一个矩阵 : : a是一个我们任意给的数(已知) : : 今天我希望v中的某行(假设第一行) 中的各元素做是否等於a的判断运算 : : 我希望能在不用for下 完成我想要的运算 : : 但是一直想不出 : : [程式范例]: : : 以下为大略范例 : : +if(v[,1]==a){执行A运算} : : +else{执行B运算}} : : Warning message: : : In if (a == v[, 1]) { : : : 条件的长度 > 1，因此只能用其第一元素 : : === : : 想请问高手能有办法步用for的状况下 : : 比对v[,1]中的元素是否等於a做逻辑运算吗? --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 61.230.45.47 ※ 编辑: Wush978 来自: 61.230.45.47 (03/21 00:36)