{{ message }}

Instantly share code, notes, and snippets.

# Hao Ye ha0ye

Last active Aug 31, 2021
generate random sequences without consecutive repeats
View generate_sequences.R
 f <- function(x) { n <- sum(x) v <- numeric(n) prev <- NULL for (i in 1:n) { idx <- x > floor(sum(x)/2) # how many counts are > n/2 if (sum(idx) >= 2) # error if more than 2 (impossible?) {
Last active Jun 26, 2020
EDM simplex calculations test (E = 1)
View main.R
 library(rEDM) # generate time series set.seed(42) ts_length <- 50 x <- rnorm(ts_length) # pre-compute distance matrix x_dists <- as.matrix(dist(x)) x_dists <- x_dists[, -ts_length] # last point can't be a nearest neighbor
Created Mar 28, 2020
pairwise overlap calculations
View gist:197de5797df1d4e2b14a8bff4fb99c95
 calc3 <- function(sets) { sets <- check_sets(sets) set_lengths <- vapply(sets, length, 0) set_order <- order(set_lengths) sets <- sets[set_order] set_lengths <- set_lengths[set_order] n_sets <- length(sets) set_names <- names(sets)
Created Mar 12, 2020
Forward forecasts using rEDM
View gist:d62b33b1d735b09610852a7113f6bda3
 library(rEDM) #### Example 1: a single forward prediction # Assume we have already done the analysis, finding that the best embedding # dimension is 2. Our objective is to make multiple 1-step ahead forecasts # using simplex and E = 2. dat <- data.frame(yr = as.numeric(time(sunspot.year)), sunspot_count = as.numeric(sunspot.year)) E <- 3
Created Feb 18, 2020
read in param output from text file
View script.R
 full_text <- readLines(file, warn = FALSE) start_par_line <- grep("^\\\$par", full_text) start_objective_line <- grep("^\\\$objective", full_text) param_lines <- seq(from = start_par_line, to = start_objective_line - 3, by = 2) params <- purrr::map_dfc(param_lines, function(line_idx) { readr::read_delim(file, delim = " ",
Last active Jul 8, 2019
View dplyr-timing-test.R
 library(dplyr) library(microbenchmark) set.seed(42) n <- 10000 df <- data.frame(a = rep(seq(n), 2), x = rnorm(2 * n)) group_summarize_mutate <- function(df) {
Created Jun 29, 2019
View change legend labels in "likert" package
 library(likert) data(pisaitems) ##### Item 24: Reading Attitudes items24 <- pisaitems[,substr(names(pisaitems), 1,5) == 'ST24Q'] head(items24); ncol(items24) names(items24) <- c( ST24Q01="I read only if I have to.", ST24Q02="Reading is one of my favorite hobbies.",
Created Jun 21, 2019
parallelization within drake targets
View drake-parallel.R
 library(drake) library(furrr) future::plan(future::multiprocess) tictoc::tic() params = data.frame(t = c(1, 2, 3, 1, 2, 3)) out = purrr::pmap(params, Sys.sleep) tictoc::toc()
Created May 3, 2019
View lavaan to visNetwork
 library(lavaan) library(visNetwork) model <- ' # measurement model ind60 =~ x1 + x2 + x3 dem60 =~ y1 + y2 + y3 + y4 dem65 =~ y5 + y6 + y7 + y8 # regressions dem60 ~ ind60
Created May 11, 2018
Travis Erroring on R package install
View log.txt
 travis_fold:start:worker_info [0K[33;1mWorker information[0m hostname: 12607e0a-d790-494c-abc0-33bec19d51ab@1.production-1-worker-org-a-4-gce version: v3.6.0 https://github.com/travis-ci/worker/tree/170b2a0bb43234479fd1911ba9e4dbcc36dadfad instance: travis-job-a2fd8785-6b65-4c56-b87e-cbd4d1016e0f travis-ci-garnet-trusty-1512502259-986baf0 (via amqp) startup: 26.186769693s travis_fold:end:worker_info [0Ktravis_fold:start:system_info [0K[33;1mBuild system information[0m Build language: ruby