[R] Advice on starting to analyze smokestack emissions?
Kevin, Sorry, was not aware of any access-restriction, but if you have problems and were interested in local AQ maybe https://www.epa.gov/outdoor-air-quality-data Karl -- [EXTERNAL] Message: 3 Date: Fri, 15 Dec 2023 13:29:36 + From: Kevin Zembower To: R-help email list Subject: Re: [R] Advice on starting to analyze smokestack emissions? Message-ID: <0100018c6dab10c0-a6da3643-a9f2-44af-8a69-9303ea534749-000...@email.amazonses.com> Content-Type: text/plain; charset="utf-8" Bert, Tim, Karl and Richard, thank you all for your suggestions and help. I will try the R-sig-ecology list. Karl, I wasn't aware of the RAQSAPI package, but it looked promising. However, when I went to the source of the data it uses, the United States Environmental Protection Agency’s (US EPA) Air Quality System (AQS) Data Mart database, it looks like interactive access to the data is restricted to those who can document a professional agency affiliation. I don't have that. I'll work with the package to see if this is true regarding obtaining the data through it. Thanks for the suggestion. Richard, the Canada study of crematoriums was very useful. Thanks. Thanks, again, all, for your help. -Kevin You might also try the R-Sig-ecology list, though I would agree that it's not clearly related. Still, air pollution effects...? -- Bert On Tue, Dec 12, 2023 at 3:15 AM Kevin Zembower via R-help < r-help@r-project.org> wrote: > Hello, all, > > [Originally sent to r-sig-geo list, with no response. Cross-posting > here, in the hope of a wider audience. Anyone with any experience in > this topic? Thanks.] > > I'm trying to get started analyzing the concentrations of smokestack > emissions. I don't have any professional background or training for > this; I'm just an old, retired guy who thinks playing with numbers is > fun. > > A local funeral home in my neighborhood (less than 1200 ft from my > home) is proposing to construct a crematorium for human remains. I > have some experience with the tidycensus package and thought it might > be interesting to construct a model for the changes in concentrations > of the pollutants from the smokestack and, using recorded wind speeds > and directions, see which US Census blocks would be affected. > > I have the US Government EPA SCREEN3 output on how concentration > varies with distance from the smokestack. > See > https://www/. > epa.gov%2Fscram%2Fair-quality-dispersion-modeling-screening-models%23s > creen3=05%7C02%7Ctebert%40ufl.edu%7C3097c182143c47a6789c08dbfb2a7 > ed2%7C0d4da0f84a314d76ace60a62331e1b84%7C0%7C0%7C638379932467260671%7C > Unknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1h > aWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=QgsYQ9w28caBmEGwJ9Kei2x0fSkH3 > 4v3%2BfAo37GdcYQ%3D=0 if curious. As a first task, I'd like > to see if I can calculate similar results in R. I'm aware of the > 'plume' steady-state Gaussian dispersion package > (https://rdr/ > r.io%2Fgithub%2Fholstius%2Fplume%2Ff%2Finst%2Fdoc%2Fplume-intro.pdf=05%7C02%7Ctebert%40ufl.edu%7C3097c182143c47a6789c08dbfb2a7ed2%7C0d4da0f84a314d76ace60a62331e1b84%7C0%7C0%7C638379932467260671%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=DN9oxiJnDFvvmY968G9t9Sagr8UfJ2ySZiGWV1%2F9AC8%3D=0), > but am a little concerned that this package was last updated 11 years ago. > > Do you have any recommendations for me on how to get started analyzing > this problem? Is 'plume' still the way to go? I'm aware that there are > many atmospheric dispersion models from the US EPA, but I was hoping > to keep my work within R, which I'm really enjoying using and learning > about. Are SCREEN3 and 'plume' comparable? Is this the best R list to > ask questions about this topic? > > Thanks for any advice or guidance you have for me. > > -Kevin > > > > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat/ > .ethz.ch%2Fmailman%2Flistinfo%2Fr-help=05%7C02%7Ctebert%40ufl.edu > %7C3097c182143c47a6789c08dbfb2a7ed2%7C0d4da0f84a314d76ace60a62331e1b84 > %7C0%7C0%7C638379932467260671%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAw > MDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C > ta=dxsuLWVRx8wNnu49SJ34AAh7oRECvDIrQh9%2Bpx48SL0%3D=0 > PLEASE do read the posting guide > http://www.r/ > -project.org%2Fposting-guide.html=05%7C02%7Ctebert%40ufl.edu%7C30 > 97c182143c47a6789c08dbfb2a7ed2%7C0d4da0f84a314d76ace60a62331e1b84%7C0% > 7C0%7C638379932467260671%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiL > CJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=QY > AiKA8xDhcPyQmRZ6Vqcr5mdszE8WSRyFmCqzQ7Rog%3D=0 > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,
Re: [R] R-help Digest, Vol 250, Issue 13
Kevin, Maybe also look at what air quality monitoring is being done in area. https://cran.r-project.org/web/packages/RAQSAPI/vignettes/RAQSAPIvignette.html Depends what and how near, but might be something relevant there? Karl Dr Karl Ropkins Transport Studies | Environment | University of Leeds -- Message: 2 Date: Tue, 12 Dec 2023 07:52:59 -0800 From: Bert Gunter To: Kevin Zembower Cc: R-help email list Subject: Re: [R] Advice on starting to analyze smokestack emissions? Message-ID: Content-Type: text/plain; charset="utf-8" You might also try the R-Sig-ecology list, though I would agree that it's not clearly related. Still, air pollution effects...? -- Bert On Tue, Dec 12, 2023 at 3:15 AM Kevin Zembower via R-help < r-help@r-project.org> wrote: > Hello, all, > > [Originally sent to r-sig-geo list, with no response. Cross-posting > here, in the hope of a wider audience. Anyone with any experience in > this topic? Thanks.] > > I'm trying to get started analyzing the concentrations of smokestack > emissions. I don't have any professional background or training for > this; I'm just an old, retired guy who thinks playing with numbers is > fun. > > A local funeral home in my neighborhood (less than 1200 ft from my > home) is proposing to construct a crematorium for human remains. I have > some experience with the tidycensus package and thought it might be > interesting to construct a model for the changes in concentrations of > the pollutants from the smokestack and, using recorded wind speeds and > directions, see which US Census blocks would be affected. > > I have the US Government EPA SCREEN3 output on how concentration varies > with distance from the smokestack. > See > https://www.epa.gov/scram/air-quality-dispersion-modeling-screening-models#screen3 > if curious. As a first task, I'd like to see if I can calculate similar > results in R. I'm aware of the 'plume' steady-state Gaussian dispersion > package > (https://rdrr.io/github/holstius/plume/f/inst/doc/plume-intro.pdf), but > am a little concerned that this package was last updated 11 years ago. > > Do you have any recommendations for me on how to get started analyzing > this problem? Is 'plume' still the way to go? I'm aware that there are > many atmospheric dispersion models from the US EPA, but I was hoping to > keep my work within R, which I'm really enjoying using and learning > about. Are SCREEN3 and 'plume' comparable? Is this the best R list to > ask questions about this topic? > > Thanks for any advice or guidance you have for me. > > -Kevin > > > > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. > [[alternative HTML version deleted]] -- Message: 3 Date: Tue, 12 Dec 2023 21:19:12 + (UTC) From: varin sacha To: "r-help@r-project.org" , Ben Bolker Subject: Re: [R] ggplot2: Get the regression line with 95% confidence bands Message-ID: <68588390.888662.1702415952...@mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Dear Ben, Dear Daniel, Dear Rui, Dear Bert, Here below my R code. I really appreciate all your comments. My R code is perfectly working but there is still something I would like to improve. The X-axis is showing 2012.5 ; 2015.0 ; 2017.5 ; 2020.0 I would like to see on X-axis only the year (2012 ; 2015 ; 2017 ; 2020). How to do? # library(ggplot2) df=data.frame(year= c(2012,2015,2018,2022), score=c(495,493, 495, 474)) ggplot(df, aes(x = year, y = score)) + geom_point() + geom_smooth(method = "lm", formula = y ~ x) + labs(title = "Standard linear regression for France", x = "Year", y = "PISA score in mathematics") + scale_y_continuous(limits=c(470,500),oob=scales::squish) # Le lundi 11 décembre 2023 à 23:38:06 UTC+1, Ben Bolker a écrit : On 2023-12-11 5:27 p.m., Daniel Nordlund wrote: > On 12/10/2023 2:50 PM, Rui Barradas wrote: >> Às 22:35 de 10/12/2023, varin sacha via R-help escreveu: >>> >>> Dear R-experts, >>> >>> Here below my R code, as my X-axis is "year", I must be missing one >>> or more steps! I am trying to get the regression line with the 95% >>> confidence bands around the regression line. Any help would be >>> appreciated. >>> >>> Best, >>> S. >>> >>> >>>
Re: [R] readline in function call with space in prompt.
Not on a machine with latest R at moment so not ruling out something there, but it is working fine for me. Karl R version 4.0.2 (2020-06-22) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 10 x64 (build 19041) -- Message: 9 Date: Tue, 9 Feb 2021 10:21:58 +1300 From: Rolf Turner To: Jeremie Juste Cc: r-help@r-project.org Subject: Re: [R] readline in function call with space in prompt. Message-ID: <20210209102158.5b04117d@rolf-Latitude-E7470> Content-Type: text/plain; charset="windows-1252" On Mon, 08 Feb 2021 14:28:33 +0100 Jeremie Juste wrote: > Hello, > > I have noticed a behavior that I don't understand. When I call the > following function from the prompt. > test <- function(){ > a <- readline("selection: ") > a > } > > > test() > > selection: | > I can only type one character and the readline function exits before > I can press enter. > > however > > test1 <- function(){ > a <- readline("selection:") > a > } > > test1() > > selection:| > works as expected. > > selection: abc[Ret] > > However calling directly readline with a space in the prompt does > what I would expect. > > > a <- readline("selection: ") > > selection: abc[Ret] > > a > > "abc" > > It is the expected behavior or am I missing something? Works fine for me, with or without the space in the prompt string. There *must* be something flaky in your system, but I'm damned if I can come up with any useful suggestions for tracking down just where that flakiness lies. Sorry. Perhaps try re-installing R??? cheers, Rolf Turner -- Honorary Research Fellow Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.
Re: [R] using edit to extract codes from vignette failed
FYI: I get the same for Windows 7 and 8, no RStudio: edit(vignette(grobs,package = grid)) works for me edit(file=vignette(grobs,package = grid)) does not edit(file=vignette(grobs,package = grid)) Error in tempfile(name$topic, fileext = .R) : argument name is missing, with no default edit(vignette(grobs,package = grid)) Re RStudio: https://support.rstudio.com/hc/en-us?community_id=public? From: PO SU [rhelpmaill...@163.com] Sent: 09 September 2014 03:07 To: William Dunlap Cc: Karl Ropkins; R. Help Subject: Re:Re: [R] using edit to extract codes from vignette failed OK, i get it, i should set the editor argument , i don't know how to report a bug to Rstudio, may you do that ? -- PO SU mail: desolato...@163.com Majored in Statistics from SJTU At 2014-09-09 00:41:33, William Dunlap wdun...@tibco.com wrote: Complain to the RStudio people - RStudio defines its own options(editor) which is not completely compatible with R's option(editor=internal). If you set options(editor=internal) in RStudio then you can look at the code in the vignette. (I tried with last year's RStudio 0.98.501 and this may have been fixed by now.) Bill Dunlap TIBCO Software wdunlap tibco.com On Mon, Sep 8, 2014 at 9:21 AM, PO SU rhelpmaill...@163.com wrote: Tks for correcting me not using the file argument, but the codes you supply seem still not work. edit(vignette(grobs,package = grid)) can't work. I am using win7, the latest version of Rstudio which using R.3.1.1.The error is: Error in editor(file = file, title = title) : argument name is missing, with no default -- PO SU mail: desolato...@163.com Majored in Statistics from SJTU At 2014-09-08 05:46:39, Karl Ropkins k.ropk...@its.leeds.ac.uk wrote: Try: edit(vignette(grobs,package = grid)) (edit is a method. It looks at the class of the first entry, name, to identify which method to use. See ?edit. You want it to use edit.vignette, so you need to drop 'file=' so you pass the vignette to edit as the first argument or name=. Then edit will pass it to edit.vignette and it'll work. Or go direct: edit.vignette(vignette(grobs,package = grid)). See ?vignette. Maybe the use of name as the first argument of a method is a little misleading? But you can work out what is going if you work through the help documentation.) Karl Message: 9 Date: Sun, 7 Sep 2014 17:06:44 +0800 (CST) From: PO SU rhelpmaill...@163.com To: R. Help r-help@r-project.org Subject: [R] using edit to extract codes from vignette failed Message-ID: 4d3c1c8a.1c96.1484f5d8d31.coremail.rhelpmaill...@163.com Content-Type: text/plain; charset=UTF-8 Dear expeRts, ? ? When i using the following code, i get a error as follows: ?edit(file=vignette(grobs,package = grid)) Error in edit.vignette(file = vignette(grobs, package = grid)) :? ? argument name is missing, with no default I investigated edit function, but still can't ?get codes from a vignette, May you help me? -- PO SU mail: desolato...@163.com Majored in Statistics from SJTU __ R-help@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. __ R-help@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. __ R-help@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.
[R] using edit to extract codes from vignette failed
Try: edit(vignette(grobs,package = grid)) (edit is a method. It looks at the class of the first entry, name, to identify which method to use. See ?edit. You want it to use edit.vignette, so you need to drop 'file=' so you pass the vignette to edit as the first argument or name=. Then edit will pass it to edit.vignette and it'll work. Or go direct: edit.vignette(vignette(grobs,package = grid)). See ?vignette. Maybe the use of name as the first argument of a method is a little misleading? But you can work out what is going if you work through the help documentation.) Karl Message: 9 Date: Sun, 7 Sep 2014 17:06:44 +0800 (CST) From: PO SU rhelpmaill...@163.com To: R. Help r-help@r-project.org Subject: [R] using edit to extract codes from vignette failed Message-ID: 4d3c1c8a.1c96.1484f5d8d31.coremail.rhelpmaill...@163.com Content-Type: text/plain; charset=UTF-8 Dear expeRts, ? ? When i using the following code, i get a error as follows: ?edit(file=vignette(grobs,package = grid)) Error in edit.vignette(file = vignette(grobs, package = grid)) :? ? argument name is missing, with no default I investigated edit function, but still can't ?get codes from a vignette, May you help me? -- PO SU mail: desolato...@163.com Majored in Statistics from SJTU __ R-help@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.
Re: [R] How to draw Bubble chart with mini pie charts as bubbles
Dear Amir, There are a couple of panel functions in loa, a package of lattice add-ins, that might be of help: require(loa) # Use a subsample of lat.lon.meuse dataset in loa temp - lat.lon.meuse[sample(1:155, 15),] # plot Cu/Pb/Zn pie plots at sampling locations loaPlot(copper+lead+zinc~longitude*latitude, panel=panel.zcasePiePlot, data=temp) # then rescale smaller pie segments on the fly loaPlot(copper*10+lead*4+zinc~longitude*latitude, panel=panel.zcasePiePlot, data=temp) Best Wishes, Karl -- Message: 5 Date: Sat, 14 Jun 2014 13:03:21 -0700 (PDT) From: Agony agony_...@yahoo.com To: r-help@r-project.org Subject: [R] How to draw Bubble chart with mini pie charts as bubbles in R Message-ID: 1402776201.7516.yahoomailba...@web120401.mail.ne1.yahoo.com Content-Type: text/plain; charset=us-ascii Dear all, Good day! Could anybody help me how to draw a bubble chart with mini pie charts as bubbles in R ? Introducing any experiences, books, booklet or source code will appreciated. Bunch of thanks. Best, Amir __ R-help@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.
Re: [R] How do I perform conditional annotation of lattice
Not sure if this adds much to Ken Knoblauch previous suggestion. But: Subject-c(1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,4,4,4,4,4,4) Day-c(1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6) Activity-c(2,3,4,3,7,4,5,8,2,8,4,6,2,5,3,8,9,5,6,3,4,5,6,7) EventA-c(Yes,NA,Yes,NA,NA,NA,Yes,NA,NA,NA,NA,NA, Yes,NA,NA,NA,NA,NA,NA,NA,Yes,NA,NA,NA) require(loa) panel.temp = function(x=x, y=y, z=z, ...){ panel.xyplot(x, y, z =1, ...) panel.abline(v=x[!is.na(z)], col=red) } loaPlot(EventA ~ Day*Activity|Subject, panel=panel.temp, type=b) The loa package is still a little messy, not yet on CRAN and the help is only half written but you are welcome to have a play: https://sites.google.com/site/karlropkins/rpackages/loa The basic plot is z ~ x * y, where z is information to be displayed at (x,y) and the package includes a very inelegant alternative to subscripting. Regards, Karl Ropkins ITS Leeds -- Date: Tue, 7 Jan 2014 15:50:39 + From: Ken Knoblauch ken.knobla...@inserm.frmailto:ken.knobla...@inserm.fr To: r-mailto:r-h...@stat.math.ethz.chh...@stat.math.ethz.ch Subject: Re: [R] How do I perform conditional annotation of lattice panel plots? Message-ID: loom.20140107t164813-...@post.gmane.orgmailto:loom.20140107t164813-...@post.gmane.org Content-Type: text/plain; charset=utf-8 Chen, George George.Chen at roswellpark.org writes: My apologies for asking this question that may have been asked before. I am trying to plot activity dependent on time conditioned by the subject. Code for sample data below. So I have something like this xyplot(Activity~Time|Subject). This works fine, but now I want to show where on these activity curves Event A occurs. This is to explore the relationship between A and activity. I tried this: xyplot(EventA+Activity~Time|Subject) but then the numerical ordering gets out of order. Ideally Event A could be denoted by a vertical line cutting through the curve at the day Event A occurred, but some other way to denote it would be great! -snip - Any ideas on how to do this would be appreciated! Thanks very much in advance! This email message may contain legally privileged and/or confidential information. If you are not the intended recipient(s), or the employee or agent responsible for the delivery of this message to the intended recipient(s), you are hereby notified that any disclosure, copying, distribution, or use of this email message is prohibited. If you have received this message in error, please notify the sender immediately by e-mail and delete this email message from your computer. Thank you. [[alternative HTML version deleted]] You used Time in your formula but Day in your data so I changed the formula to Day assuming that's what you meant Perhaps, not the most elegant solution, but how about something like this library(lattice) Subject-c(1,1,1,1,1,1,2,2,2,2,2,2,3,3,3,3,3,3,4,4,4,4,4,4) Day-c(1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6) Activity-c(2,3,4,3,7,4,5,8,2,8,4,6,2,5,3,8,9,5,6,3,4,5,6,7) EventA-c(Yes,NA,Yes,NA,NA,NA,Yes,NA,NA,NA,NA,NA, Yes,NA,NA,NA,NA,NA,NA,NA,Yes,NA,NA,NA) Data-data.frame(cbind(Subject,Day,Activity,EventA)) xyplot(Activity~ Day | Subject) xyplot(Activity ~ Day | Subject, data = Data, subscripts = TRUE, panel = function(x, y, subscripts, ...){ panel.xyplot(x, y) wh - Data[subscripts, ] panel.abline(v = wh$Day[!is.na(wh$EventA)]) }) -- Kenneth Knoblauch Inserm U846 Stem-cell and Brain Research Institute Department of Integrative Neurosciences 18 avenue du Doyen L?pine 69500 Bron France tel: +33 (0)4 72 91 34 77 fax: +33 (0)4 72 91 34 61 portable: +33 (0)6 84 10 64 10 http://www.sbri.fr/members/kenneth-knoblauch.html __ R-help@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.
Re: [R] plotting wind rose data (Karl Ropkins)
David, Following on from Jim Lemon's suggest that polar plots might be more what you need, there are both wind rose and polar plot functions in the openair package that might be of use (particularly windRose and polarFreq). They will not do everthing you are after without some careful conditioning or extra work with latticeExtra, but could get you some of the way there. Karl Ropkins, Hi List, I am trying to create a spatial representation of some wind data. I have the season, frequency, strength and direction of the wind from 10 different locations, the coverage of the area that I am interested in is not 100% there are small gaps in my coverage due to the location of the weather stations. I am trying to create a series of wind maps e.g. the Prevailing Winds, the maximum seasonal wind, etc. Could any body recommend any R-packages that would cover this type problem/issue? Hi David, While there are several packages that include plotting routines for wind roses, it looks to me as though you want to define a small number of vectors representing prevailing wind, etc., possible overlaying these on a plot. That might be a job for a circular plotting routine (e.g. polar.plot) rather than a wind rose. Jim __ R-help@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.