Saya telah melakukan tindakan berulang ANOVA di R, sebagai berikut:
aov_velocity = aov(Velocity ~ Material + Error(Subject/(Material)), data=scrd)
summary(aov_velocity)
- Sintaks apa dalam R yang dapat digunakan untuk melakukan tes post hoc setelah ANOVA dengan tindakan berulang?
- Apakah tes Tukey dengan koreksi Bonferroni sesuai? Jika demikian, bagaimana ini bisa dilakukan dalam R?
Jawaban:
Yang dapat Anda lakukan adalah menentukan model dengan
lme
dan kemudian menggunakanglht
darimultcomp
paket untuk melakukan apa yang Anda inginkan. Namun, lme memberikan nilai F yang sedikit berbeda dari ANOVA standar (lihat juga pertanyaan terakhir saya di sini ).Untuk kontras lain dari bonferroni, lihat misalnya, buku
multcomp
dari penulis paket .Anda juga mungkin ingin melihat pos ini pada daftar R-mailing , dan posting blog ini untuk menentukan tindakan berulang ANOVA di R .
Namun, seperti yang ditunjukkan dalam pertanyaan ini dari saya, saya tidak yakin apakah pendekatan ini identik dengan ANOVA. Lebih jauh, glht hanya melaporkan nilai- z alih-alih nilai t atau F biasa . Ini tampaknya tidak biasa juga.
Sejauh ini, saya belum menemukan cara lain untuk melakukan ini.
sumber
Jika Anda ingin tetap menggunakan
aov()
fungsi ini, Anda dapat menggunakanemmeans
paket yang dapat menanganiaovlist
(dan banyak objek lainnya ).Setelah membuat
emmGrid
objek sebagai berikutsangat mudah untuk mendapatkan semua perbandingan (post hoc) berpasangan menggunakan
pairs()
fungsi atau kontras yang diinginkan menggunakancontrast()
fungsiemmeans
paket. Penyesuaian beberapa pengujian dapat dicapai melaluiadjust
argumen dari fungsi-fungsi ini:Untuk informasi lebih lanjut tentang ini, saya menemukan sketsa emmeans rinci dan dokumentasi sangat membantu.
Juga, Anda dapat menemukan contoh lengkap (dapat direproduksi) termasuk deskripsi tentang cara mendapatkan bobot kontras yang benar dalam jawaban saya di sini .
Perhatikan, bagaimanapun, bahwa menggunakan model univariat untuk tes post hoc dapat menghasilkan nilai- p anti-konservatif jika kebulatan dilanggar.
sumber
Jika kebulatan terpenuhi maka Anda dapat menjalankan ANOVA dua arah:
sumber