文章标题

R学习tips1

今天打算画出非参数统计一书中P13几内亚猪的经验生存函数图
数据是pig.rar中的6个txt文件，分别代表每个实验组72头猪的死亡时间
文本文件中第一行没有标题。每个文本文件只有一列

1.导入数据

y=read.table("C:\\Users\\Administrator\\Desktop\\pig\\gpigs1.txt",header=FALSE)其中目录一定要打两个反斜杠
将存在桌面的txt文件导入了R中，getwd（）可以查看当前目录，setwd（）可以用来设置当前路径
导入之后，需要查看导入的变量在R中储存为什么类型，可以用
str（对象名）或者typeof（对象名）发现默认地储存为了data frame类型，且每个data frame中包含了一个数据列，数据列的名称分别为V1，V2。。。

2.将数据合并

v1<-y1$V1

以上的命令就是将y1这个data.frame中的域v1转为向量v1，通过类似的命令，将y1-y6都转为向量，最后，我还需要一个自然增长的序列代表2年的天数。这个由v7=1:800完成

probability	v1,2,3,4,5
1	排序之后的死亡天数
2	排序用sort(向量名称)
3	…
…	…

3.计算猪的生存比例

利用
p6=1:65
p6=(65-p6)/65计算probability，现在得到了每组猪的（死亡率，死亡天数）=（v1，p1）的数据，现在要将这六组数据画在同一张图中。

4.画图

plot(v1,p1,type="l")
lines(v2,p2)
lines(v3,p3)
...
lines(v6,p6)

先使用plot函数创建了一个折线图，其中type的参数值“l”代表图像是线，如果是b则代表是点和线。

lines函数则是向该折线图中添加其他线条。

初步画出了几内亚猪的生存概率图。