не может извлечь p_value из функции lqmm в цикле R - PullRequest
Новая
       3

не может извлечь p_value из функции lqmm в цикле R

0 голосов
/ 01 марта 2019

Я запускаю несколько моделей через цикл, используя функцию lqmm в R, и я хочу получить значение P для каждой модели.Вот простой набор данных, содержащий две зависимые переменные.

data <- structure(list(y1 = c(1, 1, 4, 4, 5, 5, 7, 7), y2 = c(1, 2, 1, 2, 7, 8, 7, 8), x = c(1, 2, 3, 4, 5, 6, 7, 8), id = c(1L, 1L, 2L, 2L, 2L, 3L, 4L, 5L), weight = c(0.9401069, 1.2730856, 0.985858, 4.1602805, 1.6042408, 1.0452717, 0.9784276, 1.5199646)), class = "data.frame", row.names = c(NA,-8L))

Сначала я запустил одну модель и извлек P-значение, используя следующий код: 

model <- lqmm(fixed= y1 ~ x, random = ~1, tau=0.5, group=id,data=data,weights=data$weight,na.action=na.omit)
summary(model)$tTable[2,5]

Это сработало хорошо.

Поскольку у меня есть несколько зависимых переменных, мне нужно использовать функцию цикла для их запуска.Вот что я пытался запустить их автоматически. 

for (j in 1:2) {
  form <- formula(paste(names(data)[j], "~ x" ))
  model <- lqmm(form, random=~1, group=id, weights=data$weight, data=data, tau=0.5, na.action=na.omit)
  print(summary(model)$tTable[2,5])
}

К сожалению, это не сработало.Ошибка «Ошибка: объект типа« символ »не может быть поднабором».

Но если я просто выведу два коэффициента, сработает приведенный ниже код: 

for (j in 1:2) {
  form <- formula(paste(names(data)[j], "~ x" ))
  model <- lqmm(form, random=~1, group=id, weights=data$weight, data=data, tau=0.5, na.action=na.omit)
  print(coef(model))
}

Но мне также нужна информация о P-значениях.Я не смог найти ответ в Google.Может ли кто-нибудь оказать мне услугу?

Большое спасибо за вашу помощь.

1 Ответ

0 голосов
/ 01 марта 2019

Добавить eval() вокруг form: 

set.seed(1234)
for(j in 1:2) {
  form <- formula(paste(names(data)[j], "~ x" ))
  model <- lqmm(eval(form), random=~1, group=id, weights=data$weight, data=data, tau=0.5, na.action=na.omit)
  print(summary(model)$tTable[2,5])
}
#[1] 7.055214e-07
#[1] 0.06919567
...