可视化学习.day3

第三天

使用Pygal模拟掷骰子

0.安装Pygal

在Linux和OS X系统中，应执行的命令类似于下面这样：

pip install --user pygal

在Windows系统中，命令类似于下面这样：

python -m pip install --user pygal

同样也可以：

pip install pygal --user

--user可以不加，但是如果系统提示没有权限，则一定需要加上

Pygal画廊

要了解使用Pygal可创建什么样的图表，请查看图表类型画廊：访问http://www.pygal.org/ ，单击Documentation，再单击Chart types。每个示例都包含源代码，让你知道这些图表是如何 生成的。

创建Die类

下面的类是在模拟掷骰(tou)子：
die.py

from random import randint

class Die():
	"""表示一个骰子的类"""
	
	def __init__(self,num_sides=6):
		"""骰子默认为六面"""
		self.num_sides = num_sides

	def roll(self):
		"""返回一个位于1和骰子面数之间的随机值"""
		return randint(1,self.num_sides)

方法__init__()接受一个可选参数，创建这个类的实例时，如果没有指定任何实参时，面数则会默认为6，如果指定了实参，这个值用于设置骰子的面数。骰子是根据面数命名的，6面的骰子名为D6，8面的骰子名为D8，以此类推。
方法roll()使用函数randint()来返回一个1和面数之间的随机数。这个函数可能返回起始值1、终止值num_sides或这两个数之间的任何整数。

掷骰子

使用这个类来创建图表前，先来掷骰子，将结果打印出来，并检查结果是否合理：
die_visual.py

from die import Die

#创建一个D6
die = Die()

#掷几次骰子，并将结果存储在一个列表中
results = []
for roll_num in range(100): #这里的roll_num可以换成任意一个为使用过的变量，就是是 _ 也可以。（_表示占位符）
	result = die.roll()
	results.append(result)

print(results)

我们创建了实例die（而且没有给予参数，所以num_sides默认为6），我们用使用for roll_num in range(100),使骰子掷100次。并且每一次的结果都存储在了列表results中。
结果大致如下（运算的结果相同的概率是1/6**100那么多）

[4, 6, 5, 6, 1, 5, 6, 3, 5, 3, 5, 3, 2, 2, 1,
3, 1, 5, 3, 6, 3, 6, 5, 4, 1, 1, 4, 2, 3, 6,
4, 2, 6, 4, 1, 3, 2, 5, 6, 3, 6, 2, 1, 1, 3,
4, 1, 4, 3, 5, 1, 4, 5, 5, 2, 3, 3, 1, 2, 3, 
5, 6, 2, 5, 6, 1, 3, 2, 1, 1, 1, 6, 5, 5, 2, 
2, 6, 4, 1, 4, 5, 1, 1, 1, 4, 5, 3, 3, 1, 3, 
5, 4, 5, 6, 5, 4, 1, 5, 1, 2]

从结果我们可以知道，每个出现的数字都在1~6之间，而且每个数字都出现过，这表面每个数字的出现都是随机的。

分析几个

为了分析掷一个D6骰子的结果，我们计算每个点出现的次数：
die_visual.py

from die import Die

#创建一个D6
die = Die()

#掷几次骰子，并将结果存储在一个列表中
results = []
for roll_num in range(1000): #这里的roll_num可以换成任意一个为使用过的变量，就是是 _ 也可以。（_表示占位符）
	result = die.roll()
	results.append(result)

#分析结果
frequencies = []
for value in range(1,die.num_sides+1):
	frequency = results.count(value)
	frequencies.append(frequency)

print(frequencies)

这样的话，我们没有显示每一次掷出的结果，但是我们根据统计，把每一个点数出现的次数进行了统计。为了使统计更精确我们这次运行1000次，而且为了存储每个点数出现的次数，我们创建了另一个空列表frequencies，用于存储它们。 通过下面这一部分：

for value in range(1,die.num_sides+1):
	frequency = results.count(value)
	frequencies.append(frequency)

我们遍历可能的点数（这里为1~6），计算每个点数在results中出现了多少次，将其传递给变量frequency，然后运用for循环和append()函数，将每个值都加到列表frequencies的末尾，最后把这个列表打印出来。

绘制直方图

