search for: basin

...to assess a variable who's name is constructed with paste. The problem is following: In a data frame, there are 12 columns whose title vary only by a digit, all other part being equal. I would like to compute a operation on a subset these variables and parse them in turn. the data frame "basin.param" contains columns called ratio1, ratio2, ratio3... in each of these ratios I want to retain only values above 0. This could be put like this crit<- which(basin.param$ratio1 > 0) basin.param$ratio1[crit] + basin.param$val[crit] OK, so what if I want to increment from ratio1 to rat...

Dear all, I wonder if anyone can help me with specifying a right model for my analysis. I am a beginner to lme methods. I was unfortunately not able to find a solution to my problem on my own. Data structure: I have sampled monthly 6 basins during two hydrological cycles, and I have taken several (2 to 4) samples (“replicate”) for each basin and month. I’m trying to relate Shannon diversity to some environmental variables, taking away pseudoreplication. Thus, I have tried an lme model with “time” and “basin” as random factors. I have...

...rying to plot multiple figures on one page through a loop, but the code just produce one graph on one page. Can someone show some light on what's wrong? Here is my code: library("quantreg") tcdata<-read.table("mydata.txt",header=TRUE) postscript("myfigure.ps") basins<- paste(c("b1","b2","b3","b4","b5","b6"),sep="") par(mfrow = c(2,3),mar=c(0,4,0,2),oma=c(5,0,3,0)) for (k in 1:6) { data0=subset(tcdata,basin==basins[k]) attach(data0) model=rq(rain~year,tau=seq(0.1,0.9,0.1)) plot(summary(...

...error is: > source('bicarb.R') Error in source("bicarb.R") : bicarb.R:15:1: unexpected symbol 14: 15: hco33 ^ The 'h' is in column 0 so the caret would be column -1, but it's presented as column 1. The file, bicarb.R is: hco31 <- qqmath(~ HCO3 | factor(basin), data = surfchem.cast, main = 'Bicarbonate (Raw)', prepanel = prepanel.qqmathline, panel = function(x, ...) { panel.qqmathline(x, ...) panel.qqmath(x, ...) }) hco32 <- qqmath(~ log10(HCO3 | factor(basin), data = surfchem.cast, main = 'Bicarbonate (Log10)',...

...'ve attempted using glm.nb from the MASS package and the negbin function from the AOD package, but both still cause some models to experience errors and warnings. Examples of my 2 largest models, run from the 2 different functions negBin.glm1 <- glm.nb(Count ~ offset(log(Tow.Area)) + Basin + Bathy + Hypoxia + Period + Depth + Basin*Depth + Bathy*Depth + Hypoxia*Depth + Period*Depth + Basin*Period + Bathy*Period + Hypoxia*Period + Hypoxia:Period:Depth + Bathy:Period:Depth + Basin:Period:Depth, control=glm.control(maxit=1000), method="glm.fit", s...

...I''ve written with comments that include the final error I get. ... initial.dir <- getwd() library(systemfit) library(multcomp) basdata <- read.table("data_into7_test.txt", header=TRUE,sep="") #show basdata basdata #output from above looks like this ... # basin Order lnLcl lnArea #1 SFK 6 3.46322000 5.19766000 #2 SFK 6 3.51767000 5.11809000 #3 SFK 6 3.79962000 5.79118000 #4 SFK 6 2.90242000 4.09768000 #5 SFK 4 1.92408000 2.31617000 # .... #21 NFK 6 3.88342000 5.50927000 #22 NFK 6 3.845...

Dear R users, It is a pleasure for me to announce the availability of the new package kwb.hantush (0.2.1)? on CRAN. Its objective is the calculation of groundwater mounding beneath an (stormwater) infiltration basin by solving the Hantush (1967) equation. For checking the correct implementation of the algorithm the R modelling results were cross-checked against alternative models assessed in Carleton (2010) by using the same model parameterisation. References: [1] Carleton, G.B., 2010, Simulation of groundw...

Dear R users, It is a pleasure for me to announce the availability of the new package kwb.hantush (0.2.1)? on CRAN. Its objective is the calculation of groundwater mounding beneath an (stormwater) infiltration basin by solving the Hantush (1967) equation. For checking the correct implementation of the algorithm the R modelling results were cross-checked against alternative models assessed in Carleton (2010) by using the same model parameterisation. References: [1] Carleton, G.B., 2010, Simulation of groundw...

