Menentukan istilah Error () dalam tindakan berulang ANOVA di R


10

Saya mengalami masalah dengan mendefinisikan istilah kesalahan untuk pengukuran dua arah berulang ANOVA di R. Data saya terdiri dari perkiraan kepadatan kayu untuk tiga posisi radial (dalam, tengah, dan luar) di sepanjang inti yang diekstraksi dari pohon. Ada total 20 spesies pohon, 6 individu dari setiap spesies, dan dua core dari masing-masing pohon.

Untuk menguji efek posisi radial pada kepadatan kayu, saya menggunakan model ANOVA dua arah berikut dengan istilah kesalahan yang menjelaskan variabilitas lintas individu:

radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual), data=Radpos)

Namun, saya tidak yakin apakah spesifikasi saya tentang istilah kesalahan memadai. Apakah saya juga harus memperhitungkan variabilitas dalam inti? Bagi saya variabilitas ini sama dengan yang disebabkan oleh posisi radial yang merupakan faktor utama yang saya minati.

Meskipun saya telah meluangkan waktu untuk membaca tentang menentukan istilah kesalahan dalam tindakan Berulang ANOVA, saya masih memiliki masalah dengan secara praktis menentukan istilah kesalahan. Saya akan menghargai bantuan dengan ini.

Jawaban:


8

Itu akan

 radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual/(Radialposition)), data=Radpos)
 summary(radpos.aov, type=3)

Itu menjelaskan kesalahan subjek dalam Radialposition. Jika Anda memiliki faktor dalam-subjek lainnya, masukkan (dalam interaksi) dengan Radialpositiondi penyebut kesalahan, misalnya + Error(Individual/(Radialpostion*wiFactorA)). Itu pemahaman saya tentang itu. Itu cocok dengan tindakan berulang SPSS GLMjika Anda tidak memiliki data yang hilang.


1
+1. Lihat lembar contekan terperinci di sini: dwoll.de/rexrepos/posts/anovaMixed.html . Dan juga di sini: personality-project.org/r/r.guide/r.anova.html
amoeba
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.