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個樣品的差異情況，方差越小，說明樣品差異越小。