kevinushey’s gists

kevinushey / position_jitterdodge.R

Last active August 29, 2015 13:57

Combine position_jitter and position_dodge so points can align with boxplots in ggplot2

	#' Jitter-dodge points to align them with a boxplot including fill aesthetic
	#'
	#' @family position adjustments
	#' @param width degree of jitter in x direction. Defaults to 40\% of the
	#' resolution of the data.
	#' @param height degree of jitter in y direction. Defaults to 40\% of the
	#' resolution of the data
	#' @export
	#' @examples
	#' dsub <- diamonds[ sample(1:nrow(diamonds), 1000), ]

kevinushey / gsl_ran_exp.cpp

Created March 5, 2014 00:48

gsl vs. R

	// [[Rcpp::depends(RcppGSL)]]

	#include <gsl/gsl_randist.h>

	#include <Rcpp.h>
	using namespace Rcpp;

	// [[Rcpp::export]]
	NumericVector gsl_dexp(NumericVector x) {
	int n = x.size();

kevinushey / magrittr.R

Last active August 29, 2015 13:56

Overloading `|` to act like a pipe operator for data.frames

	`\|` <- function(x, y) {
	if (is.data.frame(x)) {
	return( eval(call("%.%", substitute(x), substitute(y)), envir=parent.frame()) )
	} else {
	return( base::"\|"(x, y) )
	}
	}

	library(dplyr)
	mtcars \|

kevinushey / R_IsNA.cpp

Last active March 3, 2019 01:29

Comparing methods of checking whether a value is NA

	// run me by installing Rcpp and calling sourceCpp on this file

	#include <Rcpp.h>
	using namespace Rcpp;

	// some initial stuff borrowed from R sources

	#ifdef WORDS_BIGENDIAN
	static const int hw = 0;
	static const int lw = 1;

kevinushey / just_for_fun.R

Created January 31, 2014 03:18

Re: Python and R: Is Python really faster than R?

	case2 <- function(n = 10, mu = 3, sigma = sqrt(5), p = 0.025, rep = 100){

	xbar <- vapply(1:rep, FUN.VALUE=numeric(1), function(i) {
	norm <- rnorm(mean = mu, sd = sigma, n = n)
	return(mean(norm))
	})

	low <- xbar - qnorm(1-p) * (sigma / sqrt(n))
	up <- xbar + qnorm(1-p) * (sigma / sqrt(n))
	rem <- mu > low & mu < up

kevinushey / git_status.R

Created January 29, 2014 18:49

A clever use of R's print mechanism for getting the git status for the current directory

	git_status <- ''
	class(git_status) <- "__git_status__"
	print.__git_status__ <- function(x, ...) system("git status")

	git_status ## prints 'git status' for the current dir

kevinushey / enumerate.c

Last active January 3, 2016 06:59

A simple extension of `lapply` to allow passing an index as second argument of the function

	// This code copied and modified from the R sources (src/main/apply.c) for do_lapply.
	// Essentially, we remove some of the op / arity checking, and then
	// construct a call to FUN that also passes the index.

	#define USE_RINTERNALS

	#include <R.h>
	#include <Rinternals.h>

	// [[export]]

kevinushey / NA_punning.cpp

Last active January 1, 2016 16:09

Faster NA checking with punning

	#include <Rcpp.h>
	using namespace Rcpp;

	union DoublePunner {
	double d;
	uint64_t p;
	};

	static DoublePunner NA_PUNNED = { NA_REAL };

kevinushey / Rcpp_ListOf.cpp

Last active March 22, 2016 20:16

Rcpp::ListOf<T> -- An example implementation of a typed List container.

kevinushey / server.R

Created September 17, 2013 21:11

shinyGridster test

	library(shiny)

	shinyServer( function(input, output, session) {

	})

Kevin Ushey kevinushey