JD Long CerebralMastication

☠️

brilliantly executing terrible plans

Data. And Pain.

CerebralMastication / pythontest.py

Created January 6, 2010 17:28

	drawBoard <- function() {
	require(plotrix) # for the draw circle function only

	plot(-260:260, -260:260, type="n", xlab="", ylab="", asp = 1, main="Dart Board")
	draw.circle(0, 0, 12.7/2, border="purple", lty=1, lwd=1) # bull
	draw.circle(0, 0, 31.8/2, border="purple", lty=1, lwd=1) # outer bull
	draw.circle(0, 0, 107, border="purple", lty=1, lwd=1)
	draw.circle(0, 0, 99, border="purple", lty=1, lwd=1)
	draw.circle(0, 0, 162, border="purple", lty=1, lwd=1)
	draw.circle(0, 0, 170, border="purple", lty=1, lwd=1)

CerebralMastication / darts.R

Created January 15, 2010 16:01

	#draw.circle is from the plotrix library
	require(plyr)

	draw.circle <- function (x, y, radius, nv = 100, border = NULL, col = NA, lty = 1,
	lwd = 1)
	{
	xylim <- par("usr")
	plotdim <- par("pin")
	ymult <- (xylim[4] - xylim[3])/(xylim[2] - xylim[1]) * plotdim[1]/plotdim[2]
	angle.inc <- 2 * pi/nv

CerebralMastication / 3d_lattice.R

Created January 19, 2010 15:32

CerebralMastication / color in ggplot.R

Created January 20, 2010 23:16

	#fast example of hex color + transparency in ggplot2

	df <- data.frame(x=rnorm(100))
	p <- ggplot(df, aes(x=x))
	p + geom_density(fill = alpha("#335785", .6))

CerebralMastication / rmvdc.new.R

Created January 29, 2010 22:33

	require(copula)

	rmvdc.new <- function (mvdc, n)
	{
	dim <- mvdc@copula@dimension
	u <- rcopula(mvdc@copula, n)
	x <- u
	for (i in 1:dim) {
	if (mvdc@margins[i]=="Johnson") {
	qdf.expr <- copula:::asCall(copula:::P0("q", mvdc@margins[i]), list(mvdc@paramMargins[[i]])) } else {

CerebralMastication / ulimit.txt

Created February 11, 2010 23:36

	root@domU-12-31-39-0B-8D-81:~# ulimit -a
	core file size (blocks, -c) 0
	data seg size (kbytes, -d) unlimited
	scheduling priority (-e) 20
	file size (blocks, -f) unlimited
	pending signals (-i) 16382
	max locked memory (kbytes, -l) 64
	max memory size (kbytes, -m) unlimited
	open files (-n) 1024
	pipe size (512 bytes, -p) 8

CerebralMastication / AMZNloadPackage.R

Created February 14, 2010 21:04

	libPath <- paste("~/R", pid, "/", sep='')

	options(repos=c(CRAN="http://streaming.stat.iastate.edu/CRAN/"))
	dir.create(libPath)
	install.packages("Hmisc", lib=libPath)
	require(Hmisc, lib=libPath)

CerebralMastication / deprecatedShit.R

Created February 23, 2010 20:28

CerebralMastication / random_draws_peterflom.R

Created February 24, 2010 22:54

	makeDraws <- function(numberOfDraws) {
	i <- 1
	myOutput <- rep(NA, numberOfDraws) #preallocate the vector so it all goes faster
	while (i <= numberOfDraws) {
	a <- runif(7, 0, 100) #random draws from a uniform dist 0,100
	b <- runif(7, 0, 100)
	absDiff <- abs(a-b) #calculate the absolute diff
	myOutput[i] <- sum(absDiff) #sum up the diffs and put them in the output vector
	i <- i +1
	}

CerebralMastication / test_standardDev.R

Created April 9, 2010 22:12

	################################################################################################
	# AUTHOR: JD Long
	# PURPOSE: put confidence bands around a sample standard deviation
	################################################################################################


	#based on example code from here:
	#http://www.boost.org/doc/libs/1_36_0/libs/math/doc/sf_and_dist/html/math_toolkit/dist/stat_tut/weg/cs_eg/chi_sq_intervals.html

	#thanks to Vince Buffallo for sending me the link and pointing me in the right direction.