Интерпретация полученных результатов Опции

[aspir_h]

Всем здравствуйте!
Прошу строго не судить, если излагаю банальные вопросы, тем не менее, они у меня возникли, и я очень рад, что нашел ресурс, где их можно обсудить.
Читаю книгу О.Ю.Ребровой "Статистический анализ медицинских данных", раздел 10 "сравнение групп по качественному признаку", где на стр.148 приводится в качестве примера для сравнения двух групп и более использование таблиц сопряженности и метода Пирсона хи-квадрат.
Решаемая нами задача аналогична изложенной в книге, мы получили определенные результаты, и я очень прошу уважаемых форумчан подсказать, насколько правильно я интерпретировал полученные результаты.
Итак, необходимо сравнить распространенность больных с I, II,III стадиями гипертонической болезни в группах больных инфарктом миокарда и инсультом (в нашем случае, заболевания другие, данный пример придуман).
Нулевая гипотеза: - группы однородны (между группами отсутствуют различия);
- распределение по одному признаку не влияет на распределение по другому признаку.
После анализа данных получены следующие результаты (в прикрепленном файле).
Вывод: нулевая гипотеза отклоняется, изучаемые группы значимо различаются и признаки умеренно ассоциированы.

Сообщение отредактировал aspir_h - 2.02.2013 - 22:41

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

 для_сайта.rar ( 134,62 килобайт ) Кол-во скачиваний: 445

 

[DrgLena]

Матрицу я создала для этой таблицы, и задаю g=2, но
> ridit(mice,2)
Ошибка в factor(d1, labels = rownames(crosstab)) :
invalid labels; length 0 should be 1 or 2


> mice <- matrix(c(16, 20, 133, 18,61,45), nrow = 2, byrow = TRUE)
> mice # просмотр содержимого матрицы
[,1] [,2] [,3]
[1,] 16 20 133
[2,] 18 61 45

Я прицепила ridit, надеюсь на Вашу заинтересованность в R
Мои желания и возможности в данном случае не совпадают

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

 PIDIT.zip ( 930,2 килобайт ) Кол-во скачиваний: 333

 

[p2004r]

Матрицу я создала для этой таблицы, и задаю g=2, но
> ridit(mice,2)
Ошибка в factor(d1, labels = rownames(crosstab)) :
invalid labels; length 0 should be 1 or 2


> mice <- matrix(c(16, 20, 133, 18,61,45), nrow = 2, byrow = TRUE)
> mice # просмотр содержимого матрицы
[,1] [,2] [,3]
[1,] 16 20 133
[2,] 18 61 45

Я прицепила ridit, надеюсь на Вашу заинтересованность в R
Мои желания и возможности в данном случае не совпадают

Код

> library(Ridit)
> x=airquality$Ozone
>      g=airquality$Month
>      ridit(x,g)

Ridit Analysis:

Group    Label    Mean Ridit
-----    -----    ----------
1    5    0.312
2    6    0.4157
3    7    0.6673
4    8    0.6442
5    9    0.4154

Reference: Total of all groups
chi-squared = 29.2666, df = 4, p-value = 6.901e-06
> x
  [1]  41  36  12  18  NA  28  23  19   8  NA   7  16  11  14  18  14  34   6
[19]  30  11   1  11   4  32  NA  NA  NA  23  45 115  37  NA  NA  NA  NA  NA
[37]  NA  29  NA  71  39  NA  NA  23  NA  NA  21  37  20  12  13  NA  NA  NA
[55]  NA  NA  NA  NA  NA  NA  NA 135  49  32  NA  64  40  77  97  97  85  NA
[73]  10  27  NA   7  48  35  61  79  63  16  NA  NA  80 108  20  52  82  50
[91]  64  59  39   9  16  78  35  66 122  89 110  NA  NA  44  28  65  NA  22
[109]  59  23  31  44  21   9  NA  45 168  73  NA  76 118  84  85  96  78  73
[127]  91  47  32  20  23  21  24  44  21  28   9  13  46  18  13  24  16  13
[145]  23  36   7  14  30  NA  14  18  20
> g
  [1] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 6 6 6 6 6 6
