Hi,

mods <- lapply(lapply(df[,which(sapply(df, is.factor))],
function(reg)lm(df$value~reg)), summary)
res <- lapply(mods, "[", c(4,10))

And for each model adjusted:
1-pf(res[[2]][[2]][1], res[[2]][[2]][2], res[[2]][[2]][3])
1-pf(res[[1]][[2]][1], res[[1]][[2]][2], res[[1]][[2]][3])




On 09/10/2007, Jiong Zhang, PhD <jizhang@chori.org> wrote:
>
> Hi,
>
> May I ask how I can save the coefficients and the p values into a table?
>
> thanks.
>
> jiong
> The email message (and any attachments) is for the sol...{{dropped:20}}

