Dear All,
I'me having (much) trouble understanding why it happened and answering a referee's comment to part of a submitted manuscript. I've tried to google for help but... I'm really confident that although this is a R-Help list someone can help me!

I used R to do an ANCOVA w/ RNA/DNA as the dep var, sl as the indep var and gut (a factor w/ levels: prey and empty) as the covariate:

RNADNA.sl.gut<-lm(sqrt(RNADNA)~gut*sl,subset=gut!="Yolk-sac",data=cond)
summary(RNADNA.sl.gut)

The results from this are:

(...)
Coefficients:
             Estimate Std. Error t value Pr(>|t|)
(Intercept)  0.856266   0.052252  16.387  < 2e-16 ***
gutPrey     -0.009568   0.092170  -0.104    0.917
sl           0.030575   0.004648   6.578 6.35e-11 ***
gutPrey:sl   0.002285   0.007313   0.313    0.755
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 0.3312 on 1692 degrees of freedom
Multiple R-Squared: 0.05847,    Adjusted R-squared: 0.0568
F-statistic: 35.02 on 3 and 1692 DF,  p-value: < 2.2e-16

(...)

The question raised by referee is related to the "incompatibility" of the low r2 (0.057) and the high significance (p<<0.0001) of the model. I've interpreted/used this result in the following way: although there's a significant relationship between RNA/DNA and sl, it's very weak; besides, no gut effect on the relationship as been found!

Sorry for the off-topic question but...

Sincerely, Eduardo Esteves

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to