| 1 | prms <- function(name, idx) |
| 2 | { |
| 3 | load(name) |
| 4 | d <- nrow(mr[[1]][[1]])-2 |
| 5 | p <- colMeans(do.call(rbind, lapply(mr[[idx]], function(m) m[1,]))) |
| 6 | b <- colMeans(do.call(rbind, lapply(mr[[idx]], function(m) m[2+d,]))) |
| 7 | L <- length(mr[[1]]) |
| 8 | beta <- (1/L) * Reduce("+", lapply(mr[[idx]], function(m) m[2:(d+1),])) |
| 9 | list(p, beta, b, mr_params) |
| 10 | } |
| 11 | |
| 12 | pprms <- function(link) |
| 13 | { |
| 14 | for (n in c("5000", "10000", "100000", "500000", "1000000")) |
| 15 | { |
| 16 | method =1 |
| 17 | #for (method in 1:2) |
| 18 | #{ |
| 19 | toprint <- c() |
| 20 | for (d in c(2,5,10)) |
| 21 | { |
| 22 | name <- paste0("res_", n, "_", d, "_", link, ".RData") |
| 23 | params <- prms(name, method) |
| 24 | toprint <- c(toprint, c( |
| 25 | sum(abs(params[[1]] - params[[4]]$p)), |
| 26 | colSums(abs(params[[2]] - params[[4]]$beta)), |
| 27 | sum(abs(params[[3]] - params[[4]]$b)) |
| 28 | )) |
| 29 | } |
| 30 | print(toprint, digits=2) |
| 31 | #} |
| 32 | } |
| 33 | } |