[R] how to parse out fitting statistics and write them into a data frame?

weix1 weix1_2007 at hotmail.com
Thu Apr 29 17:36:21 CEST 2010


hello, everyone: 
I am conducting t test between drug and control for about 50,000 gene using
the following syntax (treatment is factor): 

result<- lapply(split(data, data$gene),function(x) lm(value~treatment,x) 
  
however, the result is a list and i do not know whether more model fitting
statistics (like p value of t test) is included in "result" or not. If i
print the first element of resut i got the followings: 

> result[1] 
$`1007_s_at` 

Call: 
lm(formula = logvalue ~ treatment, data = x) 

Coefficients: 
 (Intercept)  treatmentveh   
      8.9403        0.3232   

> summary(result[1]) 
          Length Class Mode 
1007_s_at 13     lm    list 
> 

So my question is whether more fitting statistics (other than coefficient
estimation, like p value) are included in the "result". If yes, how can I
parse them into a data frame so that i can output those statistics into a
.csv file that can be shared with my clients. If not, how can I modify the
code so that more stat can be computed and stored? 
any constructive suggestions are welcome. 

-- 
View this message in context: http://r.789695.n4.nabble.com/how-to-parse-out-fitting-statistics-and-write-them-into-a-data-frame-tp2075707p2075707.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list