[R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread V. Coudrain
Hi, I would like to test the impact of a treatment of some variable using regression (e.g. lm(var ~ trt + cov)).  However I only have four observations per factor level. Is it still possible to apply a regression with such a small sample size. I think that i should be difficult to correctly

Re: [R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread V. Coudrain
Thank you very much. If I get it right, the CI get wider, my test has less power and the probability of getting a significant relation decreases. What about the significant coefficients, are they reliable? Message du 20/10/14 à 11h30 De : Roman Luštrik A : V. Coudrain Copie à :

Re: [R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread Martin Weiser
Hi, coefficients and their p-values are reliable if your data are OK and you do know enough about the process that generated them, so you can choose appropriate model. With 4 points per line, it may be really difficult to identify bad fit or outliers. For example: simple linear regression needs

[R-sig-eco] Logistic regression with 2 categorical predictors

2014-10-20 Thread Andrew Halford
Hi Listers, I am trying to run a logistic regression to look at the effects of experiment type and age on the behavior of fish in a choice chamber experiment. I am using the glm approach and would like some advice on how or whether to perform contrasts to work out what levels of Factor1 (Age)

Re: [R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread stephen sefick
You are more or less preforming an ANOVA/ANCOVA on your data? As pointed out earlier, all of the normal theory regression assumptions apply. Assuming all of those things are satisfied then if you have large confidence intervals and there are significant differences between groups I don't see why

Re: [R-sig-eco] Logistic regression with 2 categorical predictors

2014-10-20 Thread ONKELINX, Thierry
Dear Andrew, anova() and summary() test different hypotheses. anova() tests is at least one level is different from the others. summary() tests if the coefficient is different from zero. Multiple comparison of different interaction levels is probably the most relevant in this case. The

Re: [R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread V. Coudrain
Thank you for this helpful thought. So if I get it correctly it is hopeless to try testing an interaction, but we neverless may assess if a covariate has an impact, providing it is the same in all treatments. Message du 20/10/14 à 16h46 De : Elgin Perry A : v_coudr...@voila.fr Copie à :

Re: [R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread V. Coudrain
Yes, but as I fear, the residuals behave badly as soon as the model get a little bit more complex (e.g., with two covariables or an interactions). The scope for performing an ANCOVA is thus very limited. That's why I was thinking about a potential non-parametric model. But I do not want to

Re: [R-sig-eco] Regression with few observations per factor level

2014-10-20 Thread Baldwin, Jim -FS
Yes, the analysis with a small sample size would be valid (under the assumption that the model - both fixed and random effects are correctly specified) but at some point there must be a practical assessment as to the desired precision and the costs of the consequences of either inadequate