R/EMMax.sd.const.R

EMMax.sd.const <-
function(y,z,mu) {
  sd <- sqrt( weighted.mean((rep(y,length(mu)) - rep(mu, each=length(y)))^2, as.vector(z)))
  sigma <- rep(sd,length(mu))
  sigma[sigma < 0.01] <- 0.01
  sigma
}

Try the fitTetra package in your browser

Any scripts or data that you put into this service are public.

fitTetra documentation built on May 2, 2019, 3:52 p.m.