天问

R语言笔记(1)

1.  (1) 写出元素为3, -1.5, 3E-10的向量。

c(3, -1.5, 3e-10)

   (2) 写出从3开始每次增加3,长度为100的向量。

seq(3,by=3,length=100)

   (3) 写出(0, 2)重复10次的向量。

rep(c(0,2),10)

   (4) 对向量x,写出其元素大于等于0小于1的条件。

if( all( x >= 0 && x < 1) )

   (5) 对向量x,写出其元素都等于0的条件。

if( all(x == 0 ) )

   (6) 写出包含12个月份名称的向量。

month.name[1:12]

2. x为一个长100的整数向量。比如,x <- floor(100*runif(100))

floor(100*runif(100))

(1) 显示x2130号元素。

x[c(21:30)]

(2) x313539号元素赋值为0

x[c(31,35,39)]=c(0,0,0)

(3) 显示x中除了第1号和第50号的元素之外的子集。

x[-c(1,50)]

(4) 列出x中个位数等于3的元素。

i=1

while(i<=100){

if(x[i]%%10!=3) {x[-c(i)] }

i=i+1

}

(5) 列出x中个位数等于3的元素的下标位置。Hintuse which( )

for( i in 1 : length(x) ){ if( x[ i ] %% 10 == 3 ){print(i)} }

3. 绘制y=x3 +sin(x)*cos(x)的函数图

x=1:100

y=x*x*x+sin(x)*cos(x)

plot(x,y,type='s')

#plot(x,y,xlab=’x’,ylab=’y’,main=’绘制y=x3 +sin(x)*cos(x)的函数图’,type=’s’)

rplot

4. 随机产生两个矩阵A,B,其维数都为5*6,矩阵中每一个元素服从标准正态分布。(1) A,B按行合并成一个矩阵C;(2)计算A’*B; (3) 计算矩阵C的行和及列和

AA=rnorm(30,0,1)

A=array(AA,c(5,6))

B=array(rnorm(30,0,1),c(5,6))

(1)C=rbind(A,B)

(2)t(A)%*%B

(3) 

行:sum(C[1,])   #计算矩阵第一行和,其他行同理

列:sum(C[,1])

博客地址:http://blog.yoqi.me/?p=831
扫我捐助哦
喜欢 39

这篇文章还没有评论

发表评论