is, I think, the one you can find at sacred-texts.org. curl -O 'http://www.sacred-texts.com/dna/hgp011k.htm' #get it #BORING DATA JANITORSHIP tail -n +15 hgp011k.htm __ hgp011k #remove the HTML head stuff .. up to head -n -3 hgp011k | sponge hgp011k #remove the ... [Read more...]

Clay tablet YBC 7289 with an overlay to emphasize the markings on its surface. The cuneiform inscriptions demonstrate the derivation of the square root of 2 … in sexagesimal (base-60) via Paul Murrell the tablet has been dated to sometime between 1800 BC and 1600 BC. Little is known of its rediscovery, except that it ... [Read more...]

Catalan numbers for triangulations gambler’s ruin binary trees trees …from the Flajolet/Sedgewick coursera analysis of algorithms in R: catalan % unlist %__% sum } (This is slow above catalan(10) though, so here’s a memoised version. Then if you want to run above catalan(250) you’ll need to adjust options(expressions=5... [Read more...]

[Read more...]

2006–2014 using FRBData 2006–2014 using FRBData packageTreasury yield curve from the Volcker era through Greenspan, Bernanke, and Yellen. require(YieldCurve) data(FedYieldCurve) maturities I can’t upload a larger view, but just run that code with a bigger width & height, you might get a ~10MB gif. I shrank it further with gifsicle. ... [Read more...]

List of R programmers: Hello R people. In December of 2013 I posted a cheap-o wiki-editable (thank you github) contact list which recruiters can use to find you, if they’re looking for R programmers. In what I consider a resounding success, within a few weeks it got onto the first ... [Read more...]

√(x²−1)(x²−k²), k∈ℂ x²√(x²−1)(x²−k²).      x,k∈ℂ (actually just going over the unit circle, not all of ℂ) edit: hey, are these showing up as moving gif’s for you? code: require(animation) source(wegert.R) #where I define "plat" and "Z", standard for all Wegert plots i = ... [Read more...]

@tdhopper posted his self-measurements of weight loss a few months back. I recently decided also that I wanted to lose fat-weight—the infamous “I could stand to be a few kilos lighter”—and I think I came up with a more productive way of thinking about my progress: I’m ... [Read more...]

29 39 59 79 99 119 "139"__ plot( polyroot(choose(131,14:29)) ,pch=19,col='red') __ plot( polyroot(choose(131,14:39)) ,pch=19,col='red') __ plot( polyroot(choose(131,14:59)) ,pch=19,col='red') __ plot( polyroot(choose(131,14:79)) ,pch=19,col='red') __ plot( polyroot(choose(131,14:99)) ,pch=19,col='red') __ plot( polyroot(choose(131,14:119)) ,pch=19,col='red') __ plot( polyroot(choose(131,14:139)) ,pch=19,col='red') [Read more...]

[Read more...]

A question I’ve googled before without success. Hopefully this answer will show up for someone who needs it. I’ll also go over the better-known uses of ? just in case. To get help in R about a function like subset you type ?subset . That’s like man subset from ... [Read more...]

A question I’ve googled before without success. Hopefully this answer will show up for someone who needs it. I’ll also go over the better-known uses of ? just in case. To get help in R about a function like subset you type ?subset . That’s like man subset from ... [Read more...]

tl,dr: If you want to be contacted for freelance R work, edit this list https://github.com/isomorphisms/hire-an-r-programmer/blob/master/README.md.   Background/Problem: I was looking for a list of R-programming freelancers and realised there is no such list. Other than famous people and people I already ... [Read more...]

Two interesting ideas here: "trading time" price impact of a trade proportional to exp( √size ) Code follows: require(quantmod) getSymbols("MER") #Merrill Lynch #Gatheral's model HiLo Op(symbol) #munging mer names(mer)[7] = "UpDay"names(mer)[8] = "HiLo" mer 0) #data cleaning #plot layers require(ggplot2) base 0) losers base + points + line + labels + colour + ... [Read more...]

sine mis-calculated the angles on sine mid-res, small cex, pch=46, sine low-res, small cex, pch=46 z^17th power sine in HCL rather than HSV cosh z+zz+zzz+zzzz+zzzzz+zzzzzz+zzzzzzz+zzzzzzzz+zzzzzzzzz+zzzzzzzzzzplaying along with Elias Wegert in R: X * complex(imaginary=.05) + t(X)/20 #twist & shout X hcl(... [Read more...]

“I like this concept of “low volatility, interrupted by occasional periods of high volatility”. I think I will call it “volatility”.” - Daniel Davies via nonergodic   (PS: If you didn’t see it before: try plotting this in R: vol.of.vol … and so on, to your heart’s content. ... [Read more...]

only the most active trading days $HLF (HerbaLife weight-loss supplements / MLM) $HLF regular history big loss days and big volume days for $HLF. "Ackman" should instead read "Einhorn".Whilst reading John Hempton’s post on shorting $HLF I decided to follow along in quantmod. Bronte Capital: It was the night ... [Read more...]

I was re-reading Michael Murray’s explanation of cointegration: and marvelling at the calculus. Calculus blows my mind sometimes. Like, hey guess how much we can do with subtraction. — protëa(@isomorphisms) March 28, 2013 Of course it’s not any subtraction. It’s subtracting a function from a shifted version of ... [Read more...]

The Cauchy distribution (?dcauchy in R) nails a flashlight over the number line and swings it at a constant speed from 9 o’clock down to 6 o’clock over to 3 o’clock. (Or the other direction, from 3→6→9.) Then counts how much light shone on each number. In other words we want ... [Read more...]

SETUP (CAN BE SKIPPED) We start with data (how was it collected?) and the hope that we can compare them. We also start with a question which is of the form: how much tax increase is associated with how much tax avoidance/tax evasion/country fleeing by the top 1%? how ... [Read more...]