similar to: different data series on one graph

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

Hello, I'm doing some linear regression: >lm<-lm(osas~alp,data) >summary(lm) However, the Rsquared in the output of summary() is not the same as the "standard" Rsquared calculated by spreadsheets, and outlined in statistical guidebooks, being SSR/SSTO. The output says "multiple Rsquared", but it is no multiple regression... What's the difference? Thanks,

hello, maybe a silly question, but how do you set a default directory of R, i.e. the place where R looks for source() files etc. at startup? R 1.5.0 on Mandrake linux 8.2 BTW, it seems that my doc-directory is crippled (a lot of files are lacking). A flaw of the mandrake 8.2-rpm's or a bad install? thanks! Wouter ---------------------------------------------------------------

Hi list, there always seem to be a small margin between the actual plot and the axes/box around the plot. For instance, in plot(c(0,1),c(0,1)) the 2 points do not lay on the box around the plot but slightly more to the center. Which parameter controls this margin and is it possible to eliminate it? thanks, Wouter

I'm trying to make a transition from S-plus under Windows to R under Linux. My immediate aim is to produce a scatter plot with two y-axes with different scales. That can be done in S-plus with a command of the following form: guiPlot( PlotType="Scatter", DataSet="execon", Columns="years,hstart,ship", AxisType="Multiple Y" ), where years, hstart, and

Hello, I'd like to ask, how to adjust the space between symbol/line and text in a legend. My legend() is build like that legend(xmax-netzdiff/11, # shouldn't matter here ifelse(ydiff>=0, # ymin-6.5*(netzdiff/11), # ymax-netzdiff/11), # c("stehender Stamm","Windwurf",

Having to collect hourly electricity loads and quarter-of-an-hour electricity production data for some years I think that the tidiest way of doing it is to resort to ts but I don't know how to define such a frequency starting from a set date. Leafing through r-help mail archives I've found this *ALMOST* satisfactory message: ==========================================================

Hello, (maybe a quite basic statistical question, but I'm just struggling with it) I'm doing an anova: > Res1<-aov(H2O~location+topo+horizon+pF+Error(location:topo:horizon)) (water retention of soils) which gives a significant difference at factor "location". Which function should I use to now which locations (there are 3: A, B and C) differ significantly and which do

When I try to install the tinc RPM I am missing [root@relay buytaert]# rpm -vih tinc-1.0pre2-1.i386.rpm error: failed dependencies: libgmp.so.2 is needed by tinc-1.0pre2-1 [root@relay buytaert]# ls /usr/lib/libgmp.so.3 /usr/lib/libgmp.so.3 [root@relay buytaert]# rpm -qf /usr/lib/libgmp.so.3 gmp-3.1.1-3 This seems to be a too recent version. So I tried compiling it and got the following

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

In 2.4 there used to be a megaide driver for certain ata raid controllers. In 2.6 I can't access more than my /dev/hda , I can see my other 4, I can create filesystems on them, but I can't mount them Anyone knows how to use these on a 2.6 kernel (Centos 4.2) greetings Kris 00:02.0 RAID bus controller: Silicon Image, Inc. SiI 0649 Ultra ATA/100 PCI to ATA Host Controller (rev 02)

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

Hi, I'm a newbie in using R and I would like to have a few clues as to how I could compute and plot a distribution of daily rainfall intensity in different categories. I have daily values (mm/day) for several years and I need to show the frequency of 0-1, 1-2.5, 2.5-5, 5-10, 10-20, 20+ mm/day. Can this be done easily? Thanks, Etienne

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