Re: [R] GIS in R

2022-09-01 Thread Thomas Adams
Hi Nick, First off, I would suggest going to: https://cran.r-project.org/web/views/Spatial.html also have a close look at https://journal.r-project.org/archive/2017/RJ-2017-067/RJ-2017-067.pdf And you might also consider the use of GRASS GIS, which interfaces with R pretty seamlessly:

Re: [R] forecast accuracy

2021-02-05 Thread Thomas Adams
Look at the 'verification' contributed package Tom On Fri, Feb 5, 2021 at 12:36 PM Faheem Jan via R-help wrote: > I am working in the functional time series, I obtain the one year ahead > forecast in the functional format, Know i want to forecast accuracy for > example mean absolute percentage

Re: [R] National Weather Service Data

2020-07-03 Thread Thomas Adams
Hi Philip! I'm a little familiar with rNOMADS... I tried following the example for 'ArchiveGribGrab' using a more recent date #An example for the Global Forecast System #Get data for January 1 2014 #Temperature at 2 m above ground #3 hour prediction # using GRIB abbrev <- "gfsanl" model.date <-

Re: [R] CRAN Task View: Optimization and Mathematical Programming

2019-02-06 Thread Thomas Adams
Did you Google "R stat water resources optimization" or "r stats optimization"? It seems that applying some of the references below with R packages in the optimization task view probably gets you where you want to go. https://cran.r-project.org/web/packages/reservoir/reservoir.pdf

[R] Problem with my function using as.POSIXct

