from __future__ import division # no es necesario comprender esta línea por ahora
import numpy as np # importamos la libreria numpy para poder usarla en cálculos.

%load_ext rmagic 
# Esta permite usar R en esta consola con el prefijo "%R "

# en python los vectores se encierran entre []
# para asignar valor a una variable en python usamos el "="
pedades = [44, 43, 48, 39, 33, 41, 50, 33]
print pedades # mostrar en pantalla el valor de esta variable. Podemos no usar el print, solo lo usaremos esta vez.

# la definicion anterior es de python "puro", en adelante vamos a usar numpy
pedades = np.array([44, 43, 48, 39, 33, 41, 50,33]) # descarto el valor anterior de la variable "pedades"
pedades # imprimime en pantalla el valor de pedades.

# en R la funcion c() permite construir un vector
# para asignar valor a una variable en R usamos el "<-"
%R redades <- c(44, 43, 48, 39, 33, 41, 50, 33)
# no usamos print o comando similar, la consola muestra siempre en resultado de lo que procesamos

# en R la media se calcula con la funcion mean()
%R mean(redades)

%R sum(redades)/length(redades) # sin la funcion mean se calcularia asi en R

# en numpy tambien
np.mean(pedades)

# sin la funcion mean en python y con numpy la media se calcularía:
np.sum(pedades)/pedades.size

np.median(pedades)

%R median(redades)

# al parecer numpy no incluye una funcion para calcular la moda y si la librería scipy
import scipy.stats as sci
sci.mode(pedades)

# el primer valor es la "moda" y el segundo es la cantidad de repeticiones.

# en R pasa algo similar, la funcion para calcular la moda no esta entre las funciones estándar
%R library(modeest) # carga la biblioteca modeest para dispones de la funcion que calcula la moda mfv.
%R mfv(redades)

#en R sería
%R sum((redades-mean(redades))^2)/(length(redades)-1)

# pero se puede calcular como la funcion var()
%R var(redades)

# con numpy el calculo es 
np.sum( (pedades-np.mean(pedades)) ** 2 / (pedades.size))

# pero no es necesario semejante cálculo, con en R numpy tiene una funciona para la varianza
np.var(pedades) # tambien funcionaria con pedades.var()

%R sqrt(var(redades)) #forma simple en R

%R sd(redades) # la funciona que ya hace este calculo es sd()

#en numpy se calcula con la funcion std()
np.std(pedades)