统计软件分析与应用 实验报告.doc

PAGE\*MERGEFORMAT9

实验报告一

课程___统计软件分析与应用
学生姓名_________
学号___________
学院__数学与统计学院______
专业_____统计学__________
指导教师___________

二Ｏ一四年三月十日

1、打开帮助文档首页，并查阅其中的“introductiontoR”
>help.start()

2、看一下你的软件中有哪些程序包
>search()
[1]".GlobalEnv""package:stats""package:graphics"
[4]"package:grDevices""package:utils""package:datasets"
[7]"package:methods""Autoloads""package:base"

3、载入程序包“emplik”
>install.packages("emplik")
然后选择镜像

4、查看当前能访问的数据集列表
>data()

5、找到数据包“datasets”，并打开数据Orange，编辑数据，把修改后的数据存在onew中
>fix(Orange)

>sink("D:/R/onew")
>Orange
>sink()

6、显示数据集Orange的内容
>Orange
Treeagecircumference
1111830
2148458
3166487
411004115
511231120
611372142
711582145
8211833
9248469
102664111
1121004156
1221231172
1321372203
1421582203
15311830
16348451
17366475
1831004108
1931231115
2031372139
2131582140
22411832
23448462
244664112
2541004167
2641231179
2741372209
2841582214
29511830
30548449
31566481
3251004125
3351231142
3451372174
3551582177

7、运行Orange自带的示例
>?Orangestartinghttpdhelpserver...done
>require(stats);require(graphics)
>coplot(circumference~age|Tree,data=Orange,show.given=FALSE)
>fm1<-nls(circumference~SSlogis(age,Asym,xmid,scal),
+data=Orange,subset=Tree==3)
>plot(circumference~age,data=Orange,subset=Tree==3,
+xlab="Treeage(dayssince1968/12/31)",
+ylab="Treecircumference(mm)",las=1,
+main="Orangetreedataandfittedmodel(Tree3only)")
>age<-seq(0,1600,length.out=101)
>lines(age,predict(fm1,list(age=age)))

8、显示命令”lm”自带的示例
>?Lm

9、建立一个变量X，赋值为向量分量为1到10
>x<-1:10
>x
[1]12345678910

10、建立一个变量Y，赋值为向量2-20，且全为偶数。
y<-seq(2,20,2)
>y
[1]2468101214161820

11、求向量X中所有数据的均值
>mean(x)
[1]5.5

12、求向量Y中所有数据的均值
>mean(y)
[1]11
13、求这个向量X，Y中所有数据的均值
>a<-c(x,y)
>a
[1]123456789102468101214161820
>mean(a)
[1]8.25

14、画关于数据X，Y的散点图，并在当前图上加上连接线
>plot(x,y)
>plot(x,y,type="l")

15、显示你所定义的变量
>ls()
[1]"a""ctl""group""lm.D9""lm.D90""opar""trt""weight"
[9]"x""y""a"
16、删除变量Y
>rm(y)
17、建立一个变量Z，赋值为向量3.2，4.2，5.2，6.2
>z<-c(3.2,4.2,5.2,6.2)
>z
[1]3.24.25.26.2
18、定义z为X的前四个数即在z=(1,2,,3,4)
>z<-x[1:4]
>z
[1]123