similar to: Ancova in R?

Dear R friends, I was wondering if there were factorial and modulus functions out there that I've somehow overlooked? -P. -- Peter L. Hurd, Ph.D. phurd at uts.cc.utexas.edu http://www.zo.utexas.edu/research/phurd fax 512.471-3878 Section of Integrative Biology, University of Texas, Austin TX 78712 USA

Hello R friends, I'm wondering why I get funny qqnorm() results. It seems that they should all be reflected in the normal qqline(). For instance: if I qqnorm() bimodal or uniform data I get a sigmoidal in which the qqnorm() points lie above the qqline() at -ve theoretical quantiles, and the qqnorm() points lie below the qqline() at +ve theoretical quantiles. Yet I expect such platykurtic

Hello world, I'm trying to do an anova on data in data.set, dependent variable is a column named "dep.var", grouping variable is in a column called "indep.var", and is.factor(indep.var) is TRUE... why can't I just do aov(dep.var ~ indep.var, data = data.set)? What have I done to deserve this?! What gives? Am I missing something totlly obvious? R-base-1.0.0-1,

R-users - I must preface this question by saying that I'm a relative newbie to both R and mixed-modeling. I'm using lme fit an ANCOVA-like model. My data consist of bone length measurements for a developmental series of two capuchin monkey species. I'm interested in whether the rate of bone length scaling to body mass (i.e., growth) differs between species. My call to lme

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

Dear R Users, Recently I began to use R. I`m interested about comparing several regression curves. REcently I found the package sm and the function sm.ancova which I understand allows me to do this. I applied this function to my data (seven regression curves) and I found that there are significant differences. Nevertheless, when I want to find out where are the differences, the package tells that

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

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

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

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

I've written a g.test() aka log-likelihood ratio test function for my opwn use. It's something I've seen requested (and looked to find myself) on this list a few times. It has the same basic syntax as chisq.test(). It does both goodness of fit tests and tests of independence. Yates' and Williams' corrections are implemented. I've put some examples from Sokal & Rohlf

Hi all, I have data from an experiment with 3 independent variables, 2 are within and 1 is between. In addition to the dependent variable, I have a covariate that is a single measure per subject. Below I provide an example generated data set and my approach to implementing the ANCOVA. However the output confuses me; why does the covariate only appear in the first strata? Presumably it should

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

Could you please help me on the following ANCOVA issue? This is a part of my dataset: sampling dist h 1 wi 200 0.8687212 2 wi 200 0.8812909 3 wi 200 0.8267464 4 wi 0 0.8554508 5 wi 0 0.9506721 6 wi 0 0.8112781 7 wi 400 0.8687212 8 wi 400 0.8414646 9 wi 400 0.7601675 10 wi 900 0.6577048 11 wi 900

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