search for: carslaw

...re reports/analyses are produced automatically through submission of data to a web-sever. I am mostly interested in situations where pdf reports have been produced rather than, say, a plot/table etc shown on a web page. I've had limited success finding examples on this. Many thanks. David Carslaw Environmental Research Group MRC-HPA Centre for Environment and Health King's College London Franklin Wilkins Building Stamford Street London SE1 9NH david.carslaw at kcl.ac.uk -- View this message in context: http://r.789695.n4.nabble.com/Examples-of-web-based-Sweave-use-tp3425324p3425324...

...United Kingdom.1252 [2] LC_CTYPE=English_United Kingdom.1252 [3] LC_MONETARY=English_United Kingdom.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United Kingdom.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base Dr David Carslaw King's College London Environmental Research Group Franklin Wilkins Building 150 Stamford Street London SE1 9NH -- View this message in context: http://r.789695.n4.nabble.com/difference-in-sort-order-linux-Windows-R-2-11-0-tp2234251p2234251.html Sent from the R help mailing list archive at Na...

...y difficulty is putting them into a data frame along with the other columns "fac" and "other". Note that quantile will return different length vectors due to different numbers of NAs for a and b. There's something I'm just not seeing - can you help? Many thanks. David Carslaw ----- Institute for Transport Studies University of Leeds -- View this message in context: http://www.nabble.com/subset-grouped-data-with-quantile-and-NA%27s-tp19102795p19102795.html Sent from the R help mailing list archive at Nabble.com.

...s desired. My question is, given a vector of column name(s) to split by, how can I supply split with these? I was thinking something like: splits <- c("fac1", "fac2") split(dat, list(splits)) But that is clearly wrong....and I can't see the solution Many thanks David Carslaw Science Policy Group Environmental Research Group MRC-HPA Centre for Environment and Health King's College London Room 4.129 Franklin Wilkins Building Stamford Street London SE1 9NH -- View this message in context: http://r.789695.n4.nabble.com/question-about-split-tp3605281p3605281.ht...

...th functions against the dependent variable (just as in Carlslaw et al. 2007, fig. 1*). Otherwise, I'd like to know how to plot the smoothed predictor variable values against the partial residuals. Thank you for taking the time to read my message. Any?help is gratefully appreciated, Oscar * Carslaw, D., Beevers, S.,?Tate, J. 2007. Modelling and assessing trends in traffic-related emissions using a generalised additive modelling approach. Atmospheric Environment 41 (2007) 5289?5299

Hi all, Given a data frame: my.df <- data.frame(a = c(1:5, 1:10, 1:20), b = runif(35)) I want to split it by "a" such that I end up with a list containing 3 components i.e. the first containing a = 1 to 5, the second a = 1 to 10 etc. In other words, sets of sequences of a. I can't seem to find the right form using the split function - can you help? Much appreciated. David

...process it. I am finding however that lines of code that try and import data in this way do not seem to be evaluated and subsequent functions using the data fail. I can't find any references to similar problems in R-help, hence my enquiry to you. Any tips would be welcome! David Dr David Carslaw Science Policy Group Environmental Research Group MRC-HPA Centre for Environment and Health King's College London Room 4.129 Franklin Wilkins Building Stamford Street London SE1 9NH david.carslaw at kcl.ac.uk -- View this message in context: http://r.789695.n4.nabble.com/remote-data-acce...

I seem to have a Friday afternoon block and can't see the easiest way of doing this. Given a data frame like: dat <- data.frame(x = runif(100), y = runif(100), group = rep(letters[1:10], each = 10)) > head(dat) x y group 1 0.876751503 0.6518345 a 2 0.627067150 0.8801790 a 3 0.632465192 0.1768305 a 4 0.060359554 0.8835652 a 5 0.675868776 0.7721177

.../m3, mg/m3, ng/m3...). My question is this - is there a way to separately process the pollutant and units and supply the correct overall expression to plot? I'm not sure of the best strategy for this having read through previous posts etc., and would appreciate your help! Many thanks. David Carslaw ----- Institute for Transport Studies University of Leeds -- View this message in context: http://www.nabble.com/text-processing-for-plots-tp19331219p19331219.html Sent from the R help mailing list archive at Nabble.com.

...t(tick.number = 10))) But this has no effect. While I can set ticks at specific points (and labels) this is difficult because a wide range of time intervals are used. Any suggestions? Thanks. David ----- Institute for Transport Studies University of Leeds Leeds LS2 9JT Tel: 0113 343 7522 d.c.carslaw at its.leeds.ac.uk -- View this message in context: http://www.nabble.com/xyplot-axis-scaling-with-dates-times-tp24620903p24620903.html Sent from the R help mailing list archive at Nabble.com.

...With thanks, John =================================================================== John Twyman | Computer Systems Officer | Faculty of Science | College of Sciences & Technology | Room 207/210 | Email: J.Twyman@cst.usyd.edu.au Carslaw Building, F07 | Phone: +61 2 9351 5453 University of Sydney NSW 2006 | Fax: +61 2 9351 5467 ===================================================================

It seems to be working for us now. I have the following relevant settings in the [global] section of smb.conf ... remote announce = 199.170.245.255/ENGINEER 199.170.224.255/ENGINEER 199.170.232.255/ENGINEER workgroup = ENGINEER wins support = no wins server = 199.170.224.6 domain master = yes local master = yes preferred master = yes ; os level = 0 don't want Samba to

Hi R-users, Since R.2.11 aggregate can now deal with non-scalar functions, which is very useful to me. However, I have a question about how best to process the output. test <- data.frame(a = rep(c("g1", "g2"), each = 50), b = runif(100)) res <- aggregate(test$b, list(group = test$a), function(x) quantile(x, probs = c(0.05, 0.95))) > res group x.5%

I frequently work with hourly data in GMT format, and ensure that any data I read into R are expressed in GMT through as.POSIXct(mydate, tz = "GMT"). I am interested in processing air pollution data in GMT where where peaks in emissions (say from road traffic sources) tend to occur at the same LOCAL time each day. I am interested in calculating mean concentrations by hour of day i.e.