Boyut 1 ve 3'ün iki örneğini karşılaştırmak için istatistiksel bir test var mı?


19

Bir ekoloji projesi için laboratuvar grubum, eşit hacimde havuz suyu içeren 4 tanka sirke, elodisiz 1 kontrol (su bitkisi) ve her birinde aynı miktarda elodea ile 3 tedavi sirke ekledi. Sirke eklemenin amacı pH'ı azaltmaktı. Hipotez, elodealı tankların normal pH'larına daha hızlı geri dönecekleriydi. Gerçekten de böyleydi. Her bir tankın pH değerini yaklaşık iki hafta boyunca günlük olarak ölçtük. Tüm tanklar sonunda doğal pH'larına döndü, ancak bunun süresi elodealı tanklar için çok daha kısaydı.

Profesörümüze deneysel tasarımımızı anlattığımızda, kontrolü tedaviyle karşılaştırmak için veriler üzerinde yapılabilecek hiçbir istatistiksel test olmadığını söyledi. Bu, kontrol için replikat olmadığından (sadece bir kontrol tankı kullandık) varyansı hesaplayamayız ve bu nedenle kontrol ve tedavinin örnek araçlarını karşılaştıramayız. Benim sorum şu, bu doğru mu? Ne demek istediğini kesinlikle anlıyorum. Örneğin, bir erkeğin ve bir kadının boyunu aldıysanız, ilgili popülasyonları hakkında sonuç çıkaramazsınız. Ama 3 tedavi yaptık ve varyans küçüktü. Varyansın kontrolde benzer olacağını varsaymak makul görünüyor mu?

Güncelleme:

Mükemmel cevap için teşekkürler. Sulak alandan daha fazla su ve elodea aldık ve deneyi daha küçük tanklarla tekrar yapmaya karar verdik, ancak bu sefer 5 kontrol ve 5 tedavi ile. Bunu orijinal verilerimizle birleştirecektik, ancak tankların başlangıç ​​pH'ları, yeni denemenin orijinal deneyle aynı popülasyondan örneklendiğini düşünmek için geçerli görünmeyecek kadar farklıydı.

Farklı miktarlarda elodea eklemeyi ve pH düzeltme hızını (pH orijinal değerine dönene kadar geçen süre olarak ölçülür) elodea miktarıyla ilişkilendirmeye çalıştık, ancak bunun gerekli olmadığına karar verdik. Amacımız sadece elodea'nın, pH'ın farklı miktarlardaki elodea'ya tam olarak nasıl tepki verdiğine dair bir tür tahmin modeli oluşturmak değil, olumlu bir fark yarattığını göstermektir. İdeal elodea miktarını belirlemek ilginç olurdu, ancak muhtemelen hayatta kalabilecek maksimum miktardır. Verilere bir regresyon eğrisi uydurmaya çalışmak, büyük miktarda eklenirken toplulukta meydana gelen çeşitli karmaşık değişiklikler nedeniyle özellikle aydınlatıcı olmaz. Elodea ölür, ayrışır, yeni organizmalar hakim olmaya başlar, vb.


4
3 'tedavi' tankının her birine aynı miktarda Elodea eklediniz mi?
gung - Monica'yı eski durumuna döndürün

2
Evet, her tedaviye aynı miktarda Elodea ekledik.
Simon Hunt

Yanıtlar:


21

Gung'un sorusuna dikkat edin; fark eder, önemi var. Tedavinin, tedavi grubundaki her tank için aynı olduğunu varsayacağım.

Varyansın iki grup için eşit olacağını iddia ederseniz (ki yine de iki örnekli t testi için varsayarsınız), bir test yapabilirsiniz. Ne kadar kötü ihlal edildiğine bakılmaksızın, bu varsayımı kontrol edemezsiniz .

İlgili bir sorunun bu cevabında ifade edilen endişeler durumunuzla daha da ilgilidir, ancak bu konuda yapabileceğiniz daha az şey vardır.

[Varyansların eşit olduğunu varsaymanın makul olmasını istiyorsunuz. Size cevap veremeyiz, bu konu uzmanlarını (yani ekologlar) makul bir varsayım olduğuna ikna etmeniz gerekir. Bu düzeylerin hem tedavi hem de kontrol altında ölçüldüğü başka çalışmalar var mı? Benzer testlerin yapıldığı diğerleri ( t-testleri veya özellikle anova - daha iyi bir emsal bulabilirsin) veya benzer varsayımlar yapıldı mı? Başvurmak için görebileceğiniz bir tür genel muhakeme?

x¯y¯σ2x¯-y¯μx-μyσ2(1/nx+1/ny)n

ny

(x¯-y¯)sx1/nx+1

sxtnx-1

σsxspny

Düzenle:

İşte bu test için simüle edilmiş bir güç eğrisi. Null'daki örnek boyutu 10000, diğer noktalarda 1000 idi. Gördüğünüz gibi, null'daki reddetme oranı 0.05'tir ve güç eğrisi, iyi bir güce sahip olmak için nüfus araçlarında büyük bir fark gerektirirken, doğru şekil. Yani, bu test ne gerekiyorsa yapar.

güç eğrisi

(Düzenlemeyi bitir)

Bununla birlikte örnek boyutları çok küçük olduğunda, bu durum dağıtım varsayımlarına karşı biraz hassas olacaktır.

Farklı varsayımlar yapmaya hazırsanız veya başka bir nüfus miktarının eşitliğini test etmek istiyorsanız, bazı testler yine de mümkün olabilir.

Yani her şey kaybolmaz ... ama mümkün olduğunda, her iki grupta da en azından biraz çoğaltılması daha iyidir.


1
@Glen_b ana hatlarıyla verilen formülleri izlemeniz gerekeceğini unutmayın. Hem Excel hem de Minitab bunu hesaplamaz.
zbicyclist

10
(+1) Eşdeğer bir yaklaşım (aynı formülü kullanarak) - ve dolayısıyla bu cevabı haklı çıkarmak için daha fazla mühimmat - tedavi grubundan gelecekteki bir değer için bir tahmin aralığı hesaplayabilmenizdir . Kontrol değeri bu tahmin aralığına girmezse, iki grup arasındaki farkın önemli kanıtlarına sahipsiniz. Fark, ortalamadaki bir farkın veya varyanslardaki bir farkın bir kombinasyonu olabilir , ancak (muhtemelen) bir fark vardır.
whuber

3
σ(x¯-12)sx1/nxt

2
@Glen_b: Bunun geçen Kasım ayında değişip değişmediğinden emin değilim, ancak örnek boyutlarından biri bir olduğunda R 3.0 toplanmış bir t testi yapacak ve bir anova ile aynı cevabı verecektir.
Aaron - Monica'yı eski

2
R denemek isteyen herkes için: t.test(x=c(4.5,4.8,4.6),y=5.2, var.equal=TRUE) - Bu R2.15.2 ve R3.0.0 (kullanışlı ben sadece iki sürümü) hem de çalışır gibi görünüyor.
Glen_b -Monica'yı geri yükleAğustos
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.