search for: indec

...e?")){ from <- length(startVec) shortVec<-numeric() tics <- from*to for(j in 1:to){ interval <- ((j/to)*tics - (1/to)*tics + 1):((j/to)*tics) benchmarks <- interval/to #FIRST RUN ASSUMES FINAL BENCHMARK/TO IS AN INTEGER... positions <- which(round(benchmarks) == benchmarks) indeces <- benchmarks[positions] fracs <- numeric() #SINCE MUCH OF THE TIME THIS WILL NOT BE THE CASE, THIS SCRIPT DEALS WITH THE REMAINDER... for(i in 1:length(positions)){ if(i == 1) fracs[i] <- positions[i]/length(benchmarks) else{ fracs[i] <- (positions[i] - sum(positions[1:(i-1)]))/leng...

...9;m trying to extract a histogram over the results from a bootstrap. However I keep receiving the error message "Error in hist.default(boot.lrtest$ll, breaks = "scott") : 'x' must be numeric". The bootstrap I'm running looks like: > boot.test <- function(data, indeces, maxit=20) { + y1 <- fit1+e1[indeces] + mod1 <- glm(y1 ~ X1-1, maxit=maxit) + y2 <- fit2+e2[indeces] + mod2 <- glm(y2~1, maxit=maxit) + ll <- 2*(logLik(mod1)-logLik(mod2)) + ll + } > boot.lrtest <- boot(data=M1, statistic=boot.test, R=2000, maxit=100); > hist(boot.lrtest$...

Dear Contributors, I have an easy question for you which is puzzling me instead. I am running loops similar to the following: for (i in c(100,1000,10000)){ print((mean(i))) #var<-var(rnorm(i,0,1)) } This is what I obtain: [1] 100 [1] 1000 [1] 10000 In this case I ask the software to print out the result, but I would like to store it in an object. I have tried a second loop, because if I

Dear Barry, Thanks for your stimulating reply. As you see I have send a CC of this reply to the R mailing list. The R mailing list will be the best place to send your queries. You can subscribe here: http://www.r-project.org/mail.html I would think there is an interest in the provision of landscape indices analysis methods in R. There are several packages for landscape indices analysis

...it would be nice to know the total number of items in a series to be shure to have all parts together (at least I note the total number of disks on the bottom of a CDR). So what about the conceptually same: -> Add a TRACKTOTAL, DISKTOTAL and INDEXMAX (not INDEXTOTAL for disambiguation, since indeces count from zero in the red book standard) 3) Or would you consider it wiser reusing (abusing) existing standard vorbis fields and note the total number of disks/tracks/indeces inside the very same field like TRACKNUMBER=3/15 or would this break more things since apps expect a number in there?...

...ermite llamar al verdadero Python. Funciona perfectamente en mi máquina, pero necesito ver qué problemas de instalación y uso aparecen en otras plataformas y con otras configuraciones. De momento, debería correr sobre UNIX o Linux, siempre y cuando tengan instalado Python, claro. Me sorprendería lo indecible que funcionase también sobre Windows: sería toda una casualidad. Si alguien tiene unos minutos y quiere echarme una mano, lo invito a: -- Instalar el paquete desde R-Forge mediante el comando install.packages("rPython", repos="http://R-Forge.R-project.org") -- Cargar el pa...

...ngth(data) sample <- matrix(0,numsim,1) for (i in 1:numsim) {indices <- runif(pool, min = 1, max = n) sample[i,] <- data[indices] } the idea behind this for loop is dat it first creates randomly data by runif, and then it should somehow, by using indeces, store this 3 times in a matrix 'sample'. Somehow, I am stuck here and I think that the sample[,i] part is just very wrong. I am however not for a very long time an R user, so I cannot find the solution to this problem. Can someone give me a hint? Thanks very much! Bye, sarah -- Vi...

Hi all, >From Monty's letter I know that ResVQTrain runs on vqd files that contain numbers separated by commas. I know that the .vqd files of this format generated by Vorbis compiled with RES_TRAIN contains partwords, i.e. codebook indeces assigned to partitions. I think that ResVQTrain should run on real data, not codebook indeces. Did I misunderstand something? Thanks Ervin <p><p>--- >8 ---- List archives: http://www.xiph.org/archives/ Ogg project homepage: http://www.xiph.org/ogg/ To unsubscribe from this list,...

...6.49 X428 543.57 #data set2 def #output: X401 X402 X404 X405 X406 X407 X409 X411 X412 X413 X414 X416 528.46 524.15 527.18 526.04 533.71 537.79 536.80 532.38 517.14 529.32 523.29 539.58 X417 X421 X422 X429 535.38 532.68 515.28 523.10 Both are numeric values and have indeces above each of the numbers which referring to X401 to X429 indices. I would like to combine both by sorting X401, X402,X403 and so on. Could somebody please help me before I waste my time using excel to do this. Thanks! [[alternative HTML version deleted]]

