Nothing
increasing <-
function(data,x,n.rep){
k <- length(n.rep)
data <- as.numeric(data)
n <- sum(n.rep)
s <- x
w <- n.rep
a <- s*w
is <- numeric(k)
uu <- 0
while(uu < k) {
v <- uu + 1
b <- cumsum(a[v:k] )
d <- cumsum(w[v:k])
b <- b/d
u <- min(b)
m <- match(u, rev(b))
m<- length(b)-m+v
is[v:m] <- u
uu <- m
}
sig <- 0;
len1 <- 1;
len2 <- 0;
for(i in 1:k){
muu <- is[i]
len2 <- len2+n.rep[i]
x1 <- data[len1:len2]
len1 <- len1+n.rep[i]
sig <- sig+sum((x1-muu)^2)
}
sig <- sig/n;
like <- -0.5*n*log(2*pi)-0.5*n*log(sig)-0.5*n
list( logelr=like, mu=is)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.