Как мне запустить тест Wald/Score/Lioneli Ratio Test в R только с одним вектором?

Я вычислил статистику теста вручную, но хочу подтвердить, что мой ответ правильный. Я новичок в R, и я думаю, что это должно быть довольно легко, но я не знаю, с чего начать, назначив векторные значения.Как мне запустить тест Wald/Score/Lioneli Ratio Test в R только с одним вектором?

С большой выборкой, статистика испытаний Вальда, статистика результатов теста и статистика вероятности отношения правдоподобия приблизительно имеют распределение chisquare df = 1. Y ~ Binomial (1, pi) с P (Y = 1) = pi и P (Y = 0) = 1-pi. Наблюдения Yi, i = 1 до m равны Y.

Для m = 10 и данных (0, 1, 0, 1, 0, 0, 0, 1, 0) используйте эти статистические данные для проверить нулевую гипотезу pi = 0,3.

Отказ от ответственности: Хотя это домашнее задание связано, R-часть не является частью домашней работы. Мне любопытно, могу ли я запустить эту проблему в R, чтобы подтвердить свои ответы.

источник

2016-10-02 FettahP

проверить '? Prop.test' и'? Glm' ...? –

Вы провели какое-либо исследование, как это сделать в R самостоятельно? Вы должны проверить это ... http://stackoverflow.com/help/how-to-ask –

Я знаю, что могу запускать 'glm', когда у меня есть зависимая переменная и, по крайней мере, одна независимая переменная, и она будет производить эти тестовая статистика. Однако, как мне получить статистику, если есть только один вектор. – FettahP

Вам нужно указать «1», чтобы представлять постоянный термин в вашем glm.

x=c(0, 1, 0, 0, 1, 0, 0, 0, 1, 0) 
glm(x~1)

Это эквивалентно установке GLM без данных, просто перехват.

источник

2016-10-02 02:11:29 nograpes

Это было бы испытанием для pi = 0,5 –

Спасибо, это имеет такой смысл! – FettahP

Вы должны установить смещение термина, который представляет вашу нулевую гипотезу о Logg форе масштаб:

glm(dat ~ offset(rep(log(0.3/0.7),10)), family=binomial) 

Call: glm(formula = dat ~ offset(rep(log(0.3/0.7), 10)), family = binomial) 

Coefficients: 
(Intercept) 
    1.17e-12 

Degrees of Freedom: 9 Total (i.e. Null); 9 Residual 
Null Deviance:  12.22 
Residual Deviance: 12.22 AIC: 14.22

Оценка перехватывает по существу 0, так как наблюдаемые шансы на Дэ выдвинули гипотезу стоимости; однако количество событий настолько невелико, что статистическая значимость на обычных уровнях вряд ли будет установлена независимо от результата. Значение p - это вероятность того, что данные не соответствуют гипотетическому значению, тогда как в этом случае это происходит.

источник

2016-10-02 07:49:36

Как мне запустить тест Wald/Score/Lioneli Ratio Test в R только с одним вектором?

ответ

Смежные вопросы