Функция Clogit в CEDesign не сходится
Я разработал эксперимент CE, используя пакет support.CEs. Я создал CE Design с 3 атрибутами и 4 уровнями для каждого атрибута. В анкете было 4 варианта и 4 блока
des1 <- rotation.design(attribute.names = list(
Qualitat = c("Aigua potable", "Cosetes.blanques.flotant", "Aigua.pou", "Aigua.marro"),
Disponibilitat.acces = c("Aixeta.24h", "Aixeta.10h", "Diposit.comunitari", "Pou.a.20"),
Preu = c("No.problemes.€", "Esforç.economic", "No.pagues.acces", "No.pagues.no.acces")),
nalternatives = 4, nblocks = 4, row.renames = FALSE,
randomize = TRUE, seed = 987)
На вопросник ответили 15 человек (ID 1-15), поэтому 60 результатов (15 человек ответили на 4 блока:
ID BLOCK q1 q2 q3 q4
1 1 1 1 2 3 3
2 1 2 1 3 3 4
3 1 3 5 1 3 5
4 1 4 5 2 2 5
5 2 1 1 2 4 3
6 2 2 1 4 3 4
7 2 3 3 1 3 2
8 2 4 1 2 2 2
9 3 1 1 2 2 2
10 3 2 1 4 3 4
11 3 3 3 1 3 4
12 3 4 3 2 1 4
13 4 1 1 5 4 3
14 4 2 1 4 5 4
15 4 3 5 5 3 2
16 4 4 5 2 5 5
17 5 1 1 2 4 2
18 5 2 3 2 3 2
19 5 3 3 1 3 4
20 5 4 3 2 1 4
21 6 1 1 5 5 5
22 6 2 1 3 3 4
23 6 3 3 1 3 4
24 6 4 1 2 2 2
25 7 1 1 2 4 3
26 7 2 4 2 3 4
27 7 3 3 1 3 3
28 7 4 3 4 5 5
29 8 1 1 3 2 3
30 8 2 1 4 3 4
31 8 3 3 1 3 4
32 8 4 1 2 2 1
33 9 1 1 2 3 3
34 9 2 1 3 3 4
35 9 3 5 1 3 5
36 9 4 5 2 2 5
37 15 1 1 5 5 5
38 15 2 4 4 5 4
39 15 3 5 5 3 5
40 15 4 4 3 5 5
41 11 1 1 5 5 5
42 11 2 4 4 5 4
43 11 3 5 5 3 5
44 11 4 5 3 5 5
45 12 1 1 2 4 3
46 12 2 4 2 3 4
47 12 3 3 1 3 3
48 12 4 3 4 5 5
49 13 1 1 2 2 2
50 13 2 1 4 3 4
51 13 3 3 1 3 2
52 13 4 1 2 2 2
53 14 1 1 1 3 3
54 14 2 1 4 1 4
55 14 3 4 1 3 2
56 14 4 3 2 1 2
57 15 1 1 1 3 2
58 15 2 5 2 1 4
59 15 3 4 4 3 1
60 15 4 3 4 1 4
Проблема в том, что когда я объединяю матрицу вопросов и ответов с формулой
dataset1 <- make.dataset(respondent.dataset = res1,
choice.indicators = c("q1","q2","q3","q4"),
design.matrix = desmat1)
R показывает предупреждающее сообщение: В установщике (X, Y, strats, offset, init, control, weights = weights,: закончились итерации и не сходились
Я должен ожидать, что сгенерированная матрица desmat1 имела 4800 наблюдений (80 возможных комбинаций и 60 выходных данных). Вместо этого у меня есть только 1200 наблюдений. Матричный набор данных1 показывает только комбинацию из 1 набора альтернатив вместо 4.
Например, для идентификатора 1, блок 1, вопрос 1 появляется только в качестве альтернативы 1. Он совпадает с ответом, выбранным человеком, но в других случаях он не совпадает, и эта информация теряется в R, поэтому результаты, когда засорение является применяется неправильно.
Я надеюсь, что проблемы понятны. С Уважением,
Издание:
Я нашел свою проблему. Когда я делаю набор данных из респондента.dataset, сгенерированного в формате.csv, r обнаруживает только ответ q1 вместо q1-q4. DataSet1
dataset1 <- make.dataset(respondent.dataset = res1,
choice.indicators = c("q1","q2","q3","q4"),
design.matrix = desmat1)
обнаруживает q1-q4 как новые столбцы. Но ключ в том, что q1-q4 должен заполнить столбцы QES в наборе данных1. Я сделал еще один CE прежде с 1 блоком, и набор данных был правильно сделан, читая responseant.dataset. Итак, ключевой момент в том, что теперь я использую 4 блока, но я не знаю, как заставить R интерпретировать, что q1-q4 являются столбцами QUES для каждого блока.
матрица res1 (repondant.dataset) (полная матрица содержит 60 строк = 15 респондентов (ID 1-15) * 4 вопроса (столбец QES в make.dataset) введите описание изображения здесь
Добрые пожелания,