...ppreciated. negbin(formula, random, data, phi.ini = NULL, warnings = FALSE, na.action = na.omit, fixpar = list(), hessian = TRUE, control = list(maxit = 2000), ...) My largest model using glm.nb looks like this: negBin.glm1 <- glm.nb(Count ~ offset(log(Tow.Area)) + Basin + Bathy + Hypoxia + Period + Depth + Basin*Depth + Bathy*Depth + Hypoxia*Depth + Period*Depth + Basin*Period + Bathy*Period + Hypoxia*Period + Hypoxia:Period:Depth + Bathy:Period:Depth + Basin:Period:Depth, control=glm.control(maxit=1000), method="glm.fit",...

...n xyplot? The desired end result is a single xy plot of two separate (but similar in values and ranges). Full code follows, xyplot code at bottom #########Determine Frequencies ##########coastal_slope #needs the maptools package to read ESRI grid require(maptools) #import the flow slope grid basin.map <- readAsciiGrid("C:/R_PLots/coastal_slp.asc", colname="slope") basin_slope <- (basin.map$slope) #read the slopes into a dataframe freqs<-as.data.frame(table(basin_slope)) #rank the frequencies based on each unique occerence, note, ranks from 1 to n r<-rank(freqs...

Is there any way to run an R script without running R? As an example, suppose I have a tcl/tk interface that asks for a number (in a GUI) and displays its factorial. Is there a way to invoke this script without invoking R? I'm using R 2.4.1 in GNU/Linux Fedora Core 4. Alberto Monteiro

...r, decrypt the payload and then feed the decrypted data to opus. Besides, opus_decoder_get_nb_samples() on 1st packet returns reasonable number 960, but then opus_decode() fails. On Fri, Nov 29, 2013 at 11:39 PM, Jean-Marc Valin <jmvalin at jmvalin.ca>wrote: > On 29/11/13 01:49 PM, Ilya Basin wrote: > > For example, the RTP payload is decrypted successfully, but > > opus_packet_get_nb_samples() or opus_decode() return > > OPUS_INVALID_PACKET. > > Are you sure you didn't forget to strip the RTP header before decoding it? > > Jean-Marc > > -...

On 29/11/13 01:49 PM, Ilya Basin wrote: > For example, the RTP payload is decrypted successfully, but > opus_packet_get_nb_samples() or opus_decode() return > OPUS_INVALID_PACKET. Are you sure you didn't forget to strip the RTP header before decoding it? Jean-Marc

Huhu, I installed MS Freelancer via PlayOnLinux, but I have some graphics lags ingame/the game was redicoulus slow, I also wasn't really able to start a new game since the game hang up (kinda. It just didn't wanted to finish loading) One idea of me would be that it still runs in DX mode, how can I change it to OpenGL? On a side note, I have enough ressources to run the game at windows.

Dear R users, I am sorry to ask you such a pathetic newbie question, but how does one append data at the end of a data frame? I am working with GRASS/R library, but the question is about R. I have a data.frame containing the following variables basinID, distoutlet, drainage_area, slope These variables are stored for all pixels of Grass Raster objects. For each drainage basin (basinID), I'd like to find the maximum of distoutlet. How can I store the pair of info (bv$ID, max(bv$distoutlet)) at each step of the for loop ? I presume something...

...lied Statistics with S and from other sources, but I still can't figure it out. I am trying to estimate species richness (sr) in streams across minnesota. My predictor variables are depth (d), habitat diversity (habdiv), drainage area (da) and an indicator variable representing the watershed/basins that the streams are in (bas: there are 10 watersheds). The variable explaining the greatest amount of variation appears to be da. I have used a log(da) to linearize the relationship, but an asymptotic relationship is more appropriate. I have used nls: B<-c(.007,1,3,2,2,2,2,2,2,2,2,2,.05,.00...

...I have also compiled grass 6.5 from source. I have taken the r.stream* from the addons svn and placed the sources in the grass source tree - in the raster file I run make distclean then navigate to the raster/r.stream* file run make and get this error message: ~/source/grass6_devel/raster/r.stream.basins$ make MODULE_TOPDIR=/usr/local/grass-6.5.svn/lib Makefile:10: /usr/local/grass-6.5.svn/lib/include/Make/Module.make: No such file or directory make: *** No rule to make target `/usr/local/grass-6.5.svn/lib/include/Make/Module.make'. Stop. I know I am missing something simple. Thanks to eve...

Hi. I made a webrtc relay with recording and dumped the SDP requests and RTP packets into files. Then I made a java decoder based on jitsi. Although the files contain all the needed info: encription keys, codec info, timestamps, etc., I could only decode one side in one of 2 conversations. For example, the RTP payload is decrypted successfully, but opus_packet_get_nb_samples() or opus_decode()

On the call where even opus_decoder_get_nb_samples() fails, what's the value of the first two bytes of the buffer you're sending? And what's the length of the buffer? Jean-Marc On 11/29/2013 03:31 PM, Ilya Basin wrote: > Yes, I parse the RTP header, decrypt the payload and then feed the > decrypted data to opus. Besides, opus_decoder_get_nb_samples() on 1st > packet returns reasonable number 960, but then opus_decode() fails. > > > On Fri, Nov 29, 2013 at 11:39 PM, Jean-Marc Valin <j...

...browseable = No hosts allow = . . . hosts deny = ALL guest ok = no Do I need to install some client libraries for this to work? Or did the compilation not work as expected? Any hints or suggestions would be a great help! Thanks, -- Jason, Systems Administrator Basin Telecommunications, Inc. --