R语言练习题.doc

统计软件实验1
每题需注意1命令代码2结果或图形3自己的错误
1 y =
y <- sin(10*pi)*exp(-+4^2)+log(23)/log(4)
2 x = sin(223/3), y = x^2,z = y*10 ;求x+2y-5z
x <- sin(223/3);y <- x^2;x+2*y-5*z
3 建立起始值=3,增量值=,终止值=44的一维数组x
x <- seq(,by=)
4 建立等差一维数组x:首项为0,末项为,项数为15
x <- 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 ,并用class()确认。
<- (numeric)
7

8 创建一个2到50的向量,形式为 2, 4, 6, 8, ..., 48, 50并名为为vector1
vector1 <- 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创建向量1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5
rep(1:5,5)
14使用rep()创建向量0 0 0 0 0 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4
rep(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自己的错误
1计算行列式的值
A=det(matrix(c(3,4,7,2,2,8,3,6,1),nrow=3))
2 矩阵,矩阵;求出A x B及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*B
3由1,2,...,16构成两个方阵,其中矩阵A按列输入,矩阵B按行输入,并计算C=A+B,D=AB
A=matrix(1:16,nrow=4);B=matrix(1:16,nrow=4,byrow=T);C=A+B;D=A%*%B
4先复制附录数据至文本文档,然后读取数据至文件data
data <- ("D:/",header = T)
5 比值的计算weight和height的平方的比值存入bmi变量
bmi <- weight/height^2;detach(data)
6 创建对象x, 其值为1:10, ;删除x, 然后再读入该文件并赋值给x,并保证 x是 numeric
x <- 1:10
(x,file = "");rm(x);x
x <- ("",header=T);x
class(x);x <- (x$x);x
7 查看mtcars数据(mtcars回车);把vs变量所在列的元素全都改成你的“学号后两位”(mtcars$vs<-7);;;
Mtcars
Mtcars$vs<-学号后两位
mtcarsnew<-()
8 把mtcars转换成矩阵mm并判断数据类型;把第六行改成你的“学号后两位”;取矩阵mm前11行存入变量mtcars11;取mtcars11的主对角线的元素,构成主对角阵mt;取mtcars11的上三角阵存入mtupper;(需要通过网络搜索学习上三角矩阵定义)
mm<-(mtcars);class(mm)
mm[6,]<-学号后两位
mtc