agstudy’s gists

agstudy / gist:b93ad037715f80848f35

Created June 15, 2015 13:33

	---
	title: "Data/Structure Validation"
	author: "agstudy"

	---

	This post is an answer to [SO question](http://stackoverflow.com/questions/30844363/data-structure-validation-for-r#comment49735091_30844363) about creating a well typed data structure in R.

	I think, that in R the only way to define a typed data structure is to `S4 class`. I should not that even S4 classes are not strongly typed since you can define your slot as `list`.

agstudy / gist:7908344

Created December 11, 2013 10:42

[Add yz- and xz-grid to scatterplot3d](http://stackoverflow.com/questions/20448539/add-yz-and-xz-grid-to-scatterplot3d)

	sactter.grid <- function (x, y = NULL, z = NULL, color = par("col"), pch = NULL,
	main = NULL, sub = NULL, xlim = NULL, ylim = NULL, zlim = NULL,
	xlab = NULL, ylab = NULL, zlab = NULL, scale.y = 1, angle = 40,
	axis = TRUE, tick.marks = TRUE, label.tick.marks = TRUE,
	x.ticklabs = NULL, y.ticklabs = NULL, z.ticklabs = NULL,
	y.margin.add = 0, grid = TRUE, box = TRUE, lab = par("lab"),
	lab.z = mean(lab[1:2]), type = "p", highlight.3d = FALSE,
	mar = c(5, 3, 4, 3) + 0.1, bg = par("bg"), col.axis = par("col.axis"),
	col.grid = "grey", col.lab = par("col.lab"), cex.symbols = par("cex"),
	cex.axis = 0.8 * par("cex.axis"), cex.lab = par("cex.lab"),

agstudy / gist:7676651

Created November 27, 2013 14:34

	<?xml version="1.0" encoding="utf-8" ?>
	<root>
	<div class="BVRRReviewTitleContainer">
	<span class="BVRRLabel BVRRReviewTitlePrefix"></span>
	<h2>
	<span itemprop="name" class="BVRRValue BVRRReviewTitle">Perfect size for the kids and durable</span>
	</h2>
	<span class="BVRRLabel BVRRReviewTitleSuffix">, </span>
	</div>
	<div class="BVRRReviewDateContainer">

agstudy / gist:6845022

Last active December 24, 2015 18:39

DFS to get the number of linked ones (islands) in a matrix of 1 and 0

	/* package whatever; // don't place package name! */

	import java.util.*;
	import java.lang.*;
	import java.awt.Point;
	import java.io.*;

	/* Name of the class has to be "Main" only if the class is public. */
	public class Maze
	{

agstudy / gist:5473600

Last active December 16, 2015 17:50

recommender job. This is a markdown. Create a Rstudio to convert it markdown it to html file or something like : library(knitr); knit2html(input='report.Rmd',output='report.html')

	Item-based Collaborative Filtering
	========================================================

	The task is to use 3003377s-offsides.tsv file of drugs and side effects (url below) and use the R recommenderlab library or similar to create code that will let us do Item-based Collaborative Filtering, which will mean :
	* we can identify drugs with similar side effects
	* possibly predict unknown side effects by finding the average relative risk of side effects of the K-nearest neighbors using some similarity metric.

	Basically this means making a matrix where cells are relative risk (“rr” or “log2rr” column) and columns are side effects and rows are drugs (use either the “stitch_id” or “drug” column). Then we can randomly omit values from cells and attempt to answer which distance metric or transformation of relative risk allows us to predict more accurately the omitted values.

	Reading data and levelplot

agstudy / gist:5377621

Created April 13, 2013 08:40

	## read data
	ll <- readLines(textConnection('28/03/13 FACTURE CARTE DU 240313 PHOTOBOX SARTROUVILLE CARTE -30
	25/03/13 2EME TIERS SUR FACTURE DE 180 EUR DU 240113 CAROLL INTL PARI -60
	25/03/13 RETRAIT DAB 24/03/13 11H15 166936 BNP PARIBAS PARIS -200.50
	18/03/13 PRELEVEMENT TRESOR PUBLIC 94 IMPOT NUM 005002 ECH 18 -500'))
	## convert data to a data.frame
	dat <- do.call(rbind,strsplit(gsub('[ ]{2,}','\|',ll),'\\\|'))
	## name columns
	colnames(dat) <- c('date','category','desc','amount')
	dat <- as.data.frame(dat)

agstudy / gist:5024781

Last active April 1, 2019 12:58

Multi calendar heat. Giving , dates , and 2 values per date.


	##########################################
	extra.calendarHeat <- function(dates,
	values,
	ncolors=99,
	color="b2w",
	pch.symbol = 15:20,
	cex.symbol =2,
	col.symbol ="#00000044",
	pvalues=values,

agstudy / gist:5013751

Last active December 14, 2015 02:29

transform data


	transformdata <- function(dates, values, date.form = "%Y-%m-%d", ...) {
	require(lattice)
	require(grid)
	require(chron)
	if (class(dates) == "character" \| class(dates) == "factor" ) {
	dates <- strptime(dates, date.form)
	}
	caldat <- data.frame(value = values, dates = dates)
	min.date <- as.Date(paste(format(min(dates), "%Y"),

agstudy / gist:5013606

Last active June 27, 2019 04:56

add division to the calendar.


	calendar.division <- function(...)
	{
	xyetc <- list(...)
	subs <- dat[xyetc$subscripts,]
	dates.fsubs <- dat[dat$yr == unique(subs$yr),]
	y.start <- dates.fsubs$dotw[1]
	y.end <- dates.fsubs$dotw[nrow(dates.fsubs)]
	dates.len <- nrow(dates.fsubs)
	adj.start <- dates.fsubs$woty[1]

agstudy agstudy