Skip to content

Instantly share code, notes, and snippets.

@cryptospectrum
cryptospectrum / jviolin_D0Wplt1.py
Last active February 13, 2017 04:34
Python +Seaborn Violin Plot
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
sns.set_style("whitegrid")
filename='QuoteAAPL_2011-9-5_T_2013-1-401'
csv_filename=filename+'.csv'
QuotesAAPLDF=pd.read_csv(csv_filename)
@cryptospectrum
cryptospectrum / jpanel_33plt.py
Created February 13, 2017 04:34
Python+Seaborn Multi-panel Plot
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
sns.set(style="ticks")
# np.random.seed(sum(map(ord, "axis_grids")))
filename='QuoteAAPL_2011-9-5_T_2013-1-401'
csv_filename=filename+'.csv'
QuotesAAPLDF=pd.read_csv(csv_filename)
@cryptospectrum
cryptospectrum / jlmreg_hist_plt1.py
Created February 13, 2017 04:32
Python+Seaborn plot Regression with histograms
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
sns.set(style="ticks")
filename='QuoteAAPL_2011-9-5_T_2013-1-401'
csv_filename=filename+'.csv'
QuotesAAPLDF=pd.read_csv(csv_filename)
# # More data for scatter plots
@cryptospectrum
cryptospectrum / q_symbol_dict.py
Created February 13, 2017 04:31
Quote Symbol Dictionary, use with jcan_scttM2.py
def sym_dict(opt):
if opt==0:
symbol_dict = {
'XOM': 'Exxon',
'CVX': 'Chevron',
'COP': 'ConocoPhillips',
'VLO': 'Valero Energy',
'MSFT': 'Microsoft',
'IBM': 'IBM',
'T': 'AT&T',
@cryptospectrum
cryptospectrum / jcan_scttM2.py
Created February 13, 2017 04:29
Candle Stick Scatter Plot
import datetime
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.cm as cm
import time
import pandas as pd
# # import json
# # Keep q_symbol_dict.py in same directory
import q_symbol_dict
try:
@cryptospectrum
cryptospectrum / dna2midi_csum444.m
Created September 29, 2013 22:02
Turn DNA sequence into a MIDI file
% %loadfile.m @http://www.ece.rutgers.edu/~orfanidi/ece525/subspace.zip
[Y,T] = loadfile('Saccharomyces_cerevisiae_GENES.txt');
[RT,CT]=size(T);
for n=1:CT
DNA((n-1)*RT +1:n*RT)=T(:,n);
end
% % NAk=size(DNA,2);
NAk=1200;
N=floor(NAk/5);
M=zeros(N,6);
@cryptospectrum
cryptospectrum / hist_eqcv.m
Created July 22, 2013 02:42
Histogram Equalization
function [imeq,pdfim,cdfim]=hist_eqcv(im, nbins)
%__________________________________________________________________
% Description: Historgram Equalization
%
% Input Variables:
% Assume
% im = double(imread('image_file_name'));
% im=im(:,:,1);
% Get image format back
% im_retrieve=uint8(round(im));
@cryptospectrum
cryptospectrum / wiki_latex_strip.py
Last active February 1, 2017 06:57
Scraping Latex Equations from Wikipedia Pages
import sys, requests
# requests is a third party module, more robust that urllib...
#http://docs.python-requests.org/en/latest/
from lxml import html
if len(sys.argv)==1:
#url_name="https://en.wikipedia.org/wiki/Fa%C3%A0_di_Bruno%27s_formula"
#url_name="https://secure.wikimedia.org/wikipedia/en/wiki/Euler%E2%80%93Maclaurin_formula"
#url_name="https://en.wikipedia.org/wiki/Generalized_normal_distribution"
url_name="https://en.wikipedia.org/wiki/Jacobi_elliptic_functions"
@cryptospectrum
cryptospectrum / midi_logistic_svm.m
Created July 14, 2013 21:21
Analyze Midi Files with Machine Leaning Algorithms
midi_fname='Scarlatti_Sonata_Longo483_K322.mid';
% MIDI scripts from: http://www.kenschutte.com/midi
midi = readmidi(midi_fname);
Notes = midiInfo(midi,0);
%diff start and stop times
dNo5=diff(Notes(:,5));
dNo6=diff(Notes(:,6));
% To Reconstruct Time Data
@cryptospectrum
cryptospectrum / memskit_03.ino
Created April 9, 2013 23:07
Arduino + accelerometer
/*
Memsic2125
Read the Memsic 2125 two-axis accelerometer. Converts the
pulses output by the 2125 into milli-g's (1/1000 of earth's
gravity) and prints them over the serial connection to the
computer.
The circuit:
* X output of accelerometer to digital pin 2
* Y output of accelerometer to digital pin 3