Displaying 20 results from an estimated 20000 matches similar to: "feffects in model.tables"
2002 Aug 06
2
AOV in MASS not the same??
I would appreciate it if someone could explain the results of the
example from the aov() help file. The output given below is different
from book
Venables and Ripley - MASS
The results In R1.5.1 (Under windows) is as follows:
> N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0)
> P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0)
> K <-
2012 Jul 04
2
Difference between two-way ANOVA and (two-way) ANCOVA
Hi!
as my subject says I am struggling with the different of a two-way ANOVA and
a (two-way) ANCOVA.
I found the following examples from this webpage:
http://www.statmethods.net/stats/anova.html
# One Way Anova (Completely Randomized Design)
fit <- aov(y ~ A, data=mydataframe)
# Randomized Block Design (B is the blocking factor)
fit <- aov(y ~ A + B, data=mydataframe)
# Two Way
2011 Aug 06
1
multcomp::glht() doesn't work for an incomplete factorial using aov()?
Hi R users,
I sent a message yesterday about NA in model estimates (
http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html).
If I use aov() instead of lm() I get no NA in model estimates and I use
gmodels::estimable() without problems. Ok!
Now I'm performing a lot of contrasts and I need correcting for
multiplicity. So, I can use multcomp::glht() for this.
2012 Jun 25
2
Fractional Factorial - Wrong values using lm-function
Hello.
I'm a new user of R, and I have a question regarding the use of aov and
lm-functions. I'm doing a fractional factorial experiment at our production
site, and I need to familiarize myself with the analysis before I conduct
the experiment. I've been working my way through the examples provided at
http://www.itl.nist.gov/div898/handbook/pri/section4/pri472.htm
2003 Sep 20
1
factorial design
Hello all,
I´m trying to study a factorial design, but I can´t understand why did Df, Sum Sq and Mean Sq of residuals alter when I Split the interaction? I think that Split the interaction must not alter the residuals. Am I doing something wrong?
Could anyone help me?
My data and functions I tried are:
Y<-c(196,213,183,
192,253,199,
251,331,276,
2005 Feb 18
1
Two-factorial Huynh-Feldt-Test
Hi,
I'm currently working on porting some SAS scripts to R, and hence need
to do the same calculation (and get the same results) as SAS in order to
make the transition easier for users of the script.
In the script, I'm dealing with a two-factorial repeated-measures anova.
I'll try to give you a short overview of the setup:
- two between-cell factors: facBetweenROI (numbering
2003 Jun 26
1
Correct contrast for unreplicated 2K factorial design
Hi all,
I have been trying to reproduce an analysis from Douglas Montgomery?s
book on design and analysis of experiments. Table 6.10 of example 6.2 on
page 246, gives a table as follows:
> NPK <- expand.grid(A=mp,B=mp,C=mp,D=mp)
> Rate <- c(45,71,48,65,68,60,80,65,43,100,45,104,75,86,70,96)
> filtration <- cbind(NPK,Rate)
> filtration
A B C D Rate
1 - - - - 45
2
2006 Aug 18
2
4^2 factorial help
To whom it may concern:
I am trying a factorial design a system of mine that has two factors.
Each factor was set at four different levels, with one replication for
each of the combinations. My data is as follows:
A B Response
1 600 2.5 0.0257
2 600 2.5 0.0254
3 600 5
2004 Aug 02
4
Standard errors from glm
Kia ora list members:
I'm having a little difficulty getting the correct standard errors from a glm.object (R 1.9.0 under Windows XP 5.1). predict() will gives standard errors of the predicted values, but I am wanting the standard errors of the mean.
To clarify:
Assume I have a 4x3x2 factorial with 2 complete replications (i.e. 48 observations, I've appended a dummy set of data at the
2007 Sep 15
1
Cannot get contrasts to work with aov.
I have been trying for hours now to perform an orthogonal contrast through an
ANOVA in R.
I have done a two-factor factorial experiment, each factor having three
levels. I converted this dataset to a dataframe with one factor with nine
treatments, as I couldn't work out what else to do.
I have set up a matrix with the eight orthogonal contrasts that I wish to
perform, but despite
2017 Nov 29
1
2^3 confounded factorial experiment
The following R commands were written:
>help.search("factorial")
>data(npk)
>npk
>coef(npk.aov)
In the output of coef command, please explain me the interpretation of coefficients of block1 to block 6 in this 2^3 confounded factorial experiment.
Thanks.
[[alternative HTML version deleted]]
2001 Mar 10
3
Problem With Model.Tables Function
I am using R for the first time in one of my classes. My students have
alerted me to a problem for which we have not found an answer. We find
that some means returned by the model.tables function are not correct when
missing data is present in analysis of variance problems. We have
duplicated the problem using R 1.2.0, 1.2.1, and 1.2.2 under Windows 98
and several distributions of Linux (Redhat
2005 Feb 20
1
Treatment-Contrast Interactions
Hello all,
(Apologies in advance if my terminology is incorrect, I'm relatively new
to R and statistics).
I have data from a factorial design with two treatments (CRF-23), and
I'm trying to compute treatment-contrast interactions through analysis
of variance. I can't figure out how to do contrasts properly, despite
reading the help for "C" and "contrasts"
2011 Sep 27
1
ANOVA define as factor or not
Hi all
This is probably a simple problem but somehow I am having much trouble with
finding a solution, so I seek your help!
I have a data-set with continuous response variables. The explanatory
variably is 4xpH treatments (so 8.08, 7.94, 7.81 and 7.71) so also
continuous and not technically factorial.
However I have decided to do Anova's (as well as regression) to explore the
effect of
2010 Feb 26
1
factorial block design with missing data
Hello!
I have read somewhere (somehow, I can't seem to find it again, it's been a couple of months) that when analyzing factorial block design, the position where you put the block factor is important, even more when there are missing values.
I understand that when using anova.lm, the order is sequential, so that if I want to check for a treatment effect, I should put my blocking factor
2006 Aug 02
1
unbalanced mixed effects models for fully factorial designs
Does anyone know of a way of dealing with unbalanced mixed effects
(fixed and random factors) for fully factorial designs.
An example of such data is given below;
The response variable is SQRTRECRUITS
SEASON is a random factor
DENSITY is a fixed factor
Thus DENSITY:SEASON is a fixed factor.
Therefore, whereas the effects of SEASON and DENSITY:SEASON should be
tested against the overall
2009 Mar 23
1
specifying repeated measures model in lmer
Dear Colleagues,
I have what Roger Kirk (Experimental Design: Procedures for the Behavioral
Sciences, 1968) refers to as a randomized block factorial design. The anova
table would look like this:
df
A 3
Subj/A 103 (error term for A)
B 23
A*B 69
B*Subj/A 2369 (error term for B and A*B)
Subjects are nested
2003 Nov 16
1
SE of ANOVA (aov) with repeated measures and a bewtween-subject factor
Hallo!
I have data of the following design:
NSubj were measured at Baseline (visit 1) and at 3
following time points (visit 2, visit 3, visit 4).
There is or is not a treatment.
Most interesting is the question if there is a
difference in treatment between the results of visit 4
and baseline. (The other time points are also of
interest.) The level of significance is alpha=0.0179
(because of an
2007 Jul 23
1
Function to separate effect in AOV
Hi,
I have a dummy question.
Suppose that I have two explanatory variable, T1 (A, B) and T2 (C, D) and one
response variable.
> attach(dados)
> tapply(Y,list(T1,T2),mean)
C D
A 2.200000 10.20000
B 2.223333 20.26667
In this case, "A" and "B" inside "C" have no difference, but have differences
inside "D"
I make this model:
> m
2005 Jan 26
1
Specification of factorial random-effects model
I want to specify two factors and their interaction as random effects using
the function lme(). This works okay when I specify these terms using the
function Error() within the function aov(), but I can't get the same model
fitted using lme(). The code below illustrates the problem.
a <- factor(rep(c(1:3), each = 27))
b <- factor(rep(rep(c(1:3), each = 9), times = 3))
c <-