[38] 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 7 7 7 7 7 7 7 7 7 7 7 7 7
[75] 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8
[112] 8 8 8 8 8 8 8 8 8 8 8 8 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9
[149] 9 9 9 9 9
> table(x,g)
     g
x     5 6 7 8 9
  1   1 0 0 0 0
  4   1 0 0 0 0
  6   1 0 0 0 0
  7   1 0 1 0 1
  8   1 0 0 0 0
  9   0 0 0 2 1
  10  0 0 1 0 0
  11  3 0 0 0 0
  12  1 1 0 0 0
  13  0 1 0 0 3
  14  2 0 0 0 2
  16  1 0 1 1 1
  18  2 0 0 0 2
  19  1 0 0 0 0
  20  0 1 1 0 2
  21  0 1 0 1 2
  22  0 0 0 1 0
  23  2 1 0 1 2
  24  0 0 0 0 2
  27  0 0 1 0 0
  28  1 0 0 1 1
  29  0 1 0 0 0
  30  1 0 0 0 1
  31  0 0 0 1 0
  32  1 0 1 0 1
  34  1 0 0 0 0
  35  0 0 1 1 0
  36  1 0 0 0 1
  37  1 1 0 0 0
  39  0 1 0 1 0
  40  0 0 1 0 0
  41  1 0 0 0 0
  44  0 0 0 2 1
  45  1 0 0 1 0
  46  0 0 0 0 1
  47  0 0 0 0 1
  48  0 0 1 0 0
  49  0 0 1 0 0
  50  0 0 1 0 0
  52  0 0 1 0 0
  59  0 0 1 1 0
  61  0 0 1 0 0
  63  0 0 1 0 0
  64  0 0 2 0 0
  65  0 0 0 1 0
  66  0 0 0 1 0
  71  0 1 0 0 0
  73  0 0 0 1 1
  76  0 0 0 1 0
  77  0 0 1 0 0
  78  0 0 0 1 1
  79  0 0 1 0 0
  80  0 0 1 0 0
  82  0 0 1 0 0
  84  0 0 0 1 0
  85  0 0 1 1 0
  89  0 0 0 1 0
  91  0 0 0 0 1
  96  0 0 0 0 1
  97  0 0 2 0 0
  108 0 0 1 0 0
  110 0 0 0 1 0
  115 1 0 0 0 0
  118 0 0 0 1 0
  122 0 0 0 1 0
  135 0 0 1 0 0
  168 0 0 0 1 0
>      ridit(table(x,g),2)

Ridit Analysis:

Group    Label    Mean Ridit
-----    -----    ----------
1    5    0.312
2    6    0.4157
3    7    0.6673
4    8    0.6442
5    9    0.4154

Reference: Total of all groups
chi-squared = 29.2666, df = 4, p-value = 6.901e-06

думаю как то так надо поправить, чтобы тип объекта был table

Код

> as.table(matrix(c(16, 20, 133, 18,61,45), nrow = 2, byrow = TRUE))
    A   B   C
A  16  20 133
B  18  61  45
> ridit(as.table(matrix(c(16, 20, 133, 18,61,45), nrow = 2, byrow = TRUE)),2)

Ridit Analysis:

Group    Label    Mean Ridit
-----    -----    ----------
1    A    0.5531
2    B    0.6649
3    C    0.4148

Reference: Total of all groups
chi-squared = 58.6525, df = 2, p-value = 1.836e-13

> ridit(as.table(matrix(c(16, 20, 133, 18,61,45), nrow = 2, byrow = TRUE)),1)

Ridit Analysis:

Group    Label    Mean Ridit
-----    -----    ----------
1    A    0.5835
2    B    0.3862

Reference: Total of all groups
chi-squared = 44.2308, df = 1, p-value = 2.919e-11

>