R语言入门与数据分析

R语言环境安装

安装文件下载地址：https://mirrors.tuna.tsinghua.edu.cn/CRAN/

这里可以根据你的操作系统选择合适的安装程序下载安装，我这里是mac版本。安装完成后的界面如下所示

HelloWorld

现有一个研究维生素C对于牙齿生长的数据。

	len	supp	dose
1	4.2	VC	0.5
2	11.5	VC	0.5
3	7.3	VC	0.5
4	5.8	VC	0.5
5	6.4	VC	0.5
6	10	VC	0.5
7	11.2	VC	0.5
8	11.2	VC	0.5
9	5.2	VC	0.5
10	7	VC	0.5
11	16.5	VC	1
12	16.5	VC	1
13	15.2	VC	1
14	17.3	VC	1
15	22.5	VC	1
16	17.3	VC	1
17	13.6	VC	1
18	14.5	VC	1
19	18.8	VC	1
20	15.5	VC	1
21	23.6	VC	2
22	18.5	VC	2
23	33.9	VC	2
24	25.5	VC	2
25	26.4	VC	2
26	32.5	VC	2
27	26.7	VC	2
28	21.5	VC	2
29	23.3	VC	2
30	29.5	VC	2
31	15.2	OJ	0.5
32	21.5	OJ	0.5
33	17.6	OJ	0.5
34	9.7	OJ	0.5
35	14.5	OJ	0.5
36	10	OJ	0.5
37	8.2	OJ	0.5
38	9.4	OJ	0.5
39	16.5	OJ	0.5
40	9.7	OJ	0.5
41	19.7	OJ	1
42	23.3	OJ	1
43	23.6	OJ	1
44	26.4	OJ	1
45	20	OJ	1
46	25.2	OJ	1
47	25.8	OJ	1
48	21.2	OJ	1
49	14.5	OJ	1
50	27.3	OJ	1
51	25.5	OJ	2
52	26.4	OJ	2
53	22.4	OJ	2
54	24.5	OJ	2
55	24.8	OJ	2
56	30.9	OJ	2
57	26.4	OJ	2
58	27.3	OJ	2
59	29.4	OJ	2
60	23	OJ	2

这里的第一列是牙齿的长度，第二列是维生素还是橙汁，第三列是剂量。将该数据保存到一个csv文件中。

> df <- read.table("/Users/admin/Downloads/RData/ToothGrowth.csv",sep=",",header=T)
> df

以上是读入并显示数据。

样品频数显示

> table(df$supp,df$dose)

显示结果

     0.5  1  2
  OJ  10 10 10
  VC  10 10 10

它数据中表示每种种类(OJ为橙汁，VC为VC药品)的不同剂量各有10个样品。

查看每种种类，每种剂量的牙齿生长平均值

> aggregate(df$len,by=list(df$supp,df$dose),FUN=mean)

显示结果

  Group.1 Group.2     x
1      OJ     0.5 13.23
2      VC     0.5  7.98
3      OJ     1.0 22.70
4      VC     1.0 16.77
5      OJ     2.0 26.06
6      VC     2.0 26.14

由以上数据，我们可以看出在0.5和1mg的剂量下，橙汁明显比VC药品的平均值要大很多，这说明橙汁更有利于牙齿生成。

查看每种种类，每种剂量的牙齿生长方差值

> aggregate(df$len,by=list(df$supp,df$dose),FUN=sd)

显示结果

  Group.1 Group.2        x
1      OJ     0.5 4.459709
2      VC     0.5 2.746634
3      OJ     1.0 3.910953
4      VC     1.0 2.515309
5      OJ     2.0 2.655058
6      VC     2.0 4.797731

方差用来反映每种种类，每种剂量的10个样品的差异情况，方差越小，说明样品差异越小。