{{ message }}

Instantly share code, notes, and snippets.

# Michael Kuhn mkuhn

Created Jun 30, 2021
Hollow circles for size legend
View guides.md
```library(tidyverse)

set.seed(23)

d <- tibble(x = rnorm(10), y = rnorm(10), a = rep(c("x", "y"),5), b = rep(1:5, 2))

d %>%
ggplot(aes(x, y, color = a, size = b)) + geom_point() +
guides(colour = guide_legend(order = 1), ```
Created Aug 24, 2015
View fitting_with_uncertainty.Rmd
 --- title: "Fitting with uncertainty" author: "Michael Kuhn" date: "24 Aug 2015" output: html_document --- In this toy example, we assume that we've independently measured values \$x\$ and \$y\$ and want to find a linear relationship between them, accounting for measurement uncertainty. Each \$x\$ and \$y\$ value is assigned a different uncertainty, and the challenge is to take this information into account. A standard linear model will treat all points equally. When we treat each measurement as a multivariate normal distribution, we can find a point along the proposed fitted line that is maximizing the probability density function (i.e. that has a maximum likelihood). Thus, given a slope and intercept, we virtually move all measurements to their most likely point along the line, and use the likelihood at this point.
Created Sep 28, 2017
View pdz
 Scribble: B7Z2Y1_HUMAN ARHG7_HUMAN ARHG8_HUMAN KCNA5_HUMAN VWCE_HUMAN DNML1_HUMAN PKP4_HUMAN NXPE2_HUMAN
Created Mar 8, 2016
Count comparison results without additional memory allocation
View inplace_comparison.R
 library(Rcpp) `%count<%` <- cppFunction(' size_t count_less(NumericVector x, NumericVector y) { const size_t nx = x.size(); const size_t ny = y.size(); if (nx > 1 & ny > 1) stop("Only one parameter can be a vector!");
Last active Dec 21, 2015 — forked from lsauer/InChi.js
View InChi.js
 // International Chemical Identifier Regex, by lo sauer - lsauer.com // Morphine InchI: var x="InChI=1S/C17H19NO3/c1-18-7-6-17-10-3-5-13(20)16(17)21-15-12(19)4-2-9(14(15)17)8-11(10)18/h2-5,10-11,13,16,19-20H,6-8H2,1H3/t10-,11+,13-,16-,17-/m0/s1" // applying an organic character-subset // we could check for the length property, but in case of 0 matches 'null' is returned -> hence !!.. \ generally equal to Boolean(..) !!x.trim().match(/^((InChI=)?[^J][0-9BCOHNSOPrIFla+\-\(\)\\\/,pqbtmsih]{6,})\$/ig) >true //generic: x.trim().match(/^((InChI=)?[^J][0-9a-z+\-\(\)\\\/,]+)\$/ig)
Created Apr 22, 2013
View plot with densities.R
 library(ggplot2) library(gridExtra) mtcars\$cyl <- ordered(mtcars\$cyl) p <- ggplot(mtcars, aes(mpg, hp, colour = cyl)) + geom_point() p1 <- p + theme(legend.position = "none") p2 <- ggplot(mtcars, aes(x=mpg, group=cyl, colour=cyl)) p2 <- p2 + stat_density(fill = NA, position="dodge")
Last active Sep 30, 2015