非參數(shù)統(tǒng)計(jì)r軟件參考答案

1、內(nèi)容：內(nèi)容：A.3A.10A.12A.3上機(jī)實(shí)踐：將MASS數(shù)據(jù)包用命令library(MASS)加載到R中，調(diào)用自帶“老忠實(shí)”噴泉數(shù)據(jù)集geyer，它有兩個(gè)變量：等待時(shí)間waiting和噴涌時(shí)間duration，其中…(1)將等待時(shí)間70min以下的數(shù)據(jù)挑選出來(2)將等待時(shí)間70min以下且等待時(shí)間不等于57min的數(shù)據(jù)挑選出來(3)將等待時(shí)間70min以下噴泉的噴涌時(shí)間挑選出來(4)將噴涌時(shí)間大于70min噴泉的等待時(shí)間挑選出來。

2、解:讀取數(shù)據(jù)的R命令：library(MASS)#加載MASS包data(geyser)#加載數(shù)據(jù)集geyserattach(geyser)#將數(shù)據(jù)集geyser的變量置為內(nèi)存變量(1)依題意編定R程序如下：sub1geyser=geyser[which(waiting70)1]#提取滿足條件（waiting70）的數(shù)據(jù)which()，讀取下標(biāo)Sub4geyser[1:5]#顯示子數(shù)據(jù)集sub1geyser的前5行[1]80718075

3、77…….A.10如光盤文件student.txt中的數(shù)據(jù)，一個(gè)班有30名學(xué)生，每名學(xué)生有5門課程的成績，編寫函數(shù)實(shí)現(xiàn)下述要求：(1)以data.frame的格式保存上述數(shù)據(jù)；(2)計(jì)算每個(gè)學(xué)生各科平均分，并將該數(shù)據(jù)加入(1)數(shù)據(jù)集的最后一列；(3)找出各科平均分的最高分所對(duì)應(yīng)的學(xué)生和他所修課程的成績；(4)找出至少兩門課程不及格的學(xué)生，輸出他們的全部成績和平均成績；(5)比較具有(4)特點(diǎn)學(xué)生的各科平均分與其余學(xué)生平均分之間是否存在差

4、異。先將數(shù)據(jù)集讀入R系統(tǒng)student=read.table(“…“header=T)class(student):#顯示數(shù)據(jù)集student的類型，(5)依題意，要?jiǎng)?chuàng)造兩個(gè)子集data4和data2用兩樣本的比較方法比較他們的平均成績是否有顯著差異。類似創(chuàng)造data1的方法，創(chuàng)造data2。并設(shè)x=data1$meany=data2$mean比較二樣本xy是否有顯著差異，由于還沒有學(xué)非參數(shù)檢驗(yàn)，試用t檢驗(yàn)檢驗(yàn)之(R的t檢驗(yàn)函數(shù)為t.t

5、est(xy)，原假設(shè)H0是兩樣本的均值相等，備擇假設(shè)H1是兩樣本不等)。如果P值pvalue0.05則拒絕原假設(shè)。data2=student[1]k=0f(iin1:30)x=as.numeric(student[i2:6])if(sum(x60)2)k=k1data2[k]=student[i]下面做t檢驗(yàn)x=data1$meany=data2$meant.test(xy)WelchTwoSamplettestdata:xyt=3.

6、0236df=9.309pvalue=0.01386alternativehypothesis:truedifferenceinmeansisnotequalto095percentconfidenceinterval:11.4932361.684037sampleestimates:meanofxmeanofy：67.1750073.76364結(jié)論：pvalue=0.013860.05拒絕原假設(shè)，即認(rèn)為兩樣本的平均成績有顯著差異。A.

7、12在一張圖上，用取值（10，10）之間間隔均等的1000個(gè)點(diǎn)，采用不同的線型一顏色給制sin()cos()sin()cos()的函數(shù)圖形，圖形要求有主標(biāo)題和副標(biāo)題，標(biāo)示出從坐標(biāo)x=seq(1010length=50)#構(gòu)造向量xx[1:5]#顯示x的前5個(gè)數(shù)據(jù)[1]10.000009.979989.959969.939949.91992sin=sin(x)#計(jì)算sin函數(shù)值cos=cos(x)sc=sin(x)cos(x)plot(s