-
Notifications
You must be signed in to change notification settings - Fork 0
/
Burt_1.R
114 lines (87 loc) · 3.03 KB
/
Burt_1.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
#Burt? MCA_day2_new_JH.R - esimerkki - aika sähläystä vielä...
# (27.9.18) -> r-tiedosto Burt_1.R
# Jatketaan 20.9.20
# maaga - muuttuja luodaan G1_4_laaj2.Rmd - skriptissä
#
# ika liian iso luokitus, jätetään pois
# myös maaga turha, sille haetaan vaihtoehtoa
#
# Burt-idea: ensin matriisi, sitten voi pinota yms. "oikein"
test1 <- select(ISSP2012esim1b.dat, Q1b, maa,sp, age_cat)
str(test1)
maagaBurt <- mjca(test1, ps="")$Burt
dim(maagaBurt) #19 19
str(maagaBurt)
# Tarkemmin
typeof(maagaBurt) # what is it?
class(maagaBurt) # what is it? (sorry) - on matriisi (21.9.20)
storage.mode(maagaBurt) # what is it? (very sorry)
length(maagaBurt) # how long is it? What about two dimensional objects?
attributes(maagaBurt)
maagaBurt
head(maagaBurt)
rownames(maagaBurt)
# Burt-matriisi symmetrinen -> rownames()
# Q1bS Q1bs Q1b? Q1be Q1bE maaBE maaBG maaDE maaDK maaFI maaHU spm spf
# age_cat1 age_cat2 age_cat3 age_cat4 age_cat5 age_cat6
# vastaukset 1:5, maat 6:11, sp 12:13, age_cat 14:19
# Rivit
maagaBurt[1,]
maagaBurt2 <- maagaBurt[6:11,1:5]
str(maagaBurt2)
glimpse(maagaBurt2)
maagaBurt2
# sukupuoli ja vastaukset
maagaBurt2 <- rbind(maagaBurt2, maagaBurt[12:13 ,1:5])
# ikäluokka ja vastaukset
maagaBurt2 <- rbind(maagaBurt2, maagaBurt[14:19 ,1:5])
# Toimii!
maagaBurt2
maagaBurt2CA1 <- ca(maagaBurt2)
plot(maagaBurt2CA1)
# Siistitään, ja käännetään y-akselin ympäri (21.9.20)
#
# Akseleiden kääntö (MG2017 harjoitukset)
#
# Ensimmäinen akseli:
# women.Burt.ca$rowcoord[,1] <- -women.Burt.ca$rowcoord[,1]
# women.Burt.ca$colcoord[,1] <- -women.Burt.ca$colcoord[,1]
# str(maagaBurt2CA1)
maagaBurt2CA1$colnames <- c("S", "s", "?", "e", "E")
maagaBurt2CA1$colnames # uudet sarakenimet
maagaBurt2CA1$rownames <- c("BE", "BG", "DE", "DK", "FI", "HU", "m", "f",
"a1", "a2", "a3", "a4", "a5", "a6")
maagaBurt2CA1$rownames
# käännetään kuva y-akselin ympäri
# testOld.rpc <- maagaBurt2CA1$rowcoord
maagaBurt2CA1$rowcoord[, 2] <- -maagaBurt2CA1$rowcoord[, 2]
maagaBurt2CA1$colcoord[, 2] <- -maagaBurt2CA1$colcoord[, 2]
# Toimii (17.10.20)
#testOld.rpc
#maagaBurt2CA1$rowcoord
#X11()
# plot(maagaBurt2CA1, map = "rowgreen") - ei toimi, rivipisteet läjässä) (21.9.20)
plot(maagaBurt2CA1)
# ei ihan helppoa...vanhaa ja ehkä tarpeetonta tavaraa (21.9.20)
testTab1 <- table(test1$maaga, test1$V6)
testTab1
testTab_sr <- table(test1$maa, test1$V6)
testTab_sr
testTab1 <- rbind(testTab1,testTab_sr)
testTab1
dim(testTab1) #78 riviä, 5 saraketta, 1-72 data ja 73-78 täydentävät rivit
spCAmaaga1 <- ca(testTab1[,1:5], suprow = 73:78)
X11()
par("cex"= 0.75, "asp" = 1)
plot(spCAmaaga1, main = "Äiti töissä: ikäluokka ja sukupuoli maittain 2")
#asymmetrinen kartta
X11()
par("cex"= 0.75, "asp" = 1)
plot(spCAmaaga1, main = "Äiti töissä: ikäluokka ja sukupuoli maittain 3 (kontribuutiot)",
map = "rowgreen",
contrib= c("absolute", "absolute"),
mass = c(TRUE,TRUE),
arrows = c(FALSE,TRUE)
)
#numeeriset tulokset
summary(spCAmaaga1)