Re: [R] Blocking and Nesting ANOVA

2019-07-02 Thread Bob O'Hara
No, it isn't correct I'm afraid. From what you describe, Population isn't nested within Sex or Maturity. Population 1 can contain both male and female trees, as well as trees of different ages. A simpler model would be lm (PAC ~ Season + Population + Sex + Maturity) So that Sex and Maturity

[R] Blocking and Nesting ANOVA

2019-07-02 Thread Eleftheria Dalmaris
Dear all, I have a blocking and nested design. I have Season (3 levels; three different collection periods). That's my Block. I have three Populations that I sampled from (Plant material) >From each population, I sampled 30 trees. Each Tree is either Female or Male or Unidentified (Sex factor)

[R] Error Message During ANOVA

2019-03-21 Thread Sara Ciancitto via R-help
I know this is several years later but... I ran into the same issue and was able to solve it by verifying the spacing between variables within original text file. For example, Item1, Item2, and Item3 had a single space between them while the numbers in each item had 5 spaces between them,

Re: [R] Question Mixed-Design Anova in R

2018-11-23 Thread Lisa van der Burgh
Mixed-Design Anova in R Dear Lisa, > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of peter > dalgaard > Sent: Friday, November 23, 2018 10:16 AM > To: Lisa van der Burgh <40760...@student.eur.nl> > Cc: r-help@R-project.org >

Re: [R] Question Mixed-Design Anova in R

2018-11-23 Thread Fox, John
Dear Lisa, > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of peter > dalgaard > Sent: Friday, November 23, 2018 10:16 AM > To: Lisa van der Burgh <40760...@student.eur.nl> > Cc: r-help@R-project.org > Subject: Re: [R] Quest

Re: [R] Question Mixed-Design Anova in R

2018-11-23 Thread peter dalgaard
> wrote: > > Hi Everyone, > > > > I have a question about Mixed-Design Anova in R. I want to obtain Mauchly�s > test of Sphericity and the Greenhouse-Geisser correction. I have managed to > do it in SPSS: > > > > GLM Measure1 Measure2 Measure3 Measure4

[R] Question Mixed-Design Anova in R

2018-11-23 Thread Lisa van der Burgh
Hi Everyone, I have a question about Mixed-Design Anova in R. I want to obtain Mauchly�s test of Sphericity and the Greenhouse-Geisser correction. I have managed to do it in SPSS: GLM Measure1 Measure2 Measure3 Measure4 Measure5 Measure6 BY Grouping /WSFACTOR=Measure 6 Polynomial

Re: [R] Formatting multi-way ANOVA output for spectra analysis

2018-07-25 Thread Fox, John
Dear Robert, Although you don't say so, it sounds as if you may be using the Anova() function in the car package, which is what the R Commander uses for ANOVA. If so, in most cases, Anova() returns an object of class c("anova", "data.frame"), which can be manipulated as

Re: [R] Formatting multi-way ANOVA output for spectra analysis

2018-07-25 Thread Jeff Newmiller
h site.  The output of the >spectrometer per sample is values for 2048 individual wavelengths, in a > >spreadsheet with the wavelength as the first column.  Since I'm doing >the analysis wavelength-by-wavelength, I've transposed the data and >broke the data for the project down into smal

[R] Formatting multi-way ANOVA output for spectra analysis

2018-07-25 Thread Robert D. Bowers M.A.
the data for the project down into smaller spreadsheets (so that R can perform ANOVA on each wavelength). The problem is, I can do ANOVA now on each wavelength, but I don't need a full output table for each... I just need to know if there is significant variation between any of the sites

Re: [R] question about the anova() function for deviance analysis

2017-04-17 Thread Duncan Murdoch
On 16/04/2017 9:46 AM, Sophie Dubois wrote: Dear Maintener, I have recently had a bad experience with the anova() function. Indeed, I wanted to process a deviance analysis between 2 mixed linear models and I was really surprise to see that depending on the ordre in which I gave my models, the

Re: [R] question about the anova() function for deviance analysis

2017-04-16 Thread Jeff Newmiller
You are sending your email to a whole mailing list of volunteers, not a specific "maintainer" (and I am not one). However, your assertions convey unfamiliarity with statistics rather than deficiencies in R, and this mailing list is not a stats tutoring list. I did a quick Google search and

Re: [R] question about the anova() function for deviance analysis

