0% found this document useful (0 votes)
7 views15 pages

syntax multi

Uploaded by

cfpuny
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
7 views15 pages

syntax multi

Uploaded by

cfpuny
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 15

#memanggil data

> #dat1 <- read.csv2 ("D:/Kuliah/S2 PEP/Thesis Multidimensi/Thesis/data


multi/Paket_1062.csv", sep = ",")
> dat1 <- read.csv2 ("C:/Users/ACER1/Downloads/data thesis/Paket_1062.csv", sep =
",")
> head(dat1)
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
1 0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0
2 1 0 0 0 0 0 0 1 0 0
1 0 0 0 0 0
3 0 0 0 1 1 0 0 1 0 1
1 1 0 1 1 0
4 0 0 0 1 0 0 1 0 0 1
1 0 0 0 1 1
5 1 0 0 0 1 0 0 0 0 0
0 1 0 0 1 0
6 0 0 0 0 0 0 0 0 1 0
0 0 0 0 0 1
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
1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0
2 0 0 0 0 0 0 1 0 0
1 0 0 0 1 1
3 0 0 1 1 0 0 1 0 0
0 1 1 1 0 0
4 0 0 0 0 0 0 0 0 0
0 0 0 0 0 1
5 0 0 1 0 0 0 0 0 0
1 0 0 0 0 0
6 0 0 0 0 1 0 0 1 0
0 0 0 0 0 0
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
1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0
2 0 1 0 1 0 1 0 0 1
0 1 0 0 0 0
3 1 0 0 1 1 0 0 1 1
0 1 0 0 0 0
4 1 1 0 0 0 1 1 0 1
0 0 0 0 0 0
5 0 0 0 0 0 1 0 0 1
0 0 1 0 1 0
6 0 0 0 0 0 1 0 1 1
0 0 0 0 0 0
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
1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0
2 0 0 1 0 1 0 1 0 0
1 1 0 0 0 0
3 1 0 0 0 0 0 0 1 1
1 1 0 0 0 1
4 0 0 0 0 1 1 1 0 0
0 1 0 1 0 1
5 1 0 0 1 1 1 0 0 1
1 0 0 0 1 0
6 1 0 1 0 0 0 1 1 0
1 0 0 1 0 1
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
1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0
2 1 0 0 0 1 0 0 0 1
1 0 0 1 0 1
3 0 0 0 1 1 0 1 1 0
0 1 0 0 0 0
4 1 0 0 1 1 0 1 0 0
0 0 0 0 0 1
5 0 0 0 1 1 1 0 0 1
0 0 0 0 1 0
6 0 0 0 1 0 0 0 1 0
0 0 0 0 0 0
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 skor_B91
1 0 0 0 0 0 0 0 0 0
0 0 1 0 1 1
2 0 0 0 0 0 0 1 1 0
0 0 0 0 0 1
3 1 0 0 0 0 0 0 0 0
0 1 0 0 1 1
4 0 0 1 0 0 0 0 0 1
1 1 0 0 1 1
5 0 1 0 0 0 0 0 0 0
0 1 1 0 0 1
6 0 0 1 0 0 0 0 0 0
0 0 0 0 0 1
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
1 0 1 1 0 0 1 1 1 1
0 0 1 0 0
2 1 1 0 0 1 0 1 0 1
0 1 1 1 1
3 0 1 1 1 0 1 1 0 1
1 1 1 1 1
4 0 0 1 1 1 1 1 1 0
0 0 0 1 0
5 0 1 1 0 0 1 1 1 1
1 1 1 0 0
6 1 0 1 0 0 1 1 0 0
1 0 1 0 1
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
1 0 0 0 1 1 0 1 0
0 0 1 0 0 1
2 0 0 1 1 1 1 1 1
0 0 0 1 0 1
3 1 0 0 1 0 1 0 1
1 1 0 0 0 0
4 0 0 0 1 0 1 1 0
1 0 1 1 1 1
5 0 0 0 1 1 0 0 0
1 1 1 0 0 1
6 1 0 0 1 1 1 0 0
0 0 1 0 1 0
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
1 0 0 1 0 0 1 1 1
0 0 0 0 1 1
2 1 0 1 1 0 1 0 0
0 0 1 1 0 0
3 0 1 1 0 1 0 0 1
0 0 1 0 1 1
4 1 1 1 0 1 1 1 1
0 1 1 0 0 1
5 1 1 0 0 0 0 0 1
0 0 0 0 0 0
6 1 1 1 0 1 0 0 1
0 1 0 0 0 0
skor_B134 skor_B135 skor_B136 skor_B137 skor_B138 skor_B139 skor_B140
1 1 1 0 0 0 0 0
2 0 1 0 1 0 0 1
3 1 1 0 1 0 1 0
4 0 1 1 1 0 0 0
5 1 1 1 0 0 0 0
6 0 1 0 1 1 1 1
> str(dat1)
'data.frame': 3309 obs. of 140 variables:
$ skor_B1 : int 0 1 0 0 1 0 0 0 0 1 ...
$ skor_B2 : int 0 0 0 0 0 0 1 1 1 1 ...
$ skor_B3 : int 0 0 0 0 0 0 0 0 1 1 ...
$ skor_B4 : int 0 0 1 1 0 0 0 0 0 1 ...
$ skor_B5 : int 0 0 1 0 1 0 0 1 1 1 ...
$ skor_B6 : int 0 0 0 0 0 0 0 0 1 1 ...
$ skor_B7 : int 0 0 0 1 0 0 0 1 1 1 ...
$ skor_B8 : int 0 1 1 0 0 0 0 0 1 1 ...
$ skor_B9 : int 0 0 0 0 0 1 0 0 1 1 ...
$ skor_B10 : int 0 0 1 1 0 0 1 1 1 1 ...
$ skor_B11 : int 0 1 1 1 0 0 1 0 0 1 ...
$ skor_B12 : int 1 0 1 0 1 0 1 0 1 1 ...
$ skor_B13 : int 0 0 0 0 0 0 0 0 1 1 ...
$ skor_B14 : int 0 0 1 0 0 0 0 0 1 0 ...
$ skor_B15 : int 0 0 1 1 1 0 0 1 1 0 ...
$ skor_B16 : int 0 0 0 1 0 1 0 0 1 1 ...
$ skor_B17 : int 0 0 0 0 0 0 0 1 1 1 ...
$ skor_B18 : int 0 0 0 0 0 0 1 0 1 0 ...
$ skor_B19 : int 0 0 1 0 1 0 1 0 0 0 ...
$ skor_B20 : int 0 0 1 0 0 0 0 0 0 0 ...
$ skor_B21 : int 0 0 0 0 0 1 0 0 1 1 ...
$ skor_B22 : int 0 0 0 0 0 0 0 0 0 1 ...
$ skor_B23 : int 0 1 1 0 0 0 0 0 0 0 ...
$ skor_B24 : int 0 0 0 0 0 1 1 1 0 1 ...
$ skor_B25 : int 0 0 0 0 0 0 0 0 1 1 ...
$ skor_B26 : int 0 1 0 0 1 0 0 0 0 0 ...
$ skor_B27 : int 0 0 1 0 0 0 0 0 1 0 ...
$ skor_B28 : int 0 0 1 0 0 0 0 0 0 0 ...
$ skor_B29 : int 0 0 1 0 0 0 0 0 0 0 ...
$ skor_B30 : int 0 1 0 0 0 0 0 0 0 1 ...
$ skor_B31 : int 0 1 0 1 0 0 1 0 1 1 ...
$ skor_B32 : int 0 0 1 1 0 0 0 0 1 1 ...
$ skor_B33 : int 0 1 0 1 0 0 0 0 1 1 ...
$ skor_B34 : int 0 0 0 0 0 0 0 0 0 1 ...
$ skor_B35 : int 0 1 1 0 0 0 1 0 1 1 ...
$ skor_B36 : int 0 0 1 0 0 0 0 0 1 1 ...
$ skor_B37 : int 0 1 0 1 1 1 1 1 0 1 ...
$ skor_B38 : int 0 0 0 1 0 0 1 0 1 1 ...
$ skor_B39 : int 0 0 1 0 0 1 0 1 1 1 ...
$ skor_B40 : int 0 1 1 1 1 1 1 1 1 1 ...
$ skor_B41 : int 0 0 0 0 0 0 1 0 1 1 ...
$ skor_B42 : int 0 1 1 0 0 0 0 1 0 0 ...
$ skor_B43 : int 0 0 0 0 1 0 1 0 1 1 ...
$ skor_B44 : int 0 0 0 0 0 0 0 0 1 1 ...
$ skor_B45 : int 0 0 0 0 1 0 0 0 1 1 ...
$ skor_B46 : int 0 0 0 0 0 0 0 1 1 1 ...
$ skor_B47 : int 0 0 1 0 1 1 0 0 1 1 ...
$ skor_B48 : int 0 0 0 0 0 0 0 0 1 0 ...
$ skor_B49 : int 0 1 0 0 0 1 0 0 0 0 ...
$ skor_B50 : int 0 0 0 0 1 0 0 1 1 1 ...
$ skor_B51 : int 0 1 0 1 1 0 0 0 1 1 ...
$ skor_B52 : int 0 0 0 1 1 0 0 1 1 1 ...
$ skor_B53 : int 0 1 0 1 0 1 0 0 1 1 ...
$ skor_B54 : int 0 0 1 0 0 1 0 0 1 1 ...
$ skor_B55 : int 0 0 1 0 1 0 1 0 1 0 ...
$ skor_B56 : int 0 1 1 0 1 1 1 0 1 1 ...
$ skor_B57 : int 0 1 1 1 0 0 0 1 1 1 ...
$ skor_B58 : int 0 0 0 0 0 0 0 0 0 1 ...
$ skor_B59 : int 0 0 0 1 0 1 0 0 1 1 ...
$ skor_B60 : int 0 0 0 0 1 0 0 1 0 0 ...
$ skor_B61 : int 0 0 1 1 0 1 1 1 1 1 ...
$ skor_B62 : int 0 1 0 1 0 0 0 0 1 1 ...
$ skor_B63 : int 0 0 0 0 0 0 0 0 0 1 ...
$ skor_B64 : int 0 0 0 0 0 0 0 0 0 1 ...
$ skor_B65 : int 0 0 1 1 1 1 0 0 1 1 ...
$ skor_B66 : int 0 1 1 1 1 0 0 1 1 1 ...
$ skor_B67 : int 0 0 0 0 1 0 0 1 1 1 ...
$ skor_B68 : int 0 0 1 1 0 0 0 1 1 1 ...
$ skor_B69 : int 0 0 1 0 0 1 0 0 0 1 ...
$ skor_B70 : int 0 1 0 0 1 0 1 0 0 1 ...
$ skor_B71 : int 0 1 0 0 0 0 0 0 0 1 ...
$ skor_B72 : int 0 0 1 0 0 0 0 0 0 1 ...
$ skor_B73 : int 0 0 0 0 0 0 0 0 0 0 ...
$ skor_B74 : int 0 1 0 0 0 0 0 0 0 1 ...
$ skor_B75 : int 0 0 0 0 1 0 0 0 0 1 ...
$ skor_B76 : int 0 1 0 1 0 0 0 0 0 1 ...
$ skor_B77 : int 0 0 1 0 0 0 1 0 1 0 ...
$ skor_B78 : int 0 0 0 0 1 0 0 0 0 0 ...
$ skor_B79 : int 0 0 0 1 0 1 1 1 0 1 ...
$ skor_B80 : int 0 0 0 0 0 0 0 1 0 1 ...
$ skor_B81 : int 0 0 0 0 0 0 0 1 0 0 ...
$ skor_B82 : int 0 0 0 0 0 0 0 0 1 0 ...
$ skor_B83 : int 0 1 0 0 0 0 0 0 0 0 ...
$ skor_B84 : int 0 1 0 0 0 0 0 1 0 1 ...
$ skor_B85 : int 0 0 0 1 0 0 0 1 0 0 ...
$ skor_B86 : int 0 0 0 1 0 0 0 0 0 1 ...
$ skor_B87 : int 0 0 1 1 1 0 0 1 1 0 ...
$ skor_B88 : int 1 0 0 0 1 0 0 0 0 1 ...
$ skor_B89 : int 0 0 0 0 0 0 1 0 0 0 ...
$ skor_B90 : int 1 0 1 1 0 0 0 0 0 1 ...
$ skor_B91 : int 1 1 1 1 1 1 0 0 1 0 ...
$ skor_B92 : int 0 1 0 0 0 1 1 0 1 0 ...
$ skor_B93 : int 1 1 1 0 1 0 1 0 1 1 ...
$ skor_B94 : int 1 0 1 1 1 1 1 1 1 1 ...
$ skor_B95 : int 0 0 1 1 0 0 1 1 0 1 ...
$ skor_B96 : int 0 1 0 1 0 0 0 0 1 1 ...
$ skor_B97 : int 1 0 1 1 1 1 0 0 1 1 ...
$ skor_B98 : int 1 1 1 1 1 1 0 1 1 1 ...
$ skor_B99 : int 1 0 0 1 1 0 1 1 1 1 ...
[list output truncated]
> library(FactoMineR)
Warning message:
package ‘FactoMineR’ was built under R version 4.4.2
> library(factoextra)
Loading required package: ggplot2
Welcome! Want to learn more? See two factoextra-related books at
https://goo.gl/ve3WBa
Warning message:
package ‘factoextra’ was built under R version 4.4.2
> library(psych) #fa ada di library ini

