Re: [R] GIS in R

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: https://g

Re: [R] forecast accuracy

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

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 <- 2

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

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 https://www.res

[R] Problem with my function using as.POSIXct

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 } tdate<-paste(c

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

t;, "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", > "2012-04-18 14

[R] Help needed with aggregate or other solution

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 col

Re: [R] about installing smwrGraphs package

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 tr

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

' > >yr_2000$year<-'2000' > >yr_2001$year<-'2001' > >yr_2002$year<-'2002' > >yr_2003$year<-'2003' > >yr_2004$year<-'2004' > >yr_2005$year<-'2005' > >yr_2006$year<-'20

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

#x27;2007' yr_2008$year<-'2008' yr_2009$year<-'2009' yr_2010$year<-'2010' yr_2011$year<-'2011' yr_2012$year<-'2012' yr_2013$year<-'2013' yr_2014$year<-'2014' yr_2015$year<-'2015' yr_2016$year&l

[R] Histogram plots in Lattice with spatialgrid dataframe data

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: spplot(ann_bias,xlim=

Re: [R] Help with reshape/reshape2 needed

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

[R] Help with reshape/reshape2 needed

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 1974

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

n 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 packages that are very flexible. You c

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

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 or

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

p > TIBCO Software > wdunlap tibco.com > > On Thu, Jun 2, 2016 at 8:36 AM, Thomas Adams 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 luck. >> >> I have genera

[R] Struggling trying to plot points on boxplot

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 observ

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

an 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.org] On Behalf Of Th

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

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 want

Re: [R] question about probplot in e1071 package

plot(x, line=FALSE) > #probplot(y, line=FALSE) > qq.y <- qqnorm(y, plot=F) > points(qq.y$y, qq.y$x) > > I hope this is useful for you. > > Best Regards, > > Luisfo Chiroque > > PhD Student > IMDEA Networks Institute > > http://fourier.networks.imdea.org/peo

[R] question about probplot in e1071 package

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) * *probplot

[R] NULL dev.lis()

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 S

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

nv >TRUETRUETRUETRUETRUETRUE > NLS profmem cairo ICU long.double libcurl >TRUE FALSETRUETRUETRUE FALSE > > On 10/14/2015 4:00 PM, Evan Cooch wrote: > > > >

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

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 wrote: > Tom -- > > On 10/14/2015 3:35 PM, Thomas Adams wrote: > > Evan, >

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

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 th

Re: [R] Best Mac for R

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 compared

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

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 http://cran.r-project

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

e.synnerg...@his.se > mobile: +46 (0)708 806495 > webpage: his.se/synj > > Postal Address: Visiting address: > Box 408 Kanikegränd 3A > 541 28 Skövde Skövde > > From: Thomas Adams > Date: Sunday, December 21, 2014 2:06 AM > To: Jane Synnergren > Cc: R-help >

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

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 wrote: > > Hi, > I get following error when tryi

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

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 wrote: > > This morning I was reading Jeff Leek's list of awesome things other people > did in 2014 at http://simplystatistics.org/?p=3696 (

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

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 wrote: > I wis

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

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] ubuntu 14.04

As do I... On Wednesday, September 24, 2014, Jeff Newmiller 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?

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

Simona, You need to install the dependencies: install.packages("gstat",dependencies=T) Tom On Thu, Oct 10, 2013 at 11:58 AM, Simona Augyte wrote: > Hello, > # I am able to > install.packages("gstat") > #but when I try to upload I get an error message > library(gstat) > #Error in loadNamespace

Re: [R] reshape2 issue continued

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 s

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

gt;Kevin Lang once showed how a remarkably consistent paucity of >>strong evidence in published economics studies made it unlikely >>that any of them were right)." >> >> Bill Dunlap >> Spotfire, TIBCO Software >> wdunlap tibco.com >> >> &g

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

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 func

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

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 wrote: > In case you haven't noticed, this is making the rounds in the media, > including a handful of ref

Re: [R] Odd graphic device behavior

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 >

Re: [R] Odd graphic device behavior

the other 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 wrote: > Luca, > >

Re: [R] Odd graphic device behavior

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 wrote: > Thomas, > any chance that you could have problems with your graphic driver? For > insta

Re: [R] Odd graphic device behavior

ated in the resized window, 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 wrote: > > On Mar 27, 2013, at 18:11 , David Winsemius wrote: > > > > > On Mar 27, 2013

Re: [R] Odd graphic device behavior

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 window

[R] Odd graphic device behavior

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 d

Re: [R] Running other programs from R

Use the system() command. e.g." system("Multilog") On Sat, Mar 16, 2013 at 5:09 PM, Sedat Sen 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 just want to c

Re: [R] Using grImport to create a watermark

Brian Ripley wrote: > 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… >

[R] Using grImport to create a watermark

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(…) success

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

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 wrote: > Dear R-helpers, > > I am trying to download some netcdf4 fi

