高效R语言编程

162

｜

第

7

章

练习

考虑下面的代码段：

double test1() {

double a = 1.0 / 81;

double b = 0;

for (int i = 0; i < 729; ++ i)

b = b + a;

return b;

}

1.

将

test1() 函数保存为一个单独的文件。确保它正常运行。

2.

使用

R

写一个对应函数，比较

C++

与

R

版本的速度。

3.

将 double 型变量替换成 float 型，创建一个 test2() 函数。这样做后，还

能得到正确结果吗？

4.

将 b = b + a 改写成 b += a，从而使你的代码更像

c++

。

5.

（难度较大）i++ 与 ++i 的区别是什么

?

矩阵

每个向量类型有着对应的矩阵类型：NumericMatrix、IntegerMatrix、

CharacterMatrix 与 LogicalMatrix。这些类型的使用方式与 NumericVector

的使用方式类似，主要的区别是：

•

当我们初始化时，我们需要指定行数和列数：

// 10 rows, 5 columns

NumericMatrix mat

（

10, 5

）

;

// Length 10

NumericVector v

（

10

）

;

•

使用

()

提取子集，例如，

mat(5, 4)

。

•

矩阵的第一个元素是

mat(0, 0)

，记住索引是从

0

而不是

1

开始。

•

我们使用

.nrow()

与

.ncol()

方法，确定行和列的数目。