syntax multi
syntax multi
%+%, alpha
Warning message:
package ‘psych’ was built under R version 4.4.2
> library(mirt)
> ######################Menampilkan ScreePlot ####
> datamulti<- scree(dat1,factors = T, pc = F, main = "Scree plot", hline = NULL,
add = F)
ListeningA <- rowSums(data_mod1[,1:30])
> ListeningB <- rowSums(data_mod1[,31:38])
> ListeningC <- rowSums(data_mod1[,39:50])
> GrammarA <- rowSums(data_mod1[,51:65])
> GrammarB <- rowSums(data_mod1[,66:90])
> read1 <- rowSums(data_mod1[,91:100])
> read2 <- rowSums(data_mod1[,101:110])
> read3 <- rowSums(data_mod1[,111:120])
> read4 <- rowSums(data_mod1[,121:130])
> read5 <- rowSums(data_mod1[,131:140])
> data1 <- cbind(ListeningA, ListeningB, ListeningC, GrammarA, GrammarB, read1,
read2, read3, read4, read5)
> colnames(data1) <- c('ListeningA', 'ListeningB', 'ListeningC', 'GrammarA',
'GrammarB', 'read1', 'read2', 'read3', 'read4', 'read5')
> data1
ListeningA ListeningB ListeningC GrammarA GrammarB read1 read2 read3 read4
read5
[1,] 1 0 0 0 2 7 3 3 4
4
[2,] 6 4 3 5 7 6 7 6 4
4
[3,] 14 3 4 6 7 7 7 4 5
6
[4,] 6 5 1 8 8 7 2 8 8
4
[5,] 6 1 5 6 7 7 5 5 2
3
[6,] 4 1 4 6 2 5 6 4 5
5
[7,] 7 4 3 3 4 6 3 4 3
1
[8,] 7 1 5 4 9 4 4 2 2
4
[9,] 18 6 10 11 6 8 4 6 3
2
[10,] 20 8 9 13 16 8 6 10 7
8
[11,] 2 6 7 5 6 6 3 6 1
6
[12,] 8 2 1 5 6 5 4 6 3
4
[13,] 10 3 2 5 8 6 5 6 4
4
[14,] 5 2 4 5 7 3 3 3 2
6
[15,] 8 3 5 6 5 6 4 5 2
8
[16,] 5 2 2 7 9 7 6 2 3
4
[17,] 9 0 2 4 7 3 3 3 2
2
[18,] 8 3 1 9 6 7 5 5 3
7
[19,] 19 6 3 10 12 4 7 9 6
7
[20,] 9 0 3 4 4 1 3 3 2
2
[21,] 5 4 4 5 5 8 2 3 3
4
[22,] 6 2 5 7 7 6 3 6 4
4
[23,] 9 2 5 4 9 8 7 4 5
6
[24,] 12 6 7 7 3 5 7 7 3
4
[25,] 4 2 2 5 7 3 4 7 2
4
[26,] 7 3 4 6 10 2 3 5 5
5
[27,] 5 2 5 6 9 3 2 0 3
4
[28,] 7 4 2 4 5 5 6 5 2
4
[29,] 4 3 4 7 6 3 5 3 2
2
[30,] 4 0 2 6 10 5 5 2 1
6
[31,] 11 1 3 7 6 4 4 5 6
4
[32,] 10 4 1 5 12 3 4 3 2
3
[33,] 10 2 5 7 9 5 3 7 4
4
[34,] 10 4 7 7 5 6 3 3 2
6
[35,] 11 2 8 9 15 9 6 6 5
5
[36,] 14 6 9 12 20 7 7 10 8
9
[37,] 13 5 6 8 8 6 3 7 5
6
[38,] 2 3 4 0 8 6 4 5 3
1
[39,] 7 2 1 3 8 1 1 4 4
7
[40,] 8 2 2 4 6 3 4 3 2
4
[41,] 15 6 7 12 15 7 6 7 7
8
[42,] 4 2 0 5 6 3 3 3 4
3
[43,] 6 4 6 1 6 5 4 4 0
4
[44,] 11 4 6 9 12 9 8 9 7
7
[45,] 18 4 7 7 6 6 6 6 7
5
[46,] 6 2 3 6 9 6 2 3 3
3
[47,] 8 7 4 8 7 9 5 5 4
7
[48,] 5 4 6 5 8 4 4 4 6
1
[49,] 6 4 4 5 6 4 3 3 3
2
[50,] 11 4 3 11 13 3 6 2 4
5
[51,] 4 3 2 6 8 2 6 4 2
2
[52,] 6 2 4 1 5 3 2 4 3
2
[53,] 9 3 3 6 7 8 3 4 4
5
[54,] 9 3 6 4 12 2 4 6 2
6
[55,] 5 3 3 2 8 0 2 3 1
6
[56,] 6 4 2 5 6 5 4 5 0
2
[57,] 3 3 0 3 8 2 2 5 2
5
[58,] 4 2 3 2 8 4 2 3 5
3
[59,] 8 7 5 1 8 5 5 6 5
5
[60,] 15 6 7 12 13 7 4 4 4
4
[61,] 6 6 3 7 6 3 4 4 4
5
[62,] 4 1 3 5 4 6 5 3 1
3
[63,] 8 3 6 2 10 1 4 4 3
3
[64,] 3 3 5 3 5 3 3 1 3
3
[65,] 2 1 4 4 5 6 2 4 0
4
[66,] 4 2 4 5 8 3 4 4 1
4
[67,] 4 3 1 4 7 4 3 3 3
2
[68,] 6 1 2 4 8 6 6 6 3
5
[69,] 10 2 2 5 6 3 4 4 1
1
[70,] 15 5 4 7 6 6 6 7 4
5
[71,] 5 2 3 5 6 1 1 1 1
2
[72,] 11 4 5 5 5 5 3 4 3
5
[73,] 8 1 7 6 9 2 4 3 2
2
[74,] 6 4 1 5 5 5 5 5 0
4
[75,] 15 4 8 11 11 5 4 8 6
5
[76,] 22 8 9 12 16 7 7 6 9
6
[77,] 5 1 3 3 6 1 3 3 2
2
[78,] 8 3 4 2 5 5 7 3 3
5
[79,] 24 7 8 9 13 9 9 9 7
9
[80,] 18 4 7 8 13 3 4 8 5
7
[81,] 22 7 8 12 15 7 6 8 5
8
[82,] 6 4 4 6 4 3 5 4 3
6
[83,] 7 2 5 11 15 8 5 4 5
8
[84,] 15 8 8 10 11 9 8 4 8
5
[85,] 6 3 1 4 8 5 2 6 1
2
[86,] 11 1 2 4 2 2 4 3 1
5
[87,] 24 8 9 13 21 7 8 8 9
9
[88,] 23 7 8 12 17 4 7 7 7
7
[89,] 20 7 8 14 21 7 8 8 8
7
[90,] 21 7 12 13 23 8 7 8 5
8
[91,] 28 8 10 15 24 8 8 6 8
9
[92,] 19 7 7 8 13 9 7 9 9
8
[93,] 29 7 12 15 21 9 8 9 9
9
[94,] 24 6 9 9 14 7 7 7 6
7
[95,] 5 2 4 10 5 5 4 3 1
5
[96,] 7 3 2 2 8 2 7 3 3
3
[97,] 14 5 7 6 10 9 4 4 5
4
[98,] 17 7 7 12 16 8 8 7 8
8
[99,] 6 0 3 7 9 4 3 4 4
2
[100,] 8 2 1 3 4 3 1 1 4
3
[ reached getOption("max.print") -- omitted 3209 rows ]
> #simpan data ke csv
> data.frame(data1)
ListeningA ListeningB ListeningC GrammarA GrammarB read1 read2 read3 read4
read5
1 1 0 0 0 2 7 3 3 4
4
2 6 4 3 5 7 6 7 6 4
4
3 14 3 4 6 7 7 7 4 5
6
4 6 5 1 8 8 7 2 8 8
4
5 6 1 5 6 7 7 5 5 2
3
6 4 1 4 6 2 5 6 4 5
5
7 7 4 3 3 4 6 3 4 3
1
8 7 1 5 4 9 4 4 2 2
4
9 18 6 10 11 6 8 4 6 3
2
10 20 8 9 13 16 8 6 10 7
8
11 2 6 7 5 6 6 3 6 1
6
12 8 2 1 5 6 5 4 6 3
4
13 10 3 2 5 8 6 5 6 4
4
14 5 2 4 5 7 3 3 3 2
6
15 8 3 5 6 5 6 4 5 2
8
16 5 2 2 7 9 7 6 2 3
4
17 9 0 2 4 7 3 3 3 2
2
18 8 3 1 9 6 7 5 5 3
7
19 19 6 3 10 12 4 7 9 6
7
20 9 0 3 4 4 1 3 3 2
2
21 5 4 4 5 5 8 2 3 3
4
22 6 2 5 7 7 6 3 6 4
4
23 9 2 5 4 9 8 7 4 5
6
24 12 6 7 7 3 5 7 7 3
4
25 4 2 2 5 7 3 4 7 2
4
26 7 3 4 6 10 2 3 5 5
5
27 5 2 5 6 9 3 2 0 3
4
28 7 4 2 4 5 5 6 5 2
4
29 4 3 4 7 6 3 5 3 2
2
30 4 0 2 6 10 5 5 2 1
6
31 11 1 3 7 6 4 4 5 6
4
32 10 4 1 5 12 3 4 3 2
3
33 10 2 5 7 9 5 3 7 4
4
34 10 4 7 7 5 6 3 3 2
6
35 11 2 8 9 15 9 6 6 5
5
36 14 6 9 12 20 7 7 10 8
9
37 13 5 6 8 8 6 3 7 5
6
38 2 3 4 0 8 6 4 5 3
1
39 7 2 1 3 8 1 1 4 4
7
40 8 2 2 4 6 3 4 3 2
4
41 15 6 7 12 15 7 6 7 7
8
42 4 2 0 5 6 3 3 3 4
3
43 6 4 6 1 6 5 4 4 0
4
44 11 4 6 9 12 9 8 9 7
7
45 18 4 7 7 6 6 6 6 7
5
46 6 2 3 6 9 6 2 3 3
3
47 8 7 4 8 7 9 5 5 4
7
48 5 4 6 5 8 4 4 4 6
1
49 6 4 4 5 6 4 3 3 3
2
50 11 4 3 11 13 3 6 2 4
5
51 4 3 2 6 8 2 6 4 2
2
52 6 2 4 1 5 3 2 4 3
2
53 9 3 3 6 7 8 3 4 4
5
54 9 3 6 4 12 2 4 6 2
6
55 5 3 3 2 8 0 2 3 1
6
56 6 4 2 5 6 5 4 5 0
2
57 3 3 0 3 8 2 2 5 2
5
58 4 2 3 2 8 4 2 3 5
3
59 8 7 5 1 8 5 5 6 5
5
60 15 6 7 12 13 7 4 4 4
4
61 6 6 3 7 6 3 4 4 4
5
62 4 1 3 5 4 6 5 3 1
3
63 8 3 6 2 10 1 4 4 3
3
64 3 3 5 3 5 3 3 1 3
3
65 2 1 4 4 5 6 2 4 0
4
66 4 2 4 5 8 3 4 4 1
4
67 4 3 1 4 7 4 3 3 3
2
68 6 1 2 4 8 6 6 6 3
5
69 10 2 2 5 6 3 4 4 1
1
70 15 5 4 7 6 6 6 7 4
5
71 5 2 3 5 6 1 1 1 1
2
72 11 4 5 5 5 5 3 4 3
5
73 8 1 7 6 9 2 4 3 2
2
74 6 4 1 5 5 5 5 5 0
4
75 15 4 8 11 11 5 4 8 6
5
76 22 8 9 12 16 7 7 6 9
6
77 5 1 3 3 6 1 3 3 2
2
78 8 3 4 2 5 5 7 3 3
5
79 24 7 8 9 13 9 9 9 7
9
80 18 4 7 8 13 3 4 8 5
7
81 22 7 8 12 15 7 6 8 5
8
82 6 4 4 6 4 3 5 4 3
6
83 7 2 5 11 15 8 5 4 5
8
84 15 8 8 10 11 9 8 4 8
5
85 6 3 1 4 8 5 2 6 1
2
86 11 1 2 4 2 2 4 3 1
5
87 24 8 9 13 21 7 8 8 9
9
88 23 7 8 12 17 4 7 7 7
7
89 20 7 8 14 21 7 8 8 8
7
90 21 7 12 13 23 8 7 8 5
8
91 28 8 10 15 24 8 8 6 8
9
92 19 7 7 8 13 9 7 9 9
8
93 29 7 12 15 21 9 8 9 9
9
94 24 6 9 9 14 7 7 7 6
7
95 5 2 4 10 5 5 4 3 1
5
96 7 3 2 2 8 2 7 3 3
3
97 14 5 7 6 10 9 4 4 5
4
98 17 7 7 12 16 8 8 7 8
8
99 6 0 3 7 9 4 3 4 4
2
100 8 2 1 3 4 3 1 1 4
3
[ reached 'max' / getOption("max.print") -- omitted 3209 rows ]
> fa.diagram(efa_result) # Visualisasi struktur faktor
> ###################### Analisis Kecocokan Model ####
> head(data1)
ListeningA ListeningB ListeningC GrammarA GrammarB read1 read2 read3 read4
read5
[1,] 1 0 0 0 2 7 3 3 4
4
[2,] 6 4 3 5 7 6 7 6 4
4
[3,] 14 3 4 6 7 7 7 4 5
6
[4,] 6 5 1 8 8 7 2 8 8
4
[5,] 6 1 5 6 7 7 5 5 2
3
[6,] 4 1 4 6 2 5 6 4 5
5
> #analisis data All:
> spec <-
+ "F1 = skor_B1, skor_B2, skor_B3, skor_B4, skor_B5, skor_B6, skor_B7, skor_B8,
skor_B9, skor_B10,
+ skor_B11, skor_B12, skor_B13, skor_B14, skor_B15, skor_B16, skor_B17, skor_B18,
skor_B19, skor_B20,
+ skor_B21, skor_B22, skor_B23, skor_B24, skor_B25, skor_B26, skor_B27, skor_B28,
skor_B29, skor_B30,
+ skor_B31, skor_B32, skor_B33, skor_B34, skor_B35, skor_B36, skor_B37, skor_B38,
skor_B39, skor_B40,
+ skor_B41, skor_B42, skor_B43, skor_B44, skor_B45, skor_B46, skor_B47, skor_B48,
skor_B49, skor_B50,
+ skor_B51, skor_B52, skor_B53, skor_B54, skor_B55, skor_B56, skor_B57, skor_B58,
skor_B59, skor_B60,
+ skor_B61, skor_B62, skor_B63, skor_B64, skor_B65, skor_B66, skor_B67, skor_B68,
skor_B69, skor_B70,
+ skor_B71, skor_B72, skor_B73, skor_B74, skor_B75, skor_B76, skor_B77, skor_B78,
skor_B79, skor_B80,
+ skor_B81, skor_B82, skor_B83, skor_B84, skor_B85, skor_B86, skor_B87, skor_B88,
skor_B89, skor_B90
+
+ F2 = skor_B91, skor_B92, skor_B93, skor_B94, skor_B95, skor_B96, skor_B97,
skor_B98, skor_B99, skor_B100,
+ skor_B101, skor_B102, skor_B103, skor_B104, skor_B105, skor_B106, skor_B107,
skor_B108, skor_B109, skor_B110,
+ skor_B111, skor_B112, skor_B113, skor_B114, skor_B115, skor_B116, skor_B117,
skor_B118, skor_B119, skor_B120,
+ skor_B121, skor_B122, skor_B123, skor_B124, skor_B125, skor_B126, skor_B127,
skor_B128, skor_B129, skor_B130,
+ skor_B131, skor_B132, skor_B133, skor_B134, skor_B135, skor_B136, skor_B137,
skor_B138, skor_B139, skor_B140"
> #######################Ujikecocokanmodel#######
> uji_rasch_all <- mirt(data = data_mod1, model = spec, itemtype = "Rasch", SE = T)
Iteration: 17, Log-Lik: -262335.772, Max-Change: 0.00004