Skip to content

Instantly share code, notes, and snippets.

Keybase proof

I hereby claim:

  • I am peitalin on github.
  • I am peita (https://keybase.io/peita) on keybase.
  • I have a public key ASAcLMvAKFK0RI_-2hPyaxZjAeH1p19AxYCA55O_Np2B6go

To claim this, I am signing this object:

Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@peitalin
peitalin / linreg.py
Last active December 18, 2015 03:49
import numpy as np
import statsmodels.api as sm
import pandas as pd
def FF3_Reg(y, x):
ones = np.ones(len(x[0]))
X = sm.add_constant(np.column_stack((x[0], ones)))
for factor in x[1:]:
X = sm.add_constant(np.column_stack((factor, X)))
estimates = sm.OLS(y, X).fit()
@peitalin
peitalin / lingreg.hs
Last active December 18, 2015 03:49
CSV parsing with cassava and lens -> Hmatrix types -> basic OLS regression.
{-# LANGUAGE FlexibleContexts #-}
{-# LANGUAGE ConstraintKinds #-}
import qualified Data.ByteString.Lazy.Char8 as BL
import qualified Data.Vector as V (toList)
import Data.Csv
import Data.List
import Control.Lens hiding ((<.>))
import Numeric.LinearAlgebra