R語言練習(xí)習(xí)題

1、統(tǒng)計軟件實驗1每題需注意1命令代碼2結(jié)果或圖形3自己的錯誤1 y = y <- sin(10*pi)*exp+42)+log(23)/log(4)2 x = sin(223/3)， y = x2，z = y*10 ；求x+2y-5zx <- sin(223/3);y <- x2;x+2*y-5*z3 建立起始值=3,增量值=,終止值=44的一維數(shù)組x x <- seq,by=4 建立等差一維數(shù)組x：首項為0,末項為,項數(shù)為15 x <- seq(0,pi,length=15)5 將100,200,200,200,400,600,600,600,800 輸入R中，保

2、存到numeric變量中numeric <- c(100,200,200,200,400,600,600,600,800)6 將numeric轉(zhuǎn)換為factor存入變量，并用class（）確認。 <- (numeric)7 查看的內(nèi)容8 創(chuàng)建一個2到50的向量,形式為 2, 4, 6, 8, ., 48, 50并名為為vector1vector1 <- seq(from=2,to=50,by=2)9 選取vector1中的第20個元素 vector12010 選取vector1中的第10，15，20個元素 vector1c(10,15,20)11 選取vector1中的第10

3、到20個元素vector10:2012 選取vector1中值大于40的元素vector1vector1>4013創(chuàng)建向量1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5rep(1:5,5)14使用rep()創(chuàng)建向量0 0 0 0 0 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4rep(0:4,rep(5,5)15 用函數(shù)rep()構(gòu)造一個向量x,它由3個3,4個2,5個1構(gòu)成x=c(rep(3,3),rep(2,4),rep(1,5)統(tǒng)計軟件實驗2每題需注意1命令代碼2結(jié)果或圖形3自己的錯誤1計算行列式

4、的值 A=det(matrix(c(3,4,7,2,2,8,3,6,1),nrow=3)2 矩陣,矩陣;求出A x B及A與B中對應(yīng)元素之間的乘積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構(gòu)成兩個方陣，其中矩陣按列輸入，矩陣按行輸入，并計算C=A+B,D=ABA=matrix(1:16,nrow=4);B=matrix(1:16,nrow=4,byrow=T);C=A+B;D=A%*%B4先復(fù)制附錄數(shù)據(jù)至文本文檔,然后讀取數(shù)據(jù)至文件datadata <- (&q

5、uot;D:/",header = T)5 比值的計算weight和height的平方的比值存入bmi變量 bmi <- weight/height2;detach(data)6 創(chuàng)建對象x, 其值為1:10, 使用write函數(shù)將其寫入文件；刪除x, 然后再讀入該文件并賦值給x,并保證 x是 numericx <- 1:10(x,file = "");rm(x);xx <- ("",header=T);xclass(x);x <- (x$x);x7 查看mtcars數(shù)據(jù)（mtcars回車）；把vs變量所在列的元素全都改

6、成你的“學(xué)號后兩位”（mtcars$vs<-7）；把mtcars存為文件；讀入文件存入mtcarsnew；MtcarsMtcars$vs<-學(xué)號后兩位mtcarsnew<()8 把mtcars轉(zhuǎn)換成矩陣mm并判斷數(shù)據(jù)類型；把第六行改成你的“學(xué)號后兩位”；取矩陣mm前11行存入變量mtcars11；取mtcars11的主對角線的元素，構(gòu)成主對角陣mt；取mtcars11的上三角陣存入mtupper；(需要通過網(wǎng)絡(luò)搜索學(xué)習(xí)上三角矩陣定義)mm<(mtcars);class(mm)mm6,<-學(xué)號后兩位mtcars11(mtcars11<-head(mtcars,

7、11)mt(mt<-diag(diag(mtcars11)mtupper(x(x)<-0;mtupper<-x)x(x)<-0(x,diag=FALSEx(x)<-0(x,diag=FALSE)9把mtupper的行名和列名改為NULL。rownames(mtupper)<-NULLcolnames(mtupper)<-NULLrownames(mtupper)<-NULLcolnames(mtupper)<-NULLrownames(mtupper)<-NULLcolnames(mtupper)<-NULLrownames(m

8、tupper)<-NULLcolnames(mtupper)<-NULL10. mtupper每一行求和，存入mtsum。mtsum<-apply(mtupper,MARGIN=1,FUN=sum)附錄：數(shù)據(jù)：1 身高體重數(shù)據(jù)weight height 統(tǒng)計軟件實驗31 數(shù)據(jù)的讀?。〝?shù)據(jù)見附錄） data<("路徑/",header=T)2 繪圖 plot(height, weight) # 繪圖3 更改繪圖的參數(shù)的取值例如：plot(height, weight, pch=2) #更改后的繪圖 5.根據(jù)cityrain數(shù)據(jù)作圖。令Y軸范圍（0,30

9、0），xaxt="n"，type="b",顏色為黑，X軸標題為“Month”,Y軸標題為“Tokyo Rainfall(mm)”，主標題“Monthly Rainfall in major cities”。用axis()函數(shù)添加X軸，使在X軸1到12的位置分別對應(yīng)12個月份。rain <- ("")plot(rain$Tokyo,type = "b",xaxt="n",ylim=c(0,300),col="black",xlab="Month",yl

10、ab="Rainfall(mm)",main="Monthly Rainfall in major cities")axis(1,at=1:12,labels = rain$Month)4. 在0,4pi畫sin(x),cos(x)(在同一個圖象中); 其中sin(x)和cos(x)圖象用不同的顏色和形狀表示，并在函數(shù)圖上適當?shù)奈恢脴俗?“箭頭+y=sin(x)”, “箭頭+y=cos(x)” ,標記x軸,y軸,標題為“正弦余弦函數(shù)圖象”.plot(sin,0,4*pi,main="正弦余弦函數(shù)圖像",xlab="x軸&qu

11、ot;,ylab="y軸",col="red",type="b",pch=18)curve(cos,0,4*pi,col="blue",type="b",pch=1,add=T)arrows,text,"sin(x)")arrows,text(10,"cos(x)")例如：附錄：數(shù)據(jù)：1 身高體重數(shù)據(jù)weight height 統(tǒng)計軟件實驗41使用三種循環(huán)，輸出向量1:100中所有數(shù)據(jù)。for (I in 1:100) print(I)while (i&l

12、t;=100) sum <- sum+ii <- i+1 print(sum)repeat if(i%2!=0)sum <- sum+i i <- i+1 if(i>100)break print(sum)2使用while循環(huán)求1+2+3+100的和。 i <- 1sum <- 0while(i<=100) sum<- sum+1 i <- i+1print(sum)3使用repeat循環(huán)求1至100之間的奇數(shù)和。 i <- 1sum <- 0while(i<=100) sum<- sum+1 i <-

13、i+1print(sum)4使用for循環(huán)輸出1至100之間的能夠被3或5整除的數(shù)，并求和。 sum <- 0for(i in 1:100) if(i%3=0|i%5=0) print(i) sum <- sum+i print(sum)5定義一個長度為100、mean=1，sd=1的正態(tài)分布隨機向量，并找出該向量的最小值及其出現(xiàn)的位置。 rnorm（n,mean,sd），種子設(shè)為123 (123)norm <- rnorm(100,1,1) <- 100for(i in 1:100) if(normi<= =normi =i cat,6找出100次產(chǎn)生的長度為1

14、00、mean=1，sd=1的正態(tài)分布隨機向量的最小值與出現(xiàn)位置的關(guān)系。(用plot,表示，其中為位置向量，為最小值向量)for(i in 1:100) norm <- rnorm(100,1,1) i <- 100 for(i in 1:100) if(normi<=i) i=normi i=i plot,7 Calculate the first 50 powers of 2, 2*2, 2*2*2, etc.Calculate the squares of the integer numbers from 1 to 50.Which pairs are equal, . which integer numbers fulfill the condition .How many pairs are there(Use