...val' for the large enough selection of rows. (I didn't check, how many are needed, but it seems that 10000 rows will suffice.) When trying to 'Get R val', there will be an error. (If there are too much rows, RExcel will just hang.) Some of the VBA routines seemed to have integer indeces. Could this be the source of problems? I would need this on next week (by wednesday 17.11.2004), but I'm afraid that it is impossible to get this fixed in this short time. (I don't know how to do this myself and don't have enough time to learn.) In any case, for future use, it woul...

...Look at the following R idiom: a <- array(1:30,c(3,5,2)) M <- (matrix(1:15,c(3,5)) %% 4) < 2 a[M,] <- 0 Now, I think that "a[M,]" has an unambiguous meaning (to a human). However, the last line doesn't work as desired, but I expected it to...and it recently took me an indecent amount of time to debug an analogous case. Just to be explicit, I would expect a[M,] to extract a[i,j,] where M[i,j] is TRUE. (Extract.Rd is perfectly clear here, and R is behaving as documented). The best I could cobble together was the following: ind <- which(M,arr.ind=TRUE) n <- 3...

...m: - nind[0][0..nb_subvect], or - nind[1][0..nb_subvect], or - nind[2][0..nb_subvect], or - ... Or do I have to do a little more hacking in the search function to get the N best _combinations_ of CB IDs? Wouldn't it be easier and just as effective (for N==2) to do the following: Get nb_subvect indeces to write into the stream. Write all but the last of these into the stream. If nind[0][nb_subvect-1] and nind[1][nb_subvect-1] have different LSBs and if the difference between ndist[0] and ndist[1] is small enough (smaller than the normal variation of ndist[0], which can be measured), we can deci...

...t 13:23, C Bergström <cbergstrom at pathscale.com> wrote: > Is the list PG, PG-13, R or at what level do "we" adults all consider > "ok". Even on broadcast tv (in the US) you'll hear some profanity. > (context) > https://www.fcc.gov/consumers/guides/obscene-indecent-and-profane-broadcasts Excellent context! > Some people have pointed out that they don't like the R-rated style of > the LKML. Profanity and no holds barred just isn't for some people. I > can respect that, but personally I find it more funny and raw/honest. I don't care...

...n, or copying of this message is strictly prohibited. If you have received this communication in error, please notify us immediately by e-mail, and delete the original message. (ii) Madan.Kundu at ranbaxy.com confirms that Ranbaxy shall not be responsible if this email message is used for any indecent, unsolicited or illegal purposes, which are in violation of any existing laws and the same shall solely be the responsibility of Madan.Kundu at ranbaxy.com and that Ranbaxy shall at all times be indemnified of any civil and/ or criminal liabilities or consequences there.

On 5 May 2016 at 12:42, Renato Golin <renato.golin at linaro.org> wrote: > I'm against the ownership of firearms, and go at great lengths and > poorly choosing words in a discussion, which some could consider rude, > with person X about it. I know person X for decades and have earned > the right to offend him/her personally as they know I don't mean it > (could be a

...rot++; double *psi_ln_psiq_c; psi_ln_psiq_c = REAL(psi_ln_psiq); int not_converged; int maxIter = 10000; int iter; int start_c; /* loop indeces */ int i; int j; int k; /* loop over observational units to find choice probabilities for i=1,...,N */ for (i=0;i<N_c;i++) {...

Hi, I'm having trouble understanding how to construct a random number generator for a parametric bootstrap. My aim is to bootstrap a Likelihood Ratio statistic (under the null) for a linear model. The function at this point is given by boot.test.n01 <- function(data, indeces, maxit=20) { y1 <- fit1+se(e2)*rnorm(314) mod1 <- glm(y1 ~ X1-1, maxit=maxit) y2 <- fit2+se(e2)*rnorm(314) mod2 <- glm(y2~1, maxit=maxit) t <- 2*(logLik(mod1)-logLik(mod2)) t } boot.lrtest.n01 <- boot(data=M1, statistic=boot.test.n01, R=3999, maxit=100, sim="parametric&quo...

Hello, I gave a list of 2 x 2 matrices called matlist. I have about 5000 2 x 2 matrices. I would like to count how many of each 2 x 2 unique matrix I have. So I am thinking that I need a list of the unique 2 x 2 matrices and their counts. Can anyone help. -- Thanks, Jim. [[alternative HTML version deleted]]

Hi all, There's a myriad of options these days and I haven't been keeping up to date with what's respectable any longer. I essentially need a provider that will provide me with one DID to start and let me trunk over SIP or IAX2. I'd like to obviously trunk with Asterisk on my end and have full control over the dial plan. This way I can branch out my DID into extensions and have

Hello, I want to reduce the bandwidth taken by an IAX trunk when used with a small number of voice channels. When only one call is passed through an IAX trunk, the IP overhead is indecent. I would like to increase the IAX voice packet emission interval (20ms - i'm using speex) to something much larger, in order for the packets to transport as much data as possible. I've been trying to fuddle with the "trunkfreq" parameter and push it to 80ms, but that did...