本博客所有内容是原创,未经书面许可,严禁任何形式的转载。
http://blog.csdn.net/u010255642
1)array使用
Array函数的参数有3个,第一个是需要形成数组元素的数据,第二个是dim参数提示维度
> c(1:20)->h
> mya<-array(h,dim=c(4,5))
> mya
[,1] [,2] [,3] [,4] [,5]
[1,] 1 5 9 13 17
[2,] 2 6 10 14 18
[3,] 3 7 11 15 19
[4,] 4 8 12 16 20
> mydim<-c(2,10)
> mya<-array(h,dim=c(2,10))
> mya
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 1 3 5 7 9 11 13 15 17 19
[2,] 2 4 6 8 10 12 14 16 18 20
>
> dim(mya)
[1] 2 10
>
第一个参数既可以是向量也可以是单个值
> mya<-array(1,dim=c(2,10))
> mya
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 1 1 1 1 1 1 1 1 1 1
[2,] 1 1 1 1 1 1 1 1 1 1
2)数组运算
逐元素运算
> mya
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 1 3 5 7 9 11 13 15 17 19
[2,] 2 4 6 8 10 12 14 16 18 20
> myb
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 2 2 2 2 2 2 2 2 2 2
[2,] 2 2 2 2 2 2 2 2 2 2
> mya+myb
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 3 5 7 9 11 13 15 17 19 21
[2,] 4 6 8 10 12 14 16 18 20 22
> mya*myb
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 2 6 10 14 18 22 26 30 34 38
[2,] 4 8 12 16 20 24 28 32 36 40
> 3*mya*myb
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 6 18 30 42 54 66 78 90 102 114
[2,] 12 24 36 48 60 72 84 96 108 120
>
> mya*myb+mya
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 3 9 15 21 27 33 39 45 51 57
[2,] 6 12 18 24 30 36 42 48 54 60
>
2个数组的外积
定义以下向量:
列向量 u(b1,b2,b3,b4)
行向量 v(a1,a2,a3)
它们的外积 被定义为:
> b<-array(c(1:4))
> a<-array(c(5:6))
> b%o%a
[,1] [,2]
[1,] 5 6
[2,] 10 12
[3,] 15 18
[4,] 20 24
> b
[1] 1 2 3 4
> a
[1] 5 6
>
再举一个例子
> b<-array(c(1:4))
> a<-array(c(5:8))
> a*b
[1] 5 12 21 32
> b
[1] 1 2 3 4
> a
[1] 5 6 7 8
> a%o%b
[,1] [,2] [,3] [,4]
[1,] 5 10 15 20
[2,] 6 12 18 24
[3,] 7 14 21 28
[4,] 8 16 24 32
>
生成的数组向量则由 2个数数组向量元素所有可能乘积得到