similar to: Moving Average

Hi all, I have a very large binary vector, I wish to calculate the number of 1's over sliding windows. this is my very slow function slide<-function(seq,window){ n<-length(seq)-window tot<-c() tot[1]<-sum(seq[1:window]) for (i in 2:n) { tot[i]<- tot[i-1]-seq[i-1]+seq[i] } return(tot) } this works well for for reasonably sized vectors. Does

I am trying to understand what rbinom function does. Here is some sample code. Are both the invocations of bfunc effectively doing the same or I am missing the point? Thanks, Pieter bfunc <- function(n1,p1,sims) { c<-rbinom(sims,n1,p1) c } a=c() b=c() p1=.5 for (i in 1:10000){ a[i]=bfunc(30,p1,1) } b=bfunc(30,p1,10000)

Dear R gurus, As part of finding initial values for a much more complicated fit I want to fit a function of the form y ~ a + bx + cx^d to fairly "noisy" data and have hit some problems. To demonstrate the specific R-related problem, here is an idealised data set, smaller and better fitting than reality: # idealised data set aDF <- data.frame( x= c(1.80, 9.27, 6.48, 2.61, 9.86,

Dear List and Hadley, I would like to have a boxplot with ggplot2 and have the outlier values labelled with their "name" attribute. So I did > library(ggplot2) > dat=data.frame(num=rep(1,20), val=c(runif(18),3,3.5), name=letters[1:20]) > p=ggplot(dat, aes(y=val, x=num))+geom_boxplot(outlier.size=4, outlier.colour="green") >

Hi i have a problem compiling as mantioned on object of this mail: i'm changing dovecot on my debian 6 from version 2.2.32 to 2.2.36 and pigeonhole from 0.4.20 to 0.4.24 this is output of dovecot -n for listing versions installed dovecot -n # 2.2.32 (dfbe293d4): /usr/local/etc/dovecot/dovecot.conf # Pigeonhole version 0.4.20 (7cd71ba) # OS: Linux 2.6.32-5-amd64 x86_64 Debian 6.0.10 i have

I have been very busy for the last few months with school, during which time Jack and others have taken over ogg123. Now I have a week off for Christmas ("Winter Break" for the politically correct), so I have some free time. And what better place to spend it than at my computer? (heh) So anyway I decided to take a look at (what's left of) ogg123. Yuck. Eeew. Horror. The good news:

Hi, Thanks for your response, it nearly worked! But it only wrote one coloumn of data and not the three columns I need. Using fixef(m1) doesnt give the same results as coef(m1) when you are using more than one random effect. I need the coefficients for each individual so I use coef(m1) to get this which results in an object of class lmer.coef, 3 columns by 700 rows. as.data.frame() wont work on

Hi ! I have a small problem. When switching a call (pstn -> sip user), I get the sip phone ringing - ie. everything is OK, but I do not get a ringtone in the handset on the pstn side. Can anyone help me out in how to make * play tones ? My setup: E1 IP pstn ------ Asterisk ------ sip phone Regards, Dave

Hello all, Have some time series data stored in a data.frame, and am plotting it with ggplot2 (which is totally awesome). I have explored the documentation and mailing list archives, and I can't see any way to plot a 'smoother' that is just the K-step moving average. For example, imagine I had a data.frame called 'sleep' with 'date' as the date (from as.Date()) and

hai there, today i tried to implement the prepaid application to my * box. I do the step that mantion in to voip-info. i copy the app_prepaid.c and Make file to my asteris/apps, then i run the make. but it show an error like : gcc -shared -Xlinker -x -o app_prepaid.so app_prepaid.o -lpq /usr/bin/ld: cannot find -lpq collect2: ld returned 1 exit status make[1]: *** [app_prepaid.so] Error 1

I had strange windowing behavior with R 9.1 on Windows XP when executing code from the diagram package. Perhaps it is even normal behavior that I don't understand as it was an unintentional discovery. I don't even know if it is a package issue or a windowing issue. Reproducing it. Execute the code fragment below from vignette(diagram) that uses par() to draw a four panel

Yep On 10.08.2018 11:06, davide marchi wrote: > Thanks Aki for your quick response!! I'm sorry, so i have to compile > and install dovecot version 2.2.36 before compiling and install > pigeonhole 0.4.24? If so sorry for my noob approach > > Il giorno ven 10 ago 2018 alle ore 09:55 Aki Tuomi > <aki.tuomi at dovecot.fi <mailto:aki.tuomi at dovecot.fi>> ha scritto:

Dear All, This is my first post, and probably (and hence apologies that) my question is very silly! I'm having issues with a the mvabund package (http://cran.r-project.org/web/packages/mvabund/index.html), and would be great to get some help! Here is the code (and files are attached): library(mvabund) ##visualizing data florabund <- read.csv("CPL_floristics_abund_v1d.csv",

"mustek" is all but deprecated, so the effort should be to get powermust/megatec working. Since you have the latest development code installed, can you please run "/path/to/megatec -DDDD -a the_name_for_your_ups" (this without starting NUT normally - i.e. without upsd, upsmon, etc) and send me the output? The output should help me figure out why you are getting data stale

Dear list, My name is Juan Fernandez-Tajes and I´m trying to use MALDIquant for analysing MS data from ABI 4700 series, I´ve exported a spot set to T2Dextractor in order to convert in mzXML format. I´ve imported this into R with mqReadMzXml function and have followed the suggested flowchart, however when I´m in the step for peak detections I obtained the following error: peaks <-