Openair scatterplot: Table of linear regression coefficients

[kenreid4]

Hello. thank you for the openair project.

Is it possible to obtain a table of linear coefficients , created with the scatterPlot - , linear = TRUE function.

If not that would be a great future enhancement.

Otherwise are there any other R libraries or functions that provide such a table as an output option please.

Thank you
Regards
Ken

[jack-davison]

Hi Ken,

Thanks for getting in touch, and sorry for the delay - I didn't notice your issue on the openair book repository.

At the moment one cannot extract the linear coefficients simply from a scatterPlot() function.

This may be something we can look at; its possible they could be returned in a similar way to the parameters in a theilSen() plot.

Best,
Jack

[kenreid4]

Hi Jack

Thank you very much for your response.

Yes being able to extract the coefficients from the scatterPlot/corPlot functions would help a lot in our local pm2.5 project.

---

Meanwhile I found this discussion of some help.

https://stackoverflow.com/questions/30371516/how-to-save-summarylm-to-a-file

"The suggestion above work great. Depending what you need you can use the tidy() function for the coefficients and glance() for the table.

library( broom )
a <- lm(cars$speed ~ cars$dist)
write.csv( tidy( a ) , "coefs.csv" )
write.csv( glance( a ) , "an.csv" )"

---

Using Broom::glance helped me.

Thank you again for openair.

Best regards
Ken

[jack-davison]

So for corPlot() you definitely can already extract the correlation coefficients, if you'd like:

corplot <- openair::corPlot(openair::mydata)
corplot$data

# A tibble: 81 × 6
   x     y     default                           row   col     cor
   <chr> <chr> <fct>                           <int> <int>   <dbl>
 1 pm10  pm10  01 January 1998 to 23 June 2005     1     1  1     
 2 pm25  pm10  01 January 1998 to 23 June 2005     2     1  0.836 
 3 no2   pm10  01 January 1998 to 23 June 2005     3     1  0.575 
 4 so2   pm10  01 January 1998 to 23 June 2005     4     1  0.490 
 5 nox   pm10  01 January 1998 to 23 June 2005     5     1  0.617 
 6 co    pm10  01 January 1998 to 23 June 2005     6     1  0.455 
 7 o3    pm10  01 January 1998 to 23 June 2005     7     1 -0.284 
 8 ws    pm10  01 January 1998 to 23 June 2005     8     1  0.0186
 9 wd    pm10  01 January 1998 to 23 June 2005     9     1 -0.0799
10 pm10  pm25  01 January 1998 to 23 June 2005     1     2  0.836 
# ℹ 71 more rows
# ℹ Use `print(n = ...)` to see more rows

The one thing you don't get here is the order of the x/y axes if cluster = TRUE; might be another useful feature for me to consider.

Best,
Jack

[kenreid4]

Hi Jack

Thank you very much, I will give the "corplot$data" output a try in the next few days.

Best
Ken