[R] Select components of a list

arun smartpink111 at yahoo.com
Sun Feb 17 02:55:20 CET 2013

```Hi Gustav,

Try this:

lapply(1:length(models),function(i) lapply(models[[i]],function(x) summary(x)\$coef[2,]))[[1]] #1st list component
[[1]]
#    Estimate   Std. Error      z value     Pr(>|z|) # pm10
#5.999185e-04 1.486195e-04 4.036606e+00 5.423004e-05

#[[2]]
#    Estimate   Std. Error      z value     Pr(>|z|) #ozone
#0.0010117294 0.0003792739 2.6675428048 0.0076408155

#[[3]]
#    Estimate   Std. Error      z value     Pr(>|z|) #so2
#0.0026595441 0.0009352046 2.8438097399 0.0044577658

#pm10
do.call(rbind,lapply(1:length(models),function(i) do.call(rbind,lapply(models[[i]],function(x) summary(x)\$coef[row.names(summary(x)\$coef)=="pm10",c(1:2,4)]))))
#         Estimate   Std. Error     Pr(>|z|)
#[1,] 0.0005999185 0.0001486195 5.423004e-05
#[2,] 0.0005720549 0.0001740368 1.012696e-03
#[3,] 0.0005099552 0.0001559620 1.076462e-03
#[4,] 0.0009285593 0.0001766520 1.468764e-07

#ozone

do.call(rbind,lapply(1:length(models),function(i) do.call(rbind,lapply(models[[i]],function(x) summary(x)\$coef[row.names(summary(x)\$coef)=="ozone",c(1:2,4)]))))
#        Estimate   Std. Error     Pr(>|z|)
#[1,] 0.0010117294 0.0003792739 7.640816e-03
#[2,] 0.0009128304 0.0004364390 3.647954e-02
#[3,] 0.0023896044 0.0004109854 6.087769e-09
#[4,] 0.0005455392 0.0004301502 2.047076e-01
#so2

do.call(rbind,lapply(1:length(models),function(i) do.call(rbind,lapply(models[[i]],function(x) summary(x)\$coef[row.names(summary(x)\$coef)=="so2",c(1:2,4)]))))
#        Estimate   Std. Error    Pr(>|z|)
#[1,] 0.002659544 0.0009352046 0.004457766
#[2,] 0.002825612 0.0010150314 0.005373144
#[3,] 0.002409738 0.0009563814 0.011747444
#[4,] 0.001725140 0.0011635156 0.138155175

Hope it helps.
A.K.
________________________________
From: Gustav Sigtuna <gsigtuna at gmail.com>
To: arun <smartpink111 at yahoo.com>
Sent: Saturday, February 16, 2013 7:44 PM
Subject: Re: Select components of a list

Hi Arun,

Thanks for taking your time to find a solution.

I have attached a R script that will recreate a comparable list from publicly available data. My list is longer and created by various models  than the one created here. However, the final output is similar to the one produced by script. My interest is to extract only the coefficients for  pm10., ozone and so2 (  Estimate,  Std. Error  and p value)    .

Thanks

On Fri, Feb 15, 2013 at 9:04 PM, arun <smartpink111 at yahoo.com> wrote:

Dear Gustav,
>Thank you for the data.  Could you select a smaller subset of the list and dput() that subset?  Your data is useful, but I would have to recreate list of lists from that to test and sometimes that may not accurate represent the format in your list as it is the summary().
>Arun
>
>
>
>
>
>
>
>________________________________
>From: Gustav Sigtuna <gsigtuna at gmail.com>
>To: smartpink111 at yahoo.com
>Sent: Friday, February 15, 2013 4:56 AM
>Subject: Re: Select components of a list
>
>
>
>Hi Arun,
>
>Thanks for your help. Your mail landed in my spam folder and just saw it by chance.
>
>I have attached a text file that contains the list of my model. It was extremely long, thus I took out the last part which is think is more important.
>
>
>In brief I have an output from GAM model which resulted from analysis of  ozone at three time points on 12 data sets
>
>
>
>
>
>
>On Wed, Feb 13, 2013 at 8:21 PM, <smartpink111 at yahoo.com> wrote:
>
>Dear Lungo,
>>
>>If you can email (smartpink111 at yahoo.com) me the `list` (dput(list)), I can take a look at it.  Probably, you understand that my previous solution was just guesswork. With regards to GLM, GAM, it is good to check the structure of the list (str()). It gives information about whether a `generic` tool could be applied to extract them or not.
>>Cheers.
>>Arun
>>
>><quote author='Lungo'>
>>Dear Arun,  Your code and the example works fine. However my list is quite
>>different from the one showed in your example.   As I have shown in my
>>question above I have 12 lists each having 3 lists underneath.  I get the
>>lists by different models (GLM, GAM ) but the output  I aim to have is the
>>estimates of the explanatory variable which is placed next to the intercept.
>>Thus I am looking for a “generic” tool that would extract these lists.
>></quote>
>>Quoted from:
>>http://r.789695.n4.nabble.com/Select-components-of-a-list-tp4658295p4658389.html
>>
>

```