統計分析軟件課程模擬試卷（3）(時間:年月日，共100分鐘)48學時,3學分,閉卷,總分100分,占總評成績50%注：此頁不作答題紙，請將答案寫在答題紙上一、填空題(本題20分，每小題2分)1.x=c(7,12,9,15,8,7,14,16)，則R命令median(x)返回結果是:2.在R命令行鍵入rm()的含義是3.x<-c(23,42,25,25,37)，則R命令rank(x)的輸出結果是4.R命令2*c(1,NA,2)的輸出結果是5.R語句age=c(89,20,10,66);c('A','B')[(age>65)+1];的輸出結果是6.在R循環語句中，__________語句的作用是提前結束本層循環。7.在Python中__________表示空類型。8.Python中用pip工具升級科學計算擴展庫numpy的完整命令是______________9.Python語句aa=[1,2,3,4];aa.pop();則aa.index(3)結果為10.Python語句a=2；b=14；c=b**a；則c的值為：二、單項選擇題(本題20分，每小題2分)1.R語言中，下列（）不可以作為變量名。A.1a2b B.a1b2 C.a_b D.a.b2.下列R語句中（）不能生成向量（[1]1234512345）。A.rep(1:5,2) B.sequence(rep(5,2))C.1:10%%6D.(0:9%%5)+13.已知a<-1:100，若要取出向量a中的偶數位置上的元素，表達式應為（）。A.a[2:100] B.a/50C.a[seq(2,100,by=2)] D.a/24.若a<-matrix(1:12,c(3,4))，則>apply(a,1,sum)的運行結果為()A.222630 B.6152433C.232731 D.71625345.?R語言中，為變量修改元素名稱的函數是()A.changeName B.ModifyNameC.name D.names6.有如下兩個R語言函數： g<-function(){ x<-2} f<-function(){ x<-1 g() return(x)}則f()的運行結果是（）A.1B.2C.程序報錯D.以上均不正確7.R語言中，下列命令中哪一條不能實現：A.X<-matrix(1:10,nrow=2,ncol=5);Y=X-XB.X<-matrix(1:10,nrow=2,ncol=5);Y=X%*%XC.X<-matrix(1:10,nrow=2,ncol=5);Y=X*XD.X<-matrix(1:10,nrow=2,ncol=5);Y=X+X8.?Python3.x環境創建Python腳本文件時，需要將文件編碼格式設置為（）。A.UNICODEB.ANSIC.GBKD.UTF-89.已知在?Python腳本文件N.py中有函數調用“A.B.C.d()”，則import語句的正確寫法是（）。A.importA.B.C.dB.importA.B.CC.fromA.B.CimportdD.fromA.BimportC10.關于Python中的函數，以下描述錯誤的是（）。A．函數能完成特定的功能，對函數的使用不需要了解函數內部實現原理，只要了解函數的輸入輸出方式即可。B．使用函數的主要目的是減低編程難度和代碼重用。C．Python使用del保留字定義一個函數。D．函數是一段具有特定功能的語句組。三、簡答題(本題20分，共4題，每小題5分)1.用R函數rep()構造一個向量x,它由5個3，10個2，98個1構成，寫出對應的R語句。2.R語言中表述概率分布的函數前分別加有d、p、q、r四個前綴，分別表示該分布的什么功能？3.R語言中，對lm()函數返回的對象使用plot()函數可生成評價模型擬合情況的四幅圖形，請簡述這些圖形的主要作用。4.簡述Python中元組、列表和字典的區別。四、分析題(本題20分)用Pearson擬合優度卡方檢驗驗證學生成績是否服從正態分布，R代碼如下:X<-c(25,45,50,54,55,61,64,68,72,75,75,78,79,81,83,84,84,84,85,86,86,86,87,89,89,89,90,91,91,92,100)A<-table(cut(X,br=c(0,69,79,89,100)))#代碼Ap<-pnorm(c(70,80,90,100),mean(X),sd(X))pnew<-c(p[1],p[2]-p[1],p[3]-p[2],1-p[3])#代碼Bchisq.test(A,p=pnew)#代碼C1.Pearson擬合優度卡方檢驗的基本思想是什么?（5分）2.代碼A的主要功能是什么？（5分）3．代碼B的主要功能是什么，為什么要這么處理？（5分）4．代碼C的運行結果如下：data:AX-squared=8.3,df=3,p-value=0.04據此，你的結論是什么？（5分）五、改錯題(本題10分)下面R程序用循環實現求x1，若x1達到允許誤差范圍內停止循環。x0<-1;tol<-1e-8while{x1<-computeEstimate()#計算x1值的函數if(abs(x1-x0)<tol){continue}else{x0<-x1}}請指出上述程序中的錯誤(1-2處)，并給出正確的寫法。六、編程題(本題10分)應用逆變換方法生成1000個隨機數服從下列概率密度函數，寫出對應的R代碼。

試卷答案一、填空題（每題2分，共20分）1.15.42.設置工作目錄31223335.1236.1234247.a[seq(2,100,2)]8.rownames(UGdata)<-NULL9.110.deldic['abc']二、單項選擇題（每題2分，共20分）題號12345678910答案CABBDCCDBA（無嚴格標準答案，以下內容可供參考）三、簡答題（20分，每題5分）1.A=matix(1:16nrow=4);B=matix(1:16,mow=4,byrow=T);2.C=A+B:D=A%*%B3.優點：免費，永遠正版，資源公開(不是黑匣子)，可以在UNIX,Windows和MacosX上運行，R有優秀的內在幫助系統，R有優秀的畫圖功能，學生能夠輕松地轉到商業支持的S-Plus程序；R語言有一個強大的，容易學習的語法，有許多內在的統計函數（每答對一點的一分，5分為止，言之有理即可）4.set.seed(123);u=runif(100);x=sqrt(u,2)（最后一步3分，前兩步各一分）四、分析題（15分）1.（5分）set.seed(1234)（1分）num<-seq(1:90)+1303160000regression<-round(runif(90,min=70,max=100))software<-round(rnorm(90,mean=75,sd=15))software[which(software>100)]<-100（4分）2．（5分）grade<-data.frame(num,regression,software)write.csv(grade,file="學號.csv")3.（5分） high<-apply(grade[c("regression","software")],2,max)s.sum<-apply(grade[c("regression"