Q026 - Quali sono le modalità di accesso al laboratorio?


In [1]:
# -*- coding: UTF-8 -*-

# Render our plots inline
%matplotlib inline 

import pandas as pd
import matplotlib
import matplotlib.pyplot as plt
import numpy as np
import seaborn
import shutil

pd.set_option('display.mpl_style', 'default') # Make the graphs a bit prettier, overridden by seaborn
pd.set_option('display.max_columns', None) # Display all the columns
plt.rcParams['font.family'] = 'sans-serif' # Sans Serif fonts for all the graphs

# Reference for color palettes: http://web.stanford.edu/~mwaskom/software/seaborn/tutorial/color_palettes.html

# Change the font
matplotlib.rcParams.update({'font.family': 'Source Sans Pro'})

In [2]:
# Load csv file first
data = pd.read_csv("data/lab-survey.csv", encoding="utf-8")

In [3]:
# Check data
#data[0:4] # Equals to data.head()

In [4]:
# Range: D26[SQ001] - D26[SQ008] - D21[other]

lab_columns = ['D26[SQ001]','D26[SQ002]','D26[SQ003]','D26[SQ004]','D26[SQ005]','D26[SQ006]','D26[SQ007]','D26[SQ008]']
lab_options = ['Su prenotazione',
               'Su invito',
               'Per partecipazione a corsi o eventi',
               'Accesso libero gratuito (senza registrazione)',
               'Accesso libero gratuito (con registrazione)',
               'Accesso libero con tessera mensile',
               'Accesso libero con tessera annuale',
               'Accesso libero con tessera (altra durata)']
lab = data[lab_columns]
lab.replace(u'Sì', 'Si', inplace=True) # Get rid of accented characters 
lab_other = data['D26[other]'].str.lower().value_counts()


-c:13: SettingWithCopyWarning: 
A value is trying to be set on a copy of a slice from a DataFrame

See the the caveats in the documentation: http://pandas.pydata.org/pandas-docs/stable/indexing.html#indexing-view-versus-copy

In [5]:
#lab[0:4]

In [6]:
%%capture output

# Save the output as a variable that can be saved to a file
# Gather data
lab_b = {}

for k,i in enumerate(lab_columns):
    lab_b[k] = lab[i].value_counts(dropna=False)
    print "Data:",lab_options[k]
    print lab_b[k]
    print
    print "Data %:",lab_options[k]
    print lab[i].value_counts(normalize=True,dropna=False)*100
    print

In [7]:
# Save+show the output to a text file
%save Q026-ModalitàAccesso01.py str(output)
shutil.move("Q026-ModalitàAccesso01.py", "text/Q026-ModalitàAccesso01.txt")


The following commands were written to file `Q026-ModalitàAccesso01.py`:
Data: Su prenotazione
No     35
Si     35
NaN     0
dtype: int64

Data %: Su prenotazione
No     50
Si     50
NaN     0
dtype: float64

Data: Su invito
No     53
Si     17
NaN     0
dtype: int64

Data %: Su invito
No     75.714286
Si     24.285714
NaN     0.000000
dtype: float64

Data: Per partecipazione a corsi o eventi
Si     42
No     28
NaN     0
dtype: int64

Data %: Per partecipazione a corsi o eventi
Si     60
No     40
NaN     0
dtype: float64

Data: Accesso libero gratuito (senza registrazione)
No     51
Si     19
NaN     0
dtype: int64

Data %: Accesso libero gratuito (senza registrazione)
No     72.857143
Si     27.142857
NaN     0.000000
dtype: float64

Data: Accesso libero gratuito (con registrazione)
No     57
Si     13
NaN     0
dtype: int64

Data %: Accesso libero gratuito (con registrazione)
No     81.428571
Si     18.571429
NaN     0.000000
dtype: float64

Data: Accesso libero con tessera mensile
No     62
Si      8
NaN     0
dtype: int64

Data %: Accesso libero con tessera mensile
No     88.571429
Si     11.428571
NaN     0.000000
dtype: float64

Data: Accesso libero con tessera annuale
No     40
Si     30
NaN     0
dtype: int64

Data %: Accesso libero con tessera annuale
No     57.142857
Si     42.857143
NaN     0.000000
dtype: float64

Data: Accesso libero con tessera (altra durata)
No     63
Si      7
NaN     0
dtype: int64