有了频率后，我们就可以绘制结果的直方图(条形统计图）了。
die_visual.py

import pygal
from die import Die

#创建一个D6
die = Die()

#掷几次骰子，并将结果存储在一个列表中
results = []
for roll_num in range(1000): #这里的roll_num可以换成任意一个为使用过的变量，就是是 _ 也可以。（_表示占位符）
	result = die.roll()
	results.append(result)

#分析结果
frequencies = []
for value in range(1,die.num_sides+1):
	frequency = results.count(value)
	frequencies.append(frequency)
#对结果进行可视化
hist = pygal.Bar()

hist.title = "Results of rolling one D6 1000 times."
hist.x_labels = ['1', '2', '3', '4', '5', '6']
hist.x_title = "Result"
hist.y_title = "Frequency of Result"

hist.add("D6",frequencies)
hist.render_to_file("die_visual.svg")

为了创建条形图，我们创建了一个pygal.Bar()实例，并将其存储在变量hist中，接下来，我们设置hist的属性:

• hist.title(用于表示直方图名称的的字符串）

• hist.x_labels（用于表示x轴的标签）在这里我们通过hist.x_labels = ['1', '2', '3', '4', '5', '6'],用掷D6骰子的可能的结果做为x轴的标签。

• hist.x_title 与 hist.y_title（分别表示x轴与y轴的标题）。
  然后，我们使用hist.add()将一系列值添加 到图表中（向它们传递要给添加的值指定的标签，还有一个列表，其中包含将出现在图表中的值）。
  最后，我们使用hist.render_to_file()将这个图表渲染成一个名为 “die_visual.svg” 的SYG文件，而且这种危机的扩展名必须为.svg的格式。

   die_visual.svg默认存储在die_visual.py所在的文件夹里。
  

要想查看生成的直方图，最简单的方式就是用Wed浏览器。

同时掷两个面数相同的骰子

同时掷两个骰子，得到的点数更多，结果分布情况也不同。所以我们需要创建两个D6骰子，且将两个骰子的点数相加。并将结果存储在列表results中：
dice_visual.py

import pygal
from die import Die

#创建一个D6
die_1 = Die()
die_2 = Die()

#掷几次骰子，并将结果存储在一个列表中
results = []
for roll_num in range(1000): #这里的roll_num可以换成任意一个为使用过的变量，就是是 _ 也可以。（_表示占位符）
	result = die_1.roll() + die_2.roll()
	results.append(result)

#分析结果
frequencies = []
max_result = die_1.num_sides + die_2.num_sides
for value in range(1,max_result+1):
	frequency = results.count(value)
	frequencies.append(frequency)
#对结果进行可视化
hist = pygal.Bar()

hist.title = "Results of rolling one D6 1000 times."
hist.x_labels = ['2', '3', '4', '5', '6', '7', '8', '9', '10', '11', '12']
hist.x_title = "Result"
hist.y_title = "Frequency of Result"

hist.add("D6 + D6",frequencies)
hist.render_to_file("die_visual.svg")

创建两个Die实例后，我们掷骰子多次，并计算每次的总点数，可能出现的最大点数位12为两个骰子的最大点数之和，我们定义了一个变量max_result，并将最大点数之和存储在其中。尽管我们计算2到max_result的各种点数出现的次数（仅仅使用range(2,13)就可以了），但是这个只适用于数据较小时，但在模拟现实世界情况，最后编写轻松地模拟各种情形的代码。因此，不管骰子有多少面，我们都可以计算出。

同时掷两个面数不同的骰子

比如我们创建一个6面和10面的骰子。看看同时运行这两个骰子50000次的结果。
different_dice.py

import pygal
from die import Die

#创建一个D6
die_1 = Die()
die_2 = Die(10)

#掷几次骰子，并将结果存储在一个列表中
results = []
for roll_num in range(50000): #这里的roll_num可以换成任意一个为使用过的变量，就是是 _ 也可以。（_表示占位符）
	result = die_1.roll() + die_2.roll()
	results.append(result)

#分析结果
frequencies = []
max_result = die_1.num_sides + die_2.num_sides
for value in range(1,max_result+1):
	frequency = results.count(value)
	frequencies.append(frequency)
#对结果进行可视化
hist = pygal.Bar()

hist.title = "Results of rolling a D6 and a D10 50000 times."
hist.x_labels = ['2', '3', '4', '5', '6', '7', '8', '9', '10', '11', '12', '13', '14', '15', '16']
hist.x_title = "Result"
hist.y_title = "Frequency of Result"

hist.add("D6 + D10",frequencies)
hist.render_to_file("die_visual.svg")