similar to: distribution of daily rainfall values in binned categories

I'm using barplot with the following call: barplot(stat_data[[5]][,],axes=TRUE,axisnames=TRUE,axis.lty=1,xlab=xlab,ylab=ylab,beside=TRUE,las=1,font.lab=2,font.axis=1,legend.text=TRUE) On some data, the vertical axis does not cover the whole plot area and the last tick mark is smaller than the maximum value. I tried setting the ylim values but even with that, some plots are still not OK,

Hi @all, I am biting my nails with the following problem: I have a data set of daily rainfall measurements for the last 20 years. What I want to do is calculate the daily cumulative sum of rainfall but only for every year which means that the cumulative sum has to be reset each year. After the calculations I want to plot each year of cumulative rainfall as a separate line in one graph preferably

I have a time series of rainfall in a dataframe. I would like to be able to aggregate this using a sliding window approach- i.e. a new 24 hourly total is calculated for each hours rainfall. I'm struggling to understand how this might be achieved - currently I've tried looping a sum function to re-calculate a new total at every stage of the loop. for (inp[[9]] in

Dear List Members, I need a bit help about fitting some theoretical distributions (such as geometric, exponential, lognormal or weibull distribution) to the following *dry spell*, *wet spell*, *cycles (Wet-Dry or Dry-Wet)* from my meteorological (daily rainfall) data http://www.angelfire.com/ab5/get5/R.rainfall.txt only for rainy seasen (july - september) of 14 years only:

How can I change the position of the ylab, after enlarging the margins with par(mar=...)? Here is the relevant code snippet ---- par(mar=c(5.1,5.1,4.1,2.1)) plot(c(1979,2003),c(40,50),ylim=c(1,73),lab=c(20,10,1),pch=21,col='blue',bg='blue',axes=FALSE,xlab="Years",ylab="Onset/Withdrawl Date",font.lab=2) box() axis(1,las=2)

Dear R users, I need to produce rainfall maps using R. I know that this is possible, I looked though the web, I found the example below reported (the author is Andrew Tredennick). I would ask you if this is the most performing way to make rainfall maps; if yes would someone be able to give me an example of how file.asc and pointfile.csv should be? If no would somebody please show me another way

Fwiw the engine behind geom_raster needs explicit observation-per-row form for input (with no structural normalization), so conversion to points is perfectly proper here, albeit confusing in context. (It's closer to what graphics devices actually use ultimately, but the expansion is laid out very early in ggplot2 because there's no standard for intermediate forms.) Cheers, Mike On Wed,

Hi, You might get more help from the R-sig-geo list, which is devoted to spatial topics. However. The *.asc file is an ArcGIS raster export format. You should use whatever the appropriate import commands are for your own gridded rainfall data. If you have a different format, you might or might not be able to import it directly with raster. ?raster will tell you more about the kinds of formats

Thank you Sarah and Mike for your explanations. My final objective is to produce maps (png image or any kind of extension I can import in LaTeX) where rainfall data are interpolated, using the Inverse Distance method or Kriging. My input file (pointfile.csv in the reported example) reports the station code, lat and long of the meteorological station and the rainfall value (which might be the

Hello, I have some daily rainfall data from rain gages. E.g.: date p1 p2 p3 20/04/2002 10.2 8.6 6.3 21/04/2002 0.4 1.6 1.4 22/04/2002 0.2 0.0 0.4 23/04/2002 5.2 1.0 0.2 I'd like to plot them in a column plot by day, to be able to compare them. I made an example in exell, that I attached, but can't find out how to do it in R. The data vector is a POSIXct. Thanks, Wouter

I need sequences that have gaps in them, such as the following: 4 5 6 | 12 13 14 | 20 21 22 a simple question, I've been scratching my head for a R function that will do this The : and seq do not allow this, and the c() can be used although not in an automatic way. I'm sure there is a way to do it without using a for() construct. Thank you

I am a student doing my MSc Research Methods, i am working on my thesis research on analysing and modelling of crop failure risks due to drought in selected districts in Malawi. The analysis and modelling will focus on two crop stages of development: just after planting and flowering stages. I have rainfall, temperature, Relative humidity, wind speed, coordinate location data for the past 30 years

I've been trying to split a space delimited string with double-quotes in R for some time but without success. An example of a string is as follows: /rainfall snowfall "Channel storage" "Rivulet storage"/ It's important for us because these are column headings that must match the subsequent data. Here is some code I've been trying: str <- 'rainfall

Hi, For my master thesis I have 24 micro-plots on which I did measurements during 3 months. The measurements were: - Rainfall and runoff events throughout 3monts (runoff being dependant on the rainfall, a coefficient (%) has been made per rainfall event and per 3 months) - Soil texture (3 different textures were differentiated) - Slope (3 classes of slopes) - Stoniness (one time measurement)

Dear All, I am a research student in environment. I have only little programming knowledge. I am currently doing the last project about rainfall impact on ground water quality in an area. It happens that I have to use R to read rainfall data (3 dimension) from ASC file (*.asc), and then write them into one NCDF file (*.nc). I have been working very hard on study R, but I

Hi I think I have been struggling to use replace correctly, I usually work my way around this using a loop, but I think this is in fact inefficient. I have a dataset with runoff from three plots and associated rainfall. However either the datarecording was sloppy, or the rainfall very patchy. So I am trying to remove data from my dataset for which the runoff is larger than the rainfall on the

Dear List Members, I need some help about programming in S language. My problem is as follows: I have meteorological data (about rainfall measurement each day from 1989-2002), say like http://www.angelfire.com/ab5/get5/data.rainfall.txt or http://www.angelfire.com/ab5/get5/R.rainfall.txt in a sequence of 0(denoting dry day)'s and 1(denoting wet day)'s. I want to construct a frequency

Hi, I would like to plot time against rainfall data (data is at the end) using xyplot. The basic code looks like this: xyplot(rainfall~time, type="a") When I do this, the graph looks ok except that the x-axis has too many values. I would just like to display the years and not the months on the x-axis. I've been fiddling around with 'scales', and read previous posts about

Hello, I'm looking for a way to plot different data series on one graph. I have a series of hourly rainfall and quarterly flow measurements (i.e. 4 times an hour) of a catchment. The rainfall should be plotted in bars, the flow as a line. Both on the same X axe (time) but with different Y axes. The problem is the plot() function does not support add=TRUE... Furthermore I'm not sure

type count 0 20 1 15 0 10 1 35 0 28 I would like to create two vectors from the data above. For example, type1=c(15, 35) and type0 = c(20, 10, 28). Can any one help Raphael