3.3 不能向量化
2016-04-12 21:33
176 查看
3.3 不能向量化
有一些操作是不能被向量化的.比如,如果当前迭代建立在上一次迭代的结果上,这样的情况下向量化没有用武之地.(但是一些情况被filter(),cumsum()掩盖了,等等.)如果你确实需要使用一个循环,那就让它尽可能简洁:
尽可能将更多地操作放在循环之外.例如:如果一个相同或者相似的的序列在每次循环都会被创建,那就在循环开始前首先创建它并且重复使用.创建一个序列是非常快速的,但是如果创建成千上万次的话那就是一笔昂贵的时间开销了.
尽可能将循环的次数变少.如果你在一个因素的元素层面进行迭代或者在这个因素的级别层面进行迭代之间可以选择的话,那就在它的级别层面进行迭代,几乎可以肯定的是,这样更好.
下面的代码返回矩阵的每一列元素之和(假设列的数量是正数):
sumxcol <- numeric(ncol(x)) for(i in 1:ncol(x)) sumxcol[i] <- sum(x[,i])
一个更普遍的方式肯能是这样的:
sumxcol <- apply(x, 2, sum)
除了这种常用的方法外,下面的代码是在R代码中没有包含循环的做此类操作的特殊的函数:
sumxcol <- colSums(x)
当然还有函数rowSums(),colMeans(),rowMeans().
另外一个方式是:
# %*%为求矩阵内积 sumxcol <- rep(1, nrow(x)) %*% x
这种方法用了矩阵的数乘,稍微动一下脑筋很多问题就会被蹂躏成矩阵数乘的形式.这确实是一个有效的替代方法.
相关文章推荐
- scoi2010
- ThinkPHP 笔记2 路由 控制器 视图
- 【面试】小东最少需要拉多少票才能让她获胜(京东实习笔试题)
- javascript copy拷贝对象
- 关于Maven-3.3.3版本的问题
- 【SSM-MyBatis框架】关联查询--多对多查询
- Xtext内部类型的设置
- 负载均衡(一)基于接口、方法、参数的路由
- java解析XML Node与Element的区别
- 技术记录---Handle收不到消息的问题分析
- ZZULI 1869: Mathematics and Geometry【数学】
- 异或运算实现两个整数的交换
- 正则表达式
- ACM:POJ-3009 Curling2.0
- 关于ereg()函数过时问题
- LoadRunner如何监控Linux系统资源
- pip install提示编码错误
- java.lang.NullPointerException
- Androd自定义控件(五)打造自己的Camera
- ThinkPHP 笔记1 起步