{ "metadata": { "name": "poking_eyeball_rps" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "# import some modules using standard naming conventions\n", "from pylab import *\n", "import pandas as pd\n", "\n", "# find the .csv files and read them into pandas data sets\n", "import glob\n", "flist = glob.glob('/usgs/data1/csherwood/poking_eyeball/*.csv')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "def Phi(mm):\n", " return -1.44296504*log(mm)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "def my_read(csv):\n", " phi = float(csv[-6:-4])/10. # extract phi size from CSV file name\n", " df = pd.read_csv(csv) # create Pandas DataFrame from each CSV file\n", " mg = df.geom_mean.mean() # mean of geometric means\n", " m = df.arith_mean \n", " ma = m.mean() # mean of arithmetic means\n", " sa = m.std() # std of arithmetic means\n", " return phi,ma,sa,mg" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "# list of stats for each file\n", "stats = [my_read(f) for f in flist]\n", "# convert to array\n", "stats = array(stats)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "# create a data frame from the stats array\n", "df_stats = pd.DataFrame(stats,columns=['Phi','Ma','Sa','Mg'])\n", "df_stats" ], "language": "python", "metadata": {}, "outputs": [ { "html": [ "
\n", " | Phi | \n", "Ma | \n", "Sa | \n", "Mg | \n", "
---|---|---|---|---|
0 | \n", "0.0 | \n", "0.557322 | \n", "0.067476 | \n", "0.468722 | \n", "
1 | \n", "0.5 | \n", "0.543682 | \n", "0.134539 | \n", "0.430470 | \n", "
2 | \n", "1.0 | \n", "0.487743 | \n", "0.082410 | \n", "0.381277 | \n", "
3 | \n", "1.5 | \n", "0.373238 | \n", "0.048243 | \n", "0.267245 | \n", "
4 | \n", "2.0 | \n", "0.326144 | \n", "0.031638 | \n", "0.209683 | \n", "
5 | \n", "2.5 | \n", "0.301667 | \n", "0.034709 | \n", "0.199168 | \n", "
6 | \n", "3.0 | \n", "0.235943 | \n", "0.032090 | \n", "0.134345 | \n", "
7 | \n", "3.5 | \n", "0.238555 | \n", "0.015709 | \n", "0.135005 | \n", "
8 | \n", "4.0 | \n", "0.204032 | \n", "0.015977 | \n", "0.106308 | \n", "