{
"metadata": {
"name": "Untitled0"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"import pandas as pd\n",
"from bs4 import BeautifulSoup as bs\n",
"import TOU_pricing\n",
"import numpy as np\n",
"#reload(TOU_pricing)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 6
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"xml_file = 'raw_data/Inland_Single_Family_Jan_1_2011_to_Jan_1_2012_RetailCustomer_9.xml'\n",
"soup = bs(open(xml_file))\n",
"\n",
"# initiate a numpy array to store values\n",
"data_array = np.zeros((len(soup.find_all('value')),1))\n",
"\n",
"for i, intervalreading in enumerate(soup.find_all('intervalreading')):\n",
" data_array[i,0]=intervalreading.value.get_text()\n",
" \n",
"data_start = '1/1/2011 08:00:00'\n",
"demand_data = pd.DataFrame(data_array*.001,columns=['USAGE'],index=pd.date_range(data_start, periods=len(soup.find_all('value')), freq='H'))\n",
"#demand_data.to_csv('raw_data/Green_Button_Sample_Inland_SingleFamily.csv')"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 7
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"demand_data[:2]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"
\n",
"
\n",
" \n",
" \n",
" | \n",
" USAGE | \n",
"
\n",
" \n",
" \n",
" \n",
" 2011-01-01 08:00:00 | \n",
" 1.002 | \n",
"
\n",
" \n",
" 2011-01-01 09:00:00 | \n",
" 0.867 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"output_type": "pyout",
"prompt_number": 8,
"text": [
" USAGE\n",
"2011-01-01 08:00:00 1.002\n",
"2011-01-01 09:00:00 0.867"
]
}
],
"prompt_number": 8
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"a = TOU_pricing.BGE_elec_cost(demand_data)\n",
"a[30:60]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"
\n",
" \n",
" \n",
" | \n",
" USAGE | \n",
" BGE-R_cost_perkWh | \n",
" BGE-RL_cost_perkWh | \n",
" BGE-EV_cost_perkWh | \n",
"
\n",
" \n",
" \n",
" \n",
" 2011-01-02 14:00:00 | \n",
" 0.926 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 15:00:00 | \n",
" 1.041 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 16:00:00 | \n",
" 1.111 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 17:00:00 | \n",
" 1.193 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 18:00:00 | \n",
" 1.182 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 19:00:00 | \n",
" 1.149 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 20:00:00 | \n",
" 1.111 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 21:00:00 | \n",
" 1.014 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 22:00:00 | \n",
" 0.967 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-02 23:00:00 | \n",
" 1.005 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 00:00:00 | \n",
" 1.090 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 01:00:00 | \n",
" 1.405 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 02:00:00 | \n",
" 1.500 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 03:00:00 | \n",
" 1.404 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 04:00:00 | \n",
" 1.314 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 05:00:00 | \n",
" 1.247 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 06:00:00 | \n",
" 1.089 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 07:00:00 | \n",
" 0.924 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-03 08:00:00 | \n",
" 0.808 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-03 09:00:00 | \n",
" 0.748 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-03 10:00:00 | \n",
" 0.714 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-03 11:00:00 | \n",
" 0.690 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 12:00:00 | \n",
" 0.711 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 13:00:00 | \n",
" 0.829 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 14:00:00 | \n",
" 0.964 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 15:00:00 | \n",
" 0.962 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 16:00:00 | \n",
" 1.050 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-03 17:00:00 | \n",
" 1.096 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-03 18:00:00 | \n",
" 1.111 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-03 19:00:00 | \n",
" 1.127 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"output_type": "pyout",
"prompt_number": 9,
"text": [
" USAGE BGE-R_cost_perkWh BGE-RL_cost_perkWh \\\n",
"2011-01-02 14:00:00 0.926 0.08616 0.06974 \n",
"2011-01-02 15:00:00 1.041 0.08616 0.06974 \n",
"2011-01-02 16:00:00 1.111 0.08616 0.06974 \n",
"2011-01-02 17:00:00 1.193 0.08616 0.06974 \n",
"2011-01-02 18:00:00 1.182 0.08616 0.06974 \n",
"2011-01-02 19:00:00 1.149 0.08616 0.06974 \n",
"2011-01-02 20:00:00 1.111 0.08616 0.06974 \n",
"2011-01-02 21:00:00 1.014 0.08616 0.06974 \n",
"2011-01-02 22:00:00 0.967 0.08616 0.06974 \n",
"2011-01-02 23:00:00 1.005 0.08616 0.06974 \n",
"2011-01-03 00:00:00 1.090 0.08616 0.06974 \n",
"2011-01-03 01:00:00 1.405 0.08616 0.06974 \n",
"2011-01-03 02:00:00 1.500 0.08616 0.06974 \n",
"2011-01-03 03:00:00 1.404 0.08616 0.06974 \n",
"2011-01-03 04:00:00 1.314 0.08616 0.06974 \n",
"2011-01-03 05:00:00 1.247 0.08616 0.06974 \n",
"2011-01-03 06:00:00 1.089 0.08616 0.06974 \n",
"2011-01-03 07:00:00 0.924 0.08616 0.11924 \n",
"2011-01-03 08:00:00 0.808 0.08616 0.11924 \n",
"2011-01-03 09:00:00 0.748 0.08616 0.11924 \n",
"2011-01-03 10:00:00 0.714 0.08616 0.11924 \n",
"2011-01-03 11:00:00 0.690 0.08616 0.08468 \n",
"2011-01-03 12:00:00 0.711 0.08616 0.08468 \n",
"2011-01-03 13:00:00 0.829 0.08616 0.08468 \n",
"2011-01-03 14:00:00 0.964 0.08616 0.08468 \n",
"2011-01-03 15:00:00 0.962 0.08616 0.08468 \n",
"2011-01-03 16:00:00 1.050 0.08616 0.08468 \n",
"2011-01-03 17:00:00 1.096 0.08616 0.11924 \n",
"2011-01-03 18:00:00 1.111 0.08616 0.11924 \n",
"2011-01-03 19:00:00 1.127 0.08616 0.11924 \n",
"\n",
" BGE-EV_cost_perkWh \n",
"2011-01-02 14:00:00 0.05209 \n",
"2011-01-02 15:00:00 0.05209 \n",
"2011-01-02 16:00:00 0.05209 \n",
"2011-01-02 17:00:00 0.05209 \n",
"2011-01-02 18:00:00 0.05209 \n",
"2011-01-02 19:00:00 0.05209 \n",
"2011-01-02 20:00:00 0.05209 \n",
"2011-01-02 21:00:00 0.05209 \n",
"2011-01-02 22:00:00 0.05209 \n",
"2011-01-02 23:00:00 0.05209 \n",
"2011-01-03 00:00:00 0.05209 \n",
"2011-01-03 01:00:00 0.05209 \n",
"2011-01-03 02:00:00 0.05209 \n",
"2011-01-03 03:00:00 0.05209 \n",
"2011-01-03 04:00:00 0.05209 \n",
"2011-01-03 05:00:00 0.05209 \n",
"2011-01-03 06:00:00 0.05209 \n",
"2011-01-03 07:00:00 0.18266 \n",
"2011-01-03 08:00:00 0.18266 \n",
"2011-01-03 09:00:00 0.18266 \n",
"2011-01-03 10:00:00 0.18266 \n",
"2011-01-03 11:00:00 0.05209 \n",
"2011-01-03 12:00:00 0.05209 \n",
"2011-01-03 13:00:00 0.05209 \n",
"2011-01-03 14:00:00 0.05209 \n",
"2011-01-03 15:00:00 0.05209 \n",
"2011-01-03 16:00:00 0.05209 \n",
"2011-01-03 17:00:00 0.18266 \n",
"2011-01-03 18:00:00 0.18266 \n",
"2011-01-03 19:00:00 0.18266 "
]
}
],
"prompt_number": 9
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"trial = '2011-01-06'\n",
"print a[trial].index.weekday\n",
"a.ix[trial]"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"[3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3]\n"
]
},
{
"html": [
"\n",
"
\n",
" \n",
" \n",
" | \n",
" USAGE | \n",
" BGE-R_cost_perkWh | \n",
" BGE-RL_cost_perkWh | \n",
" BGE-EV_cost_perkWh | \n",
"
\n",
" \n",
" \n",
" \n",
" 2011-01-06 00:00:00 | \n",
" 0.906 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 01:00:00 | \n",
" 1.311 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 02:00:00 | \n",
" 1.429 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 03:00:00 | \n",
" 1.383 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 04:00:00 | \n",
" 1.360 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 05:00:00 | \n",
" 1.161 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 06:00:00 | \n",
" 1.025 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 07:00:00 | \n",
" 0.821 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 08:00:00 | \n",
" 0.705 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 09:00:00 | \n",
" 0.652 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 10:00:00 | \n",
" 0.651 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 11:00:00 | \n",
" 0.650 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 12:00:00 | \n",
" 0.699 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 13:00:00 | \n",
" 0.906 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 14:00:00 | \n",
" 1.036 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 15:00:00 | \n",
" 1.051 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 16:00:00 | \n",
" 0.992 | \n",
" 0.08616 | \n",
" 0.08468 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 17:00:00 | \n",
" 0.917 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 18:00:00 | \n",
" 0.893 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 19:00:00 | \n",
" 0.900 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 20:00:00 | \n",
" 0.885 | \n",
" 0.08616 | \n",
" 0.11924 | \n",
" 0.18266 | \n",
"
\n",
" \n",
" 2011-01-06 21:00:00 | \n",
" 0.832 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 22:00:00 | \n",
" 0.809 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
" 2011-01-06 23:00:00 | \n",
" 0.834 | \n",
" 0.08616 | \n",
" 0.06974 | \n",
" 0.05209 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"output_type": "pyout",
"prompt_number": 5,
"text": [
" USAGE BGE-R_cost_perkWh BGE-RL_cost_perkWh BGE-EV_cost_perkWh\n",
"2011-01-06 00:00:00 0.906 0.08616 0.06974 0.05209\n",
"2011-01-06 01:00:00 1.311 0.08616 0.06974 0.05209\n",
"2011-01-06 02:00:00 1.429 0.08616 0.06974 0.05209\n",
"2011-01-06 03:00:00 1.383 0.08616 0.06974 0.05209\n",
"2011-01-06 04:00:00 1.360 0.08616 0.06974 0.05209\n",
"2011-01-06 05:00:00 1.161 0.08616 0.06974 0.05209\n",
"2011-01-06 06:00:00 1.025 0.08616 0.06974 0.05209\n",
"2011-01-06 07:00:00 0.821 0.08616 0.11924 0.18266\n",
"2011-01-06 08:00:00 0.705 0.08616 0.11924 0.18266\n",
"2011-01-06 09:00:00 0.652 0.08616 0.11924 0.18266\n",
"2011-01-06 10:00:00 0.651 0.08616 0.11924 0.18266\n",
"2011-01-06 11:00:00 0.650 0.08616 0.08468 0.05209\n",
"2011-01-06 12:00:00 0.699 0.08616 0.08468 0.05209\n",
"2011-01-06 13:00:00 0.906 0.08616 0.08468 0.05209\n",
"2011-01-06 14:00:00 1.036 0.08616 0.08468 0.05209\n",
"2011-01-06 15:00:00 1.051 0.08616 0.08468 0.05209\n",
"2011-01-06 16:00:00 0.992 0.08616 0.08468 0.05209\n",
"2011-01-06 17:00:00 0.917 0.08616 0.11924 0.18266\n",
"2011-01-06 18:00:00 0.893 0.08616 0.11924 0.18266\n",
"2011-01-06 19:00:00 0.900 0.08616 0.11924 0.18266\n",
"2011-01-06 20:00:00 0.885 0.08616 0.11924 0.18266\n",
"2011-01-06 21:00:00 0.832 0.08616 0.06974 0.05209\n",
"2011-01-06 22:00:00 0.809 0.08616 0.06974 0.05209\n",
"2011-01-06 23:00:00 0.834 0.08616 0.06974 0.05209"
]
}
],
"prompt_number": 5
},
{
"cell_type": "code",
"collapsed": false,
"input": [
" "
],
"language": "python",
"metadata": {},
"outputs": []
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n"
],
"language": "python",
"metadata": {},
"outputs": []
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}