search for: lcss

...+0.5*wt[i+1])/n} r[1]<-0.5*wt[1]/n wmean<-sum(y*wt)/sum(wt) ci<-2*sum(wt*(y-wmean)*(r-0.5))/sum(wt*y) devci<-abs(ci)*(1-abs(ci)) if(continous){ lss<-cumsum(wt*y)[-n] rss<-sum(wt*y)-lss lw<-cumsum(wt)[-n] rw<-sum(wt)-lw lm<-lss/lw rm<-rss/rw lcss<-cumsum(wt[1:length(lm)]*(y[1:length(lm)]-lm)*(r[1:length(lm)]-0.5)) rcss<-sum(wt*(y-wmean)*(r-0.5))-lcss lci<-2*lcss/lss #concentration index for left side rci<-2*rcss/rss #concentration index for right side devlci<-abs(lci)*(1-abs(lci)) #deviance for left...