
R语言练习试题.docx
8页R语言练习试题精选文档统计软件实验1每题需注意1命令代码2结果或图形3自己的错误21y=sin(10)e(0.34)log423y<-sin(10*pi)*exp(-0.3+4^2)+log(23)/log(4)2 x=sin(223/3),y=x^2,z=y*10;求x+2y-5zx<-sin(223/3);y<-x^2;x+2*y-5*z3 建立初步值=3,增量值=5.5,停止值=44的一维数组xx<-seq(3.47,by=5.5)4建立等差一维数组x:首项为0,末项为,项数为15x<-seq(0,pi,length=15)5将100,200,200,200,400,600,600,600,800输入R中,保存到numeric变量中numeric<-c(100,200,200,200,400,600,600,600,800)6将numeric变换为factor存入变量,并用class()确认factor.numeric<-as.factor(numeric)7查察的内容8创立一个2到50的向量,形式为2,4,6,8,...,48,50并名为为vector1vector1<-seq(from=2,to=50,by=2)9采用vector1中的第20个元素vector1[20]10采用vector1中的第10,15,20个元素vector1[c(10,15,20)]11采用vector1中的第10到20个元素vector[10:20]12采用vector1中值大于40的元素vector1[vector1>40]13 创立向量1234512345123451234512345rep(1:5,5)14 使用rep( )创立向量0000011111222223333344444rep(0:4,rep(5,5))15用函数rep( )构造一个向量x,它由3个3,4个2,5个1构成x=c(rep(3,3),rep(2,4),rep(1,5))统计软件实验2每题需注意1命令代码2结果或图形3自己的错误/.精选文档3231计算行列式的值A426781A=det(matrix(c(3,4,7,2,2,8,3,6,1),nrow=3))3231112矩阵A426,矩阵B222781333;求出AxB及A与B中对应元素之间的乘积A=matrix(c(3,4,7,2,2,8,3,6,1),nrow=3);B=matrix(rep(1:3,3),nrow=3);A%*%B;A*B3由1,2,...,16构成两个方阵,此中矩阵A按列输入,矩阵B按行输入,并计算C=A+B,D=ABA=matrix(1:16,nrow=4);B=matrix(1:16,nrow=4,byrow=T);C=A+B;D=A%*%B4先复制附录数据至文本文档,而后读取数据至文件datadata<-read.table("D:/t.data.txt",header=T)5比值的计算weight和height的平方的比值存入bmi变量bmi<-weight/height^2;detach(data)6创立对象x,其值为1:10,使用write函数将其写入文件;删除x,而后再读入该文件并赋值给x,并保证x是numericx<-1:10write.table(x,file="x.txt");rm(x);xx<-read.table("x.txt",header=T);xclass(x);x<-as.numeric(x$x);x7查察mtcars数据(mtcars回车);把vs变量所在列的元素全都改成你的“学号后两位”(mtcars$vs<-7);把mtcars存为文件;读入文件存入mtcarsnew;MtcarsMtcars$vs<-学号后两位mtcarsnew<-read.csv( )8把mtcars变换成矩阵mm并判断数据种类;把第六行改成你的“学号后两位”;取矩阵mm前11行存入变量mtcars11;取mtcars11的主对角线的元素,构成主对角阵mt;取mtcars11的上三角阵存入mtupper;(需要经过网络搜寻学习上三角矩阵定义)mm<-as.matrix(mtcars);class(mm)mm[6,]<-学号后两位mtcars11(mtcars11<-head(mtcars,11))mt(mt<-diag(diag(mtcars11)))mtupper(x[lower.tri(x)]<-0;mtupper<-x)x[lower.tri(x)]<-0.精选文档lower.tri(x,diag=FALSEx[upper.tri(x)]<-0upper.tri(x,diag=FALSE)9把mtupper的行名和列名改为NULL。
rownames(mtupper)<-NULLcolnames(mtupper)<-NULL10.mtupper每一行乞降,存入mtsummtsum<-apply(mtupper,MARGIN=1,FUN=sum)附录:数据:1 身高体重数据weightheight统计软件实验3rownames(mtupper)<-NULL1数据的读取(数据见附录)colnames(mtupper)<-NULL.精选文档data<-read.table("路径/t.data.txt",header=T)2 绘图plot(height,weight)#绘图3 改正绘图的参数的取值比方:plot(height,weight,pch=2)#改正后的绘图5.依据cityrain数据作图令Y轴范围(0,300),xaxt="n",type="b",颜色为黑,X轴标题为“Month”,Y轴标题为“TokyoRainfall(mm)”,主标题“MonthlyRainfallinmajorcities”用axis( )函数增加X轴,使在X轴1到12的地点分别对应12个月份rain<-read.csv("cityrain.csv")plot(rain$Tokyo,type="b",xaxt="n",ylim=c(0,300),col="black",xlab="Month",ylab="Rainfall(mm)",main="MonthlyRainfallinmajorcities")axis(1,at=1:12,labels=rain$Month)4.在[0,4pi]画sin(x),cos(x)(在同一个图象中);此中sin(x)和cos(x)图象用不一样的颜色和形状表示,并在函数图上合适的地点注明“箭头+y=sin(x)”,“箭头+y=cos(x)”,标志x轴,y轴,标题为“正弦余弦函数图象”.plot(sin,0,4*pi,main="正弦余弦函数图像",xlab="x轴",ylab="y轴",col="red",type="b",pch=18)curve(cos,0,4*pi,col="blue",type="b",pch=1,add=T)arrows(3.3,0.8,2.3,0.8)text(3.8,0.8,"sin(x)")arrows(10.8,0.8,11.8,0.8)text(10,0.8,"cos(x)")比方:.精选文档附录:数据:1 身高体重数据weightheight统计软件实验41使用三种循环,输出向量1:100中全部数据。
①for(Iin1:100)print(I)②while(i<=100){sum<-sum+ii<-i+1}print(sum).精选文档③repeat{if(i%%2!=0)sum<-sum+ii<-i+1if(i>100)break}print(sum)2使用while循求1+2+3+⋯+100的和i<-1sum<-0while(i<=100){sum<-sum+1i<-i+1}print(sum)3使用repeat循求1至100之的奇数和i<-1sum<-0while(i<=100){sum<-sum+1i<-i+1}print(sum)4使用for循出1至100之的能被3或5整除的数,并乞降sum<-0for(iin1:100){if(i%%3==0|i%%5==0){print(i)sum<-sum+i}}print(sum)5定一个度100、mean=1,sd=1的正分布随机向量,并找出向量的最小及其出的地点rnorm(n,mean,sd),种子123set.seed(123)norm<-rnorm(100,1,1)min.norm<-100for(iin1:100){if(norm[i]<=min.norm){min.norm=norm[i]min.count=i}}cat(min.norm,min.count)。












