search for: catchment

...#####################################3 Dear all, I am exporting grass map into R thanks to the very useful spgrass6 package. library(spgrass6) # I have 3 map I am working with a MASK map of a specific area. # 1) a landuse map landuse<-readRAST6("landuse_mapname_in_grass") # 2) a catchment map which divide the area in several catchements catchment<-readRAST6("catchement_mapname_in_grass") # 3) a precipitation map precipitation<-readRAST6("precipitation_mapname_in_grass") # then I would like to sum the precipitation spatialy over each catchment and landuse...

...rying 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: http://r.789695.n4.nabble.com/cdplot-error-tp3714454p3714454.html Sent from the R help mailing list archive at Nabble.com.

...zaar_(software) Some projects that have moved to Mercurial (hg) are Mozilla Octave Python ... http://en.wikipedia.org/wiki/Mercurial_(software) Joel Spolky's take on it: http://www.joelonsoftware.com/items/2010/03/17.html Regards -Felix -- Felix Andrews / ??? Postdoctoral Fellow Integrated Catchment Assessment and Management (iCAM) Centre Fenner School of Environment and Society [Bldg 48a] The Australian National University Canberra ACT 0200 Australia M: +61 410 400 963 T: + 61 2 6125 4670 E: felix.andrews at anu.edu.au CRICOS Provider No. 00120C -- http://www.neurofractal.org/felix/

...n 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     stations",col="yellowgreen", cex.axis=1,     xlab="rainfall",ylab="Relative frequency", ylim= c(0,.003), xlim=c(0,1200))     lines(density(dd), lwd=3,col="red")     #legend("topright",c("observed","generated"...

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 what's the best format for input data. Now it are two data frames with a time (POS...

...LATE=English_Australia.1252 LC_CTYPE=English_Australia.1252 [3] LC_MONETARY=English_Australia.1252 LC_NUMERIC=C [5] LC_TIME=English_Australia.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base -- Felix Andrews / ??? Postdoctoral Fellow Integrated Catchment Assessment and Management (iCAM) Centre Fenner School of Environment and Society [Bldg 48a] The Australian National University Canberra ACT 0200 Australia M: +61 410 400 963 T: + 61 2 6125 4670 E: felix.andrews at anu.edu.au CRICOS Provider No. 00120C -- http://www.neurofractal.org/felix/

....org/posting-guide.html > >> > and provide commented, minimal, self-contained, reproducible code. > >> > > >> > > >> > >> > >> > >> -- > >> Felix Andrews / 安福立 > >> Postdoctoral Fellow > >> Integrated Catchment Assessment and Management (iCAM) Centre > >> Fenner School of Environment and Society [Bldg 48a] > >> The Australian National University > >> Canberra ACT 0200 Australia > >> M: +61 410 400 963 > >> T: + 61 2 6125 4670 > >> E: felix.andrews@anu....

...requency starting from a set date. Leafing through r-help mail archives I've found this *ALMOST* satisfactory message: ========================================================== ......... > I have a series of hourly rainfall and quarterly flow > measurements (i.e. 4 times an hour) of a catchment ........ > Maybe time series are easier, but in > > ts(data = NA, start = X,... > > X should be a number or a vector. how does this coresponds to a > data and hour (e.g. april 26,2002, 15:00:00)? If your observations are equidistant, e.g. you've got 24 hourly measurem...

...ses for larger ones), which violates one assumption of the OLS method. Which of the numerous options should I choose? glm, robust lm, ... The problem is kept simple for now. I only try to explain the log of local topographic slope (dependent variable) with regard to the distance to the outlet of a catchment (independent variable) for a fixed drained area. Both variables are continuous. I ordered Venables and Ripley 2002, which I suspect is a sound reading for advanced stats with R, but it has not arrived yet and I need to move on asap. Any advice or pointer to the appropriate literature is greatly ap...

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