Attaching package: ‘psych’

The following objects are masked from ‘package:ggplot2’:

%+%, 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

M method may be unstable


3: Could not invert information matrix; model may not be (empirically) identified.
> ###################### #MencarMOdelFIT######
> item_rasch <- as.data.frame(itemfit(uji_rasch_all,fit_stats = "Zh"))
> uji_3PL_all <- mirt(data = data_mod1, model = spec, itemtype = "3PL", SE = T)
Iteration: 500, Log-Lik: -258145.227, Max-Change: 0.00066

Calculating information matrix...


Warning message:
EM cycles terminated after 500 iterations.
> uji_4PL_all <- mirt(data = data_mod1, model = spec, itemtype = "4PL", SE = T)
Iteration: 500, Log-Lik: -258029.864, Max-Change: 0.00929

Calculating information matrix...


Warning messages:
1: EM cycles terminated after 500 iterations.
2: Log-likelihood was decreasing near the ML solution. EM method may be unstable
3: Could not invert information matrix; model may not be (empirically) identified.
> ###################### #MencarMOdelFIT######
> item_rasch <- as.data.frame(itemfit(uji_rasch_all,fit_stats = "Zh"))
> item_2PL <- as.data.frame(itemfit(uji_2PL_all,fit_stats = "Zh"))
> item_3PL <- as.data.frame(itemfit(uji_3PL_all,fit_stats = "Zh"))
> item_4PL <- as.data.frame(itemfit(uji_4PL_all,fit_stats = "Zh"))
> #######################Keputusan####
> item_rasch$keputusan <- ifelse(abs(item_rasch$Zh)> 1.96, "Tidak Cocok", "Cocok")
> item_2PL$keputusan <- ifelse(abs(item_2PL$Zh) > 1.96, "Tidak Cocok", "Cocok")
> item_3PL$keputusan <- ifelse(abs(item_3PL$Zh) > 1.96, "Tidak Cocok", "Cocok")
> item_4PL$keputusan <- ifelse(abs(item_4PL$Zh) > 1.96, "Tidak Cocok", "Cocok")
> #BuattabelbiarkebacanyagampaangBOZZZZ
> tabel_rasch <- table(item_rasch$keputusan)
> tabel_2PL <- table(item_2PL$keputusan)
> tabel_3PL <- table(item_3PL$keputusan)
> tabel_4PL <- table(item_4PL$keputusan)
> ####################### rangkuman item fit data_all:#####
> rangkuman_fit_all <- rbind(table(item_rasch$keputusan, exclude = 'Tidak Cocok'),
+ table(item_2PL$keputusan, exclude = 'Tidak Cocok'),
+ table(item_3PL$keputusan, exclude = 'Tidak Cocok'),
+ table(item_4PL$keputusan, exclude = 'Tidak Cocok'))
> rownames(rangkuman_fit_all) <- c('Rasch', '2-PL', '3-PL', '4-PL')
> #######################AIC, BIC, Log-likelihood data_all:#####
> anova(uji_rasch_all, uji_2PL_all, uji_3PL_all, uji_4PL_all)
AIC SABIC HQ BIC logLik X2 df p
uji_rasch_all 524955.5 525371.2 525265.8 525822.4 -262335.8
uji_2PL_all 518703.5 519523.1 519315.3 520412.8 -259071.8 6528.009 138 0
uji_3PL_all 517130.5 518359.8 518048.1 519694.3 -258145.2 1853.082 140 0
uji_4PL_all 517179.7 518818.8 518403.2 520598.2 -258029.9 230.728 140 0
> #menyimpan rangkuman item fit data all, listen, grammar, reading:
> sink('Rangkuman item fit.txt')
> cat('Kecocokan Model Rasch:\n')
> item_rasch
> cat('Kecocokan Model 2PL:\n')
> item_2PL
> #######################Uji Asumsi#####
> #INVARIANSIPARAMETER
> ganjil0 <- seq(1, nrow(dat1), by = 2)
> genap0 <- seq(2, nrow(dat1), by = 2)
> uji_ganjil <- mirt(data_mod1[ganjil0, ], model = spec, itemtype = "3PL", SE =
TRUE)
Warning messages:
1: EM cycles terminated after 500 iterations.
2: Could not invert information matrix; model may not be (empirically) identified.
> uji_genap <- mirt(data_mod1[genap0, ], model = spec, itemtype = "3PL", SE = TRUE)
Warning messages:
1: EM cycles terminated after 500 iterations.
2: Could not invert information matrix; model may not be (empirically) identified.
> coef_ganjil <- coef(uji_ganjil, IRTpars = TRUE, simplify = TRUE) # Mengambil
parameter a1 dari masing-masing model
> coef_genap <- coef(uji_genap, IRTpars = TRUE, simplify = TRUE)
> # Invariansi parameter butir a1
> a1_ganjil <- coef_ganjil$items[, "a1"]
> a1_genap <- coef_genap$items[, "a1"]
> plot(a1_ganjil, a1_genap,
+ xlab = "a1 Ganjil",
+ ylab = "a1 Genap",
+ main = "Scatterplot Parameter a1 (Ganjil vs Genap)",
+ pch = 19,
+ col = "blue",
+ cex = 0.7)
> abline(a = 0, b = 1, col = "red", lty = 2) # Garis invariansi ideal
> grid(col = "gray", lty = "dotted")
> # Invariansi parameter butir a2
> a2_ganjil <- coef_ganjil$items[, "a2"]
> a2_genap <- coef_genap$items[, "a2"]
> plot(a2_ganjil, a2_genap,
+ xlab = "a2 Ganjil",
+ ylab = "a2 Genap",
+ main = "Scatterplot Parameter a2 (Ganjil vs Genap)",
+ pch = 19,
+ col = "blue",
+ cex = 0.7)
> abline(a = 0, b = 1, col = "red", lty = 2) # Garis invariansi ideal
> grid(col = "gray", lty = "dotted")
> # Invariansi parameter butir b
> b_ganjil <- coef_ganjil$items[, "b"]
> b_genap <- coef_genap$items[, "b"]
> plot(b_ganjil, b_genap,
+ xlab = "b Ganjil",
+ ylab = "b Genap",
+ main = "Scatterplot Parameter b (Ganjil vs Genap)",
+ pch = 19,
+ col = "blue",
+ cex = 0.7)
> abline(a = 0, b = 1, col = "red", lty = 2) # Garis invariansi ideal
> grid(col = "gray", lty = "dotted")
> g_ganjil <- coef_ganjil$items[, "g"]
> g_genap <- coef_genap$items[, "g"]
> plot(g_ganjil, g_genap,
+ xlab = "g Ganjil",
+ ylab = "g Genap",
+ main = "Scatterplot Parameter g (Ganjil vs Genap)",
+ pch = 19,
+ col = "blue",
+ cex = 0.7)
> abline(a = 0, b = 1, col = "red", lty = 2) # Garis invariansi ideal
> grid(col = "gray", lty = "dotted")
> #INVARIANSIPARAMETERKEMAMPUAN
> index_ganjil <- seq(1, dim(data1)[2], 2)
> itemplot(uji_3PL_all, type = "trace", item = 54, shiny = FALSE) # Item 54
> itemplot(uji_3PL_all, type = "trace", item = 55, shiny = FALSE) # Item 55
> itemplot(uji_3PL_all, type = "trace", item = 56, shiny = FALSE) # Item 56
> itemplot(uji_3PL_all, type = "trace", item = 57, shiny = FALSE) # Item 57
> itemplot(uji_3PL_all, type = "trace", item = 58, shiny = FALSE) # Item 58
> itemplot(uji_3PL_all, type = "trace", item = 63, shiny = FALSE) # Item 63
> itemplot(uji_3PL_all, type = "trace", item = 64, shiny = FALSE) # Item 64
> itemplot(uji_3PL_all, type = "trace", item = 65, shiny = FALSE) # Item 65
> itemplot(uji_3PL_all, type = "trace", item = 99, shiny = FALSE) # Item 99
> itemplot(uji_3PL_all, type = "trace", item = 100, shiny = FALSE) # Item 100
> itemplot(uji_3PL_all, type = "trace", item = 101, shiny = FALSE) # Item 101
> itemplot(uji_3PL_all, type = "trace", item = 102, shiny = FALSE) # Item 102
> itemplot(uji_3PL_all, type = "trace", item = 103, shiny = FALSE) # Item 103
> itemplot(uji_3PL_all, type = "trace", item = 104, shiny = FALSE) # Item 104
> itemplot(uji_3PL_all, type = "trace", item = 105, shiny = FALSE) # Item 105
> itemplot(uji_3PL_all, type = "trace", item = 106, shiny = FALSE) # Item 106
> itemplot(uji_3PL_all, type = "trace", item = 117, shiny = FALSE) # Item 117
> itemplot(uji_3PL_all, type = "trace", item = 118, shiny = FALSE) # Item 118
> itemplot(uji_3PL_all, type = "trace", item = 119, shiny = FALSE) # Item 119
> itemplot(uji_3PL_all, type = "trace", item = 120, shiny = FALSE) # Item 120
> itemplot(uji_3PL_all, type = "trace", item = 121, shiny = FALSE) # Item 121

You might also like