Можно ли изменить порог ошибки типа I с помощью функции t.test ()? - PullRequest
0 голосов
/ 10 февраля 2019

Меня просят вычислить статистику теста, используя функцию t.test(), но мне нужно уменьшить ошибку типа I.Мой профессор показал нам, как изменить уровень достоверности для этой функции, но не допустимую ошибку типа I для проверки нулевых гипотез.Цель состоит в том, чтобы аргумент автоматически вычислял значение p на основе коэффициента ошибок 0,01, а не нормального значения 0,05.

Приведенный ниже код r содержит загруженный мной набор данных.

t.test(mid$log_radius_area, mu=8.456)

1 Ответ

0 голосов
/ 10 февраля 2019

Мне кажется, что я где-то ответил на это, но не могу найти его в SO или CrossValidated.

Аналогично этому вопросу ,ответ таков: t.test() не указывает любой порог отклонения / отказа от отклонения нулевой гипотезы;он сообщает значение p, и вы сами решаете, отклонить или нет.(Аргумент conf.level предназначен для настройки того, какой доверительный интервал выводит отчет.)

С ?t.test:

t.test(1:10, y = c(7:20))

    Welch Two Sample t-test

data:  1:10 and c(7:20)
t = -5.4349, df = 21.982, p-value = 1.855e-05
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -11.052802  -4.947198
sample estimates:
mean of x mean of y 
      5.5      13.5 

Здесь сообщается значение pкак 1.855e-05, поэтому нулевая гипотеза будет отклонена для любого ( предварительно определенного ) альфа-уровня> 1.855e-05.Обратите внимание, что в выходных данных нигде не говорится «нулевая гипотеза отклоняется при альфа = 0,05» или что-то в этом роде.Вы могли бы написать свою собственную функцию для этого, используя элемент $p.value, который сохраняется как часть результатов теста:

report_test <- function(tt, alpha=0.05) {
    cat("the null hypothesis is ")
    if (tt$p.value > alpha) {
       cat("**NOT** ")
    }
    cat("rejected at alpha=",alpha,"\n")
}
tt <- t.test(1:10, y = c(7:20))
report_test(tt)
## the null hypothesis is rejected at alpha= 0.05 

Большинство авторов пакетов / функций R этого не делаютпотрудитесь сделать это, потому что они считают, что это должно быть достаточно просто для пользователей, чтобы сделать для себя.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...