fixed R script, updated README, polls.dat now has 1 row per survey

Author: mitzimorris

ARM/Ch.2/2.3_CI.R

@@ -41,8 +41,12 @@ m + geom_point() + scale_y_continuous("Percentage Support for the Death Penalty"
 
 # Weighted averages
 
+N <- c(65633200,80523700,59685200) # population sizes FR, DE, IT
+p <- c(0.55,0.61,0.38)  # estimated proportions of Yes responses
+se <- c(0.02,0.03,0.03)
+
 w.avg <- sum(N*p)/sum(N)
-se.w.av <- sqrt (sum ((N*se/sum(N))^2))
+se.w.avg <- sqrt (sum ((N*se/sum(N))^2))
 int.95 <- w.avg + c(-2,2)*se.w.avg
 
 # CI using simulations

ARM/Ch.2/README

@@ -1,7 +1,13 @@
 Data
 ====
 
+1 polls.dat - Gallup poll data support for death penalty by year - fig 2.3
+  5 columns per row
+  c1 : year
+  c2 : month
+  c3 : supports death penalty
+  c4 : opposed to death penalty
+  c5 : no opinion
 
 Models
 ======
-

ARM/Ch.2/polls.dat

@@ -1,128 +1,32 @@
-2002 10
- 70
- 25
- 05
-2002 5
- 72
- 25
- 03
-2001 10
- 68
- 26
- 06
-2001 5
- 65
- 27
- 08
-2001 2
- 67
- 25
- 08
-2000 8
- 67
- 28
- 05
-2000 6
- 66
- 26
- 08
-2000 2
- 66
- 28
- 06
-1999 5
- 71
- 22
- 07
-1995 9
- 77
- 13
- 10
-1994 9
- 80
- 16
- 04
-1991 6
- 76
- 18
- 06
-1988 9.5
- 79
- 16
- 05
-1988 9
- 79
- 16
- 05
-1986 1
- 70
- 22
- 08
-1985 9
- 75
- 17
- 08
-1985 1
- 72
- 20
- 08
-1981 1.5
- 66
- 25
- 09
-1978 3
- 62
- 27
- 11
-1976 4
- 66
- 26
- 08
-1972 11
- 57
- 32
- 11
-1972 3
- 50
- 41
- 09
-1971 10.5
- 49
- 40
- 11
-1969 1
- 51
- 40
- 09
-1967 6
- 54
- 38
- 08
-1966 5
- 42
- 47
- 11
-1965 1
- 45
- 43
- 12
-1960 3
- 53
- 36
- 11
-1957 8.5
- 47
- 34
- 18
-1956 3.5
- 53
- 34
- 13
-1953 9
- 68
- 25
- 07
-1937 12
- 60
- 33
- 07
+2002 10 70 25 05
+2002 5 72 25 03
+2001 10 68 26 06
+2001 5 65 27 08
+2001 2 67 25 08
+2000 8 67 28 05
+2000 6 66 26 08
+2000 2 66 28 06
+1999 5 71 22 07
+1995 9 77 13 10
+1994 9 80 16 04
+1991 6 76 18 06
+1988 9.5 79 16 05
+1988 9 79 16 05
+1986 1 70 22 08
+1985 9 75 17 08
+1985 1 72 20 08
+1981 1.5 66 25 09
+1978 3 62 27 11
+1976 4 66 26 08
+1972 11 57 32 11
+1972 3 50 41 09
+1971 10.5 49 40 11
+1969 1 51 40 09
+1967 6 54 38 08
+1966 5 42 47 11
+1965 1 45 43 12
+1960 3 53 36 11
+1957 8.5 47 34 18
+1956 3.5 53 34 13
+1953 9 68 25 07
+1937 12 60 33 07