porotypes/rateweight.R

In RGLab/flowReMix: Generalized Linear Modeling of FCM Data With Clustered Coefficients

iterations <- 1000
weights <- matrix(0, nrow = iterations, ncol = iterations)
weights[1, 1] <- 1
rate <- 0.55
for(i in 2:iterations) {
  weights[i, 1:(i - 1)] <- weights[i - 1, 1:(i-1)]  - weights[i - 1, 1:(i-1)] * 1/i^rate
  weights[i, i] <- 1/i^rate
}

sums <- t(apply(weights, 1, function(x) cumsum(x[iterations:1])))
discardAt <- 0.9
nsamps <- apply(sums, 1, function(x) sum(x != 0 & x < 0.9))



a <- list(c(1, 1, 1), c(2, 2))
attr(a[[1]], "iter") <- 1
attr(a[[2]], "iter") <- 2
weightList <- list()
for(i in 1:length(a)) {
  weightList[[i]] <- rep(attr(a[[i]], "iter"), length(a[[i]]))
}