similar to: Aggregating on Water Year Rather Than Calendar Year

A data set I obtained has the hours running from 01 through 24 rather than the conventional 00 through 23. My favorite, strptime, balks at hour 24. I thought it would be easy to correct but it must be too late on Friday for my brain and caffeine isn't helping. TIA for a hint, Clint -- Clint Bowman INTERNET: clint at ecy.wa.gov Air Quality Modeler INTERNET: clint at math.utah.edu

I have a 2 by 226200 table, conveniently read in by read.table(), which exhibits some strange behavior when plotted by plot(V1,V2). The general pattern for the range of windspeeds, [0<V1<50] is as expected -- the wind gust falls in the interval [V1<V2<65] except for certain values of V2. For V2 == c(15,26,37,48,59), the V2 values are positioned at one-tenth of the V1 (i.e., as if I

I have x, y, z data. The x, y fields dont change but Z does. How do I add a very small number onto the end of each x, y data point. For example: Original (X) Original (Y) Original (Z) 15 20 30 15 20 40 New (X) New (Y)

Hello everyone I have a dataframe with 4 colums (year, month, day, runoff) for 1993-2009. Now I like to calculate the average runoff for each day. Finally I like to plot the median runoff for all this years. I tried with some loops, but it didn't work. Do you have any Tips for my problem? Any help is greatly appreciated! Thank you very much. Dominique -- View this message in context:

I have a large data structure that looks like: > strsplit(st,",")[14395] [1] "KGEG" [2] "SA => KGEG" [3] "72785" [4] "47.62139" [5] "-117.52778" [6] "723" [7] "WA" [8] "US" [9] "2" [10] "SPOKANE SPOKANE INTERNATIONAL AIRPORT" [11] "1" I'd like to be able to

I have shapefiles for the state climatic divisions for the United States and read.shape brings them in wonderfully. Now I wish to run through a list of several thousand observation sites to find out in which division each is located. I figure that I can compute the winding number for each site and be done. However a search doesn't find any references and I can't find a winding number

Folks: I thought the following excerpt from Bruce McCullough's post would be a good candidate for the R fortunes package -- except that it's about Excel, not R! So I nominate it... but leave it to others to say whether it's really "qualified" to be nominated. ---- "The idea that the Excel solver "has a good reputation for being fast and accurate" does not

Dear all,   I have the following written example > coords <- "51°30'48.58\"N" > > as.integer(strsplit(coords, "°")[[1]][1]) [1] 51 > as.integer(strsplit(strsplit(coords, "°")[[1]][2], "'")[[1]][1]) [1] 30 > as.numeric(strsplit(strsplit(strsplit(coords, "°")[[1]][2], "'")[[1]][2],

Does anyone have information about using the National Practitioner Data Bank (NPDB) in R? It's public use data available from the U.S. Department of Health and Human Services that can be downloaded as a .dat or .por file. I can't even figure out how to open the file in R. [[alternative HTML version deleted]]

How can I get the outlier in this boxplot of "Score" to be represented by the corresponding value in "SubNo"? score=c(6,6,7,14,5,7,6,8) SubNo=1:8 mydata=data.frame(SubNo, score) boxplot(mydata$score) Thanks! Kevin [[alternative HTML version deleted]]

Dear All, I have a query : what is the command to count number of repeated words in a column. for ex: a = oranges oranges apples apples grape oranges apple pine the result should be oranges 3 apples 3 grape 1 pine 1 is there an easy way for this. Thanks, Nataraju GM R & D Bangalore -- "No relationship is Static .. You either Step up or Step down" [[alternative HTML

I have a file with 3 cols as: x y value 1 1 0.5000 1 2 1.8333 1 3 4.3333 2 1 0.7500 2 2 0.7500 2 3 4 3 1 0.7500 3 2 1.4167 3 3 5.1667 I read them in using p <- read.table("data.txt", header=TRUE) My questoin: how to draw a image to show the values at different positions (x, y)? I would like to see different colors for

Hello there, would you please look into my codes? Here I have following: > set.seed(100) > samp <- sample(c(1,-1,0), 20, replace=T); samp [1] 1 1 -1 1 -1 -1 0 -1 -1 1 -1 0 1 -1 0 0 1 -1 -1 0 Here I want to calculate the length of each unique number for above vector. How can I do that? Thanks in advance [[alternative HTML version deleted]]

I'm trying to plot smooth trend using smoothTrend in OpenAir but I'm having problems. I used the following code. --------------------------------------------------------------------------------- #Set my working dir to the dir with my files setwd("c:/R") #Load the openair library library(openair) #Load the data mydata <- read.table("MCNP-pH.csv", header=TRUE,

Hi All: This follows my earlier post on webized slides on lecture presentation on introducing R. I learned that in Mozilla (Firefox) browsers, the slides did not show up. Sorry for the no show. As a reluctant windows user, I kind of carelessly clicked through Powerpoint to convert the presentation file to its html form, unwittingly leading to the mess. See if it got corrected now (I do not have

Hello, I'm currently trying to find a method to interpolate or smooth data that represent a trajectory in space. For example, I have an ordered (=time) set of (x,y) tuples which constitute a path in a 2D space. Is there a way using R to interpolate between these points in a way similar to spline interpolation so that I get a smooth path in space? Greetings, Dieter -- Dieter Vanderelst

Dear List I hope you can help me: I?ve got a dataframe (df) within which I am looking for Peak Over Threshold values as well as the length of the events. An event starts when walevel equals 5.8 and it should end when walevel equals the lower threshold value (5.35). I tried ?clusters (?)? from ?evd package?, and varied r (see example) but it did not work for all events (again

Do you know of any references that verify the accuracy of R for basic statistical calculations and tests. The results of these studies should indicate that R results are the same as the results of other statistical packages to a certain number of decimal places on some benchmark calculations. Thanks, Corey Powell Clinical Data Analyst Broncus Technologies cpowell@broncus.com [[alternative

Attention R users, especially those that are experienced enough to be opinionated, I need your input. Consider the following simple plot: x <- rnorm(100) y <- rnorm(100) plot(x, y, bty='n') A colleague (and dreaded SAS user) commented that she thought that my plots could be "cleaned up" by connecting the X and Y axes. I know that I can do that with bty='l' but I

I have a matrix where one column has a date and another column has a time. I need to delete all times before 6am. I had combined the Date and Time column into DateTime. Mat1: Weight Date Time 7.6 04/28/11 09:03 8.4 04/29/11 03:11 8.6 04/29/11 05:32 8.6 04/29/11 09:53 1.4 05/01/11 19:52 I tried just picking up the time: as.POSIXct(Mat1$Time, format = "%H:%M") but