R語(yǔ)言案例分析：多元數(shù)據(jù)的基本統(tǒng)計(jì)分析

我們利用該數(shù)據(jù)集中的Case1來(lái)完成下面的R語(yǔ)言操作：

> options(digits = 4) #輸出結(jié)果位數(shù)
> par(mar=c(4, 4, 2, 1) + 0.1, cex=0.8) # 圖形修飾
> case1 <- read.csv("clipboard", header=T, sep = "\t") #復(fù)制表中的數(shù)據(jù)，直接創(chuàng)建case1
> head(case1)

地區(qū) 性別 教育程度   觀點(diǎn) 年齡 月收入 月支出
1    A   女       中 不支持   55   2299   1423
2    A   女       低 不支持   39   3378   2022
3    A   女       中   支持   33   3460   1868
4    B   男       高  支持   41   4564   1918
5    B   女       高 不支持   55   3206   1906
6    A   女       中 不支持   48   4043   2233
> summary(case1)

地區(qū)    性別     教育程度     觀點(diǎn)     年齡          月收入         月支出    
 A:204   男:603   低:319   不支持:628   Min.   : 6.0   Min.   : 637    Min.   : 797  
 B:401   女:597   高:303   支持  :568  1st Qu.:34.0    1st Qu.:2388     1st Qu.:1722  
 C:384          中:578   NA's  : 4   Median :40.0    Median :2978     Median :1993  
 D:211                          Mean   :40.1   Mean   :3006    Mean   :1997  
                              3rd Qu.:47.0    3rd Qu.:3624     3rd Qu.:2262  
                              Max.   :72.0   Max.   :6239    Max.   :3385
# 定性分析
> attach(case1) #綁定數(shù)據(jù)
> T1 <- table(地區(qū))> T1

地區(qū)
A   B   C   D
204 401 384 211
> barplot(T1) #繪制條形圖

# 定量分析 > f <- hist(月收入) #直方圖

# 定性定量分析 > boxplot(月收入~性別) #箱線圖

> t.test(月收入~性別) #t檢驗(yàn)
Welch Two Sample t-test

data:  月收入 by 性別
t = 0.51, df = 1200, p-value = 0.6
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -75.43 128.49
sample estimates:
mean in group 男 mean in group 女
         3019         2993
# 接受男女的月收入無(wú)顯著差異的假設(shè)（p>0.5）

# 二維列聯(lián)表分析
> T2 <- table(性別, 觀點(diǎn))
> T2
     觀點(diǎn)
性別 不支持 支持
  男    319  282
  女    309  286
> barplot(T2, beside = TRUE) # 條形圖

> barplot(T2, beside = F) # 條形圖

# beside=T表示繪制分組條形圖，beside=F表示繪制堆疊條形圖

# 多維列聯(lián)表分析
> T3 <- ftable(性別, 教育程度, 觀點(diǎn)) # 創(chuàng)建一個(gè)緊湊的"平鋪"式列聯(lián)表
> T3
              觀點(diǎn)     不支持 支持
性別 教育程度                 
男   低                81   88
    高                78   66
    中                 160  128
女   低                82   68
    高                86   72
    中                 141  146
> barplot(T3, beside = TRUE, col = 3:4) #條形圖

> T4 <- ftable(教育程度, 性別, 觀點(diǎn))
> T4
                       觀點(diǎn) 不支持 支持
教育程度   性別                 
低       男           81   88
        女           82   68
高       男           78   66
        女           86   72
中       男           160  128
        女           141  146
> barplot(T4, beside = TRUE, col = 3:4) #條形圖

> detach(case1) #解除綁定