2018-03-28 Thread Thomas Adams
Hello all: I wrote a function: my.bastimeToSynoptic <- function(x) { f<-unlist(strsplit(as.character(x), " ")) hr<-unlist(strsplit(f[2], ":")) if(as.numeric(hr[1])<6) { synoptic<-"00" } else { synoptic<-as.integer(as.numeric(hr[1])/6)*6 }

Re: [R] Help needed with aggregate or other solution

2017-10-26 Thread Thomas Adams
:13:00", "2012-04-08 12:59:00", > "2012-04-09 13:53:00", "2012-04-10 14:15:00", "2012-04-11 14:26:00", > "2012-04-12 14:05:00", "2012-04-13 13:37:00", "2012-04-14 13:52:00", > "2012-04-15 15:00:00", "2012-04-16 14:42:00", "2012-04-17 14:21:00",

[R] Help needed with aggregate or other solution

2017-10-25 Thread Thomas Adams
Hello all! I've been struggling with is for many hours today; I'm close to getting what I want, but not close enough... I have a dataframe consisting of two date-time columns followed by two numeric columns. what I need is the max value (in the first numeric column) based on the 2nd date-time

Re: [R] about installing smwrGraphs package

2017-07-23 Thread Thomas Adams
Lily, This: https://cran.r-project.org/web/packages/githubinstall/vignettes/githubinstall.html will probably be useful to you. Best, Tom On Sun, Jul 23, 2017 at 12:31 PM, Uwe Ligges < lig...@statistik.tu-dortmund.de> wrote: > > > On 23.07.2017 05:28, lily li wrote: > >> Hi R users, >> >> I'm

Re: [R] Histogram plots in Lattice with spatialgrid dataframe data

2017-07-09 Thread Thomas Adams
gt; [2] http://adv-r.had.co.nz/Reproducibility.html > > [3] https://cran.r-project.org/web/packages/reprex/index.html > > -- > Sent from my phone. Please excuse my brevity. > > On July 9, 2017 12:32:32 PM PDT, Thomas Adams <tea...@gmail.com> wrote: > >Hello all, >

Re: [R] Histogram plots in Lattice with spatialgrid dataframe data

2017-07-09 Thread Thomas Adams
' yr_2012$year<-'2012' yr_2013$year<-'2013' yr_2014$year<-'2014' yr_2015$year<-'2015' yr_2016$year<-'2016' bias<-rbind(yr_1997,yr_1998,yr_1999,yr_2000,yr_2001,yr_2002,yr_2003,yr_2004,yr_2005,yr_2006,yr_2007,yr_2008,yr_2009,yr_2010,yr_2011,yr_2012,yr_2013,yr_2014,yr_2015,yr_2016

[R] Histogram plots in Lattice with spatialgrid dataframe data

2017-07-09 Thread Thomas Adams
Hi all, I can not seem to get what I want using the Lattice package to generate an array of histograms of spatialgrid dataframe data. I can use the sp package and spplot to generate an array of maps that display an array of spatialgrid dataframe data -- that's good. I have:

Re: [R] Help with reshape/reshape2 needed

2017-07-05 Thread Thomas Adams
9231NA NANA > > - > David L Carlson > Department of Anthropology > Texas A University > College Station, TX 77840-4352 > > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of

[R] Help with reshape/reshape2 needed

2017-07-05 Thread Thomas Adams
Hi all: I'm struggling with getting my data re-formatted using functions in reshape/reshape2 to get from: 1957 0.86250 1958 0.75000 1959 0.3 1960 0.28750 1963 0.67500 1964 0.93750 1965 0.02500 1966 0.38750 1969 0.08750 1970 0.27500 1973 0.5

Re: [R] can we visualize water flows with 3d in R?

2016-10-13 Thread Thomas Adams
..@gmail.com> wrote: > On 13/10/2016 8:35 AM, Thomas Adams wrote: > >> All, >> >> Very respectfully, there are no R packages that can do what Marna desires. >> > > I would guess that's not literally true, in that there are several > graphics pack

Re: [R] can we visualize water flows with 3d in R?

2016-10-13 Thread Thomas Adams
All, Very respectfully, there are no R packages that can do what Marna desires. His/Her data, undoubtably, comes from a 1-D hydraulic model simulation -- where output is generated at channel cross-sections -- representing the sloping water surface elevation of the centerline of flow in a stream

Re: [R] Struggling trying to plot points on boxplot

2016-06-02 Thread Thomas Adams
> > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Thu, Jun 2, 2016 at 8:36 AM, Thomas Adams <tea...@gmail.com> wrote: > >> Hello all: >> >> I've been beating my head on this for over a day and I have done a lot of >> Google'ing with no lu

[R] Struggling trying to plot points on boxplot

2016-06-02 Thread Thomas Adams
Hello all: I've been beating my head on this for over a day and I have done a lot of Google'ing with no luck. I have generated a 'time-series' of boxplots (227), covering more than a 30-day period at 6-hour intervals, which summarize an ensemble forecast. What I want to do is over-plot the

Re: [R] Help needed to format data for boxplot time-series

2016-06-01 Thread Thomas Adams
gt; so that we can use your data. > > From your description maybe > > boxplot(split(yourdata$value, yourdata$valid_time)) > > can give you what you want. > > Regards > Petr > > > -Original Message- > > From: R-help [mailto:r-help-boun...@r-project.o

[R] Help needed to format data for boxplot time-series

2016-06-01 Thread Thomas Adams
All: I have used R in combination with GRASS GIS spatial data (using spgrass) many times in the past to generate a 'time series' of boxplots, to show variations over time. But I have a new problem, not involving spatial data, but rather, true time-series data (snippet shown below). So, what I

Re: [R] question about probplot in e1071 package

2016-04-04 Thread Thomas Adams
.networks.imdea.org/people/~luis_nunez/ > > > El 4 abr 2016, a las 18:19, Thomas Adams <tea...@gmail.com> escribió: > > Hello! > > I am using probplot in the e1071 package and want to do something like the > following, only with the the 2nd plot overlaying the fir

[R] question about probplot in e1071 package

2016-04-04 Thread Thomas Adams
Hello! I am using probplot in the e1071 package and want to do something like the following, only with the the 2nd plot overlaying the first. I can't seem to make it work. Any suggestions? *library(e1071) **x <- rnorm(100, mean=5)* *y <- rnorm(100, mean=3)* *probplot(x, line=FALSE) *

[R] NULL dev.lis()

2015-11-08 Thread Thomas Adams
All, I have previous built R from source many times, generally, without problems. However on my new Ubuntu 15.04 Linux system with R 3.2.2 when I run the command dev.list() I get: > dev.list() NULL At the completion of running ./configure, I have R is now configured for x86_64-pc-linux-gnu

Re: [R] 'strange' R graphics problem | Linux...

2015-10-14 Thread Thomas Adams
Evan, I have Ubuntu 14.04 and 15.10 at home and have not had problems, but I don't think I've been using R 3.2.2 — I'll try this evening. Tom On Wed, Oct 14, 2015 at 2:47 PM, Evan Cooch <evan.co...@gmail.com> wrote: > Tom -- > > On 10/14/2015 3:35 PM, Thomas Adams

Re: [R] 'strange' R graphics problem | Linux...

2015-10-14 Thread Thomas Adams
Evan, Not that this helps you, but I am using a very similar platform and I am having the identical problem. My test simply comes from the first help(plot) example. I tried doing some things to 'correct' the problem and ended up mucking-up my Gnome environment. In the process, I was able to get

Re: [R] 'strange' R graphics problem | Linux...

2015-10-14 Thread Thomas Adams
iconv >TRUETRUETRUETRUETRUETRUE > NLS profmem cairo ICU long.double libcurl >TRUE FALSETRUETRUETRUE FALSE > > On 10/14/2015 4:00 PM, Evan Cooch wrote: > > &g

Re: [R] Best Mac for R

2015-02-26 Thread Thomas Adams
Dan, FWIW, I have basically the system you describe, except a larger HD — I'm quite happy, but I'm a biased Mac user, although I love my Ubuntu Linux machine as well… One can bring any machine to its knees, so there is the element of expectations. A MacBook Pro stacks up as well or better

Re: [R] plot circles where z values are circles radius

2015-01-23 Thread Thomas Adams
I think what you should look at are these web sites I found with a Google search: http://flowingdata.com/2010/11/23/how-to-make-bubble-charts/ http://www.r-bloggers.com/bubble-plots-ggplot2/ http://cran.r-project.org/web/packages/HSAUR/vignettes/Ch_logistic_regression_glm.pdf

Re: [R] how update to latest version of R on mac

2014-12-21 Thread Thomas Adams
Postal Address: Visiting address: Box 408 Kanikegränd 3A 541 28 Skövde Skövde From: Thomas Adams tea...@gmail.com Date: Sunday, December 21, 2014 2:06 AM To: Jane Synnergren jane.synnerg...@his.se Cc: R-help r-help@r-project.org Subject: Re: [R] how update to latest version of R

Re: [R] how update to latest version of R on mac

2014-12-20 Thread Thomas Adams
Jane, It's possible, since you are using Mac OS X Lion (10.7.5), you need the R version for that, which can be found here: http://cran.r-project.org/bin/macosx/old/R-3.1.1-mavericks.pkg Cheers! Tom On Fri, Dec 19, 2014 at 2:56 PM, Jane Synnergren jane.synnerg...@his.se wrote: Hi, I get

Re: [R] We won an award, and didn't even notice...

2014-12-19 Thread Thomas Adams
Duncan, Congratulations to you and all the founders and contributors — very much deserved; thank you!! Tom On Fri, Dec 19, 2014 at 9:24 AM, Duncan Murdoch murdoch.dun...@gmail.com wrote: This morning I was reading Jeff Leek's list of awesome things other people did in 2014 at

Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread Thomas Adams
Matteo, I tried your example code using R 3.1.1 on an iMac (24-inch, Early 2009), 3.06 GHz Intel Core 2 Duo, 8 GB 1333 MHz DDR3, NVIDIA GeForce GT 130 512 MB running Mac OS X 10.10 (Yosemite). After entering your code, the elapsed time from the time I hit return to when the graphics appeared was

Re: [R] speed issue in simulating a stochastic process

2014-11-06 Thread Thomas Adams
Matteo, Ah — OK, N=20, I did not catch that. You have nested for loops, which R is known to be exceedingly slow at handling — if you can reorganize the code to eliminate the loops, your performance will increase significantly. Tom On Thu, Nov 6, 2014 at 7:47 AM, Matteo Richiardi

Re: [R] ubuntu 14.04

2014-09-24 Thread Thomas Adams
As do I... On Wednesday, September 24, 2014, Jeff Newmiller jdnew...@dcn.davis.ca.us wrote: Any incompatibility is a high standard, but I run it just fine on that platform. --- Jeff NewmillerThe

Re: [R] Best way to merge 300+ .5MB dataframes?

2014-08-11 Thread Thomas Adams
Grant, Assuming all your filenames are something like file1.txt, file2.txt,file3.txt... And using the Mac OSX terminal app (after you cd to the directory where your files are located... This will strip off the 1st lines, that is, your header lines: for file in *.txt;do sed -i '1d'${file}; done

Re: [R] installing package gstat

2013-10-10 Thread Thomas Adams
Simona, You need to install the dependencies: install.packages(gstat,dependencies=T) Tom On Thu, Oct 10, 2013 at 11:58 AM, Simona Augyte simona.aug...@uconn.eduwrote: Hello, # I am able to install.packages(gstat) #but when I try to upload I get an error message library(gstat) #Error in

Re: [R] reshape2 issue continued

2013-06-05 Thread Thomas Adams
Bruce, I'm not sure what's going on since I tried this on my Linux system running R 3.0.0 and just did: library(reshape2) help(cast) and the help for 'cast' came up. There was no indication to me that 'reshape' was needed and I can not see a dependency in CRAN for 'reshape'. But I built R from

Re: [R] the joy of spreadsheets (off-topic)

2013-04-27 Thread Thomas Adams
wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-** project.org r-help-boun...@r-project.org] On Behalf Of S Ellison Sent: Friday, April 26, 2013 9:08 AM To: Thomas Adams; peter dalgaard Cc: r-help Subject: Re: [R] the joy

Re: [R] How to export graph value in R

2013-04-26 Thread Thomas Adams
Anup, You should have provided some additional information, such as that the function 'hypsometric' is found in the hydroTSM contributed package. Nevertheless, here's what I did (maybe not elegant, but it works) : (1) at the R command prompt simply type hypsometric -- the source code for the

Re: [R] the joy of spreadsheets (off-topic)

2013-04-24 Thread Thomas Adams
One might wonder if the Excel error was indeed THAT or perhaps a way to get the desired results, give the other issues in their analysis? On Wed, Apr 24, 2013 at 11:58 AM, peter dalgaard pda...@gmail.com wrote: In case you haven't noticed, this is making the rounds in the media, including a

Re: [R] Odd graphic device behavior

2013-03-29 Thread Thomas Adams
Thank you Peter, that sounds pretty reasonable. Best, Tom On Friday, March 29, 2013, peter dalgaard wrote: On Mar 28, 2013, at 22:27 , Thomas Adams wrote: All, Well to my relief and embarrassment, I discovered my problem. About 5 weeks ago, I shutdown my computer and moved it. When

Re: [R] Odd graphic device behavior

2013-03-28 Thread Thomas Adams
Luca, Thank you for the suggestion; I do have an Nvidia graphics card and I am using the Nvida driver; still searching for a solution, quite odd... Tom On Thu, Mar 28, 2013 at 12:28 PM, Luca Nanetti luca.nane...@gmail.comwrote: Thomas, any chance that you could have problems with your

Re: [R] Odd graphic device behavior

2013-03-28 Thread Thomas Adams
port and presto changeo -- problem solved!! The odd thing was, besides my R graphics device window problems, everything else was fine... Out of curiosity, can anyone explain this? Thanks for all the help... Tom On Thu, Mar 28, 2013 at 4:06 PM, Thomas Adams tea...@gmail.com wrote: Luca, Thank

Re: [R] Odd graphic device behavior

2013-03-27 Thread Thomas Adams
John, Thanks for the suggestion, but no. I have even gone so far as to rebuild R from source, re-booted my computer, and tried the 'experiment': require(stats) plot(cars) immediately after starting R. Still the same result. I think it must be related to some default Ubuntu Unity

Re: [R] Odd graphic device behavior

2013-03-27 Thread Thomas Adams
, but the lettering was large as before and the lines were thick as they were previously. Regards, Tom On Wed, Mar 27, 2013 at 1:43 PM, peter dalgaard pda...@gmail.com wrote: On Mar 27, 2013, at 18:11 , David Winsemius wrote: On Mar 27, 2013, at 9:07 AM, Thomas Adams wrote: John, Thanks

[R] Odd graphic device behavior

2013-03-26 Thread Thomas Adams
I'm experiencing odd graphics device behavior running R 2.15.3 on Ubuntu. Regardless of what I try like: require(stats) plot(cars) lines(lowess(cars)) plot(sin, -pi, 2*pi) for example, the graphics device fills the entire screen with the graphic and a very large font. When I shrink the graphics

Re: [R] Running other programs from R

2013-03-17 Thread Thomas Adams
Use the system() command. e.g. system(Multilog) On Sat, Mar 16, 2013 at 5:09 PM, Sedat Sen sedatse...@gmail.com wrote: Dear list, I want to run a statistical program (using its .exe file) from R by writing a script. I know there are some packages that call WinBUGS, Mplus etc. form R. I

[R] Using grImport to create a watermark

2012-12-28 Thread Thomas Adams - NOAA Federal
Hi… I want to use grImport to create a watermark on a plot() using the methods Paul Murrell describes here: http://cran.r-project.org/web/packages/grImport/vignettes/import.pdf (page 28). I can essentially reproduce this manually at the R prompt, and independently I can use grid.picture(…)

Re: [R] Using grImport to create a watermark

2012-12-28 Thread Thomas Adams - NOAA Federal
Brian Ripley rip...@stats.ox.ac.ukwrote: I suspect you need a device supporting translucency: PostScript does not and hence postscript() cannot. Try the pdf() device (and convert the output if you need it). On 28/12/2012 17:23, Thomas Adams - NOAA Federal wrote: Hi… I want to use grImport

Re: [R] Auomatic Download of climate netcdf4 files from the Climate Explorer Website

2012-09-14 Thread Thomas Adams
I'm not sure about R code, but wget (http://www.gnu.org/software/wget/) should work. Here are some examples: http://www.editcorp.com/Personal/Lars_Appel/wget/v1/wget_7.html Cheers! Tom On Thu, Sep 13, 2012 at 1:37 PM, ucakmde maris_...@hotmail.com wrote: Dear R-helpers, I am trying to

Re: [R] On Reproducible Code

2012-07-30 Thread Thomas Adams
Hadley, Thank you for posting this. I think the danger is that novices (and I'm not far removed from that category) can be intimidated by R, but more so by R experts that make people AFRAID to ask questions. The danger is that these intimidating R experts could turn people away from using R; at

Re: [R] On Reproducible Code

2012-07-30 Thread Thomas Adams
a reproducible example, etc., but they don't deserve to be treated disrespectfully — I think we are agreeing on this… Tom On Mon, Jul 30, 2012 at 8:02 AM, peter dalgaard pda...@gmail.com wrote: On Jul 30, 2012, at 13:05 , Thomas Adams wrote: Nothing is gained by punishing people over their internet

Re: [R] is it possible to insert a figure into into another new figure by r script

2012-07-10 Thread Thomas Adams
Jie, I think the R contributed package, grImport, by Paul Murrell does what you want. See this: http://www.jstatsoft.org/v30/i04/paper/ Tom On Mon, Jul 9, 2012 at 11:14 PM, Mikhail Titov m...@gmx.us wrote: Jie Tang totang...@gmail.com writes: hi R-users Now I have a figure in emf or

Re: [R] Plotting tidal speed and direction in R

2012-06-20 Thread Thomas Adams
Zoe, If you use a Wind Rose diagram (e.g., http://rgm2.lab.nig.ac.jp/RGM2/func.php?rd_id=circular:wind rose), which are available in several packages, you could plot tidal speed and direction Wind Rose diagrams in a x-y plot to plot location. If you used the lattice package, this could give you

Re: [R] On the Design of the R Language

2012-04-25 Thread Thomas Adams
Hmmm… an 'objective' assessment? Maybe. But it looks to me that some commercial enterprise paid for this study as a means to argue against the use of R in favor of a commercial package… On Wed, Apr 25, 2012 at 1:48 PM, R. Michael Weylandt michael.weyla...@gmail.com michael.weyla...@gmail.com

Re: [R] Organizations where IT has approved the use of R software

2012-04-12 Thread Thomas Adams
All: R is used by the NOAA/U.S. National Weather Service to generate graphics representing real-time hydrologic ensemble (probabilistic) forecasts. Go to: http://www.erh.noaa.gov/mmefs/ to see. It is also used in research and development for forecast verification and analyses for the calibration

Re: [R] time-series features skewness kurtosis periodicity trend seasonality

2012-04-03 Thread Thomas Adams
John, If the paper is the one by Xiaozhe Wang, Kate A. Smith, Rob Hyndman, and Damminda Alahakoon, the way I read it, it does not say they created an R package (that may be on CRAN). They may have just used R in their analyses; just a guess... Tom On Tue, Apr 3, 2012 at 12:34 PM, John Kohr

[R] RSqlite UPDATE command problem

2012-03-27 Thread Thomas Adams
All: I am using RSqlite and want to be able to update individual values in a record, such as with this simple example: library(RSQLite) drv-dbDriver(SQLite) con-dbConnect(drv,test.db) my.data-data.frame(countries=c(US,UK,Canada,Australia,NewZealand),vals=c(52,36,74,10,98))

Re: [R] RSqlite UPDATE command problem

2012-03-27 Thread Thomas Adams
2012 14:18, Thomas Adams thomas.ad...@noaa.gov wrote: All: I am using RSqlite and want to be able to update individual values in a record, such as with this simple example: library(RSQLite) drv-dbDriver(SQLite) con-dbConnect(drv,test.db) my.data-data.frame(countries=c(US,UK,Canada,Australia

Re: [R] Struggling with zoo and aggregate

2012-03-26 Thread Thomas Adams
Gabor, That does it! I can't thank you enough… Many thanks, Tom On Mon, Mar 26, 2012 at 7:22 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Sun, Mar 25, 2012 at 10:20 PM, Thomas Adams thomas.ad...@noaa.gov wrote: Gabor, Thank you for your help -- it did help me a lot

Re: [R] Struggling with zoo and aggregate

2012-03-25 Thread Thomas Adams
ggrothendi...@gmail.com wrote: On Sat, Mar 24, 2012 at 10:44 PM, Thomas Adams thomas.ad...@noaa.gov wrote: All: I have a SQlite database where I have stored some verification data by date time (cycle Z/UTC), lead_time as well as type, duration, etc. I would like to analyze plot

[R] Struggling with zoo and aggregate

2012-03-24 Thread Thomas Adams
All: I have a SQlite database where I have stored some verification data by date time (cycle Z/UTC), lead_time as well as type, duration, etc. I would like to analyze plot the data as monthly averages. I have looked at a bunch of examples which use some combination of zoo and aggregate, but I

Re: [R] a somewhat related cartoon

2012-03-07 Thread Thomas Adams
Sarah, Thanks! Great stuff... Tom On Wed, Mar 7, 2012 at 10:07 AM, Sarah Goslee sarah.gos...@gmail.comwrote: Quoting from today's PhD Comics, available at: http://www.phdcomics.com/comics.php?f=1476 What the methodology section says: Analysis was performed using a commercially available

Re: [R] The Future of R | API to Public Databases

2012-01-13 Thread Thomas Adams
Sarah, I agree; I think it would be the exception rather than the rule that one would access these public data sources given the range of needs of R users, who are generally analyzing their own data. Plus, IMO, it just is not very difficult to reformat the data to a suitable format, if need be,

Re: [R] Boxplot of multiple vectors with different lengths

2011-12-12 Thread Thomas Adams
Ryan, I think you could do what you want by having the vector data written to separate files; then create a file containing the individual file names. In R, read the file containing the list of file names and loop through this reading in the individual vector files. Maybe this is an inelegant,

Re: [R] Boxplot of multiple vectors with different lengths

2011-12-12 Thread Thomas Adams
, Dec 12, 2011 at 1:43 PM, Thomas Adams thomas.ad...@noaa.gov wrote: Ryan, I think you could do what you want by having the vector data written to separate files; then create a file containing the individual file names. In R, read the file containing the list of file names and loop through

Re: [R] Question about ggplot2 and stat_smooth

2011-10-06 Thread Thomas Adams
Dennis Hadley, This does exactly what I need — thank you so much! Regards, Tom On 10/4/11 5:34 PM, Dennis Murphy wrote: Hi Hadley: When I tried your function on the example data, I got the following: dd- data.frame(year = rep(2000:2008, each = 500), y = rnorm(4500)) g- function(df, qs =

Re: [R] Question about ggplot2 and stat_smooth

2011-10-04 Thread Thomas . Adams
Adams thomas.ad...@noaa.gov Cc: R-help forum r-help@r-project.org On Mon, Oct 3, 2011 at 12:24 PM, Thomas Adams thomas.ad...@noaa.gov wrote:  I'm interested in creating a graphic -like- this: c - ggplot(mtcars, aes(qsec, wt)) c + geom_point() + stat_smooth(fill=blue, colour=darkblue

Re: [R] Question about ggplot2 and stat_smooth

2011-10-04 Thread Thomas Adams
Hadley: Below is an example of what I am trying to do, I just don't understand how to supply the limits to the blue and pink shaded regions and the values of the black line, which are meant to represent from bottom to top, the 5%, 25%, 50%, 75%, 95% limits that I get from quantile(): h +

[R] Question about ggplot2 and stat_smooth

2011-10-03 Thread Thomas Adams
I'm interested in creating a graphic -like- this: c - ggplot(mtcars, aes(qsec, wt)) c + geom_point() + stat_smooth(fill=blue, colour=darkblue, size=2, alpha = 0.2) but I need to show 2 sets of bands (with different shading) using 5%, 25%, 75%, 95% limits that I specify and where the heavy

Re: [R] Question about ggplot2 and stat_smooth

2011-10-03 Thread Thomas Adams
Andrés, Thank you for your help, but that does not capture what I'm looking for. I need to be able to control the shaded bound limits and they need to be coincident. Tom On 10/3/11 3:37 PM, Andrés Aragón wrote: Hi, Try some like this: c- ggplot(mtcars, aes(qsec, mpg, colour=factor(cyl)))

Re: [R] Looking for Filliben (correlation test)

2011-06-30 Thread Thomas . Adams
- Original Message - From: gaiarrido gaiarr...@usal.es Date: Thursday, June 30, 2011 1:11 pm Subject: Re: [R] Looking for Filliben (correlation test) To: r-help@r-project.org Mario, I did a google search and found this: http://genepi.qimr.edu.au/staff/davidD/R/filliben.R Cheers! Tom

Re: [R] Thiessen Method

2011-05-26 Thread Thomas Adams
Federico, I understand what you are after — you want time-series estimates based on the Thiessen polygon estimates taken from the station time-series data. My recommendation is that the process of doing this would be far easier using something like GRASS GIS, possibly in conjunction with R

Re: [R] recommendation on r scripting tutorial?

2011-04-02 Thread Thomas . Adams
Wensui, use google and search for r stats scripts and find among others: http://cran.r-project.org/doc/contrib/Lemon-kickstart/kr_scrpt.html Tom - Original Message - From: Wensui Liu liuwen...@gmail.com Date: Saturday, April 2, 2011 11:22 am Subject: [R] recommendation on r scripting

[R] 2 questions about probplot in package e1071

2011-03-26 Thread Thomas Adams
The contributed package e1071 does exactly what I want except that I need to have (1) the abscissa and ordinate axes swapped, with the probability scale on the bottom and the quantiles scale on the LHS. Using the following example: library(e1071) x - rnorm(100, mean=5) probplot(x,

Re: [R] 2 questions about probplot in package e1071

2011-03-26 Thread Thomas Adams
David, Thanks! This is very helpful! I'm still very much a novice… Tom On 3/26/11 4:11 PM, David Winsemius wrote: On Mar 26, 2011, at 4:06 PM, Thomas Adams wrote: The contributed package e1071 does exactly what I want except that I need to have (1) the abscissa and ordinate axes swapped

[R] Help needed with plot axis labeling

2011-03-24 Thread Thomas Adams
I have looked at many examples and tried many different combinations of doing this, but with no luck. I have something like this: plot(1:10, xaxt = n) axis(1, xaxp=c(2, 9, 7)) axis(4) but, what I need is to have different labels for axis-4 than those for axis-2 (the vertical axes) — that is,

Re: [R] Timeseries Data Plotted as Monthly Boxplots

2011-02-16 Thread Thomas Adams
Katrina, What I have done, if I understand what you are after, was to create a list for each month of data - in order. Then, create a boxplot - in order - by month/year. I do this for our ensemble streamflow forecasts. The key us to create the list of values by month. Regards, Tom Sent from

Re: [R] Layout of mulitpage conditioned lattice plots

2010-12-19 Thread Thomas Adams
Dennis, Thank you; this helps me, too! Tom On 12/19/10 11:45 AM, Dennis Murphy wrote: Hi Dieter: If I read your intention correctly, you need a third element in layout = . Here's a little example: df- data.frame(month = rep(month.abb, each = 20), time = rep(1:20, 12),

Re: [R] Why software fails in scientific research

2010-07-01 Thread Thomas Adams
OK… My Grandfather, who was a farmer, was outstanding in his field… Cheers… Murray M Cooper, PhD wrote: For what its worth! A good friend who also happens to be an ecologist told me An ecologist is a statistician who likes to be outside. Murray M Cooper, Phd Richland Statistics -

Re: [R] two questions for R beginners

2010-02-26 Thread Thomas Adams
Paul, I think your point you need [to] spend at least a few hours a week on it is key. Since I am not doing statistics daily, more in fits starts as my latest project -may- require, my approach has been more task oriented. A less-than-ideal approach. So, I think your suggestion is

[R] help needed using t.test with factors

2010-02-04 Thread Thomas Adams
I am trying to use t.test on the following data: datetypeINTERVALnCASESMTFSDFMTOSDO nFSTMFnOBSMOMBBIASCVBIASEVMEMAE RMSECRCF 2001-06-15avnGE1.0043850.2460.3001.502 0.55613671.3734385

Re: [R] help needed using t.test with factors

2010-02-04 Thread Thomas Adams
Dennis, Thank you for the suggestion, but I get this error: t.test(MAE ~ type,data=data) Error in t.test.formula(MAE ~ type, data = data) : grouping factor must have exactly 2 levels Tom Dennis Murphy wrote: Hi: On Thu, Feb 4, 2010 at 11:07 AM, Thomas Adams thomas.ad...@noaa.gov

Re: [R] help needed using t.test with factors

2010-02-04 Thread Thomas Adams
Thomas Adams wrote: I am trying to use t.test on the following data: datetypeINTERVALnCASESMTFSDFMTOSDO nFSTMFnOBSMOMBBIASCVBIASEVMEMAE RMSECRCF 2001-06-15avnGE1.0043850.2460.3001.502 0.5561367

Re: [R] help needed using t.test with factors

2010-02-04 Thread Thomas Adams
Peter, Thank you very much! That did the trick… Regards, Tom Peter Ehlers wrote: Tom, t.test(MAE ~ type, data=data, subset=type %in% c('hpc','rfc')) -Peter Ehlers Thomas Adams wrote: Dennis, Thank you for the suggestion, but I get this error: t.test(MAE ~ type,data=data) Error

[R] ggplot2/qplot question regarding reducing the no. of x-axis labels

2010-02-03 Thread Thomas Adams
All: I am using the command: qplot(date,MAE,data=data,facets=INTERVAL~type) which works fine except that the dates for my date axes are crunched together so much that they are unreadable. I can not find an option that I can set that will automatically reduce the x-axis labels to fit the

Re: [R] ggplot2/qplot question regarding reducing the no. of x-axis labels

2010-02-03 Thread Thomas Adams
Dennis, Thanks; that's just what I needed! Regards, Tom Dennis Murphy wrote: Hi: ggplot2 works with dates of class date and times of class POSIXct. See scale_x_date(), which is discussed in section 6.4.2 of the ggplot2 book. HTH, Dennis On Wed, Feb 3, 2010 at 11:15 AM, Thomas Adams

Re: [R] Are there free R webinar recordings somewhere ?

2009-12-07 Thread Thomas Adams
Tal, You might look here: http://www.fort.usgs.gov/brdscience/LearnR.htm I just googled r stats webinar Tom Tal Galili wrote: Hi all, A friend just sent me this: http://www.mathworks.com/company/events/webinars/index.html?id=language=en

Re: [R] From R to LaTeX to pdf?

2009-11-24 Thread Thomas Adams
Joel, You should consider using Sweave: http://www.stat.umn.edu/~charlie/Sweave/ -or- http://en.wikipedia.org/wiki/Sweave Regards, Tom Joel Fürstenberg-Hägg wrote: Hi all, Anyone experienced in the LaTeX format? I'm trying to use the xtable package to create nice anova tables,

Re: [R] Problem building R 2.10.0 - Matrix package

2009-11-13 Thread Thomas Adams
, but ... RPH == R P Herrold herr...@owlriver.com on Thu, 12 Nov 2009 14:13:26 -0500 (EST) writes: TA == Thomas Adams thomas.ad...@noaa.gov TA Attached is the output file from building R 2.10.0 on RedHat Linux. I TA have never previously experienced any problems when

[R] Problem building R 2.10.0 - Matrix package

2009-11-02 Thread Thomas Adams
All: Attached is the output file from building R 2.10.0 on RedHat Linux. I have never previously experienced any problems when building R from source with new releases. But, now I get a compile error with the Matrix package: CHOLMOD/Include/cholmod.h:87:22: error: UFconfig.h: No such file

Re: [R] Temperature Prediction Model

2009-10-22 Thread Thomas Adams
Aneeta, You will have to have a seasonal component built into your model, because the seasonal variation does matter, particularly -where- you are geographically (San Diego, Chicago, Denver, Miami are very different). Generally, there is a sinusoidal daily temperature variation, but frontal

Re: [R] Google's R Style Guide

2009-08-28 Thread Thomas . Adams
Esmail, Very nice; thanks! Tom - Original Message - From: Esmail esmail...@gmail.com Date: Friday, August 28, 2009 8:59 am Subject: [R] Google's R Style Guide Perhaps most of you have already seen this? http://google-styleguide.googlecode.com/svn/trunk/google-r- style.html

Re: [R] Wind-data analysis with R?

2009-08-20 Thread Thomas Adams
Oliver, You may consider looking at the climatol: Some Tools for Climatology (http://cran.r-project.org/) contributed package. This may come close to what you're looking for. I would suggest using R with GRASS GIS (http://grass.osgeo.org/), which is a powerful combination. Tom Oliver

Re: [R] USGS stream flow data automatic download R

2009-07-29 Thread Thomas . Adams
Stephen, I don't have anything or know of anything in R to do this. But, we download USGS streamflow data routinely. I have a Perl script that will reformat the downloaded data into a R-importable format (basically two columns date/time flow value). Are you interested in mean daily or

[R] Plotting a time series

2009-04-01 Thread Thomas Adams
I have data that I read in using: data-read.table(RAVK2.obs.data,sep=\t) 'data' looks like this: V1V2 1 2009-03-25 06:00:00 12.86 2 2009-03-25 12:00:00 12.80 3 2009-03-25 18:00:00 12.76 4 2009-03-26 00:00:00 12.68 5 2009-03-26 06:00:00 12.66 6 2009-03-26 12:00:00

Re: [R] popular R packages

2009-03-07 Thread Thomas Adams
I don't think At least one of the participants in the 2004 thread suggested that it would be a good thing to track the numbers of downloads by package. is reasonable because I download R packages for 2 home computers (laptop desktop) and 2 at work (1 Linux 1 Mac). There must be many such

Re: [R] The Origins of R

2009-02-04 Thread Thomas Adams
Rolf, Yes, that's what I was referring to as well… Cheers! Tom Rolf Turner wrote: On 4/02/2009, at 8:15 PM, Mark Difford wrote: Indeed. The postings exuded a tabloid-esque level of slimy nastiness. Indeed, indeed. But I do not feel that that is necessarily the case. Credit should be

Re: [R] The Origins of R

2009-02-03 Thread Thomas Adams
John, I certainly had that same impression of mischief making — I would call it trolling with the intent of trying to discredit R, its developers contributors. Mischief making indeed! Regards, Tom John Maindonald wrote: In another thread on this list, various wild allegations have been

Re: [R] R in the NY Times

2009-01-07 Thread Thomas Adams
Wacek, One would hope that if someone were to use software to build engines for aircraft, that said person would sufficiently test the software to have confidence in it, whether it had a Warranty or not — at least that's my mode of operation… Cheers! Tom Wacek Kusnierczyk wrote: Kevin E.

Re: [R] Boxplot 5% and 95% quantile instead of 25% and 75%

2008-08-21 Thread Thomas Adams
Joris, I found this (http://ceae.colorado.edu/~balajir/r-session-files/) on the web. It will do exactly what you want. Get the files: myboxplot-stats.r myboxplot.r Leesferry-mon-data.txt = example data The usage is: #Boxplots #Source the ‘myboxplot’ codes from Balaji’s directory.

  1   2   >