search for: vervoort

...month 04 day 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

...s > "nntik"<-function(A,L,data,lambda) > { > 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...

Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 > -51^B [1] -3.190824e+81 # seems fine # now this: > x <- seq(-51,-49,length=100) > x^B [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN <snip> > is.numeric(x^B) [1] TRUE > is.real(x^B) [1]

Hi, Found the answer to the background myself: 'trellis.par.get' and 'trellis.par.set' background<-trellis.par.get("background") background$col<-"white" trellis.par.set("background", background) 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...

Dear R-help list, I am trying to construct a lattice histogram using 3 factors. My dataframe looks like this: (simulating a waterbalance over groundwater with different salinities) s days model EC EC_max 0.4 1 "A" 10 9 0.42 2 "A" 10 9 0.44 3 "A" 10 9 : : : :

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