2017-04-16 Thread Bert Gunter
This list is about R programming; your question seems mostly about statistics, and is therefore off topic here. I suggest you consult a local statistical expert who *is* comfortable with such statistical analyses. In general, partitions of sums of squares in statistical models can depend on the

[R] question about the anova() function for deviance analysis

2017-04-16 Thread Sophie Dubois
Dear Maintener, > I have recently had a bad experience with the anova() function. > Indeed, I wanted to process a deviance analysis between 2 mixed linear > models and I was really surprise to see that depending on the ordre in > which I gave my models, the function did not the same thing: once it

Re: [R] use one way ANOVA to select genes

2016-04-05 Thread PIKAL Petr
etr <petr.pi...@precheza.cz> Subject: Re: [R] use one way ANOVA to select genes thanks Petr Actually , i have wrote code for One way ANOVA ,but still i am sure is it true or not for analysis data , because I have 5 groups i want to select the significant gene from this gro

Re: [R] use one way ANOVA to select genes

2016-04-04 Thread PIKAL Petr
Sent: Sunday, April 3, 2016 4:45 PM To: R-help Mailing List <r-help@r-project.org> Subject: [R] use one way ANOVA to select genes i want to select the significant genes form 5 clusters (groups) by one way ANOVA

[R] use one way ANOVA to select genes

2016-04-03 Thread hehsham alpukhity via R-help
i want to select the significant genes form  5 clusters (groups) by one way ANOVA  in r### # i want use One way ANOVA to select the siginificant from the clusters above  selectgene <- function(GropuData,pva

Re: [R] R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-07 Thread Meyners, Michael
...@virgilio.it Cc: Meyners, Michael; r-help@r-project.org Subject: Re: [R] R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design Dear Angelo, The Bonferroni p-value is just the ordinary p-value times the number of tests, so, since R

[R] R: Re: R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-07 Thread angelo.arc...@virgilio.it
originale Da: j...@mcmaster.ca Data: 7-lug-2015 3.23 A: angelo.arc...@virgilio.itangelo.arc...@virgilio.it Cc: meyner...@pg.com, r-help@r-project.orgr-help@r-project.org Ogg: Re: [R] R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design Dear

[R] R: RE: R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-07 Thread angelo.arc...@virgilio.it
-help@r-project.orgr-help@r-project.org, John Foxj...@mcmaster.ca Ogg: RE: [R] R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design Angelo, the conservatism of Bonferroni aside for a moment, I don't really see what you are after. mcp

[R] R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-06 Thread angelo.arc...@virgilio.it
Angelo Messaggio originale Da: meyner...@pg.com Data: 6-lug-2015 17.52 A: angelo.arc...@virgilio.itangelo.arc...@virgilio.it, r-help@r-project.orgr-help@r-project.org Ogg: RE: [R] Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design Untested

Re: [R] R: RE: Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-06 Thread John Fox
Angelo Messaggio originale Da: meyner...@pg.com Data: 6-lug-2015 17.52 A: angelo.arc...@virgilio.itangelo.arc...@virgilio.it, r-help@r-project.orgr-help@r-project.org Ogg: RE: [R] Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between

[R] Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-06 Thread angelo.arc...@virgilio.it
Dear List Members, I need to perform a Bonferroni post hoc test in R on a table with three within subjects factors (Emotion, having 5 levels, Material, having 4 levels, Shoes, having 2 levels) and one between subject factor (Musician, having 2 levels). I normally use the Tukey method

Re: [R] Bonferroni post hoc test in R for repeated measure ANOVA with mixed within and between subjects design

2015-07-06 Thread Meyners, Michael
test in R for repeated measure ANOVA with mixed within and between subjects design Dear List Members, I need to perform a Bonferroni post hoc test in R on a table with three within subjects factors (Emotion, having 5 levels, Material, having 4 levels, Shoes, having 2 levels) and one

[R] Help on phylogenetic ANOVA functions

