Últimamente parece que me persiguen los Fractales. Es una especie de atracción que no sé de dónde me está viniendo, seguro que significa algo pero aún no sé el qué.
El caso es que, después de leer el primer libro de la colección El Mundo es Matemático (a la que me suscribí inmediatamente), me llamó la atención el capítulo de los fractales, de los que no había oído hablar siquiera.
Ahora estoy leyendo la obra "germen" de este campo, Los objetos Fractales, de Benoît Mandelbrot, que tenía en casa en una edición de Círculo de Lectores y todavía no me había leído. Otro día recomendaré este libro, cuando lo acabe de leer.
El caso es que compré la revista que os recomiendo, buscando el artículo "Fractales, ¡QUÉ ARTE!". El artículo en sí no es gran cosa una vez que los has estudiado un poco, pero es muy ameno para quien no los conozca. Sin embargo la revista en conjunto me ha gustado tanto, que he osado recomendarla en este foro, porque trata temas muy variados, con entrevistas muy interesantes, reportajes, relatos de ficción, ...
En fin, os aconsejo humildemente que la leáis, a mí me ha gustado. De paso he conocido la web de Punset, que es un genio más por cómo explica las cosas que por lo que sabe.
La revista se puede hojear en la siguiente dirección:
http://www.youkioske.com/ciencia/redes-para-la-ciencia-mayo-2010/
Bitácora personal de Emilio López Cano, Estadístico (Data Scientist dicen ahora). Enseñando Estadística y Matemáticas en la Universidad de Castilla-La Mancha, y Data Science en la Universidad Rey Juan Carlos. Investigando en proyectos de la Universidad Rey Juan Carlos. Contribuyendo a estándares en AENOR. Formando en la Asociación Estadística de la Calidad. Escribo sobre estos temas.
martes, 11 de mayo de 2010
miércoles, 5 de mayo de 2010
Control de Calidad con R: Muestreos de Aceptación
Para el diseño y verificación de muestreos de aceptación, disponemos en R del paquete AcceptanceSampling.
Uso del paquete AcceptanceSampling
Muestreo simple
Primero tenemos que cargar la librería, en el intérprete de comandos, o alguna de las interfaces gráficas (RCommander o RKWard):
> library(AcceptanceSampling)
Crear un plan de muestreo (ejemplo: n=300, c=5, distribución binomial):
> x <- OC2c(300, 5, type="b")
Mostrar el plan de muestreo:
> x
Mostrar el plan de muestreo detallado (muchos puntos OC):
> summary(x, full=TRUE)
Mostrar gráfico de la curva OC:
> plot(x, xlim = c(0, 0.3))
Encontrar un plan, dada OC (alfa=0.05, p1=0.05, beta=0.15, p2=0.10, distribución binomial):
> find.plan(PRP = c(0.05, 0.95), CRP = c(0.15, 0.10), type = "binom")
Muestreo doble y múltiple
Se procede igual que para el muestreo simple, pero indicando los valores de n, c y r como vectores:
> x <- OC2c(n=c(50, 100), c=c(0, 1), r=(2, 2))
Otros recursos
En la siguiente página se pueden calcular automáticamente muestreos de aceptación utilizando las normas MIL-STD y otras:
www.sqconline.com/
Gracias por el aporte María.
Paquete AcceptanceSampling
Disponible para instalar de los repositorios
Artículo en Journal of Statistical Software
Uso del paquete AcceptanceSampling
Muestreo simple
Primero tenemos que cargar la librería, en el intérprete de comandos, o alguna de las interfaces gráficas (RCommander o RKWard):
> library(AcceptanceSampling)
Crear un plan de muestreo (ejemplo: n=300, c=5, distribución binomial):
> x <- OC2c(300, 5, type="b")
Mostrar el plan de muestreo:
> x
Mostrar el plan de muestreo detallado (muchos puntos OC):
> summary(x, full=TRUE)
Mostrar gráfico de la curva OC:
> plot(x, xlim = c(0, 0.3))
Encontrar un plan, dada OC (alfa=0.05, p1=0.05, beta=0.15, p2=0.10, distribución binomial):
> find.plan(PRP = c(0.05, 0.95), CRP = c(0.15, 0.10), type = "binom")
Muestreo doble y múltiple
Se procede igual que para el muestreo simple, pero indicando los valores de n, c y r como vectores:
> x <- OC2c(n=c(50, 100), c=c(0, 1), r=(2, 2))
Otros recursos
En la siguiente página se pueden calcular automáticamente muestreos de aceptación utilizando las normas MIL-STD y otras:
www.sqconline.com/
Gracias por el aporte María.
martes, 4 de mayo de 2010
Redes Neuronales con R
Introducción
Para empezar, existen distintos modelos de redes neuronales, siendo uno de los más utilizados el del "Perceptron".
Esta red se basa en una "caja negra", donde lo importante es la predicción, y no cómo se hace.
El proceso incluye una fase de entrenamiento (training) para la optimización de las predicciones.
Los elementos de la red son:
Una descripción "grosso modo" podría ser: Los nodos (neuronas) de la capa de entrada, se combinan con los nodos de la capa oculta mediante la función de combinación, que suele ser una combinación lineal de los nodos de entrada mediante los pesos. A las neuronas de las capas ocultas, se aplica una función de activación, que suele ser la tangente hiperbólica de la anterior combinación más un parámetro por nodo oculto, con la que estimamos las neuronas de la capa de salida, y sus errores.
Recursos:
En Wikipedia: RNA, Perceptrón.
En este documento de la Universidad Carlos III de Madrid hay una introducción bastante bien explicada, con ejemplos en R:
Introducción a las redes neuronales aplicadas uc3m)
Paquete nnet
Se puede instalar desde los repositorios
Documentación
Paquete RWeka
Se puede instalar desde los repositorios
Documentación
Paquete neural
Documentación
neural no está en los repositorios, para instalar:
Paquete AMORE
Se puede instalar desde los repositorios. Es un paquete desarrollado en la Universidad de León, con abundante documentación.
Documentación
Wiki
Ejemplo
Próximamente
Para empezar, existen distintos modelos de redes neuronales, siendo uno de los más utilizados el del "Perceptron".
Esta red se basa en una "caja negra", donde lo importante es la predicción, y no cómo se hace.
El proceso incluye una fase de entrenamiento (training) para la optimización de las predicciones.
Los elementos de la red son:
- Las neuronas o nodos
- Las capas
- De entrada
- De salida
- Oculta (puede tener a su vez varias capas)
- Los pesos
- La función de combinación
- La función de activación
- El objetivo (target)
Una descripción "grosso modo" podría ser: Los nodos (neuronas) de la capa de entrada, se combinan con los nodos de la capa oculta mediante la función de combinación, que suele ser una combinación lineal de los nodos de entrada mediante los pesos. A las neuronas de las capas ocultas, se aplica una función de activación, que suele ser la tangente hiperbólica de la anterior combinación más un parámetro por nodo oculto, con la que estimamos las neuronas de la capa de salida, y sus errores.
Recursos:
En Wikipedia: RNA, Perceptrón.
En este documento de la Universidad Carlos III de Madrid hay una introducción bastante bien explicada, con ejemplos en R:
Introducción a las redes neuronales aplicadas uc3m)
Paquete nnet
Se puede instalar desde los repositorios
Documentación
Paquete RWeka
Se puede instalar desde los repositorios
Documentación
Paquete neural
Documentación
neural no está en los repositorios, para instalar:
- Descargar de http://cran.r-project.org/src/contrib/Archive/neural/
- Instalar desde ubicación local
Paquete AMORE
Se puede instalar desde los repositorios. Es un paquete desarrollado en la Universidad de León, con abundante documentación.
Documentación
Wiki
Ejemplo
Próximamente
Paquete estadístico "R": Uso de Rcommander
Notas
">" : Ejecutar desde el intérprete de comandos de R
"~$": Ejecutar en la consola de Linux
Para instalar:
~$ sudo R
> install.packages("Rcmdr", dependencies=TRUE)
Para iniciar desde la consola:
~$ R
> library(Rcmdr)
Para iniciar desde la consola cuando ya hemos cerrado Rcmdr:
> Commander()
En Windows
Se puede arrancar también desde la consola de R, pero para evitar errores tenemos que abrir R en modoSDI . Si no elegimos esta opción al instalar R, podemos forzar este tipo de ejecución, añadiendo al comando del acceso directo que abre la aplicación la opción --r (Propiedades del acceso directo, caja de texto "Destino", quedaría algo así: "C:\Program Files\R\R-2.11.0\bin\Rgui.exe" --sdi
">" : Ejecutar desde el intérprete de comandos de R
"~$": Ejecutar en la consola de Linux
Para instalar:
~$ sudo R
> install.packages("Rcmdr", dependencies=TRUE)
Para iniciar desde la consola:
~$ R
> library(Rcmdr)
Para iniciar desde la consola cuando ya hemos cerrado Rcmdr:
> Commander()
En Windows
Se puede arrancar también desde la consola de R, pero para evitar errores tenemos que abrir R en modo
sábado, 1 de mayo de 2010
La Estadística, "sexy job" para los próximos 10 años, según Google
Esta es una de las afirmaciones que se hacen en un artículo del New York Times, nada menos que por el economista jefe de Google:
La propia IBM ha creado un grupo de análisis y optimización, contratanto a más de 200 analistas (matemáticos, estadísticos, ...) y tienen previsto contar con más de 4.000 en toda la compañía.
Pues nada, a ver si es verdad.
Enlaces:
http://www.nytimes.com/2009/08/06/technology/06stats.html?_r=2
http://www.newser.com/story/66223/todays-sexy-job-statistician.html
http://www.conversion-rate-experts.com/the-datarati/
NOTA: el artículo del NYT es de agosto del 2009, pero como los americanos nos llevan varios años de adelanto en estos temas, tiene absoluta vigencia.
"I keep saying that the sexy job in the next 10 years will be statisticians,” said Hal Varian, chief economist at Google. “And I’m not kidding".Según otro economista, del MIT, el reto es la capacidad de los humanos para usar, analizar, y dar sentido a los datos, en un mundo en que todo puede ser medido y monitorizado.
Los Datos son simplemente el material en bruto del conocimiento.
La propia IBM ha creado un grupo de análisis y optimización, contratanto a más de 200 analistas (matemáticos, estadísticos, ...) y tienen previsto contar con más de 4.000 en toda la compañía.
Pues nada, a ver si es verdad.
Enlaces:
http://www.nytimes.com/2009/08/06/technology/06stats.html?_r=2
http://www.newser.com/story/66223/todays-sexy-job-statistician.html
http://www.conversion-rate-experts.com/the-datarati/
NOTA: el artículo del NYT es de agosto del 2009, pero como los americanos nos llevan varios años de adelanto en estos temas, tiene absoluta vigencia.
Suscribirse a:
Entradas (Atom)