Wald & LRT Test - PullRequest
Новая
       6

Wald & LRT Test

0 голосов
/ 20 января 2019

Я разработал 3 кодовых сценария на R, которые описаны следующим образом:

Сценарий I: Я импортировал набор данных из Национального института статистики и информатики из Перу, где я получил параметры ссылок покаждая страта после моделирует популяцию со структурой модели линейной регрессии.

Сценарий II: Я моделировал нормальную многомерную популяцию с параметрами, полученными в пяти стратах (тест Мардии был выполнен для каждой страты)

Сценарий III: После того, как я уже смоделировал популяцию, я извлекаю образцы двумя способами.

Сначала я извлек образец со сложным дизайном (настройка комплекса образцовдизайн) в то время как с другой стороны я извлек простую случайную выборку.Идея состоит в том, что я должен сделать тест Вальда и тест LRT для обоих образцов, в отличие от параметров, полученных для популяции.

Итак, наличие вектора параметра для контраста называется coeficientesfinal.

1) После того, как объявить мой дизайн, когда я выполняю функцию "regTermTest" 

regTermTest(reg_compleja_diseno,~IngresoHogar+MiembrosHogar+NivelInstruccion, null= coeficientesfinal[c(1,3,5)],
            method=c("Wald"), lrt.approximation="saddlepoint")

regTermTest(reg_compleja_diseno,~IngresoHogar+MiembrosHogar+NivelInstruccion, null= coeficientesfinal[c(1,3,5)],
            method=c("LRT"), lrt.approximation="saddlepoint")

ПервыйЧасть результатов теста Вальда логична, но вторая (для LRT) приводит к обратному. Может быть, я неправильно использую функцию?

2) Для того, чтобы повторить те же тесты, но для простой случайной выборки, Какие функции я должен использовать? (IE:Тест Вальда и Тест LRT.)

...