prop.test как тут быть Опции

[kont]

Коллеги, смотрите, как пользоваться проп.тестом, если у нас всего двупольная таблица ,я знаю! Но что делать, если возможных градаций 3?

Код

    полный ответ    Стабилизация    Частичный ответ
Группа    n    n    n
R    0    2    28
T    1    1    28

Т.е. в каждой группе по 30, всего 60.
И тут я впал в ступор, я не могу понять, а как составить таблицу для трех градаций, чтобы посмотреть межгрупповое различие. Мне просто нужно ещё и 95% ДИ рассчитать.
Может кто подскажет

[kont]

Коллеги, спасибо! Помогло! DrgLena prop.test это функция сравнения долей на языке R. Правда когда у нас 3 градации в факторе, 95%ДИ не считает ни хи-квадрат , ни по - пролетарски Фишер.

> x=matrix(c(0,2,28,1,1,28),ncol=3,byrow=T)
> x
[,1] [,2] [,3]
[1,] 0 2 28
[2,] 1 1 28
> fisher.test(x)

Fisher's Exact Test for Count Data

data: x
p-value = 1
alternative hypothesis: two.sided

Если 2х2, то ДИ высчитывается.

[100$]

Коллеги, спасибо! Помогло! DrgLena prop.test это функция сравнения долей на языке R. Правда когда у нас 3 градации в факторе, 95%ДИ не считает ни хи-квадрат , ни по - пролетарски Фишер.

> x=matrix(c(0,2,28,1,1,28),ncol=3,byrow=T)
> x
[,1] [,2] [,3]
[1,] 0 2 28
[2,] 1 1 28
> fisher.test(x)

Fisher's Exact Test for Count Data

data: x
p-value = 1
alternative hypothesis: two.sided

Если 2х2, то ДИ высчитывается.

/С искренним любопытством/
А какой, стесняюсь спросить, доверительный интервал может выдать точный метод Фишера, если он по определению выдает вероятность (к тому же одностороннюю!) заполнения таблицы сопряженности? А для вычисления двусторонней вероятности надо генерировать таблицы по определенным правилам. Н-р, для расширения точного теста Фишера на таблицы RxC (известного как критерий Фримана-Холтона), их надо генерить, н-р, по алгоритму Пэтфилда (Patefield, algorithm AS 159). Жду ответа, аки соловей лета.

Сообщение отредактировал 100$ - 19.02.2017 - 14:42