similar to: Multiple comparisons for a two-factor ANCOVA

Dear R users, i'm using a custom function to fit ancova models to a dataset. The data are divided into 12 groups, with one dependent variable and one covariate. When plotting the data, i'd like to add separate regression lines for each group (so, 12 lines, each with their respective individual slopes). My 'model1' uses the group*covariate interaction term, and so the coefficients

Please forgive me if this has been asked before but I could not readily find an answer. First, from the example I was able to determine that graphics device commands such as plot can be redirected using 'png' and 'bmp' type commands. I tried this and it works as I understand it. The question is how this graphics device works. If I issue another 'plot' command will it also

Hi R-Users, I am working with sm.ancova (in the package sm) and I have two problems with the graph, which is automatically generated when sm.ancova() is run. 1-Besides of the fitted lines, the observed data appeared automatically in the graph. I prefer that only fitted lines appear. I check the sm.options, but I could not find the way that the observed data do not appear in the graph. 2-I

Hi, what would be a good way of enhancing the ancova function from the HH package when using a 2 factor ANCOVA? The current behaviour for the "ancova" function from package HH is: ---------------------------------------------- | P1 || P1 || P3 || P4 | | PS | <- the lattice strip ---------------------------------------------- | x|| x ||x ||x | |

I've run an ancova, edadysexo is a factor with 3 levels,and log(lcc) is the covariate (continous variable) I get this results > ancova<-aov(log(peso)~edadysexo*log(lcc)) > summary(ancova) Df Sum Sq Mean Sq F value Pr(>F) edadysexo 2 31.859 15.9294 803.9843 <2e-16 *** log(lcc) 1 11.389 11.3887 574.8081 <2e-16 ***

I have some data with two categorises plus/minus (p53) and a particular time (Time) and the outcome is a continuous vairable (Result). I set up a maximum model. ancova <- lm(Result~Time*p53) > summary(ancova) .. Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 0.05919 0.55646 0.106 0.916 Time -0.02134 0.01785 -1.195 0.241 p53plus

Dear R user: I have a question about doing ANCOVA in S-plus or R. I know that many users use lm to do the regression and check the ANCOVA. But is there a way to get the traditional Table form of the ANCOVA test through S-plus (like what we would get from SPSS or SAS)? The problem I’m interested in is whether or not there is a treatment effect on some medical measurement. I will

How to Ancova in R? I know this has got to be an FAQ, because I see it asked in the lists, but I haven't seen an answer to it. I see the R-sm has the ancova thing happening, but I kind of doubt that what I'm trying to do is "smoothing"... -- Pete Hurd phurd at uts.cc.utexas.edu http://www.zo.utexas.edu/research/phurd Section of Integrative Biology, University of Texas, Austin

Dear list members: I searched the R-help for packages to calculate power for an ANCOVA problem I have. I have found power.t.test, power.anova.test. But it seems that I can not found one with ANCOVA. I have two datasetsets with variables: univariate response(one data with continous response and one with 0 and 1), treatment(two levels), covariates(x1,x2,x3). I would appreciate your help. Tony

Hello, I am VERY new to R, just picking it up infact. I have got my head around the basics of ANOVA with post hoc tests but I am struggling with regression, especially with ANCOVAs. I have two sets of data, one of type A, one of type B. Both have been placed in a wind tunnel and sampled every week. The co variate is of course the days since the start. An example is day A B 0 10.0 10.0 7 9.0

Hello,   Could you please help me in the following question: I have 16 persons 6 take 0.5 mg, 6 take 0.75 mg and 4 take placebo! Can I use the ANCOVA and t-test in this case? Is it possible in R?   Thank you in advance, Samuel [[alternative HTML version deleted]]

Hi! Can anyone tell me how to convert the anova/ancova summary output into a data frame? Thanks! Shane Phillips [[alternative HTML version deleted]]

Hi, I have a modelo like this: Yvar <- c(0, 0, 0, 0, 1, 0, 0, 0, 1, 2, 1, 1, 2, 3, 6, 6, 3, 3, 4) TIME <- 4:22 ID <- rep("PlotA",19) m <- lmer(Yvar~TIME+(TIME|ID),family=poisson) anova(m) summary(m) How to get the p-value for this case? Thanks Ronaldo -- Just because you're paranoid doesn't mean they AREN'T after you. -- > Prof. Ronaldo Reis J?nior |

Hello, Despite my search, I didn't find a post-hoc test for an ANCOVA. I used the functions aov() and lm() to run the ANCOVA then I tried TukeyHSD() but it didn't work (because of the covariable is a continuous variable?). Furthermore, I would like to plot the adjusted values (i.e. the values of the tested variable taking into account the covariable). Thanks for your help! N. Poulet

Dear All, I am a researcher doing research in plant growth and I have a statistical problem that seems to not be able to handle. Recently, I conducted an experiment about plant growing in three different nutrient-level sediments. I harvested these every three week (three harvests in all). Some growth traits of these plants were recorded (e.g. total biomass, leaf biomass and stem biomass). In

Hello, I have a problem with using the following design with ANCOVA in R. There are two groups (control + treatment), each with ten subjects. The subjects show a response that is monitored over time (four time points). For a single given subject, the response can be analysed with linear regression with time as the independent variable. The question is, how does the response differ between the

Hallo! Does anyone know a possibility to perform a sample size estimation for an ANCOVA? Would be great! Thanks Karl

I initially posted this question to one of the StackExchange sites, and they suggested that I repost my problem here. After using ezANOVA as my primary way of specifying mixed ANOVAs, I've hit a stumbling block when it come to adding a covariate to the model. I am using an ANCOVA in order to determine if there is a developmental trajectory in my data; namely, I need to be able to see the

Greetings! I have one quick question: How do you do a Tukey test on an ANCOVA? Thanks for any tips! -Paul

Please help me interpret the following results. The full model (Schwa~Dialect*Prediction*Reduction) was reduced via both update() and step(). The minimal adequate model is: ancova<-lm(Schwa~Dialect+Prediction+Reduction+Dialect:Prediction) Schwa is response variable Dialect is factor, two levels ("QF","SF") Prediction is factor, two levels ("High","Low")