Re: [R] On Reproducible Code

not have supplied 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 wrote: > > On Jul 30, 2012, at 13:05 , Thomas Adams wrote: > > > > > Nothing is

Re: [R] On Reproducible Code

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 th

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

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 wrote: > Jie Tang writes: > > > hi R-users > >Now I have a figure in emf or png or tiff format that

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

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 th

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

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> wrote: > http://www.cs.purd

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

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

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 wrote

Re: [R] RSqlite UPDATE command problem

dbGetQuery(con, sql) > > instead... > > b > > On 27 March 2012 14:18, Thomas Adams 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)

[R] RSqlite UPDATE command problem

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)) dbWrit

Re: [R] Struggling with zoo and aggregate

Gabor, That does it! I can't thank you enough… Many thanks, Tom On Mon, Mar 26, 2012 at 7:22 AM, Gabor Grothendieck wrote: > On Sun, Mar 25, 2012 at 10:20 PM, Thomas Adams > wrote: > > Gabor, > > > > Thank you for your help -- it did help me a lot. However, with

Re: [R] Struggling with zoo and aggregate

Tom On Sat, Mar 24, 2012 at 10:58 PM, Gabor Grothendieck < ggrothendi...@gmail.com> wrote: > On Sat, Mar 24, 2012 at 10:44 PM, Thomas Adams > wrote: > > All: > > > > I have a SQlite database where I have stored some verification data by > date > > & time (cy

[R] Struggling with zoo and aggregate

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

Sarah, Thanks! Great stuff... Tom On Wed, Mar 7, 2012 at 10:07 AM, Sarah Goslee wrote: > 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 software package."

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

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, to

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

the lattice ?bwplot function for th > eplot. > > Cheers, > Bert > > > > > result > > > > On Mon, Dec 12, 2011 at 1:43 PM, Thomas Adams > wrote: > >> Ryan, > >> > >> I think you could do what you want by having the vector data wr

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

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, bru

Re: [R] Question about ggplot2 and stat_smooth

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

previous email where I can specify the limits calculated for the 'boxplots' using 5%, 25%,75%, 95% limits as we have with the 'boxplots'. Tom - Original Message - From: Hadley Wickham Date: Tuesday, October 4, 2011 10:23 am Subject: Re: [R] Question about ggplot2 and

Re: [R] Question about ggplot2 and stat_smooth

ious email where I can specify the limits calculated for the 'boxplots' using 5%, 25%,75%, 95% limits as we have with the 'boxplots'. Tom - Original Message - From: Hadley Wickham Date: Tuesday, October 4, 2011 10:23 am Subject: Re: [R] Question about ggplot2 and stat_

Re: [R] Question about ggplot2 and stat_smooth

ctor(cyl))) c + stat_smooth(aes(group=cyl))+stat_smooth(aes(fill=factor(cyl)))+geom_point() Andrés AM 2011/10/3, 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="darkblu

[R] Question about ggplot2 and stat_smooth

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 hea

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

- Original Message - From: gaiarrido 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 > Thanks very much

Re: [R] Thiessen Method

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 (s

Re: [R] recommendation on r scripting tutorial?

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 Date: Saturday, April 2, 2011 11:22 am Subject: [R] recommendation on r scripting tutorial? To: r-hel

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

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 sw

[R] 2 questions about probplot in package e1071

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, line=FALSE

[R] Help needed with plot axis labeling

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

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 m

Re: [R] Layout of mulitpage conditioned lattice plots

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

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 - Orig

Re: [R] two questions for R beginners

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 on-the-

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

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 ~ t

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

quot; or even close to minimal.) -Peter Ehlers Thomas Adams wrote: I am trying to use t.test on the following data: datetypeINTERVALnCASESMTFSDFMTOSDO nFSTMFnOBSMOMBBIASCVBIASEVMEMAE RMSECRCF 2001-06-15avnGE1.00

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

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 <mailto:tho

[R] help needed using t.test with factors

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] ggplot2/qplot question regarding reducing the no. of x-axis labels

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

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

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 avail

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

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?

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, bu

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

Russ, but ... "RPH" == R P Herrold on Thu, 12 Nov 2009 14:13:26 -0500 (EST) writes: >>>>>>> "TA" == Thomas Adams TA> Attached is the output file from building R 2.10.0 on RedHat Linux. I TA> have never previously experien

[R] Problem building R 2.10.0 - Matrix package

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 or

Re: [R] Temperature Prediction Model

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 p

Re: [R] Google's R Style Guide

Esmail, Very nice; thanks! Tom - Original Message - From: Esmail 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 > Comments/Critiques?

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

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 Bandel

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

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 instantan

[R] Plotting a time series

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:

Re: [R] popular R packages

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 suc

Re: [R] The Origins of R

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 gi

Re: [R] The Origins of R

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

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

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

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. source("

  1   2   >