【R】数值型向量及运算

Python018

【R】数值型向量及运算,第1张

向量 是将若干个基础类型相同的值存储在一起, 各个元素可以按序号访问。 如果将若干个数值存储在一起可以用序号访问, 就叫做一个数值型向量

用c()函数把多个元素或向量组合成一个向量。如

10:13这样的写法表示从10到13的整数组成的向量。

用print()函数显示向量或在命令行中显示向量时, 每行显示的行首会有方括号和数字序号, 代表该行显示的第一个向量元素的下标。如

length(x)可以求x的长度。 长度为零的向量表示为numeric(0)。 numeric()函数可以用来初始化一个指定元素个数而元素都等于零的数值型向量, 如numeric(10)会生成元素为10个零的向量。

单个数值称为标量, R没有单独的标量类型, 标量实际是长度为1的向量。

R中四则运算用+ - * / ^表示(加、减、乘、除、乘方),如

R中四则运算仍遵从通常的优先级规则, 可以用圆括号()改变运算的先后次序。 如

除了加、减、乘、除、乘方, R还支持整除运算和求余运算。 用%/%表示整除,用%%表示求余。如

向量与标量运算

向量与标量的运算为每个元素与标量的运算, 如

一个向量乘以一个标量, 就是线性代数中的数乘运算。

四则运算时如果有缺失值,缺失元素参加的运算相应结果元素仍缺失。 如

等长向量的运算为对应元素两两运算。 如

两个等长向量的加、减运算就是线性代数中两个向量的加、减运算。

两个不等长向量的四则运算, 如果其长度为倍数关系,规则是每次从头重复利用短的一个。 如

不仅是四则运算,R中有两个或多个向量按照元素一一对应参与某种运算或函数调用时, 如果向量长度不同,一般都采用这样的规则。

如果两个向量的长度不是倍数关系,会给出警告信息。如

R中的函数一般都是向量化的: 在R中, 如果普通的一元函数以向量为自变量,一般会对每个元素计算。 这样的函数包括sqrt, log10, log, exp, sin, cos, tan等许多。 如

为了查看这些基础的数学函数的列表,运行命令help.start(), 点击链接“Search Engine and Keywords”, 找到“Mathematics”栏目, 浏览其中的“arith”和“math”链接中的说明。 常用的数学函数有:

有一些不太常用的数学函数:

如果自己编写的函数没有考虑向量化问题, 可以用Vectorize()函数将其转换成向量化版本。

sort(x)返回排序结果。 rev(x)返回把各元素排列次序反转后的结果。 order(x)返回排序用的下标。如

例子中, order(x)结果中3是x的最小元素11所在的位置下标, 1是x的第二小元素33所在的位置下标, 2是x的最大元素55所在的位置下标。

sum(求和), mean(求平均值), var(求样本方差), sd(求样本标准差), min(求最小值), max(求最大值), range(求最小值和最大值)等函数称为统计函数, 把输入向量看作样本,计算样本统计量。 prod求所有元素的乘积。

cumsum和cumprod计算累加和累乘积。如

其它一些类似函数有pmax, pmin, cummax, cummin等。

1. 显示1到100的整数的平方根和立方根(提示:立方根就是三分之一次方)。

2. 设有10个人的小测验成绩为:

解答:

你可以先按照名字给他们排序一下

a <- a[order(names(a))]

b <- b[order(names(b))]

这样他们俩个就是一样的顺序了

R语言中两个数组(或向量)的外积如何计算

所谓数组(或向量)a和b的外积,指的是a的每一个元素和b的每一个元素搭配在一起相乘得到的新元素.当然运算规则也可自定义.外积运算符为 %o%(注意:百分号中间的字母是小写的字母o).例如:

[plain] view plain copy

>a <- 1:2

>b <- 3:5

>d <- a %o% b

>d

[,1] [,2] [,3]

[1,]345

[2,]68 10

注意维数公式为:

[plain] view plain copy

dim(d) = c( dim(a) , dim(b) )

实际上R语言提供了一个更为一般化得外积函数outer(X,Y,Fun),比较简单,看个示例吧.

下面的例子是为了绘制 z = x**2 +50cos(y)的3维表面图:

[plain] view plain copy

>x <- seq(from = -5, to = 5, by = 0.25)

>y <- seq(from = -2, to = 2, by = 0.1)

>f <- function(x,y){x**2+50*cos(y)}

>z <- outer(x,y,f)

>persp(x,y,z)

绘制图像如下:

其它看帮助文档吧.