...ev(ylim0) ylim <- xl tmp <- xaxt xaxt <- yaxt yaxt <- tmp } if (missing(xlim)) xlim <- xlim0 if (missing(ylim)) ylim <- ylim0 plot(0, xlim = xlim, ylim = ylim, ...... Regards, Felix -- Felix Andrews / ??? PhD candidate Integrated Catchment Assessment and Management Centre The Fenner School of Environment and Society The Australian National University (Building 48A), ACT 0200 Beijing Bag, Locked Bag 40, Kingston ACT 2604 http://www.neurofractal.org/felix/ 3358 543D AAC6 22C2 D336 80D9 360B 72DD 3E4C F5D8

...ambda) > { > H<-rbind(A,lambda*L) > i<-1:(nrow(L)+length(data)) > q<-ifelse(i<=length(data),data[i],0) > nntik<-nnls.fit(H,q) > return(nntik) > } I think this is the same as what pcls states: Willem -- Dr R.W. Vervoort McCaughey Senior Lecturer Hydrology and Catchment Management Faculty of Agriculture, Food and Natural Resources Rm 503, Watt Bldg. http://www.agric.usyd.edu.au/mccaughey Postal: Bldg A03 The University of Sydney, NSW, 2006 phone: +61 (2) 9351 8744 fax: +61 (2) 9351 5108 e-mail: w.vervoort at acss.usyd.edu.au

...ay 17 svn rev 48333 language R version.string R version 2.9.0 (2009-04-17) ----------------------------------------------------------------- Dr. R.W. Vervoort McCaughey Senior Lecturer Hydrology and Catchment Management Faculty of Agriculture, Food and Natural Resources, Bldg A04, The University of Sydney, NSW 2006 http://tinyurl.com/mccaughey http://blogs.usyd.edu.au/waterhydrosu ph: +61 2 9351 8744 fax: +61 2 9351?4953

...thods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: lattice "0.13-8" Regards Karl _________________________________ Dr Karl J Sommer, Department of Primary Industries, Catchment & Agriculture Services, PO Box 905 Mildura, VIC, 3502 Australia Tel: +61 (0)3 5051 4390 Fax +61 (0)3 5051 4534 Email: karl.sommer at dpi.vic.gov.au

...w 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: Precipitation_18o: 2 columns - column 1 is a date (Excel number format) and column 2 is an isotope ratio (i.e. -8.12) Runoff_18o: same as above Daily_Precip: 2 columns - column 1 is the same date format but column 2 is a weekly bulk precipitation value...

...r column. However I would like to use zoo instead, because of its many utility functions and easy conversion to ts. Is there any prospect of zoo handling such numeric + factor data? Other suggestions on elegant ways to do it are also welcome. Felix -- Felix Andrews / ??? PhD candidate Integrated Catchment Assessment and Management Centre The Fenner School of Environment and Society The Australian National University (Building 48A), ACT 0200 Beijing Bag, Locked Bag 40, Kingston ACT 2604 http://www.neurofractal.org/felix/ xmpp:foolish.android at gmail.com 3358 543D AAC6 22C2 D336 80D9 360B 72DD 3E4C...

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 to plot graphs of all the chemical variables from 1985-2006 to analyze the trends in the concentration of water in the area. I would like to request...

...Devices utils datasets [8] methods base other attached packages: [1] ggplot2_0.5.7 colorspace_0.95 RColorBrewer_1.0-2 MASS_7.2-38 [5] proto_0.3-8 reshape_0.8.0 loaded via a namespace (and not attached): [1] tools_2.6.1 -- Felix Andrews / $B0BJ!N)(B PhD candidate Integrated Catchment Assessment and Management Centre The Fenner School of Environment and Society The Australian National University (Building 48A), ACT 0200 Beijing Bag, Locked Bag 40, Kingston ACT 2604 http://www.neurofractal.org/felix/ 3358 543D AAC6 22C2 D336 80D9 360B 72DD 3E4C F5D8

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

...mostattributes(timelim) <- attributes(from) axis(side = 3, at = prettyDate(timelim), labels = prettyDate(timelim, do.format=TRUE)) }) devAskNewPage(FALSE) Is it appropriate / desirable for this to be incorporated into R? -- Felix Andrews / ??? Postdoctoral Fellow Integrated Catchment Assessment and Management (iCAM) Centre Fenner School of Environment and Society [Bldg 48a] The Australian National University Canberra ACT 0200 Australia M: +61 410 400 963 T: + 61 2 6125 4670 E: felix.andrews at anu.edu.au CRICOS Provider No. 00120C -- http://www.neurofractal.org/felix/