scripts/ISPRS/whit_lambda/check_result.R

In kongdd/gee_whittaker: Whittaker smoother

source("test/load_pkgs.R")
source('R/smooth_whit_lambda.R')
# source('test/07_whit/whit_lambda/smooth_whit_lambda.R')

lst <- get_sbatch("result/whit_lambda/whit2_grp1/")
# lst <- get_sbatch("result/whit_lambda/wBisquare_grp1/")

I <- sapply(lst, length) %>% {which(. <= 1 | names(.) == "")} %>%
    sort() #%T>% print
lst_good <- lst[-I]

a <- lst[I]

res <- lambda_optim_FUN(I[1])
res <- par_sbatch(as.numeric(names(I))[51:100],
                  lambda_optim_FUN, wFUN = wSELF,
                  Save = F, return.res = T,
                  outdir = paste0("result/whit_lambda/wSELF", subfix))

sapply(res, length) %>% {which(. == 1)} %>% as.numeric() %>%
    sort() #%T>% print

res <- par_sbatch(sites[134], lambda_optim_FUN, wFUN = wSELF,
                  Save = F, return.res = T,
                  outdir = paste0("result/whit_lambda/wSELF", subfix))