similar to: [HH] extending ancova function for 2 factors

Hi, I would like to add confidence intervales to an ANCOVA with 2 covariates when using xyplot. What would be a good way of accomplishing this? --8<---------------cut here---------------start------------->8--- rm(list = ls(all = TRUE)) rm(list = c(ls())) library(lattice) ## 1. generate data random <- rnorm(200) y <- abs(random) x1.cont <- abs(random) x2.fac <-

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 ***

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

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

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

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

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

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, 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 am trying to test some results I have for significance. It has been recommended that I use R and I am completely new to this. Set-up: Groups: two groups of 8 subjects (16 total) Two conditions: alert and passive Measurements: responses for three different stimuli (A, B, and C) measured in each condition Experiment: Testing the order of conditions Group one: Alert A, B

Hi list, does anybody know a easy to understand example in R which shows how to perfom an ANCOVA? I already tried to understand the example which you get when you type ?aov, but I did not really understand the output. Best, syrvn -- View this message in context: http://n4.nabble.com/R-example-for-ancova-tp949324p949324.html Sent from the R help mailing list archive at Nabble.com.

I am trying to run an ANCOVA with defined error terms. Thus I have to use AOV and not lm. my response variable is proportion of mice paw prints on track plates. These plates were placed on plots that had vegetation and fruit manipulated to two levels each (present or absent), and were sampled monthly for 14 months (repeated measures). The fully crossed factor design was also blocked. My sample

Hi, I am having a hard time getting what the summary.lm-output for an ANOVA / ANCOVA means. Examples I find always seem to deal with simpler cases than what I meet in my data. My main problem is understanding the output when getting significant INTERACTION TERMS (what never occurs in examples :(). The following is the output after summary.lm(ancova) where "week" is continuous,

Hi everyone, I've only been using R for a week or so now, and am now required to work with some pretty mixed data at my hospital. Basically, we're looking to see if hypertension status (nominal) affects cardiac function (multiple & continuous). In order to do this, we've collected multiple echocardiograph data reports on each patient. 1. Each patient is coded either as 0 or 1

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")

Hello R experts, I have having a difficult time figuring out how to perform and interpret an ANCOVA of my nested experimental data and would love any suggestions that you might have. Here is the deal: 1) I have twelve tanks of fish (1-12), each with a bunch of fish in them 2) I have three treatments (1-3); 4 tanks per treatment. (each tank only has one treatment applied to it) 3) I sampled

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

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