高雄市長罷免電訪
M096020001 林庭如
2020/12/02
研究動機
2018舉行高雄市長選舉，最終由韓國瑜當選。多家媒體爭相報導韓流帶起的現象，有分析指出韓國瑜是草根型候選人，與之前國民黨候選人不同，其支持者大多屬於農工階級。此外，台灣媒體日漸有政黨色彩，報導的內容逐漸有了特定立場，因此形成了不同政黨的支持者有特定喜好的新聞頻道。因此想透過分析高雄市長罷免電訪，來找出受訪者罷免韓國瑜與否，與其職業和平常看的新聞台之間的關聯性。
研究假設
假設一:支持罷免韓國瑜的受訪者，其職業非農工階級。
假設二:支持罷免韓國瑜的受訪者，平常收看不偏藍的電視新聞台。
變數說明
本次資料分析選擇以投票支持罷免韓國市長與否為依變數，自變數則以受訪者的職業，以及受訪者平常喜愛看的新聞台。
資料說明
這次實作分析資料是使用政大選舉研究中心的研究資料，受訪對象為設籍在臺灣地區戶籍在高雄市且年滿二十歲以上的成年人。
抽樣方法有二，第一是以電話簿抽樣法，以「中華電信住宅部107-108年版電話號碼簿」為母體清冊，依據高雄各區電話簿刊登的電話數佔高雄市所刊電話總數比例，決定各區抽出之電話數比例，以等距抽樣法抽出各區電話樣本後，為求涵蓋的完整性，再以隨機亂數修正電話號碼的最後二碼或四碼，以求接觸到未登錄電話的住宅戶。其二為利用利用政治大學選舉研究中心歷年所累積的電話號碼資料庫，先按照各縣市比例以機率與大小成比例抽樣法（Probabilities Proportional to Size，簡稱PPS）抽出局碼，再以隨機亂數產生電話號碼的最後四碼來製作電話樣本。
訪問期間自109年6月10日（星期三）至6月14日（星期日），本次訪問預定完成市話1100個樣本，經實際訪問完成1,106個有效樣本，以95%之信心水準估計，最大可能隨機抽樣誤差為：±2.95%。
讀取sav資料檔、描述與移除個資欄位
使用sjlabelled::read_spss()讀取sav檔案
library(sjlabelled)
ED206462 <- read_spss("ED206462.sav")
str(ED206462, list.len=5) # 為節省版面只顯示其中的五筆## 'data.frame': 1106 obs. of 124 variables:
## $ ID : num 10001 10014 10030 10051 10063 ...
## ..- attr(*, "label")= chr "樣本編號"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ NAME : chr "唯一合格女" "最年輕的男" "最年輕的男" "最年輕的女" ...
## ..- attr(*, "label")= chr "受訪者姓名"
## ..- attr(*, "format.spss")= chr "A10"
## ..- attr(*, "display_width")= int 12
## $ TEL_A : num 7 7 7 7 7 7 7 7 7 7 ...
## ..- attr(*, "label")= chr "電話區碼"
## ..- attr(*, "format.spss")= chr "F5.0"
## ..- attr(*, "display_width")= int 7
## $ TEL : num 1e+07 1e+07 1e+07 1e+07 1e+07 ...
## ..- attr(*, "label")= chr "電話號碼"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## $ TEL_NAME : chr "04" "04" "04" "04" ...
## ..- attr(*, "label")= chr "電話套數"
## ..- attr(*, "format.spss")= chr "A6"
## ..- attr(*, "display_width")= int 10
## [list output truncated]ED206462$TEL<-NULL存成R專用的rda檔
save(ED206462, file="ED206462_raw.rda", compress=TRUE)讀入R資料檔
load("ED206462_raw.rda")觀察值
nrow(ED206462)## [1] 1106變數及名稱
names(ED206462)## [1] "ID" "NAME" "TEL_A" "TEL_NAME" "DATE"
## [6] "BTIME" "ETIME" "RESULT" "CALLER" "GENDER"
## [11] "LEVEL" "HOURS" "VERSION" "HOME_A" "HOME_M"
## [16] "HOME_O" "SAMPLE" "DATE_1" "BTIME_1" "ETIME_1"
## [21] "RESULT_1" "DATE_2" "BTIME_2" "ETIME_2" "RESULT_2"
## [26] "DATE_3" "BTIME_3" "ETIME_3" "RESULT_3" "DATE_4"
## [31] "BTIME_4" "ETIME_4" "RESULT_4" "DATE_5" "BTIME_5"
## [36] "ETIME_5" "RESULT_5" "DATE_6" "BTIME_6" "ETIME_6"
## [41] "RESULT_6" "DATE_7" "BTIME_7" "ETIME_7" "RESULT_7"
## [46] "Q1" "Q2" "Q3" "Q4" "Q5"
## [51] "Q6" "Q7" "Q8" "Q9" "Q10"
## [56] "Q11" "Q12" "Q13" "Q14" "Q15"
## [61] "Q16" "Q17" "Q18" "Q19" "Q20"
## [66] "Q21" "Q22" "Q23" "Q24" "Q25"
## [71] "Q26" "Q27" "Q28" "Q29" "Q30"
## [76] "Q31" "Q32" "Q33" "Q34" "Q35"
## [81] "Q36" "Q37" "Q38" "Q39" "Q40"
## [86] "NQ1" "NQ2" "NQ3" "NQ4" "NQ5"
## [91] "NQ6" "NQ7" "NQ8" "NQ9" "NQ10"
## [96] "NQ11" "NQ12" "NQ13" "NQ14" "NQ15"
## [101] "NQ16" "NQ17" "NQ18" "NQ19" "NQ20"
## [106] "NQ21" "NQ22" "NQ23" "NQ24" "NQ25"
## [111] "NQ26" "SEX" "age" "edu" "sengi"
## [116] "MASENGI" "T_Cidentity" "partyid" "PARTY" "CAREER"
## [121] "CAREER8" "tondu" "w"資料結構
str(ED206462)## 'data.frame': 1106 obs. of 123 variables:
## $ ID : num 10001 10014 10030 10051 10063 ...
## ..- attr(*, "label")= chr "樣本編號"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ NAME : chr "唯一合格女" "最年輕的男" "最年輕的男" "最年輕的女" ...
## ..- attr(*, "label")= chr "受訪者姓名"
## ..- attr(*, "format.spss")= chr "A10"
## ..- attr(*, "display_width")= int 12
## $ TEL_A : num 7 7 7 7 7 7 7 7 7 7 ...
## ..- attr(*, "label")= chr "電話區碼"
## ..- attr(*, "format.spss")= chr "F5.0"
## ..- attr(*, "display_width")= int 7
## $ TEL_NAME : chr "04" "04" "04" "04" ...
## ..- attr(*, "label")= chr "電話套數"
## ..- attr(*, "format.spss")= chr "A6"
## ..- attr(*, "display_width")= int 10
## $ DATE : num 200610 200613 200610 200612 200612 ...
## ..- attr(*, "label")= chr "訪問日期"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME : num 202059 190500 201324 203250 215214 ...
## ..- attr(*, "label")= chr "開始時間"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ ETIME : num 202331 191357 202354 204502 220151 ...
## ..- attr(*, "label")= chr "結束時間"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ RESULT : num 0 0 0 0 0 0 0 0 0 0 ...
## ..- attr(*, "label")= chr "訪問結果"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "labels")= Named num [1:38] 0 10 20 30 40 50 60 70 80 81 ...
## .. ..- attr(*, "names")= chr [1:38] "訪問成功" "受訪者暫時不在 或 不便接聽 (非當日約訪者)" "受訪者因臨時有事而中途拒訪 (非當日約訪者)" "受訪者暫時不在 或 不便接聽 (可當日約訪者)" ...
## $ CALLER : num 1.07e+08 1.00e+05 1.00e+05 1.00e+05 1.00e+05 ...
## ..- attr(*, "label")= chr "訪員編號"
## ..- attr(*, "format.spss")= chr "F10.0"
## ..- attr(*, "display_width")= int 12
## $ GENDER : num 2 2 2 2 2 2 2 2 2 1 ...
## ..- attr(*, "label")= chr "訪員性別"
## ..- attr(*, "format.spss")= chr "F1.0"
## ..- attr(*, "labels")= Named num [1:2] 1 2
## .. ..- attr(*, "names")= chr [1:2] "男" "女"
## $ LEVEL : num 1 1 1 3 1 1 1 1 1 1 ...
## ..- attr(*, "label")= chr "訪員級數"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 7
## ..- attr(*, "labels")= Named num [1:3] 1 2 3
## .. ..- attr(*, "names")= chr [1:3] "第一級" "第二級" "第三級"
## $ HOURS : num 156 984 1420 144 252 213 120 992 984 528 ...
## ..- attr(*, "label")= chr "累計時數"
## ..- attr(*, "format.spss")= chr "F7.0"
## ..- attr(*, "display_width")= int 9
## $ VERSION : num 1 1 1 1 1 1 1 1 1 1 ...
## ..- attr(*, "label")= chr "版本"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 9
## $ HOME_A : num 1 4 3 4 2 2 4 5 2 7 ...
## ..- attr(*, "label")= chr "戶中全部合格人數"
## ..- attr(*, "format.spss")= chr "F2.0"
## $ HOME_M : num 0 2 2 2 1 1 2 2 1 4 ...
## ..- attr(*, "label")= chr "戶中合格男性人數"
## ..- attr(*, "format.spss")= chr "F2.0"
## $ HOME_O : num 1 4 3 4 1 1 1 5 1 1 ...
## ..- attr(*, "label")= chr "戶中排行"
## ..- attr(*, "format.spss")= chr "F2.0"
## $ SAMPLE : num 1 1 1 1 2 1 2 1 1 2 ...
## ..- attr(*, "label")= chr "受訪樣本"
## ..- attr(*, "format.spss")= chr "F1.0"
## ..- attr(*, "labels")= Named num [1:2] 1 2
## .. ..- attr(*, "names")= chr [1:2] "正取樣本" "替代樣本"
## $ DATE_1 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_1"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_1 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_1"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_1 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_1"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_1 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_1"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ DATE_2 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_2"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_2 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_2"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_2 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_2"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_2 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_2"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ DATE_3 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_3"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_3 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_3"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_3 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_3"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_3 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_3"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ DATE_4 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_4"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_4 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_4"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_4 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_4"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_4 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_4"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ DATE_5 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_5"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_5 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_5"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_5 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_5"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_5 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_5"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ DATE_6 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_6"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_6 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_6"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_6 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_6"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_6 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_6"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ DATE_7 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問日期_7"
## ..- attr(*, "format.spss")= chr "F6.0"
## $ BTIME_7 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "開始時間_7"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ ETIME_7 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "結束時間_7"
## ..- attr(*, "format.spss")= chr "F6.0"
## ..- attr(*, "display_width")= int 9
## $ RESULT_7 : num NA NA NA NA NA NA NA NA NA NA ...
## ..- attr(*, "label")= chr "訪問結果_7"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 10
## $ Q1 : num 14 21 11 9 12 5 4 19 2 5 ...
## ..- attr(*, "label")= chr "首先，想請問您的戶籍是在高雄市_________區"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:40] 1 2 3 4 5 6 7 8 9 10 ...
## .. ..- attr(*, "names")= chr [1:40] "鹽埕區" "鼓山區" "左營區" "楠梓區" ...
## $ Q2 : num 92 92 96 92 1 8 5 6 92 4 ...
## ..- attr(*, "label")= chr "請問您平常最常看哪一臺的電視新聞？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:20] 1 2 3 4 5 6 7 8 9 10 ...
## .. ..- attr(*, "names")= chr [1:20] "臺視" "中視" "華視" "民視" ...
## $ Q3 : num 98 2 98 3 2 1 1 3 2 4 ...
## ..- attr(*, "label")= chr "請問您對韓國瑜擔任高雄市長以來的整體施政表現滿不滿意？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常滿意" "滿意" "不滿意" "非常不滿意" ...
## $ Q4 : num 98 2 3 2 3 96 1 2 2 1 ...
## ..- attr(*, "label")= chr "有人說：「政治人物當選後，很快就會跟民意脫節。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q5 : num 98 1 3 3 3 2 4 2 3 1 ...
## ..- attr(*, "label")= chr "有人說：「政治人物並不會關心像我這樣一般百姓的想法。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q6 : num 98 4 3 2 2 96 4 3 2 4 ...
## ..- attr(*, "label")= chr "有人說：「政治人物會用心為一般百姓做事。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q7 : num 98 2 2 2 2 1 1 3 3 4 ...
## ..- attr(*, "label")= chr "有人說：「市政府所做的事大多數是正確的。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q8 : num 98 1 2 2 2 4 4 2 3 1 ...
## ..- attr(*, "label")= chr "有人說：「市政府官員時常浪費一般民眾所繳納的稅金。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q9 : num 98 4 3 2 2 1 3 3 2 3 ...
## ..- attr(*, "label")= chr "請問您認為市政府決定重大政策時，會不會把「民眾的福利」放在第一優先考慮的地位？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 4
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "經常會" "有時會" "不太會" "絕對不會" ...
## $ Q10 : num 98 1 98 95 2 98 3 3 2 98 ...
## ..- attr(*, "label")= chr "整體來說，請問您認為我們高雄市的公務人員貪污的情形嚴不嚴重？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常嚴重" "嚴重" "不嚴重" "非常不嚴重" ...
## $ Q11 : num 2 2 2 1 2 2 2 1 2 3 ...
## ..- attr(*, "label")= chr "在這次高雄市長罷免案中，請問您有沒有參與連署？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:7] 1 2 3 4 5 95 98
## .. ..- attr(*, "names")= chr [1:7] "有" "沒有" "只參加第一階段" "只參加第二階段" ...
## $ Q12 : num 98 1 2 3 2 1 1 3 3 3 ...
## ..- attr(*, "label")= chr "請問您認為韓國瑜擔任高雄市長的表現與前任的陳菊市長相比，是比較好、比較差還是差不多？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:7] 1 2 3 95 96 97 98
## .. ..- attr(*, "names")= chr [1:7] "比較好" "差不多" "比較差" "拒答" ...
## $ Q13 : num 98 4 2 1 2 4 4 1 2 1 ...
## ..- attr(*, "label")= chr "有人說：「高雄市長韓國瑜講話常常引起爭議」。請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q14 : num 98 1 2 4 2 1 1 3 2 4 ...
## ..- attr(*, "label")= chr "有人說：「韓國瑜就任高雄市長後，很認真替高雄市打拼」。請問您同不同意個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q15 : num 98 4 2 1 3 4 3 1 1 1 ...
## ..- attr(*, "label")= chr "有人說：「韓國瑜的高雄市長任期未滿一年就請假參選總統，違背對高雄市民的承諾」。請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q16 : num 98 1 3 2 2 2 1 3 3 4 ...
## ..- attr(*, "label")= chr "有人說：「這次韓國瑜被罷免，主要原因是政黨惡鬥（臺：鬥爭）」。請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常同意" "同意" "不同意" "非常不同意" ...
## $ Q17 : num 2 2 1 1 2 2 2 1 2 1 ...
## ..- attr(*, "label")= chr "請問6月6日高雄市長罷免案，您有沒有去投票？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:4] 1 2 92 95
## .. ..- attr(*, "names")= chr [1:4] "有" "沒有" "無投票權" "拒答"
## $ Q18 : num 99 99 1 1 99 99 99 95 99 1 ...
## ..- attr(*, "label")= chr "請問您是投同意票還是不同意票？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:6] 1 2 91 94 95 99
## .. ..- attr(*, "names")= chr [1:6] "投同意票" "投不同意票" "忘記了" "投廢票" ...
## $ Q19 : num 99 99 3 4 99 99 99 3 99 4 ...
## ..- attr(*, "label")= chr "請問您是什麼時候決定要去投票，是一年以前、半年以前、三個月以前還是一週以前？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:9] 1 2 3 4 95 96 97 98 99
## .. ..- attr(*, "names")= chr [1:9] "一年以前" "半年以前" "三個月以前" "一週以前" ...
## $ Q20 : num 7 1 99 99 1 3 3 99 11 99 ...
## ..- attr(*, "label")= chr "請問您沒有去投票的原因是什麼？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:21] 1 2 3 4 5 6 7 8 9 10 ...
## .. ..- attr(*, "names")= chr [1:21] "不支持罷韓" "避免感染新冠肺炎" "支持韓國瑜" "上班、工作" ...
## $ Q21 : num 98 92 1 3 92 92 92 1 92 92 ...
## ..- attr(*, "label")= chr "在這次高雄市市長罷免案中，誰對您的投票決定影響最大？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:22] 1 2 3 4 5 6 7 8 9 10 ...
## .. ..- attr(*, "names")= chr [1:22] "韓國瑜本人" "罷韓團體" "家人" "親戚" ...
## $ Q22 : num 99 99 1 1 99 99 99 95 99 99 ...
## ..- attr(*, "label")= chr "請問〔021〕是影響您去投同意票、要去投反對票，還是不要去投票？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 95 99
## .. ..- attr(*, "names")= chr [1:5] "投同意票" "投反對票" "不要去投票" "拒答" ...
## $ Q23 : num 98 4 3 1 4 4 4 1 4 3 ...
## ..- attr(*, "label")= chr "請問您會不會擔心（臺：咁唔擔心）別人知道您有沒去投票？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 95 96 97 98
## .. ..- attr(*, "names")= chr [1:8] "非常擔心" "擔心" "不擔心" "非常不擔心" ...
## $ Q24 : num 93 1 93 2 4 92 1 95 1 2 ...
## ..- attr(*, "label")= chr "請問前年（2018）11月高雄市長選舉時，您把票投給哪一位候選人？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:10] 1 2 3 4 89 91 92 93 94 95
## .. ..- attr(*, "names")= chr [1:10] "韓國瑜" "陳其邁" "璩美鳳" "蘇盈貴" ...
## $ Q25 : num 95 93 92 2 2 92 2 95 3 2 ...
## ..- attr(*, "label")= chr "請問2014年11月高雄市長選舉時，您把票投給哪一位候選人？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:9] 1 2 3 89 91 92 93 94 95
## .. ..- attr(*, "names")= chr [1:9] "周可盛" "陳菊" "楊秋興" "沒領票" ...
## $ Q26 : num 95 2 93 3 2 2 2 95 94 3 ...
## ..- attr(*, "label")= chr "請問2020年1月的總統選舉中，您把票投給哪一組候選人？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:9] 1 2 3 89 91 92 93 94 95
## .. ..- attr(*, "names")= chr [1:9] "宋楚瑜與余湘" "韓國瑜與張善政" "蔡英文與賴清德" "沒領票" ...
## $ Q27 : num 98 2 6 5 3 6 5 5 6 6 ...
## ..- attr(*, "label")= chr "關於臺灣和大陸的關係，有下面幾種不同的看法："
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:12] 1 2 3 4 5 6 7 8 95 96 ...
## .. ..- attr(*, "names")= chr [1:12] "儘快統一" "儘快宣布獨立" "維持現狀，以後走向統一" "維持現狀，以後走向獨立" ...
## $ Q28 : num 7 38 7 2 1 1 38 7 7 2 ...
## ..- attr(*, "label")= chr "在國內的政黨之中，請問您認為您比較支持哪一個政黨？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:21] 1 2 3 4 5 6 7 16 18 19 ...
## .. ..- attr(*, "names")= chr [1:21] "國民黨" "民進黨" "新黨" "親民黨" ...
## $ Q29 : num 99 2 99 2 2 2 2 99 99 1 ...
## ..- attr(*, "label")= chr "請問您支持〔028〕的程度是非常支持、還是普普通通？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:3] 1 2 99
## .. ..- attr(*, "names")= chr [1:3] "非常支持" "普普通通" "跳題"
## $ Q30 : num 8 99 8 99 99 99 99 96 8 99 ...
## ..- attr(*, "label")= chr "那請問您有沒有比較偏向哪一個政黨？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:22] 3 4 5 6 7 8 18 19 20 21 ...
## .. ..- attr(*, "names")= chr [1:22] "偏國民黨" "偏民進黨" "偏新黨" "偏親民黨" ...
## $ Q31 : num 98 1 1 2 2 2 2 1 2 1 ...
## ..- attr(*, "label")= chr "我們社會上，有人說自己是「臺灣人」，也有人說自己是「中國人」，也有人說都是。請問您認為自己是「臺灣人」、「中國人」，或者都是？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:7] 1 2 3 95 96 97 98
## .. ..- attr(*, "names")= chr [1:7] "臺灣人" "都是" "中國人" "拒答" ...
## $ Q32 : num 95 80 87 78 57 54 49 95 48 30 ...
## ..- attr(*, "label")= chr "請問您是民國那一年出生的？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num 95
## .. ..- attr(*, "names")= chr "拒答"
## $ Q33 : num 95 6 5 4 7 3 5 6 4 3 ...
## ..- attr(*, "label")= chr "請問您的最高學歷是什麼？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:8] 1 2 3 4 5 6 7 95
## .. ..- attr(*, "names")= chr [1:8] "不識字及未入學" "小學" "國、初中" "高中、職" ...
## $ Q34 : num 98 2 2 2 3 95 1 2 2 2 ...
## ..- attr(*, "label")= chr "請問您的父親是本省客家人、本省閩南人、大陸各省市人、原住民，還是新住民？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:9] 1 2 3 4 5 6 7 95 98
## .. ..- attr(*, "names")= chr [1:9] "本省客家人" "本省閩南人" "大陸各省市人" "原住民" ...
## $ Q35 : num 98 2 2 2 2 95 1 2 2 2 ...
## ..- attr(*, "label")= chr "請問您的母親是本省客家人、本省閩南人、大陸各省市人、原住民，還是新住民？"
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:9] 1 2 3 4 5 6 7 95 98
## .. ..- attr(*, "names")= chr [1:9] "本省客家人" "本省閩南人" "大陸各省市人" "原住民" ...
## $ Q36 : num 911 302 801 905 103 401 105 602 901 906 ...
## ..- attr(*, "label")= chr "請問您的職業是什麼？"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:39] 101 102 103 104 105 106 201 202 203 204 ...
## .. ..- attr(*, "names")= chr [1:39] "民意代表" "政府行政主管" "公營事業主管" "民營事業主管" ...
## $ Q37 : num 999 999 999 999 999 999 999 999 104 999 ...
## ..- attr(*, "label")= chr "請問您配偶的職業是什麼？"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:36] 101 102 103 104 105 106 201 202 203 204 ...
## .. ..- attr(*, "names")= chr [1:36] "民代" "政府行政主管" "公營事業主管" "民營事業主管" ...
## $ Q38 : num 999 999 999 602 999 999 999 999 999 601 ...
## ..- attr(*, "label")= chr "請問您失業或退休前的職業是什麼？"
## ..- attr(*, "format.spss")= chr "F3.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:33] 101 102 103 104 105 106 201 202 203 204 ...
## .. ..- attr(*, "names")= chr [1:33] "民代" "政府行政主管" "公營事業主管" "民營事業主管" ...
## $ Q39 : num 2 1 1 2 1 1 1 2 2 1 ...
## ..- attr(*, "label")= chr "性別："
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:2] 1 2
## .. ..- attr(*, "names")= chr [1:2] "男性" "女性"
## $ Q40 : num 1 1 1 1 1 1 1 1 1 1 ...
## ..- attr(*, "label")= chr "訪問時使用語言："
## ..- attr(*, "format.spss")= chr "F2.0"
## ..- attr(*, "display_width")= int 5
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 5
## .. ..- attr(*, "names")= chr [1:5] "國語" "臺語" "客語" "國、臺語" ...
## $ NQ1 : num 14 21 11 9 12 5 4 19 2 5 ...
## ..- attr(*, "label")= chr "首先，想請問您的戶籍是在高雄市_________區"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:39] 1 2 3 4 5 6 7 8 9 10 ...
## .. ..- attr(*, "names")= chr [1:39] "鹽埕區" "鼓山區" "左營區" "楠梓區" ...
## $ NQ2 : num 92 92 99 92 1 8 5 6 92 4 ...
## ..- attr(*, "label")= chr "請問您平常最常看哪一臺的電視新聞？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:19] 1 2 3 4 5 6 7 8 9 10 ...
## .. ..- attr(*, "names")= chr [1:19] "臺視" "中視" "華視" "民視" ...
## $ NQ3 : num 9 2 9 3 2 1 1 3 2 4 ...
## ..- attr(*, "label")= chr "請問您對韓國瑜擔任高雄市長以來的整體施政表現滿不滿意？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常滿意" "滿意" "不滿意" "非常不滿意" ...
## $ NQ4 : num 9 2 3 2 3 9 1 2 2 1 ...
## ..- attr(*, "label")= chr "有人說：「政治人物當選後，很快就會跟民意脫節。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## $ NQ5 : num 9 1 3 3 3 2 4 2 3 1 ...
## ..- attr(*, "label")= chr "有人說：「政治人物並不會關心像我這樣一般百姓的想法。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## $ NQ6 : num 9 4 3 2 2 9 4 3 2 4 ...
## ..- attr(*, "label")= chr "有人說：「政治人物會用心為一般百姓做事。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## $ NQ7 : num 9 2 2 2 2 1 1 3 3 4 ...
## ..- attr(*, "label")= chr "有人說：「市政府所做的事大多數是正確的。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## $ NQ8 : num 9 1 2 2 2 4 4 2 3 1 ...
## ..- attr(*, "label")= chr "有人說：「市政府官員時常浪費一般民眾所繳納的稅金。」請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## $ NQ9 : num 9 4 3 2 2 1 3 3 2 3 ...
## ..- attr(*, "label")= chr "請問您認為市政府決定重大政策時，會不會把「民眾的福利」放在第一優先考慮的地位？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "經常會" "有時會" "不太會" "絕對不會" ...
## $ NQ10 : num 9 1 9 9 2 9 3 3 2 9 ...
## ..- attr(*, "label")= chr "整體來說，請問您認為我們高雄市的公務人員貪污的情形嚴不嚴重？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常嚴重" "嚴重" "不嚴重" "非常不嚴重" ...
## $ NQ11 : num 2 2 2 1 2 2 2 1 2 3 ...
## ..- attr(*, "label")= chr "在這次高雄市長罷免案中，請問您有沒有參與連署？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:6] 1 2 3 4 5 9
## .. ..- attr(*, "names")= chr [1:6] "有" "沒有" "只參加第一階段" "只參加第二階段" ...
## $ NQ12 : num 9 1 2 3 2 1 1 3 3 3 ...
## ..- attr(*, "label")= chr "請問您認為韓國瑜擔任高雄市長的表現與前任的陳菊市長相比，是比較好、比較差還是差不多？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:4] 1 2 3 9
## .. ..- attr(*, "names")= chr [1:4] "比較好" "差不多" "比較差" "無反應"
## $ NQ13 : num 9 4 2 1 2 4 4 1 2 1 ...
## ..- attr(*, "label")= chr "有人說：「高雄市長韓國瑜講話常常引起爭議」。請問您同不同意這個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## $ NQ14 : num 9 1 2 4 2 1 1 3 2 4 ...
## ..- attr(*, "label")= chr "有人說：「韓國瑜就任高雄市長後，很認真替高雄市打拼」。請問您同不同意個說法？"
## ..- attr(*, "format.spss")= chr "F8.0"
## ..- attr(*, "display_width")= int 10
## ..- attr(*, "labels")= Named num [1:5] 1 2 3 4 9
## .. ..- attr(*, "names")= chr [1:5] "非常同意" "同意" "不同意" "非常不同意" ...
## [list output truncated]批次清理無效值 課本p97
library(sjmisc)##
## Attaching package: 'sjmisc'## The following objects are masked from 'package:sjlabelled':
##
## to_character, to_factor, to_label, to_numericED206462 <-set_na(ED206462,na=c(91:99,"NA"))要使用的變數編碼(自變數)
120、職業 (1)高、中級白領 (2)中低、低級白領 (3)農林漁牧 (4)藍領 (5)其他
table(ED206462$CAREER)##
## 1 2 3 4 5
## 430 240 41 340 55ED206462$CAREER<-rec(ED206462$CAREER, rec="1,2,5=0[Others];
3,4=1[PI worker];", as.num = F)
library(sjPlot)
plot_frq(ED206462$CAREER)
47、請問您平常最常看哪一臺的電視新聞？(自變數)
table(ED206462$Q2)##
## 1 2 3 4 5 6 7 8 9 10 11 12 13 14 16 18
## 29 36 11 112 178 123 72 124 48 2 7 7 18 1 50 7ED206462$Q2 <-rec(ED206462$Q2, rec ="1,2,3,5,7,8=1[blue];
4,6,9,10,11,12,13,14,15,16=0[others];",
as.num = F)
library(sjPlot)
plot_frq(ED206462$Q2)
63、投同意票還是不同意票(依變數) (1)同意 (2)不同意
table(ED206462$Q18)##
## 1 2
## 431 16ED206462$Q18 <-rec(ED206462$Q18,rec = "1=1[Yes];2=2[No]", as.num = F)
library(sjPlot)
plot_frq(ED206462$Q18)
將包含標籤的資料檔存為rda格式
save(ED206462, file="ED206462_raw.rda", compress=TRUE)用製表的方式描述變數
load("ED206462_raw.rda")sjmisc::frq() 有frq的指令要加上sjmisc
library(sjmisc)
frq(ED206462$CAREER)##
## 職業五分類 (x) <categorical>
## # total N=1106 valid N=1106 mean=0.34 sd=0.48
##
## Value | Label | N | Raw % | Valid % | Cum. %
## --------------------------------------------------
## 0 | Others | 725 | 65.55 | 65.55 | 65.55
## 1 | PI worker | 381 | 34.45 | 34.45 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>frq(ED206462$CAREER, weights = ED206462$w) ##
## 職業五分類 (xw) <categorical>
## # total N=1106 valid N=1106 mean=0.38 sd=0.48
##
## Value | Label | N | Raw % | Valid % | Cum. %
## --------------------------------------------------
## 0 | Others | 689 | 62.30 | 62.30 | 62.30
## 1 | PI worker | 417 | 37.70 | 37.70 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>frq(ED206462$Q2)##
## 請問您平常最常看哪一臺的電視新聞？ (x) <categorical>
## # total N=1106 valid N=818 mean=0.55 sd=0.50
##
## Value | Label | N | Raw % | Valid % | Cum. %
## -----------------------------------------------
## 0 | others | 368 | 33.27 | 44.99 | 44.99
## 1 | blue | 450 | 40.69 | 55.01 | 100.00
## <NA> | <NA> | 288 | 26.04 | <NA> | <NA>frq(ED206462$Q2, weights = ED206462$w)##
## 請問您平常最常看哪一臺的電視新聞？ (xw) <categorical>
## # total N=794 valid N=794 mean=0.53 sd=0.50
##
## Value | Label | N | Raw % | Valid % | Cum. %
## -----------------------------------------------
## 0 | others | 371 | 46.73 | 46.73 | 46.73
## 1 | blue | 423 | 53.27 | 53.27 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>frq(ED206462$Q18)##
## 請問您是投同意票還是不同意票？ (x) <categorical>
## # total N=1106 valid N=447 mean=1.04 sd=0.19
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 1 | Yes | 431 | 38.97 | 96.42 | 96.42
## 2 | No | 16 | 1.45 | 3.58 | 100.00
## <NA> | <NA> | 659 | 59.58 | <NA> | <NA>frq(ED206462$Q18,weights = ED206462$w)##
## 請問您是投同意票還是不同意票？ (xw) <categorical>
## # total N=470 valid N=470 mean=1.03 sd=0.17
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 1 | Yes | 456 | 97.02 | 97.02 | 97.02
## 2 | No | 14 | 2.98 | 2.98 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>使用sjPlot來視覺化變數之間的關係
使用plot_grpfrq()
library(sjPlot)
library(sjlabelled)
library(sjmisc)
load("ED206462_raw.rda")
frq(ED206462$CAREER)##
## 職業五分類 (x) <categorical>
## # total N=1106 valid N=1106 mean=0.34 sd=0.48
##
## Value | Label | N | Raw % | Valid % | Cum. %
## --------------------------------------------------
## 0 | Others | 725 | 65.55 | 65.55 | 65.55
## 1 | PI worker | 381 | 34.45 | 34.45 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>ED206462$CAREER <- set_label(ED206462$CAREER, label="career") # 重設變數標籤
ED206462$CAREER <- set_labels(ED206462$CAREER, labels= c("PI worker","other")) #重設選項標籤
frq(ED206462$Q2)##
## 請問您平常最常看哪一臺的電視新聞？ (x) <categorical>
## # total N=1106 valid N=818 mean=0.55 sd=0.50
##
## Value | Label | N | Raw % | Valid % | Cum. %
## -----------------------------------------------
## 0 | others | 368 | 33.27 | 44.99 | 44.99
## 1 | blue | 450 | 40.69 | 55.01 | 100.00
## <NA> | <NA> | 288 | 26.04 | <NA> | <NA>ED206462$Q2 <- set_label(ED206462$Q2, label = "TV station")
ED206462$Q2 <- set_labels(ED206462$Q2,labels =c("Blues","other"))
frq(ED206462$Q18)##
## 請問您是投同意票還是不同意票？ (x) <categorical>
## # total N=1106 valid N=447 mean=1.04 sd=0.19
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 1 | Yes | 431 | 38.97 | 96.42 | 96.42
## 2 | No | 16 | 1.45 | 3.58 | 100.00
## <NA> | <NA> | 659 | 59.58 | <NA> | <NA>ED206462$Q18 <- set_label(ED206462$Q18, label = "vote")
ED206462$Q18 <- set_labels(ED206462$Q18,labels = c("Yes","No"))製圖
frq(ED206462$CAREER)##
## career (x) <categorical>
## # total N=1106 valid N=1106 mean=0.34 sd=0.48
##
## Value | Label | N | Raw % | Valid % | Cum. %
## --------------------------------------------------
## 0 | PI worker | 725 | 65.55 | 65.55 | 65.55
## 1 | other | 381 | 34.45 | 34.45 | 100.00
## <NA> | <NA> | 0 | 0.00 | <NA> | <NA>frq(ED206462$Q18)##
## vote (x) <categorical>
## # total N=1106 valid N=447 mean=1.04 sd=0.19
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 1 | Yes | 431 | 38.97 | 96.42 | 96.42
## 2 | No | 16 | 1.45 | 3.58 | 100.00
## <NA> | <NA> | 659 | 59.58 | <NA> | <NA>plot_grpfrq(ED206462$CAREER, ED206462$Q18) # 橫向比較
frq(ED206462$Q2)##
## TV station (x) <categorical>
## # total N=1106 valid N=818 mean=0.55 sd=0.50
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 0 | Blues | 368 | 33.27 | 44.99 | 44.99
## 1 | other | 450 | 40.69 | 55.01 | 100.00
## <NA> | <NA> | 288 | 26.04 | <NA> | <NA>frq(ED206462$Q18)##
## vote (x) <categorical>
## # total N=1106 valid N=447 mean=1.04 sd=0.19
##
## Value | Label | N | Raw % | Valid % | Cum. %
## ----------------------------------------------
## 1 | Yes | 431 | 38.97 | 96.42 | 96.42
## 2 | No | 16 | 1.45 | 3.58 | 100.00
## <NA> | <NA> | 659 | 59.58 | <NA> | <NA>plot_grpfrq(ED206462$Q2, ED206462$Q18)
卡方檢定
library(sjPlot)
tab_xtab(ED206462$CAREER, ED206462$Q18,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 = "gray", # 將欄百分比顏色改為灰色 （預設為綠色）
tdcol.row = "brown" # 將列百分比顏色改為褐色 （預設為藍色）
) | career | vote | Total | |
|---|---|---|---|
| Yes | No | ||
| PI worker |
270 95.7 % 62.6 % |
12 4.3 % 75 % |
282 100 % 63.1 % |
| other |
161 97.6 % 37.4 % |
4 2.4 % 25 % |
165 100 % 36.9 % |
| Total |
431 96.4 % 100 % |
16 3.6 % 100 % |
447 100 % 100 % |
χ2=0.550 · df=1 · φ=0.048 · Fisher’s p=0.431 |
library(sjPlot)
tab_xtab(ED206462$Q2, ED206462$Q18,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 = "gray", # 將欄百分比顏色改為灰色 （預設為綠色）
tdcol.row = "brown" # 將列百分比顏色改為褐色 （預設為藍色）
) | TV station | vote | Total | |
|---|---|---|---|
| Yes | No | ||
| Blues |
256 98.5 % 73.8 % |
4 1.5 % 30.8 % |
260 100 % 72.2 % |
| other |
91 91 % 26.2 % |
9 9 % 69.2 % |
100 100 % 27.8 % |
| Total |
347 96.4 % 100 % |
13 3.6 % 100 % |
360 100 % 100 % |
χ2=9.508 · df=1 · φ=0.179 · Fisher’s p=0.002 |
模型一
library(car)## Loading required package: carData## Registered S3 methods overwritten by 'car':
## method from
## influence.merMod lme4
## cooks.distance.influence.merMod lme4
## dfbeta.influence.merMod lme4
## dfbetas.influence.merMod lme4model.1 <- glm(Q18 ~ Q2+CAREER,
data=ED206462,
family=binomial)
summary(model.1) ##
## Call:
## glm(formula = Q18 ~ Q2 + CAREER, family = binomial, data = ED206462)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -0.4891 -0.2918 -0.2021 -0.1186 3.1497
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) -3.8814 0.5216 -7.442 9.92e-14 ***
## Q21 1.8183 0.6151 2.956 0.00312 **
## CAREER1 -1.0719 0.7842 -1.367 0.17166
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 111.87 on 359 degrees of freedom
## Residual deviance: 99.54 on 357 degrees of freedom
## (746 observations deleted due to missingness)
## AIC: 105.54
##
## Number of Fisher Scoring iterations: 7vif(model.1)## Q2 CAREER
## 1.000168 1.000168簡易描述
從此次分析來看，支持罷免韓國瑜的受訪者與其職業關係，兩者相關性的檢驗，p值大於.05統計上並無顯著性，不拒絕假設一的虛無假設，也就是說，支持罷免韓國瑜的人其職業為非農工階級。
支持罷免韓國瑜的受訪者，與其平常收看不偏藍的電視新聞台，兩者相關性的檢驗，統計上的顯著性不及0.05的標準，統計上具有顯著性，拒絕假設二的虛無假設。因此，支持罷免韓國瑜的受訪者，不收看不偏藍的電視新聞台。