Как установить текстовую строку альтернативной гипотезы в объекте `htest` в R?

Question

Я создаю проверку гипотезы в R как htest объект . Мне удалось создать объект, который я хочу, с требуемой оценкой, тестовой статистикой c и p-значением. Моя единственная оставшаяся проблема заключается в том, что утверждение, которое я хочу дать для своей альтернативной гипотезы, не соответствует текстовой структуре, используемой в методе печати для объекта htest. Кажется, что настройка этих объектов предполагает, что у вас есть альтернативная гипотеза, которая представляет собой односторонний или двусторонний тест, работающий с неизвестным параметром. Кажется, что он не вмещает больше генерирующих формулировок альтернативных гипотез, например, для критериев согласия. Чтобы быть более конкретным c о моей проблеме, вот текстовая структура выходного сообщения печати для объекта htest:

alternative hypothesis: true [name(null.value)] is [less than/equal to/greater than] [null.value]

Я хотел бы получить более общий вывод на печать, например:

alternative hypothesis: [character string]

Когда вы создаете объект htest, вы можете установить name(null.value) и null.value на любую строку символов по вашему желанию, поэтому можно изменить начальную и конечную части сообщения печати на все, что вы хотите. Вы также можете установить alternative на NA, и это удалит среднюю часть. Однако промежуточные слова «истинно» и «есть» кажутся фиксированными. Это означает, что вы, кажется, застряли в сообщении со структурой true [character string] is [character string].

---

Мой вопрос: При создании объекта htest есть ли способ получить сообщение для печати альтернативной гипотезы, которое представляет собой произвольную строку символов? Если да, то как это сделать проще всего?

Answer 1

Пока вы установили x$null.value <- NULL, он будет печатать любую строку, которую вы построите для x$alternative

x <- t.test(1:10)
x$null.value <- NULL
x$alternative <- sprintf('%.2f on %s degrees of freedom, p %s',
                         x$statistic, x$parameter, format.pval(x$p.value, eps = 0.001))
x

# One Sample t-test
# 
# data:  1:10
# t = 5.7446, df = 9, p-value = 0.0002782
# alternative hypothesis: 5.74 on 9 degrees of freedom, p < 0.001
# 95 percent confidence interval:
#  3.334149 7.665851
# sample estimates:
# mean of x 
#       5.5