[R] Formulating anova for partially nested model

Yves Brostaux brostaux.y at fsagx.ac.be
Tue Nov 20 17:57:47 CET 2001


Hello,

I'm trying to analyse data from an incomplete design with four factor :
- fr : number of the batch
- op : ID of operator
- meth : method used
- mat : nature of the material used
and one variable
- mv : mesure
- trmv : transformed mesure

  str(matvol)
`data.frame':   120 obs. of  6 variables:
  $ fr  : Factor w/ 30 levels "1","2","3","4",..: 1 1 1 1 2 2 2 2 3 3 ...
  $ op  : Factor w/ 3 levels "1","2","3": 3 3 3 3 1 1 1 1 2 2 ...
  $ meth: Factor w/ 2 levels "1","2": 2 2 2 2 1 1 1 1 2 2 ...
  $ mat : Factor w/ 5 levels "1","2","3","4",..: 1 5 2 4 2 1 3 5 2 1 ...
  $ mv  : num  0.761 0.806 0.837 0.264 0.799 ...
  $ trmv: num  2.12 2.23 2.31 1.08 2.21 ...


op, meth and mat are crossed and fr is nested in op*meth interaction. op 
and fr are taken as random. For each batch, only 4 on the total 5 matrices 
are mesured (oven constraint), and op and meth are constant.

I spent time trying every kind of aov models seeming sensible to me, read 
the according chapters in MASS2 and still didn't succeed writing the 
formula to obtain the results I expected (e.a. meth tested vs meth*op, mat 
vs mat*op, etc.).

Can someone bring some light to my clouded mind ?

Thanks in advance,

Y. Brostaux
=====================================================================
	YVES BROSTAUX - Ingénieur agronome Orientation Eaux & Forêts
	Assistant - Unité de Statistique et Informatique
	Gembloux Agricultural University
	8, avenue de la Faculté B-5030 Gembloux (Belgium)
         Tél:      +32 (0)81 62 24 69
         E-mail :  brostaux.y at fsagx.ac.be
=====================================================================

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list