search for: catchments

Dear all, I have some additionale question concerning the spgrass6 package. * When you set a region in GRASS, does the readGRASS6 function in R only load data contained in the zoomed region or the whole map ? * When you have a MASK map in grass, does the readGRASS6 function in R only load data contained inside the MASK area ? Could this be the problem ? Thanks, Jessica

Hello I have lots of folders of individual Scottish river catchments on my uni One Drive. Each folder is labelled with the river name eg "Tay" and they are all in a folder named "Scotland" I want to list the folders on One Drive so that I can cross check that I have them all against a list of folders on my laptop. Can I somehow use list.files()...

Fairly new at this. Trying to create a conditional density plot. >cdplot(status~harvd.l,data=phy) Error in cdplot.formula(status~harvd.l,data=phy): dependent variable should be a factor What does this error mean? Status is a binary response of infestation (0/1) and harvd.l is the log of timber harvest density per catchment. Thanks. -- View this message in context:

Hi, Just wondering whether anyone had thought about moving the R sources to a "distributed" version control system such as Bazaar, Git or Mercurial. These new generation systems make it easier to work on feature branches, allow working offline, are very fast, etc. Some projects that have moved to Git are Linux Kernel Perl Ruby on Rails ... http://en.wikipedia.org/wiki/Git_(software)

Hi, Does anybody encounter the same problem when we overlap histogram and density     that the density line seem to shift to the right a little bit?           If you do have the same problem, what should we do to correct that?           Thank you.           par(mar=c(4,4,2,1.2),oma=c(0,0,0,0))     hist(datobs,prob=TRUE, main ="Volume of a catchment from four    

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

Hi, I am getting unexpected behaviour from aggregate.ts(). The 'ndeltat' argument is effectively being reduced by 1 in some cases, even when it is an integer, with the result that the blocks to be aggregated are not of the expected size, and also that the end() of the aggregated series is much later than the end() of the original series. rawts <- ts(rep(1:10, each = 5), start = 1) ##

I HAVE BEEN ASKED TO BE REMOVED FROM THIS LIST OVER 20 TIMES, PLEASE CAN I BE REMOVED!!!!!!!!!!!!! --- On Thu, 3/12/09, Peng Cai <pengcaimaillist@gmail.com> wrote: From: Peng Cai <pengcaimaillist@gmail.com> Subject: Re: [R] Help: barchart() {Lattice} To: "Felix Andrews" <felix@nfrac.org> Cc: r-help@r-project.org Received: Thursday, 3 December, 2009, 10:17 AM Got it

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

Dear R-users, There are tons of methods out there for fitting independant variables to a dependent variable. All stats books tell you about the assumptions behind OLS (ordinary least squares) and warn against abusive use of the method (which many of us do disregard by lack of a better knowledge). Most introductory text books stop there and don't tell you what the next best option might be. I

...999-12-01")&index(aa)<=as.Date("2000-02-28")) dd <- subset(aa, index(aa)>=as.Date("2000-12-01")&index(aa)<=as.Date("2001-02-28")) ee<- rbind(bb,cc,dd) The method above appears to do the job just fine except that I have around 30 locations (catchments) each with varying data availability and some with over 20 years worth of data. Ideally I would like to combine the second set of commands into a single command where I specify the start and end year and the months that I am interested in. Any advice on this matter would be greatly appreciated. M...

list(...), I would like to zoom in to the leaves of large trees in a dendrogram plot. The playwith package allows zooming by passing xlim and ylim arguments to the plot call (Hmisc does this too I think). But currently stats:::plot.dendrogram does not accept xlim or ylim. So I would like to enable that. In place of the existing code chunk: xlim <- c(x1 - 1/2, x2 + 1/2) ylim <- c(0,

Hi all, I am trying to do an inversion of electromagnetic data with non-negative least squares method (Tikhonov regularisation) and have got it programmed in S-Plus. However I am trying to move all my scripts from S-Plus to R. Is there an equivalent to nnls.fit in R? I think this can be done with pcls? Right? S-Plus script: A, L and data are matrices, lambda is a vector of possible lambda

Dear all, I am not sure what I am doing wrong, but I have some unexplained behaviour when saving a lattice graph including text to a pdf file. The text seems to move around. It must have something to do with the way coordinates are set in devices other than jpg. Any suggestions would be helpful Willem Here is some example code setwd("c:/willem/research/misc") today <-

Hello list, I have encountered a problem trying to save graphs using the R-graphics menu: File|Save as. The menu suggests that files may be saved as either Metafile, Postscript, pdf, png, bmp, jpeg. When I specify any of those file formats a menu comes up requesting a file name. After providing a name R invariably hangs and has to be restarted. I am able to save files under the various

Hello all, I'm relatively new to "R" and programming in general - I had previously used MatLab, but decided to make the transition to R, as the computational times are much better! Anyway, I'm trying to use R to run a gamma distribution model to estimate mean transit times of water moving through a hydrological catchment. My input are 3 .txt format files as follow:

list(...), I am working with environmental time series (eg rainfall, stream flow) that have attached quality codes for each data point. The quality codes have just a few factor levels, like "good", "suspect", "poor", "imputed". I use the quality codes in plots and summaries. They are carried through when a time series is aggregated to a longer time-step,

Good Day The attached document refers. I have recently attended an introduction course to R, and now wish to apply the program into my work analysis. The spreadsheet contains water quality data sampled from different locations within a catchment. I wish to group the variables per point (Point Name), Year (all samples per year eg 1985 etc) These should be averages of each year. Lastly, I wish

list(...) I can not get at the grid viewports in a ggplot2 plot. I know there is supposed to be a viewport called "panel_1_1" but downViewport() can't find it. Has the viewport been popped? I had a quick look at the functions involved (eg ggplot_plot) but there is no obvious problem there. > library(ggplot2) > qplot(1:10, 1:10) > current.viewport() viewport[ROOT] >

...ently printing two sig figs), could someone also explain how I can consistently get two significant figures in situation like above? Thanks as always. P. -- Dr Peter Dunn | dunn <at> usq.edu.au Faculty of Sciences, USQ; http://www.sci.usq.edu.au/staff/dunn Aust. Centre for Sustainable Catchments: www.usq.edu.au/acsc This email (including any attached files) is confidentia...{{dropped:15}}