Анализ данных по оборачиваемости, дисперсионный анализ Опции

[maxandron]

Здравствуйте!
Нужна Ваша помощь. Анализируются данные по оборачиваемости. Для 15 аптек в течение 6 месяцев (с мая по октябрь) были получены данные по оборачиваемости (приведены в прикрепленном файле). В части аптек (аптеки 5,6,9,14) были проведены мероприятия по улучшению показателей оборачиваемости (в июне). Основная задача: анализ эффективности проведенных мероприятий (и интересно (?) вообще так ли нужны они были именно в этих аптеках). Я понимаю, что мне в данном случае может помочь дисперсионный анализ (я его задал через модуль GLM). Возникли вопросы по полученным результатам.
1) Для всей модели дисперсионного анализа получил таблицу (в прикрепленном файле). Можно ли как-нибудь оценить двухфакторное взаимодействие? Как получить значение ошибки (Error)?
2) Можно ли в моем случае вообще применять дисперсионный анализ? В прикрепленном файле находятся гистограммы распределения остатков (нормальность не видна) и результаты проверки однородности дисперсии (получил высоко статистически значимую неоднородность дисперсии в ячейках дисперсионного комплекса). На данном форуме прочитал про робастность дисперсионного анализа. Или все-таки необходимо проводить преобразование данных? Но если их провести, не потеряют ли они свой смысл (как интерпретировать изменение не оборачиваемости, а квадратного корня из значения оборачиваемости)?
3) Можно ли в данном случае и как грамотно провести Post hoc сравнения? В прикрепленном файле находятся данные, полученные при использовании теста Фишера и Тьюки для данных в моем примере ("значимых различий нет").
Может, еще что-нибудь посоветуете по анализу данных в моем примере. Заранее благодарен за ответы и помощь.

Прикрепленные файлы

 Analyse.rar ( 90,6 килобайт ) Кол-во скачиваний: 303

 

[p2004r]

Еще и месяц фактически не указан

Код

library(reshape2)
library(fda.usc)

read.csv2("dannye.csv")
   аптека месяц оборач
1       1   105     62
2       2   105     37
3       3   105     32
4       4   105     54
5       5   105     86
6       6   105     58
7       7   105     43
8       8   105     32
9       9   105     70
10     10   105     50
11     11   105     47
12     12   105     27
13     13   105     26
14     14   105     68
15     15   105     35
16      1   106     57
17      2   106     36
18      3   106     29
19      4   106     57
20      5   106     87
21      6   106     56
22      7   106     37
23      8   106     30
24      9   106     70
25     10   106     52
26     11   106     43
27     12   106     25
28     13   106     25
29     14   106     61
30     15   106     35
31      1   107     47
32      2   107     37
33      3   107     33
34      4   107     56
35      5   107     79
36      6   107     63
37      7   107     43
38      8   107     30
39      9   107     59
40     10   107     55
41     11   107     46
42     12   107     25
43     13   107     30
44     14   107     59
45     15   107     39
46      1   108     54
47      2   108     39
48      3   108     34
49      4   108     57
50      5   108     78
51      6   108     60
52      7   108     40
53      8   108     30
54      9   108     54
55     10   108     51
56     11   108     44
57     12   108     32
58     13   108     31
59     14   108     55
60     15   108     38
61      1   109     34
62      2   109     31
63      3   109     26
64      4   109     43
65      5   109     60
66      6   109     40
67      7   109     30
68      8   109     26
69      9   109     37
70     10   109     34
71     11   109     32
72     12   109     23
73     13   109     23
74     14   109     39
75     15   109     29
76      1   110     39
77      2   110     31
78      3   110     31
79      4   110     42
80      5   110     65
81      6   110     41
82      7   110     30
83      8   110     26
84      9   110     40
85     10   110     40
86     11   110     45
87     12   110     24
88     13   110     24
89     14   110     41
90     15   110     31


data <- read.csv2("dannye.csv")

> dcast(data, аптека~месяц)
Using оборач as value column: use value.var to override.
   аптека 105 106 107 108 109 110
1       1  62  57  47  54  34  39
2       2  37  36  37  39  31  31
3       3  32  29  33  34  26  31
4       4  54  57  56  57  43  42
5       5  86  87  79  78  60  65
6       6  58  56  63  60  40  41
7       7  43  37  43  40  30  30
8       8  32  30  30  30  26  26
9       9  70  70  59  54  37  40
10     10  50  52  55  51  34  40
11     11  47  43  46  44  32  45
12     12  27  25  25  32  23  24
13     13  26  25  30  31  23  24
14     14  68  61  59  55  39  41
15     15  35  35  39  38  29  31


plot.fdata(fdata(data.recast[,-1]))

t(scale(t(data.recast[,-1])))
             105         106        107         108        109         110
[1,]  1.2153558  0.75382831 -0.1692268  0.47691179 -1.3691984 -0.90767082
[2,]  0.5437821  0.24717366  0.5437821  1.13699884 -1.2358683 -1.23586831
[3,]  0.3986033 -0.62637656  0.7402632  1.08192316 -1.6513564  0.05694332
[4,]  0.3539075  0.77859645  0.6370335  0.77859645 -1.2032854 -1.34484841
[5,]  0.9198187  1.01029264  0.2865009  0.19602693 -1.4325045 -0.98013465
[6,]  0.5020121  0.30120725  1.0040242  0.70281691 -1.3052314 -1.20482899
[7,]  0.9753883 -0.02786824  0.9753883  0.47376005 -1.1983342 -1.19833425
[8,]  1.2247449  0.40824829  0.4082483  0.40824829 -1.2247449 -1.22474487
[9,]  1.0522754  1.05227537  0.2806068 -0.07015169 -1.2627304 -1.05227537
[10,]  0.3681605  0.61360087  0.9817614  0.49088069 -1.5953623 -0.85904121
[11,]  0.7586214  0.03034486  0.5765523  0.21241400 -1.9724157  0.39448313
[12,]  0.3100868 -0.31008684 -0.3100868  1.86052102 -0.9302605 -0.62017367
[13,] -0.1528545 -0.45856338  1.0699812  1.37569014 -1.0699812 -0.76427230
[14,]  1.2285598  0.62150670  0.4480630  0.10117551 -1.2863743 -1.11293060
[15,]  0.1286713  0.12867125  1.1580413  0.90069877 -1.4153838 -0.90069877

plot.fdata(fdata(t(scale(t(data.recast[,-1])))))

> res.pca <- fdata2pc(fdata( t(scale(t(data.recast[,-1]))) ), ncomp=4, norm=TRUE)
> summary(res.pca)

     - SUMMARY:   fdata2pc  object   -

-With 4  components are explained  99.18 %
of the variability of explicative variables.

-Variability for each component (%):
  PC1   PC2   PC3   PC4
51.51 23.20 18.20  6.27

Теперь если всё таки появится вектор воздействий можно провести функциональный аналог anova.

Эскизы прикрепленных изображений