Блог

Главная - Вопросы по программированию - Проблемы с начальной загрузкой p-значений с псевдонимными коэффициентами

Проблемы с начальной загрузкой p-значений с псевдонимными коэффициентами

#p-value #coefficients

Вопрос:

У меня проблема с псевдонимами коэффициентов в моей начальной загрузке.

Это был мой вопрос от моего учителя:

Взаимодействие с начальной загрузкой p-значения.

Загрузите значение p для F-теста взаимодействия с 5000 образцами. Сообщите о доверительном интервале 90% BCA (максимум 2 предложения). Проверить ??загрузите информацию о том, как изменить диапазон доверительных интервалов. Важно: примечание.Вам нужно будет установить «единственное число.ok = TRUE» для функции Anova (автомобильный пакет), иначе вы получите сообщение об ошибке. (2 балла)

Это код, который я ввел:

 # Bootstrap CI for ANOVA 
require(car)
library(boot)
# function to obtain p.values from F-test from the data 
p_value_two_by4<- function(data, indices) {
  data_boot <- data[indices,] # allows boot to select sample 
  ANOVAtwobyfour_boot<-Anova(lm(AttractionCombined~ColorLaptop*Nationality, data=data_boot), contrasts = c("contr.sum", "contr.poly"), type = 3, singular.ok = TRUE)
  return(ANOVAtwobyfour_boot
 #p-value #coefficients




Вопрос:






У меня проблема с псевдонимами коэффициентов в моей начальной загрузке.

Это был мой вопрос от моего учителя:


Взаимодействие с начальной загрузкой p-значения.



Загрузите значение p для F-теста взаимодействия с 5000 образцами. Сообщите о доверительном интервале 90% BCA (максимум 2 предложения). Проверить ??загрузите информацию о том, как изменить диапазон доверительных интервалов. Важно: примечание.Вам нужно будет установить "единственное число.ok = TRUE" для функции Anova (автомобильный пакет), иначе вы получите сообщение об ошибке. (2 балла)


Это код, который я ввел:


Pr(>F)`[4])
}
 
 set.seed(1342)
results_twowayaov_4 <- boot(data= Lin_Paper_NEW, statistic = p_value_two_by4, R=5000)## 5000s samples

Это мое сообщение об ошибке:

Примечание: модель имеет сглаженные коэффициенты, суммы квадратов, вычисленные путем сравнения моделей

Прав ли я, думая, что эту проблему необходимо каким-то образом исправить? Или тот факт, что R удалил их из выходных данных, означает, что эта модель теперь лучше подходит?

Это и есть результат:

ВЫЧИСЛЕНИЯ ДОВЕРИТЕЛЬНОГО ИНТЕРВАЛА НАЧАЛЬНОЙ ЗАГРУЗКИ на основе 4948 повторений начальной загрузки

ЗВОНИТЕ : boot.ci(boot.out = results_twowayaov_4, conf = 0.9, тип = «bca»)

Интервалы : Уровень BCa
90% ( 0,0000, 0,2927 )
Расчеты и интервалы в исходной шкале

Теперь это основано на 4948 репликах, а не на запрошенных 5000

Метки: