Skip to content

Instantly share code, notes, and snippets.

View muschellij2's full-sized avatar

John Muschelli muschellij2

480 followers · 15 following
View GitHub Profile
@muschellij2
muschellij2 / install_ffmpeg.sh
Created September 9, 2019 04:06
Install `ffmpeg` on Mac (from https://gist.github.com/clayton/6196167#gistcomment-3017801)
brew tap homebrew-ffmpeg/ffmpeg
brew install homebrew-ffmpeg/ffmpeg/ffmpeg --with-chromaprint /
--with-fdk-aac /
--with-game-music-emu /
--with-libbluray /
--with-libbs2b /
--with-libcaca /
--with-libgsm /
--with-libmodplug /
--with-librsvg /
@muschellij2
muschellij2 / tidy_bootstrap_repeated.R
Created August 15, 2019 17:37
library(dplyr)
if (getRversion() >= "3.6.0") {
set.seed(234)
} else {
# RNGversion("3.5.0")
set.seed(1)
}
# creating example data
n = 20
ids = 1:n
@muschellij2
muschellij2 / fst_stuff.R
Last active June 6, 2019 18:15
Example of using FST data
library(ActivityIndex)
library(fst)
library(lubridate)
library(data.table)
filename = system.file("extdata","sample_GT3X+.csv.gz",
package = "ActivityIndex")
res = ActivityIndex::ReadGT3XPlus(filename)
data = res$Raw
data$DT = paste(data$Date, data$Time)
@muschellij2
muschellij2 / gist:70de5d3b5fc25c0ddcfc0971f50d4c71
Last active May 22, 2019 17:55 — forked from bcaffo/gist:c12e6ca0d65ca8b02e5d3bdcda61a729
## Performs MCMC P-value calculation for 0-1 tables
## from Darwin's Finch data
set.seed(1)
I = 7; J = 5
testDat = matrix(sample(0:1, I * J, replace= TRUE), I, J)
A = matrix(c(
0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1,
@muschellij2
muschellij2 / create_interactive_tweet.R
Created March 4, 2019 04:47
Creating interactive tweet content in #rstats based on https://datatitian.com/how-to-turn-your-ggplot2-visualization-into-an-interactive-tweet/
library(ggplot2)
library(plotly)
library(htmlwidgets)
make_card = function(handle = "@StrictlyStat",
title = "",
description = "",
img = NULL,
html = NULL,
height = 480,
@muschellij2
muschellij2 / script_to_gif.R
Created March 4, 2019 03:04
Take an R Script and make it into a gif.
library(rsciinema)
input = paste0("library(magrittr); # a comment is a gra", backspace(), "eat thing to do",
"\niris %>% \n dplyr::group_by(Species) %>%\n ",
"dplyr::summarise_all(mean) \n # a new line")
data = asciicast( input )
json = rsciinema::json_asciicast(data)
@muschellij2
muschellij2 / Ioannidis_publications_scopus.R
Last active January 23, 2019 22:34
Publications of John Ioannidis from Scopus
# Nature article: https://www.nature.com/articles/d41586-018-06185-8
library(rscopus)
library(lubridate)
library(dplyr)
x = author_data(first_name = "john", last_name = "Ioannidis")
df = x$df
df$date = ymd(df$`prism:coverDate`)
df$year = year(df$date)
@muschellij2
muschellij2 / yg_read.R
Created November 30, 2018 23:17
library(pdftools)
library(purrr)
res = pdftools::pdf_text("YG-Archive-DatingSocialMediaInternal-090818.pdf")
ss = strsplit(res, "\n")
names(ss) = 1:length(ss)
df = map_df(ss, function(x) {
data.frame(txt = x,
row = seq_along(x),
start = grepl("unweighted base", tolower(trimws(x))),
stringsAsFactors = FALSE)
@muschellij2
muschellij2 / rsciinema_example.R
Created November 16, 2018 16:56
rsciinema Example
library(rsciinema)
data = asciicast( "# a comment\niris %>% \n dplyr::group_by(Species) %>%\n dplyr::summarise_all(mean)" )
tfile= tempfile(fileext = ".cast")
write_asciicast(data, tfile)
asciinema( tfile )
@muschellij2
muschellij2 / scrape_rates.R
Last active November 5, 2018 03:46
Getting data from https://webapp.psc.state.md.us/ecm/web/Compare_new1.cfm
# see http://www.opc.state.md.us/RegulatoryActivities/Publications.aspx
# https://webapp.psc.state.md.us/ecm/web/Compare_new1.cfm
library(rvest)
library(dplyr)
# doc = read_html("rates.html")
# tab = html_table(doc, header = TRUE)[[1]]
# tab = tab %>%
# rename(
# name = `Company Name3`,
# price = `Price/KWH`,