Resolución del juego de modelos con R

Hace mucho planteé un juego de identificación de modelos con R y ya se me había olvidado daros la solución. Pensando en el Grupo de Usuarios de R y en hacer algo parecido en una presentación recordé que había que solucionar el ejercicio. Lo primero es la creación de los datos, se me ocurrió una función sencilla y una nube de puntos alrededor de ella: #Variable independiente indep = runif(500,100,500) #Función para crear la variable dependiente foo = function(x){ mean(x)*(1-sin(-0.006042*x)) } dep = sapply(indep,foo) dep=dep+(runif(length(dep),-100,100)) datos = data.frame(cbind(indep,dep)) plot(datos) ...

11 de mayo de 2016 · rvaquerizo

Reunión del Grupo de Usuarios de R de Madrid 12/05/2016

Nueva reunión del Grupo de Usuarios de R de Madriden MediaLab Prado. No está publicado en la web por temas de cerrar definitivamente las agendas. De momento las propuestas son: • Iñaki Úcar: «Presentación de simmer (Discrete Event Simulation in R)». Iñaki nos presentará simmer. Un paquete que ha creado para poder realizar análisis y simulaciones de eventos discretos (DES) en R. Con este paquete se pueden realizar simulaciones para el estudio de la evolución de sistemas teniendo en cuenta diferentes tipos de gestión de colas. ...

7 de mayo de 2016 · rvaquerizo

Valor atípico o pocos registros. Animación con R

¿Cómo influye un solo punto en una recta de regresión? Evidentemente cuanto menos observaciones tengo más puede «descolocar» la recta de regresión. Sin embargo, cuantos más puntos tengo más complicado es encontrar ese punto con una recta de regresión, sin analizar los residuos podríamos hasta pasarlo por alto, aunque puede ser que nos interese ese punto. El código de R que genera la animación es: ```r library(animation) saveGIF( for (i in c(100,50,25,10,5,1)){ x <- seq(-500,500, by = i ) y=sin(x)+x/100 y[10]=y[10]+10 plot (y,x,main=paste("Regresión lineal con ",1000/i," observaciones")) reg <- lm(y~x) points( fitted.values(reg),x, type="l", col="red", lwd=2)}, interval = .85, ,movie.name="/Users/raulvaquerizo/Desktop/R/animaciones/outlier.gif")

24 de abril de 2016 · rvaquerizo

Reunión del Grupo de Usuarios de R de Madrid 14/04/2016

Nueva reunión del Grupo de Usuarios de R de Madrid este jueves día 14 de abril de 2016, para más detalles: http://madrid.r-es.org/34-jueves-14-de-abril-2016/ La agenda de la reunión será: Pedro Concejero: “Análisis de datos de accidentes de tráfico – DGT“. Miguel Ángel Gómez: “Análisis de Sentimiento – Aplicación a noticias Financieras” Carlos J. Gil Bellosta: “Posibilidades de interacción entre R y Python a través del nuevo paquete feather” Si el tiempo me lo permite por allí estaré. Saludos

10 de abril de 2016 · rvaquerizo

Reunión del grupo de usuarios de R de Madrid

Nueva reunión del Grupo de Usuarios de R de Madrid este jueves 17 de marzo, para más información: http://madrid.r-es.org/33-jueves-17-de-marzo-2016/ Con Python como gran protagonista de este evento. Por allí nos veremos.

15 de marzo de 2016 · rvaquerizo

Reunión del Grupo de Usuarios de R de Madrid. 14 de enero de 2016

De nuevo reunión del Grupo de Usuarios de R de Madrid. La agenda será: Francisco Rodriguez: “Introducción a las Redes Neuronales con R”. Carlos Ortega: “Calidad del Aire de Madrid”. Karel Rodriguez: “Price Sensitivity Metr (PSM) in R. Para más información http://madrid.r-es.org/jueves-14-de-enero-2016/ Por motivos de trabajo no podré asistir a la reunión.

12 de enero de 2016 · rvaquerizo

TwitteR con R. El hashtag `#rstats`

El objetivo de la entrada es empezar a analizar tweets con R y que mejor comienzo que usar el #rstats para ver usuarios que más lo utilizan. [no me gusta conjugar el verbo retwitear] El primer paso es crear una app con Twitter, para ello nos dirigimos a <https://apps.twitter.com/> y creamos una aplicación. Crearla es muy sencillo, sólo necesitamos una descripción y un nombre. La aplicación será la que permitirá a R interaccionar con Twitter mediante OAuth y para realizar esta interacción entre la aplicación y nuestra sesión de R es imprescindible: ...

6 de enero de 2016 · rvaquerizo

Creacion de ranking con R

Hasta la fecha si necesitaba crear un ranking o un orden con R realizaba la tarea del siguiente modo: nombres <-c("grupo_1","grupo_2") grupo <-sample( nombres, 10, replace=TRUE, prob=c( 0.5, 0.5) ) dataset <- data.frame(grupo) dataset$importes <- runif(10,100,30000) #Creación del ranking de las variables agrupadas dataset$ranking = ave(dataset$importes,dataset$grupo, FUN= function(x) rank(x, ties.method = "first")) Es una agrupación de factores a la que asignamos el orden con rank, con ties.method=»first» esta agrupación se lleva a cabo desde el primer nivel del factor. El resultado se puede comprobar haciendo: ...

20 de diciembre de 2015 · rvaquerizo

Reunión del Grupo de Usuarios de R de Madrid. 12 de noviembre

Este jueves día 12 de noviembre nos volvemos a juntar el Grupo de Usuarios de R de Madrid en el MediaLab Prado. La reunión comenzará a las 19 horas y además de los temas administrativos habituales veremos: – Mariano Rico: “La Web de Datos desde R”. – Rafael Cobo, María Allanegui, Rafael Bermúdez Miguez: “Jugando con time series: de raw app”. – Pedro Concejero: “Análisis de Supervivencia para rotación de clientes”. ...

10 de noviembre de 2015 · rvaquerizo

Mapa de Argentina con R

Un lector necesita realizar un mapa de Argentina con R. El primer paso es descargar el mapa en formato R de la página web de siempre: http://www.gadm.org/country seleccionamos Argentina y el formato en R. Podéis descargar en otros formatos y trabajar con R, pero eso lo contaré otro día. Para ilustrar el ejemplo me he descargado el mapa de nivel 2, es decir, a nivel de Estado argentino. Una vez descargado el mapa empleamos el código de siempre: ...

9 de julio de 2015 · rvaquerizo