Существует ли статистический тест для сравнения двух выборок размера 1 и 3?


19

Для экологического проекта моя лабораторная группа добавила уксус в 4 резервуара, содержащих равные объемы прудовой воды, 1 контрольный образец без элодеи (водное растение) и 3 обработки с одинаковым количеством элодеи в каждом. Целью добавления уксуса было снижение рН. Гипотеза заключалась в том, что резервуары с элодеей быстрее вернутся к своему нормальному pH. Это действительно было так. Мы измеряли pH каждого резервуара ежедневно в течение примерно двух недель. Все резервуары в конечном итоге вернулись к своему естественному pH, но время, которое потребовалось для этого, было намного короче для резервуаров с элодеей.

Когда мы рассказали нашему профессору о нашем экспериментальном проекте, он сказал, что не существует статистического теста, который можно было бы выполнить на данных, чтобы сравнить контроль с лечением. Из-за того, что не было никакой копии для контроля (мы использовали только один контрольный резервуар), мы не можем рассчитать дисперсию и поэтому мы не можем сравнивать выборочные средства контроля и обработки. Итак, мой вопрос, это правда? Я определенно понимаю, что он имеет в виду. Например, если вы взяли рост одного мужчины и одной женщины, вы не сможете сделать выводы об их соответствующей популяции. Но мы сделали 3 процедуры, и разница была небольшой. Кажется разумным предположить, что дисперсия будет аналогичной в контроле?

Обновить:

Спасибо за отличный ответ. Мы получили больше воды и элодеи из водно-болотных угодий и решили снова провести эксперимент с небольшими резервуарами, но на этот раз с 5 контролями и 5 обработками. Мы собирались объединить это с нашими исходными данными, но начальный pH в резервуарах был достаточно разным, так что кажется неправильным считать новый эксперимент взятым из той же популяции, что и исходный эксперимент.

Мы рассмотрели добавление различных количеств элодеи и пытались соотнести скорость восстановления pH (измеренную по истечении времени, пока pH не вернулась к своему первоначальному значению) с количеством elodea, но решили, что в этом нет необходимости. Наша цель состоит только в том, чтобы показать, что элодея имеет положительное значение, а не в том, чтобы построить какую-то прогностическую модель для точного определения того, как pH реагирует на различные количества элодеи. Было бы интересно определить оптимальное количество elodea, но это, вероятно, только максимальное количество, которое может выжить. Попытка подогнать кривую регрессии к данным не будет особенно полезной из-за различных сложных изменений, которые происходят в сообществе при добавлении большого количества данных. Элодея умирает, разлагается, начинают доминировать новые организмы и так далее.


4
Вы добавили одинаковое количество Elodea в каждый из 3 «лечебных» резервуаров?
gung - Восстановить Монику

2
Да, мы добавили одинаковое количество Elodea для каждого лечения.
Саймон Хант

Ответы:


21

Обратите внимание на вопрос Гунга; это имеет значение. Я предполагаю, что обработка была одинаковой для каждого резервуара в группе обработки.

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

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

[Вы спрашиваете, разумно ли предположить, что различия равны. Мы не можем ответить на это для вас, это то, что вы должны были бы убедить, что эксперты в данной области (например, экологи) были разумным предположением. Существуют ли другие исследования, в которых такие уровни измерялись как при лечении, так и под контролем? Другие, где были проведены подобные тесты (особенно t-тесты или anova - держу пари, что вы можете найти лучший прецедент) или сделали подобные предположения? Какую форму общего мышления вы можете применить?]

Икс¯Y¯σ2Икс¯-Y¯μИкс-μYσ2(1/NИкс+1/NY)N

NY

(Икс¯-Y¯)sИкс1/NИкс+1

sИксTNИкс-1

σsИксsпNY

Редактировать:

Вот кривая смоделированной мощности для этого теста. Размер выборки при нулевом значении составлял 10000, в остальных точках - 1000. Как видите, коэффициент отклонения при нулевом значении равен 0,05, а кривая мощности, хотя для этого требуется большая разница в совокупности, означает, что приличная мощность имеет правильная форма. То есть этот тест делает то, что должен.

кривая мощности

(Конец редактирования)

Однако при столь малых размерах выборки это будет несколько чувствительно к предположениям о распределении.

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

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


1
Обратите внимание, что вам нужно будет следовать формулам @Glen_b. И Excel, и Minitab этого не вычисляют.
zbicyclist

10
(+1) Эквивалентный подход (с использованием той же формулы) - и, следовательно, больше боеприпасов для обоснования этого ответа - это то, что вы можете вычислить интервал прогнозирования для одного будущего значения из группы лечения. Если контрольное значение не попадает в этот интервал прогнозирования, у вас есть существенные доказательства различия между двумя группами. Разница может быть некоторой комбинацией разницы в среднем или разницы в отклонениях, но есть (вероятно) разница.
whuber

3
σ(Икс¯-12)sИкс1/NИксT

2
@Glen_b: Не уверен, что это изменилось в этом ноябре прошлого года, но R 3.0 выполнит объединенный t-тест, когда один из размеров выборки равен единице, и даст тот же ответ, что и анова.
Аарон - Восстановить Монику

2
Для тех, кто хочет попробовать это в R: t.test(x=c(4.5,4.8,4.6),y=5.2, var.equal=TRUE) - похоже, это работает как в R2.15.2, так и в R3.0.0 (единственные две версии, которые мне пригодятся).
Glen_b
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.