Data %: Accesso libero con tessera (altra durata)
No     90
Si     10
NaN     0
dtype: float64



In [8]:
yes = []
no = []
nanvalue = []

for k,i in enumerate(lab_columns):
    lab_presents = lab_b[k].index.tolist()
    
    # Convert NaN to "NaN"
    for o,h in enumerate(lab_presents):
        if type(h) is float:
            lab_presents.pop(o)
            lab_presents.append("NaN")
    
    # Reassign new list with "NaN"
    lab_b[k].index = lab_presents
    
    # Check for empty values, and put a 0 instead
    if "Si" not in lab_presents:
        yes.append(0)
    if "No" not in lab_presents:
        no.append(0)
    if "NaN" not in lab_presents:
        nanvalue.append(0)
    
    for j in lab_presents:
        if j == "Si":
            yes.append(lab_b[k].ix["Si"])
        elif j == "No":
            no.append(lab_b[k].ix["No"])
        elif j == "NaN":
            nanvalue.append(lab_b[k].ix["NaN"])

In [9]:
# Plot the data
plt.figure(figsize=(8,6))
plt.xlabel(u'Modalità di accesso', fontsize=16)
plt.ylabel(u'Lab', fontsize=16)
plt.title(u'Quali sono le modalità di accesso al laboratorio?', fontsize=18, y=1.02)
plt.xticks(range(len(lab_options)),lab_options,rotation=90)
ind = np.arange(len(lab_columns))   # the x locations for the groups
width = 0.25                              # the width of the bars

my_colors = seaborn.color_palette("Set1", 3) # Set color palette
rect1 = plt.bar(ind,yes,width,color=my_colors[1],align='center') # Plot Yes
rect2 = plt.bar(ind+width,no,width,color=my_colors[0],align='center') # Plot No 
rect3 = plt.bar(ind+width*2,nanvalue,width,color=my_colors[2],align='center') # Plot NaN 
plt.legend( (rect1, rect2, rect3), ('Si', 'No', 'Nessuna risposta') )
plt.savefig(u"svg/Q026-ModalitàAccesso01.svg")
plt.savefig(u"png/Q026-ModalitàAccesso01.png")
plt.savefig(u"pdf/Q026-ModalitàAccesso01.pdf")



In [10]:
%%capture output

# Save the output as a variable that can be saved to a file
# Get "other" data
lab_other = data["D26[other]"].str.lower().value_counts()
print "Data:"
print lab_other
print ""
print "Data %:"
print data["D26[other]"].str.lower().value_counts(normalize=True) * 100

In [11]:
# Save+show the output to a text file
%save Q026-ModalitàAccesso02.py str(output)
shutil.move("Q026-ModalitàAccesso02.py", "text/Q026-ModalitàAccesso02.txt")


The following commands were written to file `Q026-ModalitàAccesso02.py`:
Data:
fab lab pomeriggio lavoratori di giorno                         1
macchine grosse su prenotazione e gratuite in parte             1
siamo alla fase preliminare non è ancora attivo lo spazio comune ufficiale e dunque le regole non sono definite    1
6 mesi                                                          1
dtype: int64

Data %:
fab lab pomeriggio lavoratori di giorno                         1.428571
macchine grosse su prenotazione e gratuite in parte             1.428571
siamo alla fase preliminare non è ancora attivo lo spazio comune ufficiale e dunque le regole non sono definite    1.428571
6 mesi                                                          1.428571
dtype: float64


In [12]:
# Plot bar
plt.figure(figsize=(8,6))
plt.title(u'Quali sono le modalità di accesso al laboratorio? Altro', fontsize=18, y=1.02)
plt.xticks(range(len(lab_other.index)),lab_other.index,rotation=90)
plt.xlabel(u'Modalità di accesso', fontsize=16)
plt.ylabel('Lab', fontsize=16)
ind = np.arange(len(lab_other))   # the x locations for the groups
width = 0.35                       # the width of the bars

my_colors = seaborn.color_palette("husl", len(lab_other)) # Set color palette
rect1 = plt.bar(ind,lab_other,width,color=my_colors,align='center')
plt.savefig(u"svg/Q026-ModalitàAccesso02.svg")
plt.savefig(u"png/Q026-ModalitàAccesso02.png")
plt.savefig(u"pdf/Q026-ModalitàAccesso02.pdf")



In [12]: