## 讀入資料檔並將變數名稱字串轉為變數標籤
library(readr)
X1360 <- read.csv("1360.csv")
View(X1360)
nrow(X1360) # 共1006列
ncol(X1360) # 共46欄
library(sjlabelled)
X1360 <- read_csv("1360.csv")
save(X1360, file = "X1360.rda", compress = T)
names(X1360)
varlabels <- colnames(X1360)
colnames(X1360)[1:46] <- paste("V", 1:46, sep="")
sjlabelled::set_label(X1360) <- varlabels
save(X1360, file = "X1360.rda", compress = T)
X1360 <- sjmisc::set_na(X1360, na="NA")
varlabels #列出每個變數的標籤
library(sjmisc)
library(sjPlot)
## #refugeeswelcome
names(X1360)
# 他人:是指陌生人,而非親朋好友
# (V25)在網路世界中,你會更願意與人互動嗎?(1~6)
# table(X1360$V25)
library("sjlabelled")
X1360$V25 <- set_label(X1360$V25, label="在網路上是否比現實生活更願意與他人互動")
table(X1360$V25)
##
## 1 2 3 4 5 6
## 126 557 17 180 118 8
X1360$V25r <- rec(X1360$V25, rec="1:3=1[yes];4:6=0[no] ", as.num=F)
frq(X1360$V25)
##
## 在網路上是否比現實生活更願意與他人互動 (x) <numeric>
## # total N=1006 valid N=1006 mean=2.63 sd=1.28
##
## Value | N | Raw % | Valid % | Cum. %
## --------------------------------------
## 1 | 126 | 12.52 | 12.52 | 12.52
## 2 | 557 | 55.37 | 55.37 | 67.89
## 3 | 17 | 1.69 | 1.69 | 69.58
## 4 | 180 | 17.89 | 17.89 | 87.48
## 5 | 118 | 11.73 | 11.73 | 99.20
## 6 | 8 | 0.80 | 0.80 | 100.00
## <NA> | 0 | 0.00 | <NA> | <NA>
frq(X1360$V25r)
##
## 在網路上是否比現實生活更願意與他人互動 (x) <categorical>
## # total N=1006 valid N=1006 mean=0.70 sd=0.46
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | no | 306 | 30.42 | 30.42 | 30.42
## 1 | yes | 700 | 69.58 | 69.58 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>
par(family="jf-jinxuan-fresh2.1-Medium")
plot_frq(X1360$V25r)
# 年輕世代:30歲以下 1
# 中老年世代:31歲以上 0
X1360$V2<- set_label(X1360$V2, label = "世代")
table(X1360$V2)
##
## 1934 1939 1940 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962
## 1 1 1 3 1 2 1 3 7 6 4 3 1 4 4 5
## 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978
## 12 6 6 11 6 11 12 21 14 17 18 15 16 21 23 28
## 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994
## 26 17 36 43 33 39 42 30 33 55 30 54 25 37 41 30
## 1995 1996 1997 1998 1999 2000 2001 2002 2003
## 43 22 27 16 17 12 11 1 2
#hist(X1360$age)
X1360$V2r<- rec(X1360$V2, rec="1934:1989=0[old]; 1990:2020=1[youth]", as.num=F)
frq(X1360$V2r)
##
## 世代 (x) <categorical>
## # total N=1006 valid N=1006 mean=0.34 sd=0.47
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | old | 668 | 66.40 | 66.40 | 66.40
## 1 | youth | 338 | 33.60 | 33.60 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>
plot_frq(X1360$V2r)
# (V30)在高度匿名網路平台中,發言時會讓你更加暢所欲言、還是會讓你更加小心翼翼?
# 我覺得更加暢所欲言 1
# 我覺得更加小心翼翼 2
# 我覺得與現實生活不會有明顯的差別 3
X1360$V30 <- set_label(X1360$V30,label="相信網路匿名性")
table(X1360$V30)
##
## 1 2 3 90
## 274 515 210 7
X1360$V30r <- rec(X1360$V30, rec="1=1[yes]; 2:3=0[no]", as.num=F)
frq(X1360$V30r)
##
## 相信網路匿名性 (x) <categorical>
## # total N=1006 valid N=999 mean=0.27 sd=0.45
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | no | 725 | 72.07 | 72.57 | 72.57
## 1 | yes | 274 | 27.24 | 27.43 | 100.00
## <NA> | <NA> | 7 | 0.70 | <NA> | <NA>
plot_frq(X1360$V30r)
# (V32)你相不相信網路世界可以做到「完全匿名」這件事?
# 相信,現在科技進步有新技術可以加密 1
# 相信,只要大家不知道我是誰就是匿名 2
# 不相信,仍有許多知名網路公司被盜取資料 3
# 不相信,都會有被破解的一天只是早晚而已 4
X1360$V32 <- set_label(X1360$V32,label="相信網路匿名性的程度")
table(X1360$V32)
##
## 1 2 3 4 90
## 98 90 417 394 7
X1360$V32r <- rec(X1360$V32, rec="1:2=1[high]; 3:4=0[low]", as.num=F)
frq(X1360$V32r)
##
## 相信網路匿名性的程度 (x) <categorical>
## # total N=1006 valid N=999 mean=0.19 sd=0.39
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | low | 811 | 80.62 | 81.18 | 81.18
## 1 | high | 188 | 18.69 | 18.82 | 100.00
## <NA> | <NA> | 7 | 0.70 | <NA> | <NA>
plot_frq(X1360$V32r)
# (V36)你覺得社群網路可以彌補現實世界「人際關係」的不足嗎?
# 可以,它可以更讓我暸解朋友 1
# 可以,它可以幫助我交到新朋友 2
# 不行,無助於現實生活的人際關係 3
# 不行,大家都是隱姓埋名,關係是脆弱的 4
X1360$V36 <- set_label(X1360$V36,label="相信網路可以彌補現實生活中所缺乏的人際關係")
table(X1360$V36)
##
## 1 2 3 4 90
## 149 309 284 249 15
X1360$V36r <- rec(X1360$V36, rec="1:2=1[yes]; 3:4=0[no]", as.num=F)
frq(X1360$V36r)
##
## 相信網路可以彌補現實生活中所缺乏的人際關係 (x) <categorical>
## # total N=1006 valid N=991 mean=0.46 sd=0.50
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | no | 533 | 52.98 | 53.78 | 53.78
## 1 | yes | 458 | 45.53 | 46.22 | 100.00
## <NA> | <NA> | 15 | 1.49 | <NA> | <NA>
plot_frq(X1360$V36r)
save(X1360, file = "X1360.rda", compress = T)
load("X1360.rda")
library(sjmisc)
frq(X1360$V25r)
##
## 在網路上是否比現實生活更願意與他人互動 (x) <categorical>
## # total N=1006 valid N=1006 mean=0.70 sd=0.46
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | no | 306 | 30.42 | 30.42 | 30.42
## 1 | yes | 700 | 69.58 | 69.58 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>
frq(X1360$V2r)
##
## 世代 (x) <categorical>
## # total N=1006 valid N=1006 mean=0.34 sd=0.47
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | old | 668 | 66.40 | 66.40 | 66.40
## 1 | youth | 338 | 33.60 | 33.60 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>
frq(X1360$V30r)
##
## 相信網路匿名性 (x) <categorical>
## # total N=1006 valid N=999 mean=0.27 sd=0.45
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | no | 725 | 72.07 | 72.57 | 72.57
## 1 | yes | 274 | 27.24 | 27.43 | 100.00
## <NA> | <NA> | 7 | 0.70 | <NA> | <NA>
frq(X1360$V32r)
##
## 相信網路匿名性的程度 (x) <categorical>
## # total N=1006 valid N=999 mean=0.19 sd=0.39
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | low | 811 | 80.62 | 81.18 | 81.18
## 1 | high | 188 | 18.69 | 18.82 | 100.00
## <NA> | <NA> | 7 | 0.70 | <NA> | <NA>
frq(X1360$V36r)
##
## 相信網路可以彌補現實生活中所缺乏的人際關係 (x) <categorical>
## # total N=1006 valid N=991 mean=0.46 sd=0.50
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | no | 533 | 52.98 | 53.78 | 53.78
## 1 | yes | 458 | 45.53 | 46.22 | 100.00
## <NA> | <NA> | 15 | 1.49 | <NA> | <NA>
rm(list=ls()) # 清空所有暫存區的物件
load("X1360.rda")
names(X1360)
## [1] "V1" "V2" "V3" "V4" "V5" "V6" "V7" "V8" "V9" "V10"
## [11] "V11" "V12" "V13" "V14" "V15" "V16" "V17" "V18" "V19" "V20"
## [21] "V21" "V22" "V23" "V24" "V25" "V26" "V27" "V28" "V29" "V30"
## [31] "V31" "V32" "V33" "V34" "V35" "V36" "V37" "V38" "V39" "V40"
## [41] "V41" "V42" "V43" "V44" "V45" "V46" "V25r" "V2r" "V30r" "V32r"
## [51] "V36r"
library(sjPlot)
tab_xtab(X1360$V25r, X1360$V2r,encoding="utf8",
show.row.prc = TRUE,
show.col.prc = TRUE,
show.na = FALSE,
show.legend = FALSE,
show.exp = FALSE,
show.cell.prc = FALSE,
tdcol.col = "green",
tdcol.row = "blue"
)
| 在網路上是否比現實生活更願意與他人互動 | 世代 | Total | |
|---|---|---|---|
| old | youth | ||
| no |
212 69.3 % 31.7 % |
94 30.7 % 27.8 % |
306 100 % 30.4 % |
| yes |
456 65.1 % 68.3 % |
244 34.9 % 72.2 % |
700 100 % 69.6 % |
| Total |
668 66.4 % 100 % |
338 33.6 % 100 % |
1006 100 % 100 % |
χ2=1.454 · df=1 · φ=0.040 · p=0.228 |
tab_xtab(X1360$V25r, X1360$V30r,encoding="utf8",
show.row.prc = TRUE,
show.col.prc = TRUE,
show.na = FALSE,
show.legend = FALSE,
show.exp = FALSE,
show.cell.prc = FALSE,
tdcol.col = "green",
tdcol.row = "blue"
)
| 在網路上是否比現實生活更願意與他人互動 | 相信網路匿名性 | Total | |
|---|---|---|---|
| no | yes | ||
| no |
256 85 % 35.3 % |
45 15 % 16.4 % |
301 100 % 30.1 % |
| yes |
469 67.2 % 64.7 % |
229 32.8 % 83.6 % |
698 100 % 69.9 % |
| Total |
725 72.6 % 100 % |
274 27.4 % 100 % |
999 100 % 100 % |
χ2=32.803 · df=1 · φ=0.184 · p=0.000 |
tab_xtab(X1360$V25r, X1360$V32r,encoding="utf8",
show.row.prc = TRUE,
show.col.prc = TRUE,
show.na = FALSE,
show.legend = FALSE,
show.exp = FALSE,
show.cell.prc = FALSE,
tdcol.col = "green",
tdcol.row = "blue"
)
| 在網路上是否比現實生活更願意與他人互動 | 相信網路匿名性的程度 | Total | |
|---|---|---|---|
| low | high | ||
| no |
253 83.5 % 31.2 % |
50 16.5 % 26.6 % |
303 100 % 30.3 % |
| yes |
558 80.2 % 68.8 % |
138 19.8 % 73.4 % |
696 100 % 69.7 % |
| Total |
811 81.2 % 100 % |
188 18.8 % 100 % |
999 100 % 100 % |
χ2=1.319 · df=1 · φ=0.039 · p=0.251 |
tab_xtab(X1360$V25r, X1360$V36r,encoding="utf8",
show.row.prc = TRUE,
show.col.prc = TRUE,
show.na = FALSE,
show.legend = FALSE,
show.exp = FALSE,
show.cell.prc = FALSE,
tdcol.col = "green",
tdcol.row = "blue"
)
| 在網路上是否比現實生活更願意與他人互動 | 相信網路可以彌補現實生活中所缺乏的人際關係 | Total | |
|---|---|---|---|
| no | yes | ||
| no |
234 77 % 43.9 % |
70 23 % 15.3 % |
304 100 % 30.7 % |
| yes |
299 43.5 % 56.1 % |
388 56.5 % 84.7 % |
687 100 % 69.3 % |
| Total |
533 53.8 % 100 % |
458 46.2 % 100 % |
991 100 % 100 % |
χ2=93.530 · df=1 · φ=0.309 · p=0.000 |
load("X1360.rda")
library(car)
## 載入需要的套件:carData
## Registered S3 methods overwritten by 'car':
## method from
## influence.merMod lme4
## cooks.distance.influence.merMod lme4
## dfbeta.influence.merMod lme4
## dfbetas.influence.merMod lme4
mod.1 <- glm(V25r~ V2r+V30r+V32r+V36r,
data=X1360,
family=binomial)
summary(mod.1)
##
## Call:
## glm(formula = V25r ~ V2r + V30r + V32r + V36r, family = binomial,
## data = X1360)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -2.2910 -1.1865 0.6154 0.7906 1.2482
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 0.1255 0.1078 1.164 0.244
## V2r1 -0.1041 0.1631 -0.638 0.523
## V30r1 0.9813 0.1935 5.072 3.94e-07 ***
## V32r1 -0.1863 0.2025 -0.920 0.358
## V36r1 1.4424 0.1650 8.743 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 1204.5 on 982 degrees of freedom
## Residual deviance: 1077.6 on 978 degrees of freedom
## (因為不存在,23 個觀察量被刪除了)
## AIC: 1087.6
##
## Number of Fisher Scoring iterations: 4
vif(mod.1)
## V2r V30r V32r V36r
## 1.045772 1.016697 1.061723 1.061558
library(sjPlot)
plot_model(mod.1, type = "est", auto.label = F, colors = "gs")