The other packages that have the ability to conduct MaxDiff/BWS are very complex and I have difficulty understanding them.
Example 1: https://core.ac.uk/download/pdf/197514291.pdf
Example 2: https://cran.r-project.org/web/packages/support.BWS2/support.BWS2.pdf
The support.BWS2 package should offer some sort of solution, but I just can't figure my way through it.
The data isn't that small, it has 92000 rows and 4 collumns (22 sets per respondent, more than 500 resppondents).
Any suggestions would be helpful. I am attaching a versions (still isn't complete, as there are 10 sets per 2 respondents) where the items are coded, so that all is neater. Thank you for all your help so far!
head(maxdiff_germany_ls_RCOMMUNITY, 33)[, c("id", "block", "item", "result")]
datapasta::df_paste(head(maxdiff_germany_ls_RCOMMUNITY, 161)[, c("id", "block", "item", "result")])
data.frame(
id = c(18,18,18,18,18,
18,18,18,18,18,18,18,18,18,18,18,18,18,
18,18,18,18,18,18,18,18,18,18,18,18,18,
18,18,18,18,18,18,18,18,18,18,18,18,18,
18,18,18,18,18,18,18,18,18,18,18,18,18,
18,18,18,18,18,18,18,18,18,18,18,18,18,
18,18,18,18,18,18,18,18,18,18,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20,20,
20,20,20,20,20,20,20,20,20,20,20,20),
block = c(1,1,1,1,1,1,1,
1,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,4,
4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,6,6,
6,6,6,6,6,6,7,7,7,7,7,7,7,7,8,8,8,
8,8,8,8,8,9,9,9,9,9,9,9,9,10,10,10,
10,10,10,10,10,1,1,1,1,1,1,1,1,2,2,2,
2,2,2,2,2,3,3,3,3,3,3,3,3,4,4,4,4,4,
4,4,4,5,5,5,5,5,5,5,5,6,6,6,6,6,6,
6,6,7,7,7,7,7,7,7,7,8,8,8,8,8,8,8,
8,9,9,9,9,9,9,9,9,10,10,10,10,10,10,
10,10),
item = c(37,16,10,7,39,23,
8,28,35,18,9,4,21,36,25,19,31,16,14,
35,3,2,32,38,1,3,10,2,17,34,45,12,38,40,
13,41,18,26,29,5,17,9,15,30,20,23,22,
19,25,28,9,43,2,8,14,20,40,5,29,6,9,34,
25,28,42,29,33,19,6,35,20,5,31,19,36,
22,43,45,1,16,11,8,19,25,10,33,43,5,4,
26,1,13,12,33,21,17,36,38,3,16,13,24,11,
14,21,40,11,6,45,15,1,31,30,38,40,36,
41,27,4,22,28,8,9,30,14,35,45,42,16,20,
39,33,43,23,34,19,25,5,1,6,7,22,24,3,
14,38,40,37,32,27,31,29,44,32,22,29,20,
40,12,19),
result = c(1,0,0,-1,0,0,0,
0,0,0,0,0,1,0,-1,0,0,0,-1,1,0,0,0,0,
-1,0,0,0,0,0,0,1,0,0,-1,0,1,0,0,0,
0,0,0,0,-1,0,1,0,0,0,0,0,1,0,-1,0,0,
0,-1,0,0,1,0,0,0,0,0,1,0,0,0,-1,0,0,
0,-1,0,1,0,0,0,0,-1,1,0,0,0,0,0,0,
-1,0,0,0,0,1,0,-1,0,0,0,1,0,0,0,0,-1,
0,1,0,0,0,1,0,0,0,-1,0,0,0,1,0,0,
-1,0,0,0,0,1,0,0,0,0,0,-1,0,0,1,0,0,
-1,0,0,0,0,-1,0,0,0,1,0,0,1,0,0,0,0,
-1,0,0)
)data.frame(
id = c(18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,18,
18,18,18,18,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20,20,
20,20,20,20),
block = c(1,
1,1,1,1,1,1,1,
2,2,2,2,2,2,2,
2,3,3,3,3,3,
3,3,3,4,4,4,4,
4,4,4,4,5,5,5,
5,5,5,5,5,6,
6,6,6,6,6,6,6,
7,7,7,7,7,7,7,
7,8,8,8,8,8,
8,8,8,9,9,9,9,
9,9,9,9,10,10,
10,10,10,10,10,
10,1,1,1,1,1,
1,1,1,2,2,2,2,
2,2,2,2,3,3,3,
3,3,3,3,3,4,
4,4,4,4,4,4,4,
5,5,5,5,5,5,5,
5,6,6,6,6,6,
6,6,6,7,7,7,7,
7,7,7,7,8,8,8,
8,8,8,8,8,9,
9,9,9,9,9,9,9,
10,10,10,10,10,
10,10,10),
item = c(37,
16,10,7,39,23,
8,28,35,18,9,4,
21,36,25,19,31,
16,14,35,3,2,
32,38,1,3,10,2,
17,34,45,12,38,
40,13,41,18,26,
29,5,17,9,15,30,
20,23,22,19,25,
28,9,43,2,8,
14,20,40,5,29,6,
9,34,25,28,42,
29,33,19,6,35,
20,5,31,19,36,22,
43,45,1,16,11,
8,19,25,10,33,
43,5,4,26,1,13,
12,33,21,17,36,
38,3,16,13,24,
11,14,21,40,11,
6,45,15,1,31,30,
38,40,36,41,27,
4,22,28,8,9,
30,14,35,45,42,
16,20,39,33,43,
23,34,19,25,5,1,
6,7,22,24,3,14,
38,40,37,32,27,
31,29,44,32,22,
29,20,40,12,19),
result = c(1,
0,0,-1,0,0,0,0,
0,0,0,0,1,0,
-1,0,0,0,-1,1,
0,0,0,0,-1,0,0,
0,0,0,0,1,0,
0,-1,0,1,0,0,0,
0,0,0,0,-1,0,
1,0,0,0,0,0,1,
0,-1,0,0,0,-1,
0,0,1,0,0,0,
0,0,1,0,0,0,-1,
0,0,0,-1,0,1,
0,0,0,0,-1,1,
0,0,0,0,0,0,-1,
0,0,0,0,1,0,
-1,0,0,0,1,0,0,
0,0,-1,0,1,0,
0,0,1,0,0,0,
-1,0,0,0,1,0,0,
-1,0,0,0,0,1,
0,0,0,0,0,-1,0,
0,1,0,0,-1,0,
0,0,0,-1,0,0,
0,1,0,0,1,0,0,
0,0,-1,0,0)
)
datapasta::df_paste(head(maxdiff_germany_ls_RCOMMUNITY, 161)[, c('id', "block", "item", "result")]
library(bwsTools)
diffscoring(maxdiff_germany_ls, id="id", block = "block", item= "item", choice= "result")
Warning messages:
1: In get_checks(data, id, block, item, choice, nonbibd = TRUE) :
Analyzing non-BIBD data. Each pairwise comparison between
items does not occur for every id.
2: In get_checks(data, id, block, item, choice, nonbibd = TRUE) :
Analyzing non-BIBD data. Each pairwise comparison between items
does not occur the same amount of times for each id.