Bagaimana saya bisa menguji efek dalam ANOVA Split-Plot menggunakan perbandingan model yang sesuai untuk digunakan dengan X
dan M
argumen anova.mlm()
dalam R? Saya kenal ?anova.mlm
dan Dalgaard (2007) [1]. Sayangnya itu hanya menyikat Desain Split-Plot. Melakukan ini dalam desain acak lengkap dengan dua faktor dalam-mata pelajaran:
N <- 20 # 20 subjects total
P <- 3 # levels within-factor 1
Q <- 3 # levels within-factor 2
DV <- matrix(rnorm(N* P*Q), ncol=P*Q) # random data in wide format
id <- expand.grid(IVw1=gl(P, 1), IVw2=gl(Q, 1)) # intra-subjects layout of data matrix
library(car) # for Anova()
fitA <- lm(DV ~ 1) # between-subjects design: here no between factor
resA <- Anova(fitA, idata=id, idesign=~IVw1*IVw2)
summary(resA, multivariate=FALSE, univariate=TRUE) # all tests ...
Perbandingan model berikut mengarah ke hasil yang sama. Model terbatas tidak menyertakan efek yang dipermasalahkan tetapi semua efek lain dari urutan yang sama atau lebih rendah, model lengkap menambahkan efek yang dipermasalahkan.
anova(fitA, idata=id, M=~IVw1 + IVw2, X=~IVw2, test="Spherical") # IVw1
anova(fitA, idata=id, M=~IVw1 + IVw2, X=~IVw1, test="Spherical") # IVw2
anova(fitA, idata=id, M=~IVw1 + IVw2 + IVw1:IVw2,
X=~IVw1 + IVw2, test="Spherical") # IVw1:IVw2
Desain Split-Splot dengan satu faktor di dalam dan satu di antara subjek:
idB <- subset(id, IVw2==1, select="IVw1") # use only first within factor
IVb <- gl(2, 10, labels=c("A", "B")) # between-subjects factor
fitB <- lm(DV[ , 1:P] ~ IVb) # between-subjects design
resB <- Anova(fitB, idata=idB, idesign=~IVw1)
summary(resB, multivariate=FALSE, univariate=TRUE) # all tests ...
Ini adalah anova()
perintah untuk mereplikasi tes, tetapi saya tidak tahu mengapa itu bekerja. Mengapa tes perbandingan model berikut mengarah ke hasil yang sama?
anova(fitB, idata=idB, X=~1, test="Spherical") # IVw1, IVw1:IVb
anova(fitB, idata=idB, M=~1, test="Spherical") # IVb
Dua faktor dalam mata pelajaran dan satu faktor antara mata pelajaran:
fitC <- lm(DV ~ IVb) # between-subjects design
resC <- Anova(fitC, idata=id, idesign=~IVw1*IVw2)
summary(resC, multivariate=FALSE, univariate=TRUE) # all tests ...
Bagaimana cara mereplikasi hasil yang diberikan di atas dengan perbandingan model yang sesuai untuk digunakan dengan X
dan M
argumen anova.mlm()
? Apa logika di balik perbandingan model ini?
EDIT: suncoolsu menunjukkan bahwa untuk semua tujuan praktis, data dari desain ini harus dianalisis menggunakan model campuran. Namun, saya masih ingin mengerti bagaimana mereplikasi hasil summary(Anova())
dengan anova.mlm(..., X=?, M=?)
.
[1]: Dalgaard, P. 2007. Fungsi Baru untuk Analisis Multivariat. R News, 7 (2), 2-7.
lme4
paket yang sesuai dengan model DAN TIDAKlm
. Tapi ini mungkin tampilan berbasis buku yang sangat spesifik. Saya akan membiarkan komentar orang lain di atasnya. Saya dapat memberikan contoh berdasarkan bagaimana saya menafsirkannya yang berbeda dari milik Anda.Jawaban:
The
X
danM
pada dasarnya menentukan dua model yang ingin Anda bandingkan, tetapi hanya dalam hal efek dalam subyek; kemudian menunjukkan hasil untuk interaksi semua efek antar-subjek (termasuk intersepsi) dengan efek dalam-subjek yang berubah antaraX
danM
.Contoh Anda di
fitB
lebih mudah dimengerti jika kami menambahkan default untukX
danM
:Model pertama adalah perubahan dari tidak dalam efek subjek (semua memiliki rata-rata yang sama) ke rata-rata yang berbeda untuk masing-masing, jadi kami telah menambahkan
id
efek acak, yang merupakan hal yang benar untuk menguji keseluruhan intersep dan keseluruhan antara efek subjek di.Model kedua mengiklankan
id:IVw1
interaksi, yang merupakan hal yang benar untuk diujiIVw1
danIVw1:IVb
syarat - syarat yang dilawan. Karena hanya ada satu efek di dalam subjek (dengan tiga level), defaultdiag(3)
dalam model kedua akan menjelaskannya; itu akan setara dengan dijalankanUntuk Anda
fitC
, saya percaya perintah ini akan membuat ulangAnova
ringkasan.Sekarang, ketika Anda menemukan, perintah-perintah ini sangat rumit. Untungnya, tidak ada banyak alasan untuk menggunakannya lagi. Jika Anda ingin menganggap kebulatan, Anda harus menggunakan
aov
, atau untuk sintaks yang lebih mudah, gunakan sajalm
dan hitung sendiri uji-F yang benar. Jika Anda tidak mau menganggap kebulatan, menggunakanlme
adalah cara yang tepat karena Anda mendapatkan fleksibilitas yang jauh lebih banyak daripada yang Anda lakukan dengan koreksi GG dan HF.Misalnya, inilah kode
aov
danlm
untuk AndafitA
. Anda harus memiliki data dalam format panjang terlebih dahulu; inilah salah satu cara untuk melakukan itu:Dan ini
lm and
kode aov`:sumber
anova()
karena masalah dengan yangAnova()
dijelaskan di sini . Tetapi saran terakhir Anda bekerja dengan baik dan lebih sederhana. (Hal kecil: Saya pikir 2 baris terakhir masing-masing hilang 1 kurung penutup, dan seharusnya berbunyiError(id/(IVw1*IVw2))
)Desain petak-petak berasal dari pertanian, karenanya namanya. Tetapi mereka sering terjadi dan saya akan mengatakan - pekerja keras dari sebagian besar uji klinis. Plot utama diperlakukan dengan level satu faktor sedangkan level beberapa faktor lain dibiarkan bervariasi dengan subplot. Desain muncul sebagai akibat dari pembatasan pengacakan penuh. Misalnya: suatu bidang dapat dibagi menjadi empat subplot. Dimungkinkan untuk menanam varietas yang berbeda di subplot, tetapi hanya satu jenis irigasi yang dapat digunakan untuk seluruh lahan. Bukan perbedaan antara pemisahan dan blok. Blok adalah fitur dari unit eksperimental yang kami memiliki opsi untuk memanfaatkan dalam desain eksperimental, karena kami tahu mereka ada di sana. Split, di sisi lain, memaksakan pembatasan pada penugasan faktor apa yang mungkin dilakukan. Mereka memaksakan persyaratan pada desain yang mencegah pengacakan lengkap.
Mereka banyak digunakan dalam uji klinis di mana ketika satu faktor mudah berubah sementara faktor lain membutuhkan lebih banyak waktu untuk berubah. Jika eksperimen harus melakukan semua langkah untuk setiap tingkat faktor yang sulit untuk diubah secara berurutan, desain petak terbelah akan menghasilkan faktor yang sulit diubah yang mewakili seluruh faktor plot.
Berikut ini sebuah contoh: Dalam uji coba lapangan pertanian, tujuannya adalah untuk menentukan efek dari dua varietas tanaman dan empat metode irigasi yang berbeda. Delapan bidang tersedia, tetapi hanya satu jenis irigasi yang dapat diterapkan untuk setiap bidang. Bidang dapat dibagi menjadi dua bagian dengan variasi yang berbeda di setiap bagian. Faktor plot keseluruhan adalah irigasi, yang harus ditugaskan secara acak ke ladang. Dalam setiap bidang, varietas ditugaskan.
Ini adalah bagaimana Anda melakukan ini di
R
:Pada dasarnya, apa yang model ini katakan adalah, irigasi dan varietas adalah efek tetap dan varietas bersarang di dalam irigasi. Field adalah efek acak dan secara gambar akan menjadi seperti itu
Tapi ini adalah varian khusus dengan efek plot keseluruhan tetap dan efek subplot. Mungkin ada varian di mana satu atau lebih acak. Mungkin ada desain yang lebih rumit seperti split-split .. desain plot. Pada dasarnya, Anda bisa menjadi liar dan gila. Tetapi mengingat struktur dan distribusi yang mendasarinya (yaitu tetap atau acak, bersarang atau bersilangan, ..) dipahami dengan jelas, suatu
lmer-Ninja
kehendak tidak memiliki masalah dalam pemodelan. Mungkin interpretasi akan berantakan.Mengenai perbandingan, katakanlah Anda memiliki
lmer1
danlmer2
:akan memberi Anda tes yang sesuai berdasarkan statistik uji chi-sq dengan derajat kebebasan sama dengan perbedaan parameter.
cf: Faraway, J., Memperluas Model Linier dengan R.
Casella, G., Desain Statistik
sumber