[R] Is it feasible with R?

Louis Longchamps kpangolo at hotmail.com
Thu Nov 17 03:04:23 CET 2005


   My data are (for one experiment):

   3 Types of plant (1,2,3)

   4  Species  per  Type (ex:for type 1(a,b,c,d), for type 2(e,f,g,h) and
   for type 3(i,j,k,l)

   8 Repetitions of each Species

   3 Stages (10, 20, 30)

   2  Measures per Stages, for the stages 20 and 30 (ex: for Stage 10(C),
   for Stage 20 (A and B) and for Stage 30 (A and B)

   3  Types  x  4 Species(Type) x 8 Repetitions x (1 Stage x 1 Mesure + 2
   Stages  x  2  Measures= 480 data. There are 10 data missing. The final
   number of data is 470.

   My questions that I ask to my data:

   1- Is there a significant difference between Measure A and Measure B?

   2- Is there a significant difference between Stages 10, 20 and 30?

   3- Is there a significant difference between the 3 plant Types?

   I haven't found  any way to consider the whole data lot in one shot, I
   decided to remove the first stage with the Measure "C" that comes with
   it. It bring the number of data to 384.

   Is it possible to consider the whole lot in one shot?

   Thus my model is

   Source                                                d.l.

   Types                                                  2

   Species (Type)                                      9

   Error a)                                                84


   Stage                                                  1

   Measure                                               1

   Stage x Measure                                   1

   Stage x Type                                        2

   Measure x Type                                     2

   Stage x Measure x Type                         2

   Stage x Species (Type)                          9

   Measure x Species (Type)                       9

   Stage x Measure x Species (Type)            9

   Erreur b)                                             252

   Total                                                   383

   How  can I write this in R? My problem is that when I use %in% to nest
   Species in Type, R looks, for exemple, for specie "a" in Type 2 and it
   finds nothing and returns "Na".

   Type1:Plante[T.a]  -1.51967    0.82774  -1.836   0.0672 .
   Type2:Plante[T.a]        NA         NA      NA       NA

   Can  you  help  me?  Can  R  consider  Species a,b,c,d only in Type 1,
   Species e,f,g,h only in Type 2 and Species i,j,k,l only in Type 3?

   Thank you very much.



More information about the R-help mailing list