Saya memiliki kumpulan data tolok ukur dan sampel di setiap tolok ukur. Saya menjalankan benchmark ini dan subsampel mereka di mesin subjek. 'Individu' yang dipelajari oleh sub-sampel adalah sama untuk setiap mesin subjek, dan tolok ukurnya sama untuk setiap mesin subjek.
Bagaimana saya menjalankan ANOVA dalam R dalam situasi ini?
Terutama saya ingin menghitung rata-rata total dan interval kepercayaan. Saya tidak peduli tentang sarana sampel sub sama sekali, tetapi saya ingin mengenali replikasi di sana dalam keyakinan dan cara akhir. Saya mungkin peduli tentang tolok ukur berarti. Saya tidak tahu cara mengatur anova ini di R. Saya ingin dapat meniru cara dengan perhitungan manual.
Saya telah mencoba glm
, anova
, aov
, dan lme
tapi aku benar-benar bingung. Saya pikir hasil ANOVA harus setara untuk dua mesin subjek dengan rata-rata bersarang dari mesin / benchmark / pos pemeriksaan, tetapi artinya tidak keluar sama ketika saya mencobanya.
Edit:
Saya mulai mendapatkan petunjuk dari http://zoonek2.free.fr/UNIX/48_R/13.html
sumber
Jawaban:
Perbedaan utama antara desain plot split dan desain lain seperti desain acak lengkap dan variasi desain blok adalah struktur bersarang subjek, yaitu, ketika pengamatan berasal dari subjek yang sama (unit eksperimen) lebih dari sekali. Ini mengarah ke struktur korelasi dalam subjek dalam desain plot terpisah yang berbeda dari struktur korelasi dalam satu blok.
Mari kita ambil contoh gambar kumpulan data dari desain petak-petak sederhana (di bawah). Ini adalah studi tentang komposisi makanan pada kesehatan, empat diet secara acak diberikan kepada 12 subjek, semua status kesehatan yang sama. Tekanan darah awal telah ditetapkan, dan satu ukuran kesehatan adalah perubahan tekanan darah setelah dua minggu. Tekanan darah diukur pada pagi dan sore hari. (Contoh ini disalin dari buku Desain Statistik Casella contoh 5.1)
Beberapa hal penting yang perlu diperhatikan:
masing-masing di sini adalah mata pelajaran yang berbeda. Ini menggambarkan konsep bersarang. Artinya, subjek 1, 2, 3 bersarang di Diet 1. - Seluruh plot, unit eksperimental di seluruh plot (Diet) tingkat (Subjek) bertindak sebagai blok untuk perawatan plot terpisah (Pagi-Malam)x
Model untuk desain petak split ini adalah:
R
aov
sumber
Diet
danTime*Diet
dalamaov
panggilan? Itu harusTime:Diet
sesuai dengan rumus matematika Anda.