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