> Dear List,
>
> I have a problem with the ezANOVA package using for rm ANOVA or a  
> problem with my statistic design of the experiment. (or  both ;-))
>
> My file looks like this, dput(head(reg,20))

> structure(list(sample = c(13L, 14L, 15L, 13L, 14L, 15L, 13L,
14L, 15L, 13L, 14L, 15L, 13L, 14L, 15L, 13L, 14L, 15L, 1L, 2L
), days = c(7L, 7L, 7L, 14L, 14L, 14L, 21L, 21L, 21L, 28L, 28L,
28L, 35L, 35L, 35L, 42L, 42L, 42L, 7L, 7L), treat = structure(c(1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L), .Label = c("asynchron", "control", "synchron"), class =  
"factor"),
     repl = c(1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L,
     1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L), disp = structure(c(2L, 2L,
     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
     2L, 2L, 2L), .Label = c("high", "low", "medium"), class =  
"factor"),
     gbiov = c(205088809L, 236043290L, 182258079L, 356710654L,
     519036032L, 423556500L, 265916788L, 260495312L, 424304021L,
     456717288L, 412530454L, 495216491L, 360012083L, 246070292L,
     351623075L, 99075488L, 141565279L, 109303015L, 157719021L,
     275901505L), rich = c(15L, 15L, 15L, 12L, 12L, 10L, 10L,
     10L, 10L, 8L, 7L, 9L, 9L, 9L, 8L, 6L, 8L, 8L, 15L, 15L),
     evenness = c(0.810893754, 0.766130198, 0.82401103, 0.422734873,
     0.365797654, 0.406629574, 0.317890623, 0.434262493, 0.475385016,
     0.311819454, 0.307131166, 0.321269955, 0.365280751, 0.505205296,
     0.43643584, 0.484220445, 0.455471176, 0.423730082, 0.811921742,
     0.833306048)), .Names = c("sample", "days", "treat", "repl",
"disp", "gbiov", "rich", "evenness"), row.names = c(NA, 20L), class =  
"data.frame")

> I would like to compare within SSw (treatment,dispersal) over time  
> (days= between, SSb) for the listed variables gbiov, rich, evenness  
> (dv), of course a single analysis for each variable. I used the  
> repl. column as the SS identifier for the replicates, so I have for  
> each combination of dispersal x treatment 4 replicates.
>
> My problem in my dataframe is, that for the combination (asynchron x  
> low) I have only three replicates (one is missing) and by checking  
> with ezDesign, I add the missing replicates by reseting to zero for  
> each time point, but still I get error message, except when I add  
> days to the within-SSw and omit between-SSb. But that makes no sence  
> to me.

> My code is:
test<-ezANOVA(data=reg,dv=.(gbiov),wid=.(repl),within=. 
(treat,disp),between=.(days))

Error term is:
Warnung: Converting "repl" to factor for ANOVA.
Warnung: "days" will be treated as numeric.
Fehler in ezANOVA_main(data = data, dv = dv, wid = wid, within =  
within,  :
   One or more cells is missing data. Try using ezDesign() to check  
your data.

>
> Perhaps, the SS identifier makes no sence?! Or my design is wrong  
> and I have no within and between conditions Or it is a problem of  
> one missing replicate
>
>
> Thanks in advance for solving the problem
>
> Best,
> Nils






Nils Gülzow

Müllerstr. 21
26384 Wilhelmshaven
Tel.: 04421/3597682
E-mail: nils.guelzow@gmx.de







	[[alternative HTML version deleted]]

