search for: andrae

Dear All, would you please provide your thoughts on the following: let us say I have: a <-c(1,5,8,15,32,69) b <-c(8.5,33) and I would like to extract from "a" the two values that are closest to the values in "b", where the length of this vectors may change but b will allways be shorter than "a". So at the end based on this example I should have the result

...such, but this seems like overkill. I am wondering if there is a simple solution that would allow me to accomplish this. The code I am using is outlined below. Data <- c(1:10) #Data value <- seq(from = 0, to = 100 , by = 0.01) #Sequence Data - min(Data) + k # Equation Thanks, Jake Jake Andrae PhD Candidate Geology & Geophysics ? Sprigg Geobiology Centre Department of Earth Science School of Physical Sciences The University of Adelaide, AUSTRALIA 5005 Phone: 0407701565 Email: jake.andrae at adelaide.edu.au [[alternative HTML version deleted]]

...omplexities that you have omitted. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Mar 14, 2018 at 9:14 PM, Jake William Andrae < jake.andrae at adelaide.edu.au> wrote: > Hi All, > > I have a vector of data on which I am operating. The equation with which I > want to operate on the vector has a value k. I want to run the equation and > output a new vector, each time replacing k with each value from the...

So, like a lot of people, I am using --link-dest to do backups. My backup target is ext4 so with a hard link limit of 64K. I do end up with trying to create too many links at some point though and get the following sequence of events: rsync: link "/.snapshots/hourly.0/linux/usr/src/openwrt-r18617/build_dir/host/mpfr-2.4.1/tests/.deps/tui_sub.Po" =>

Dear All, I have a bunch of data points as follows: x  100 y  200 z  300 ... where 100, 200, 300 are the values. I would like to know the distribution of my data? how can I fit my data into a distribution? Thanks a lot, Andra [[alternative HTML version deleted]]

Dear All,   this is probably an easy one but I can not get a handle on it:   x <-c(1,2,3,4,5) y <-c(6,7,8,9,10) z <-15 w <-ifelse(z>14,x,y)   this will give me a value of 1 for w. What I would like to get is the whole string of x, so that w would become a numeric object of 5 characters exactly the same as x.   Apreciate the help,   Sincerely,   Andras [[alternative HTML version

Dear All,   I have a matrix with 33 columns and 5000 rows. I would like to find 2 specific columns in the set: the one that holds the highest values and the one that holds the lowest values. In this case the column's mean would be apropriate to use to try to find those specific columns because each columns mean is different and they all change together based on the same "change of rate

hi! I have a problem with rsync (current CVS version), running on a win2k machine (rsync compiled using Cygwin). c:\rsync>rsync --version rsync version 2.4.6dev protocol version 24 Copyright (C) 1996-2001 by Andrew Tridgell, Paul Mackerras and others Capabilities: 64-bit files, socketpairs, hard links, symlinks It can DOWNload >2gb files from my linux box (with rsync reporting negative

Hi All, When modeling with glm and family = binomial (link = logit) and response values of 0 and 1, I get the predicted probabilities of assigning to my class one, then I would like to compare it with my vector y which does have the original labels. How should I change the probabilities into values of zero and 1 and then compare it with my vector y to find out about the accuracy of my

Dear All,   I am working with a Sweave document to be converted into PDF using Rstudio. It seems to me that my R code will also show up after conversion, which I would like not to happen. Is there a way to specify a command that I could place on the beggining and at the end of the R code that would allow R to recognize the area where my R code is and would ignor it when the data is converted?  

Dear All,   I often have to work with certain models in which I try to "reproduce" a distribution the best I can with very little known information avaible. Is there a package or function in R that could best reproduce a probability distribution using only the mean, median and SD values availble without knowing the actual distribution type to begin with and/or the covariance matrix (for

Dear All   I have the following data (somewhat simplyfied):   TINF <-1 a <-c(500,750,1000,1250,1500,1750,2000) b <-c(8,12,18,24,36,48,60,72,96)   following function:   infcprodessa <-function (D, tin, tau, ts)   (D * (1 - exp(-0.048 * tin))/(tin * (0.048*79) * (1 - exp(-0.048 * tau)))) * exp(-0.048 * (ts - tin)) z <-sapply(1:1, function(n) infcprodessa(1000,TINF,12,12-TINF))   is

Dear All,   I was wondering if you could help me with the following: I have the code:   tin <-0.5 tau <-24 output0 <-10 TIMELOW <-tin TIMEHIGH <-1*tau TIME1 <-c(seq(TIMELOW,TIMEHIGH, by = sign(TIMEHIGH-TIMELOW)*(tau-tin)/3))   then I would like to calculate:   cp1 <-output0*exp(-0.3*TIME1[1]) cp2 <-output0*exp(-0.3*TIME1[2]) cp3 <-output0*exp(-0.3*TIME1[3]) cp4

Hi All, I am trying to fit my data with glm model, my data is a matrix of size n*100. So, I have n rows and 100 columns and my vector y is of size n which contains the labels (0 or 1) My question is: instead of manually typing the model as glm.fit = glm(y~ x[,1]+x[,2]+...+x[,100], family=binomial()) I have a for loop as follows that concatenates the x variables as follows: final_str=NULL for

Peter, thanks, very nice, this will work for me... could you also help with setting up the code to run the on liner "approx(sort(x), seq(0,1,,length(x)), q)$y" on the rows of a data frame using my example above? So if I cbind z and res,? df<-cbind(z,res) the "x" in your one liner would be the first 4 column values of each row and "q" is the last (5fth) column

Dear All,   as a follow up to my previous e-mail (I think I am getting closer...):   I am trying to apply the trapezoidal functions to a matric column by column. I have the following code:   a <-matrix(c(1:100),ncol=10) b <-matrix(c(2,4,6,8,10,12,14,16,18,20))   apply(a,2,function(b,a) sum(diff(b)*(a[-1]+a[-length(a)]))/2)   for some reason i get an error message: Error in FUN(newX[[, i],

Dear All,   sorry, got stuck again on the following: let us say we have:   a <-c(1:5) b <-c(6:10) d <-cbind(a,b)     from d I would like to remove total number of rows based on the length of f. So if:   f <-c(1)   my result is working great with the following solution:   d[-length(f),]   so I get: a b [1,] 2 7 [2,] 3 8 [3,] 4 9 [4,] 5 10 but if I do: f <-c(1,2) then I get: 

Hi, I''ve downloaded the latest stabil Hieraki. I have ruby 1.8.2 (2004-12-25) [i686-linux], the latest rails (and gems, too), Apache/1.3.33 (Unix) and mod_fastcgi/mod_fastcgi-SNAP-0404142202. I''m successfully running Typo (from SVN) as a FastCGI application. Hieraki works well with "script/server", and works well as a .cgi (just very slow). My problem is, that

Dear All,   another quick question, this one is on skipping part of my code, so let us say:   a <-5 b <-2 e <-0   d <-a+b f <-a-b   what I would like to do is to have R NOT to calculate the value for d in case the value of e equals to zero (essentially skip that "chunk"), but instead move on to calculate te value for f. In the code I am working with the value of e changes,

I'll try to make some shorter samples later, but for now here are the ones I have tried with: http://www.surfnonstop.com/~bandit/speex/1.2beta1_AEC_garble/ The original recordings are in mic.raw and spk.raw. Jean-Marc Valin wrote: > You may have triggered an instability problem. Can you upload your files > somewhere so I can have a look at them? > > Jean-Marc > > Andras