[R] In fact this is a Stats question, but...

From: <eesteves_at_ualg.pt>
Date: Thu, 29 May 2008 14:00:24 +0200


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_at_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.
Received on Thu 29 May 2008 - 13:18:24 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 29 May 2008 - 14:30:45 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive