Rで何かをしたり、読書をするブログ

政府統計の総合窓口のデータや、OECDやUCIやのデータを使って、Rの練習をしています。ときどき、読書記録も載せています。

内閣支持率のデータ分析1 - 各社で支持率は違うか、支持率の信頼区間をR言で調べる。

安倍晋三首相の在任日数が憲政史上最長になりました。

そこで今回は、各報道機関の世論調査、内閣支持率がどうなっているかを調べてみました。

f:id:cross_hyou:20191120192542p:plain

こちらは、NHKの世論調査です。

http://www.nhk.or.jp/senkyo/shijiritsu/

 

2019年11月8日から10日の調査で、1253人のうち、支持するが47%, 支持しないが35%です。

f:id:cross_hyou:20191120193219p:plain

こちらはFNNの調査です。

https://www.fnn.jp/posts/00427442CX/201911181205_CX_CX

11月16日と17日に1039人に調査して、支持が45.1%, 不支持が37.7%でした。

f:id:cross_hyou:20191120193906p:plain

f:id:cross_hyou:20191120193951p:plain

こちらは、朝日新聞の世論調査です。

https://www.asahi.com/articles/ASMCK5RS3MCKUZPS009.html

11月の16日と17日に1998人に調査して、支持が44%、支持しないが36%でした。

f:id:cross_hyou:20191120194601p:plain

f:id:cross_hyou:20191120194618p:plain

こちらはTBSの世論調査です。

https://news.tbs.co.jp/newsi_sp/yoron/backnumber/20191109/q1-1.html

11月の9日と10日に調査して、1139人のうち支持が54.3%, 不支持が42.4%でした。

f:id:cross_hyou:20191120195830p:plain

f:id:cross_hyou:20191120195805p:plain

これは日本テレビの調査です。

https://www.ntv.co.jp/yoron/

11月15日から17日の調査で、1051人の回答で、支持が49%、支持しないが36%でした。

これらをまとめて表にしてみました。

f:id:cross_hyou:20191120200630p:plain

このデータをR言語で読込みます。read.csv関数です。

f:id:cross_hyou:20191120200828p:plain

一番支持率が高かったのがTBS、低かったのがAsahiでしたので、この二つの会社のクロス表を作ってカイ二乗検定をします。

f:id:cross_hyou:20191120201428p:plain

まずは、このようにAsahiとTBSのマトリックスを作りました。

これをカイ二乗検定します。chisq.test関数です。

f:id:cross_hyou:20191120201716p:plain

p-valueが0.5637ですから、AsahiとTBSの世論調査に有意な違いは無い、ということですね。支持率の一番高いところと低いところに有意な差が無いということですから、各社で支持率が違うということは無いです。

 

それでは、支持率はどのくらいの範囲にあるのでしょうか?Binomial Testで調べてみましょう。関数は、binom.testです。

回答者の合計、支持者の合計を出します。

apply関数を使って合計を計算します。

f:id:cross_hyou:20191120203013p:plain

支持者の合計が3070人、回答者の合計が6480人です。この二つの数字をbinom.test関数にあてはめて支持率の信頼区間を計算します。

f:id:cross_hyou:20191120203257p:plain

95 percent confidence intervalの下の値を見ればいいです。

95%の信頼区間で46.2%から48.6%の間に支持率があります。

今回は以上です。