2015-05-12 Thread Guido Buezas
Hi everyone: I'm trying to do a phylogenetic ANOVA analysis, but i'm stuck with both possibilities: phylANOVA from phytools package, and aov.phylo from geiger package. Here is what i have Input data: CSV file with columns: Spp (categorical, for taxonomic species,, 29 species) Fam (categorical,

Re: [R] Wrong results from anova

2015-04-14 Thread William Dunlap
It is a printing problem - the default number of digits in print.summary.aov is max(3L, getOption(digits) - 3L). Set options(digits=7) instead of your current 6 (?) or try print(summary(fit.Y), digits=7) to see more digits. Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, Apr 14, 2015 at

Re: [R] Wrong results from anova

2015-04-14 Thread Phil Spector
Li Li - I belive it's a rounding error -- try setting options(digits=8) before displaying the output: summary(fit.Y) Df Sum Sq Mean Sq F value Pr(F) sample.Y 3 2203 734 190706 2e-16 *** Residuals8 0 0 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01

[R] Wrong results from anova

2015-04-14 Thread li li
Hi all, I have following data. When I perform an anova, the residual sum of square returns to be zero. But this is wrong, since we can hand calculate the RSS to be around 0.0308. Did anyone come across the same problem before? Any suggestions? Thanks. Hanna ydata Y sample.Y 1

Re: [R] Wrong results from anova

2015-04-14 Thread Michael Dewey
See in-line On 14/04/2015 18:23, li li wrote: Hi all, I have following data. When I perform an anova, the residual sum of square returns to be zero. But this is wrong, since we can hand calculate the RSS to be around 0.0308. Did anyone come across the same problem before? Any suggestions?

Re: [R] Wrong results from anova

2015-04-14 Thread li li
0.003725088 is the mean square error. 2015-04-14 13:44 GMT-04:00 Michael Dewey li...@dewey.myzen.co.uk: See in-line On 14/04/2015 18:23, li li wrote: Hi all, I have following data. When I perform an anova, the residual sum of square returns to be zero. But this is wrong, since we

[R] Calculating contrasts in Anova with interaction term

2014-08-21 Thread Lynn Govaert
Hi all, I have the following variables: a trait z, Treatment A and B, Population 1 and 2 (each having undergo the treatment), and their interaction. Now, I'm doing an Anova as follows aov( z ~ Treatment * Population ) I now want to compare the values of population 1 which underwent treatment A

Re: [R] Nonparametric k-way ANOVA

2013-10-25 Thread Gerrit Eichner
Hello, Vincent, you may want to take a look at Nonparametric methods in factorial designs by Edgar Brunner and Madan L. Puri in Statistical Papers 42, 1-52 (2001). There is the R-package nparcomp for one-way layouts, but the paper goes further (and mentions another software) and is maybe a

Re: [R] Nonparametric k-way ANOVA

2013-10-25 Thread Vicent Giner-Bosch
Thanks to all!! Gerrit, thank you for the reference. I was looking for something like that, but I was not able to find it. I partially agree with some of you, who said that nonparametric k-way ANOVA is a little bit of contradiction, but I just wanted to have an alternative to classical ANOVA.

[R] Nonparametric k-way ANOVA

2013-10-24 Thread Vicent Giner-Bosch
Sorry if this subject has been already dealt here. Which are some common tests for nonparametric k-way ANOVA? I have read about Kruskal-Wallis test as a kind of nonparametric one-way ANOVA, but I have not found anything about a general-setting (I mean k-way) nonparametric ANOVA. Can you

Re: [R] Nonparametric k-way ANOVA

2013-10-24 Thread peter dalgaard
On Oct 24, 2013, at 18:59 , Vicent Giner-Bosch wrote: Sorry if this subject has been already dealt here. Which are some common tests for nonparametric k-way ANOVA? I have read about Kruskal-Wallis test as a kind of nonparametric one-way ANOVA, but I have not found anything about a

Re: [R] Nonparametric k-way ANOVA

2013-10-24 Thread Cade, Brian
Depends on what the original poster really wants with a nonparametric k-way ANOVA. If what they are really looking for is an approach that has fewer distributional assumptions, then it is possible to perform permutation tests for linear models parameterized for k-way anova that would eliminated

[R] Few doubts about ANOVA

2013-08-29 Thread bala chand
Hi can you please give the brief explanation about anova? what is the purpose of null hypothesis in anova? how can we find future predictive value from existing data? [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] Few doubts about ANOVA

2013-08-29 Thread John Kane
Looks like school is starting up again. We don't usually help with homework especially at this level. Read a text book John Kane Kingston ON Canada -Original Message- From: bal.chan...@gmail.com Sent: Thu, 29 Aug 2013 15:57:29 +0530 To: r-help@r-project.org Subject: [R] Few doubts

Re: [R] Error Message During ANOVA

2013-03-21 Thread Rodinsky, Dr Harold R.
I am teaching myself R for use in Psych research. I don't understand the error message I am getting or how to fix it. Any suggestions will be greatly appreciated df1=read.table(fastfood.txt, header=TRUE); df1 c(t(as.matrix(df1))) r = c(t(as.matrix(df1))) # response data r

Re: [R] Error Message During ANOVA

2013-03-21 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Rodinsky, Dr Harold R. Sent: Thursday, March 21, 2013 11:32 AM To: 'r-help@r-project.org' Subject: Re: [R] Error Message During ANOVA I am teaching myself R for use in Psych

Re: [R] Error Message During ANOVA

2013-02-13 Thread Pascal Oettli
Hello, The first error message says: Error in model.frame.default(formula = r ~ tm, drop.unused.levels = TRUE) : variable lengths differ (found for 'tm') Please compare the size of 'r' and 'tm'. Regards, Pascal Le 14/02/2013 07:14, Craig O'Connell a écrit : Dear R-Help, I'm using

[R] Mean SD in ANOVA Table

2012-10-08 Thread Jhope
Hi R-listers, Is there a way to get the mean SD value after creating an ANOVA table? This is what I've been using to create an anova table. summary(aov(HSuccess ~ HTLIndex*Aeventexhumed, data=data.to.analyze)) Best, Jean -- View this message in context:

Re: [R] significance of character (Anova)

2012-08-13 Thread peter dalgaard
On Aug 13, 2012, at 15:38 , Jennifer Kaiser wrote: Model: poisson, link: log Response: sb_ek_ber Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev Pr(Chi) NULL1237837 4.4998e+10 ABWHALT_C 2

Re: [R] significance of character (Anova)

2012-08-13 Thread Jessica Streicher
Also, since there apparently were warnings, how about adding them? On 13.08.2012, at 16:34, peter dalgaard wrote: On Aug 13, 2012, at 15:38 , Jennifer Kaiser wrote: Model: poisson, link: log Response: sb_ek_ber Terms added sequentially (first to last) Df Deviance Resid.

[R] significance of character (Anova)

2012-08-10 Thread Jennifer Kaiser
Hi,  I have a problem with the output of my anova. Tabelle2 - data.frame(sb_ek_ber, JE, ABWHALT_C, ALLEINF, Alter_Ältester_inklVN_C, Alter_Jüngster_C_inkl_AlterNutz, ALTERKAU_C, Wuerfel, BEGKUNDE_C, Geoscore_C, GESCHL_VN, +  NUTZART, NUTZKREIS, NICHTK_C_korrigiert, TARIFDAT_C,

Re: [R] significance of character (Anova)

2012-08-10 Thread peter dalgaard
On Aug 10, 2012, at 15:47 , Jennifer Kaiser wrote: Model: poisson, link: log Response: sb_ek_ber Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev Pr(Chi) NULL 1237837 4.4998e+10 ABWHALT_C 2 244353722

Re: [R] Tukey Kramer with ANOVA (glm)

2012-06-15 Thread David Winsemius
On Jun 14, 2012, at 2:17 PM, Alaska_Man wrote: Dr. Winsemius, Really quick, a BACI is a Before-After-Control-Impact approach. I have a long time series of sea cucumber density estimates, which are taken at the same location(s) through time. Some are in areas Impacted by sea otters and

Re: [R] Tukey Kramer with ANOVA (glm)

2012-06-14 Thread David Winsemius
On Jun 13, 2012, at 7:36 PM, Alaska_Man wrote: Hello, I am performing a BACI analysis with ANOVA using the following glm: I admit I had no idea what a BACI analysis might be. Looking it up it appears to be a cross-over design and my statistical betters have sternly warned me about this

Re: [R] Tukey Kramer with ANOVA (glm)

2012-06-14 Thread Alaska_Man
Dr. Winsemius, Really quick, a BACI is a Before-After-Control-Impact approach. I have a long time series of sea cucumber density estimates, which are taken at the same location(s) through time. Some are in areas Impacted by sea otters and some are in areas Not Impacted by sea otters (two

[R] Tukey Kramer with ANOVA (glm)

2012-06-13 Thread Alaska_Man
Hello, I am performing a BACI analysis with ANOVA using the following glm: fit1-glm(log(Cucs_m+1)~(BA*Otter)+BA+Otter+ID+Primary, data=b1) The summary(aov(fit1)) shows significance in the interaction; however, now I would like to determine what combinations of BA and Otter are significantly

Re: [R] Regression Analysis or Anova?

2012-05-16 Thread Robert Latest
Hello Andrea, I don't know if I can help you (probably not, I'm a beginner myself), but you that you should make it a lot easier for those that can if you post a self-contained script in this forum that shows what you're trying to do. Use dput() to dump your dataset in text form. Good luck,

[R] Non-parametric (M)ANOVA with R ?

2012-05-15 Thread Panagiotis Varsamis
Dear R-team, As an extension of my previous work ( http://www.ncbi.nlm.nih.gov/pubmed/21377324) I am now about to handle new obtained ordinal data. Disabled students had to rank their top 5 incentives for participating in sport activities according to their personal (subjective) importance.

[R] Regression Analysis or Anova?

2012-05-15 Thread Andrea Sica
Dear all, I hope to be the clearest I can. Let's say I have a dataset with 10 variables, where 4 of them represent for me a certain phenomenon that I call Y. The other 6 represent for me another phenomenon that I call X. Each one of those variables (10) contains 37 units. Those units are just

Re: [R] How to read ANOVA output

2012-05-07 Thread Robert Baer
Hey all who have responded to this post. I am a newbie to ANOVA analysis in R, and let me tell you- resources for us learners are scant, horrible, unclear, imprecise.. in other words.. the worst ever. So advice like go look it up in your classical textbook or on google is not helpful at all. I am

Re: [R] How to read ANOVA output

2012-05-03 Thread kydaviddoyle
section on ANOVA that you can see at: https://sites.google.com/site/davidsstatistics/using-r/anova-nonparameteric Take Care David - Take Care David Doyle -- View this message in context: http://r.789695.n4.nabble.com/How-to-read-ANOVA-output-tp2329457p4605879.html Sent from the R help mailing list

Re: [R] How to read ANOVA output

2012-05-02 Thread aRghhhhhh
Hey all who have responded to this post. I am a newbie to ANOVA analysis in R, and let me tell you- resources for us learners are scant, horrible, unclear, imprecise.. in other words.. the worst ever. So advice like go look it up in your classical textbook or on google is not helpful at all. I am

Re: [R] How to read ANOVA output

2012-05-02 Thread Ista Zahn
, aRghh sydney.ver...@gmail.com wrote: Hey all who have responded to this post. I am a newbie to ANOVA analysis in R, and let me tell you- resources for us learners are scant, horrible, unclear, imprecise.. in other words.. the worst ever. So advice like go look it up in your classical

Re: [R] How to read ANOVA output

2012-05-02 Thread John Kane
, May 2, 2012 at 1:37 AM, aRghh sydney.ver...@gmail.com wrote: Hey all who have responded to this post. I am a newbie to ANOVA analysis in R, and let me tell you- resources for us learners are scant, horrible, unclear, imprecise.. in other words.. the worst ever. So advice like go look

[R] Nested? Two-way ANOVA with repeated measures

2012-02-02 Thread Vladimir Mikryukov
Dear R-users, I have 3 plant populations (fixed). Within each population there is the same number of “families” (random) – the seed progeny of the same plant. These families were exposed to 2 treatments (fixed) and their response was measured (mean values for 25 seedlings per family per

[R] multcomp two-way anova with interactions within and between

2012-01-11 Thread Schreiber, Stefan
Hi all, I'd like to compare all levels of my interaction with each other. I read the pdf 'Additional multcomp Examples' but even though there is an example with an interaction it doesn't work for me when I want to compare within and between groups. Here is an example:

Re: [R] about interpretation of anova results...

2011-12-05 Thread Roger Koenker
If you were to rtfm you'd see that your anova suggests that the slope coefficients are the same for your tau = .15 and tau = .30 models. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558

[R] about interpretation of anova results...

2011-12-04 Thread narendarreddy kalam
quantreg package is used. *fit1 results are* Call: rq(formula = op ~ inp1 + inp2 + inp3 + inp4 + inp5 + inp6 + inp7 + inp8 + inp9, tau = 0.15, data = wbc) Coefficients: (Intercept) inp1 inp2 inp3 inp4 inp5 -0.191528450 0.005276347 0.021414032

[R] 2^k*r experimental design and anova

2011-11-10 Thread Giovanni Azua
Hello, Can anyone point me to an online tutorial or book containing the easiest way to do ANOVA over the result data from a 2^k*r experiment. It is not clear to me if I can pass the raw data corresponding to each experiment or just the summarized data i.e. mean, sse, std, etc. I would like to

Re: [R] 2^k*r experimental design and anova

2011-11-10 Thread Nicole Marie Ford
, November 10, 2011 7:41:48 AM Subject: [R] 2^k*r experimental design and anova Hello, Can anyone point me to an online tutorial or book containing the easiest way to do ANOVA over the result data from a 2^k*r experiment. It is not clear to me if I can pass the raw data corresponding to each

Re: [R] 2^k*r experimental design and anova

2011-11-10 Thread Kevin Wright
Maybe of interest (uses the FrF2 package) http://www.stat.umn.edu/~gary/classes/5303/handouts/fractional.pdf Kevin On Thu, Nov 10, 2011 at 7:41 AM, Giovanni Azua brave...@gmail.com wrote: Hello, Can anyone point me to an online tutorial or book containing the easiest way to do ANOVA over

[R] r squared and anova for linear mixed-effects model

2011-10-07 Thread florencia bonatto
I have a linear mixed-effects model (from the package nlme) with a random effect; Is there something like an r squared for the whole model which I can state? I´d like to kown: How would I do anova for a linear mixed-effects model? Lic. Florencia BonattoUniversidad Nacional de Rio Cuarto,

Re: [R] r squared and anova for linear mixed-effects model

2011-10-07 Thread Ben Bolker
florencia bonatto mfbonatto at yahoo.com.ar writes: I have a linear mixed-effects model (from the package nlme) with a random effect; Is there something like an r squared for the whole model which I can state? I´d like to kown: How would I do anova for a linear mixed-effects model?

[R] Data simulation for ANOVA decomposition into sums of squares

2011-09-30 Thread Josh B
Dear listserv, Please consider the following dataset: x - matrix(nrow = 8, ncol = 2) colnames(matrix) - c(classification, soluble_fiber) x[1:4,1] - bagel x[5:8,1] - donut How would I simulate a dataset for a one-way fixed-effect ANOVA (where classification is the treatment variable and

Re: [R] When models and anova(model) disagree...

2011-06-16 Thread Ben Bolker
Greg Snow Greg.Snow at imail.org writes: The individual tests on coefficients in logistic regression are generally based on a Wald test statistic. Unfortunately there is a bit of a paradox possible in this case where the coefficient is highly significant, but due to a flattening of the

Re: [R] When models and anova(model) disagree...

2011-06-16 Thread Frank Harrell
Without seeing the model I'm not clear on the cause. In general it is not very useful to compute 1 d.f. tests for parameters that represent only part of a meaningful hypothesis. Note that anova.lrm computes Wald statistics, which are subject to the Hauck-Donner effect (which may not be the cause

[R] When models and anova(model) disagree...

2011-06-15 Thread Rob James
I have a situation where the parameter estimates from lrm identify a binary predictor variable (X) as clearly non-significant (p0.3), but the ANOVA of that same model gives X a chi^2-df rank of 200, and adjudicates X and one interaction of X and a continuous measure as highly significant.

Re: [R] When models and anova(model) disagree...

2011-06-15 Thread Greg Snow
-help@r-project.org Subject: [R] When models and anova(model) disagree... I have a situation where the parameter estimates from lrm identify a binary predictor variable (X) as clearly non-significant (p0.3), but the ANOVA of that same model gives X a chi^2-df rank of 200, and adjudicates X

[R] error term for ANOVA of generalized randomized block design

2011-05-03 Thread Kelsey Ketcheson
Does anyone know how to write the code for an anova for a generalized randomized block design? I have two blocks (random) and three treatments (fixed). Each treatment has two reps at each site. I know that for a RCBD with no replication that an anova can be run using aov(object~Block *Trt), but

Re: [R] is this an ANOVA ?

2011-04-13 Thread Lao Meng
You may try the following to perform anova: anova(lm(y~x)) or summary(aov(y~x)) 2011/4/13 Ubuntu Diego ubuntu.di...@gmail.com Hi all, I have a very easy questions (I hope). I had measure a property of plants, growing in three different substrates (A, B and C). The rest of the

[R] ordinal predictor in anova

2011-04-13 Thread Vincent Staszewski
the ordinal variable (the order is from a priori prediction) #NOW DEFINING THE ORDERED VARIABLE: sugars$treatment_ordered-ordered(sugars$treatment,c(AA, AB,AC, AD,AE)) The problem is that when I run the ANOVA or perform model comparison, R is giving me the same results if I use either

Re: [R] is this an ANOVA ?

2011-04-13 Thread Bill.Venables
are). -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Ubuntu Diego Sent: Wednesday, 13 April 2011 10:10 AM To: r-help@r-project.org Subject: [R] is this an ANOVA ? Hi all, I have a very easy questions (I hope). I had measure

[R] is this an ANOVA ?

2011-04-12 Thread Ubuntu Diego
Hi all, I have a very easy questions (I hope). I had measure a property of plants, growing in three different substrates (A, B and C). The rest of the conditions remained constant. There was very high variation on the results. I want to do address, whether there is any

Re: [R] is this an ANOVA ?

2011-04-12 Thread skfgla...@gmail.com
The response variable y does not make sense to me. What does it represent? friedman.st...@gmail.com 517-648-6290 -Original message- From: Ubuntu Diego ubuntu.di...@gmail.com To: r-help@r-project.org Sent: Wed, Apr 13, 2011 00:12:07 GMT+00:00 Subject: [R] is this an ANOVA ? Hi all

Re: [R] is this an ANOVA ?

2011-04-12 Thread Steven McKinney
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Ubuntu Diego Sent: April-12-11 5:10 PM To: r-help@r-project.org Subject: [R] is this an ANOVA ? Hi all, I have a very easy questions (I hope). I had measure a property

Re: [R] two-way unbalanced ANOVA

2011-02-27 Thread peter dalgaard
On Feb 27, 2011, at 01:12 , Pasicolan wrote: Hello Everyone, *Question: *How do you calculate the sum of squares for a two-way _unbalanced_ ANOVA? There is no unique way of doing it. You can either do a decomposition as f1 : SS(Y~f1) f2 : SS(Y~f1+f2) - SS(Y~f1) res: SStot -

Re: [R] Perform one-way ANOVA using standard deviation and mean

2011-02-08 Thread Greg Snow
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Rao, Niny Sent: Monday, February 07, 2011 7:18 PM To: David Winsemius Cc: r-help@r-project.org Subject: Re: [R] Perform one-way ANOVA using standard deviation and mean On Monday2

[R] Perform one-way ANOVA using standard deviation and mean

2011-02-07 Thread Rao, Niny
Hi! I need to perform ANOVA on a couple of data sets. The only information I have are N, Mean and Standard deviation. I am very new to R, so can someone point me to the right direction on where to go? Thank you very much. Sincerely Niny Rao, PhD Philadelphia University

Re: [R] Perform one-way ANOVA using standard deviation and mean

2011-02-07 Thread David Winsemius
On Feb 7, 2011, at 8:15 PM, Rao, Niny wrote: Hi! I need to perform ANOVA on a couple of data sets. The only information I have are N, Mean and Standard deviation. I am very new to R, so can someone point me to the right direction on where to go? You need to go back and get the data.

Re: [R] Perform one-way ANOVA using standard deviation and mean

2011-02-07 Thread Rao, Niny
On Monday2/7/11 9:09 PM, David Winsemius dwinsem...@comcast.net wrote: On Feb 7, 2011, at 8:15 PM, Rao, Niny wrote: Hi! I need to perform ANOVA on a couple of data sets. The only information I have are N, Mean and Standard deviation. I am very new to R, so can someone point me to the right

Re: [R] Perform one-way ANOVA using standard deviation and mean

2011-02-07 Thread Dennis Murphy
Hi: You need to check that you have the sufficient statistics necessary to obtain an ANOVA table corresponding to the model you intend to fit. If you have that, you can use the vectorization and matrix capabilities in R to manually get the ANOVA table, perform the F tests, etc. The only thing

Re: [R] Perform one-way ANOVA using standard deviation and mean

2011-02-07 Thread RICHARD M. HEIBERGER
Niny Rao, It is possible to get the ANOVA table from the sufficient statistics. The sufficient statistics include the sample sizes which were not included in the data you posted. I have fake sample sizes n in this example. You can use the aov.sufficient function in the HH package for this.

Re: [R] t-test or ANOVA...who wins? Help please!

2011-01-05 Thread Tal Galili
Hello Frodo, It is not clear to me from your questions some of the basics of your analysis. If you only have two levels of a factor, and one response - why in the anova do you use more factors (and their interactions)? In that sense, it is obvious that your results would differ from the t-test.

Re: [R] t-test or ANOVA...who wins? Help please!

2011-01-05 Thread Frodo Jedi
and understand. Please enlighten me. Thanks in advance Best regards From: Tal Galili tal.gal...@gmail.com To: Frodo Jedi frodo.j...@yahoo.com Cc: r-help@r-project.org Sent: Wed, January 5, 2011 10:15:41 AM Subject: Re: [R] t-test or ANOVA...who wins? Help

[R] t-test or ANOVA...who wins? Help please!

2011-01-04 Thread Frodo Jedi
Dear all, I need an help because I don´t know how to perform the analysis in the right way, as I get different beheaviors using t-test and two ways ANOVA. In what follow I post the table, my goal and the strange results I got. I kindly ask you an help because I really don´t know how to solve

Re: [R] Doing a mixed-ANOVA after accounting for a covariate

2010-12-27 Thread Dror D Lev
the ones in Faraway's Practical regression and anova using r) use lm(), but this doesn't handle repeated-measures. I thought of a solution in the form of first running a regression on the covariate: cov.accnt = lm (myMeasure ~ myCovMeasure, data=dat) and then run the aov() on the residuals

[R] Doing a mixed-ANOVA after accounting for a covariate

2010-12-26 Thread Dror D Lev
I noticed for handling covariates (i.e. ANCOVA, including the ones in Faraway's Practical regression and anova using r) use lm(), but this doesn't handle repeated-measures. I thought of a solution in the form of first running a regression on the covariate: cov.accnt = lm (myMeasure ~ myCovMeasure

Re: [R] Doing a mixed-ANOVA after accounting for a covariate

2010-12-26 Thread David Winsemius
to require treating practice as a covariate. All the examples I noticed for handling covariates (i.e. ANCOVA, including the ones in Faraway's Practical regression and anova using r) use lm(), but this doesn't handle repeated-measures. See if Dalgaard's piece in R-News offers better guidance: http

Re: [R] Doing a mixed-ANOVA after accounting for a covariate

2010-12-26 Thread Dror D Lev
, after accounting for any variance due to practice (31 trials in each of two blocks) in the task. It seems to require treating practice as a covariate. All the examples I noticed for handling covariates (i.e. ANCOVA, including the ones in Faraway's Practical regression and anova using r) use

Re: [R] Doing a mixed-ANOVA after accounting for a covariate

2010-12-26 Thread David Winsemius
groups? -- David. All the examples I noticed for handling covariates (i.e. ANCOVA, including the ones in Faraway's Practical regression and anova using r) use lm(), but this doesn't handle repeated-measures. See if Dalgaard's piece in R-News offers better guidance: http://www.r

Re: [R] Doing a mixed-ANOVA after accounting for a covariate

2010-12-26 Thread RICHARD M. HEIBERGER
for handling covariates (i.e. ANCOVA, including the ones in Faraway's Practical regression and anova using r) use lm(), but this doesn't handle repeated-measures. I thought of a solution in the form of first running a regression on the covariate: cov.accnt = lm (myMeasure ~ myCovMeasure, data=dat

Re: [R] repeatedrepeated measures in ANOVA or mixed model

2010-11-10 Thread Mike Marchywka
Date: Tue, 9 Nov 2010 18:25:18 -0800 From: djmu...@gmail.com To: paul.rhee...@up.ac.za CC: r-help@r-project.org Subject: Re: [R] repeatedrepeated measures in ANOVA or mixed model Hi: This sounds like a 'doubly repeated measures problem

[R] repeatedrepeated measures in ANOVA or mixed model

2010-11-09 Thread Paul Rheeder
dear List I have a dataset with blood measurements at 5 points in TIME (0,30,60,90,120) taken on 3 VISITS (same subjects). the interest is to compare these measurements between Visits, overall and at the different time points. I have problems setting up repeated measures ANOVA with 2 repeated

Re: [R] repeatedrepeated measures in ANOVA or mixed model

2010-11-09 Thread Tal Galili
Hi Paul, Here are a bunch of tutorials on the topic: http://www.r-statistics.com/2010/04/repeated-measures-anova-with-r-tutorials/ I suggest a more specific question for people to help :) Best, Tal Contact Details:--- Contact

Re: [R] repeatedrepeated measures in ANOVA or mixed model

2010-11-09 Thread Dennis Murphy
Hi: This sounds like a 'doubly repeated measures problem'. Are any treatments assigned to individuals or is this a purely observational study? Is the time horizon of the between-visit factor (much?) longer than that of the within-visit factor? You could try to assess the strength of correlation

[R] unbalanced repeated measurements Anova with mixed effects

2010-10-17 Thread Thijs Muizelaar
Dear R-list members, I've been struggling with the proper setup for analysing my data. I've performed a route choice experiment, in which participants had to make a choice at each junction for the next road. During the experiment they received traffic information, but also encountered two

  1   2   3   >