similar to: stat_density2d de ggplot2

Marta, Puedes ver en el código source de ggplot stat_density si es que lo tienes descargado como sorce o puedes ir a GitHub https://github.com/hadley/ggplot2/blob/master/R/stat-density-2d.r y ver como calcula los datos. Por lo que veo usa la función de la librería MASS kde2d, para ver que hace esta escribe kde2d en la consola. Espero te sirva. Daniel El 18 de enero de 2014, 22:35, marta

Hello, I'd be very grateful for help with some ggplot2's stat_density2d issues. First issue is with data limits. xlim() and ylim() doesn't seem to work; instead, estimates (and plotting) seems to be constrained to range(x), range(y) no matter what i do. The documentation says i can pass in kde2d's parameters to ... but pussing kde2d's "lims" parameter achieves

Hola a tod@s, Hace poco en la lista se discutieron algunas aproximaciones para determinar la distribución de probabilidad que potencialmente podrían haber generado los datos (ver [1]) y una de ellas fué el AIC. Haciendo uso del programa enviado por Pablo Verde (ver [1]), estos son los resultados para mis datos: weibull , AIC = 69839.44 exponential , AIC = 79488.77 gaussian , AIC = 69413.03

Carlos, mi intencion es entender el codigo: https://github.com/woobe/blenditbayes/blob/master/2013-06-street-level-crime/analysis.R Para luego hacer mi trabajo, algo similar a esto, pero adaptado a mi región. Si te fijas en las lineas 141-161, esta esto: ## Convert the base map into a ggplot object ## All added Cartesian coordinates to enable more geom options later on map.ggmap <-

Cordial saludo. Estoy trabajando con una densidad conjunta algo elaborada (suma y cociente de dos variables aleatorias positivas usando la cópula Frank) y requiero encontrar la distribución marginal de cada una de las variables. A manera de ejemplo (no es el caso específico, que tiene mucho más código), si partimos de la distribución exponencial bivariada #----- theta<-0.5

Estimados, tengo que hacer un calculo muy simple, pero con una BD mas o menos grande (250mil filas x 500 columnas) ... esta BD esta en formato de SPSS y la importo asi: library(foreign) bdr <- read.spss("CASEN_2013_MN_B_Principal.sav", use.value.labels=FALSE, to.data.frame=TRUE) luego, quiero transformar el DF en un data.tale pues el calculo requiere obtener promedios de

Muchísimas gracias Carlos. La he probado simulando datos con solapamiento conocido y va como un tiro! Saludos Angel Sent from Samsung tabletCarlos Pérez Carmona <carlos.perez@uam.es> wrote:Hola Ángel, Nosotros creamos recientemente una función para calcular el overlap entre distintos caracteres de especies de plantas o animales. La función (trova) puedes encontrarla aquí:

Hola a todos. Estoy interesado en obtener el grado de solapamiento de dos parametros obtenidos a partir de una MCMC. La gente que he consultado realiza inspecciones visuales (graficando ambas funciones de densidad  de los valores obtenidos y viendo si hay mucho o poco solapamiento), pero me gustaría incluirlo en una serie de simulaciones con lo que no podría realizarlas. Conoceia alguna manera

Estimados amigos y expertos del R, Les escribo para hacerles una pregunta que parece un poco sencilla pero me ha costado mucho. Estoy tratando de graficar los datos correspondientes a tasas de consumo de algunos organismos cuando están expuestos o no expuestos a un químico sobre tres tipos de presa. Por definición, debía ajustar los datos a un glm con distribución gama. Las gráficas pueden ser

Hola a tod en s, tengo una duda que se relaciona con alternativas al uso de variables globales. En principio, si se quiere usar un generador de v.a con la librería Runuran sólo se permite definir las funciones de densidad (o el núcleo de las mismas) con funciones con un único argumento en (x). Sin embargo, necesito pasar a las funciones más argumentos que van cambiando en las iteraciones de

Buenos dias R-help-es, Estoy interesado en estimar una curva dosis-respuesta para un conjunto de datos y para ello, estoy utilizando la libreria "drm". Hasta ahi todo bien. Me gustaria automatizar algunas cosas y el primer paso para ello es la estimacion del modelo. Si la estimacion funciona, todo lo demas funciona; de lo contrario, todo fallara. Tengo algunas lineas que mitigan un

Hola Alguien me podría decir como hacer una grafica del tipo persp() de la densidad de una distribucion normal bivariante estandarzada con correlacion 0.5?? gracias [[alternative HTML version deleted]]

Hola, Estoy intentando realizar un ajuste a un modelo setar con R y me surgen problemas a la hora de representar el modelo setar sin la constante. El código es el siguiente: # Estimacion TAR(2;2,2) con delta 1: mod.setar <- setar(x, m = 2, mL = 2, mH = 2, thDelay = 1) mod.setar summary(mod.setar) mod.setarc<-setar(x, m=2, mL=2, mH=2, thDelay=1, include=c("none")) ##eliminamos

Para este ejemplo especifico: Deseo saber si tmp genera un objecto ggplot > tmp <- ggplot(mtcars, aes(x = hp, y=mpg)) + geom_point() El objecto o variable tmp, devuelve una lista y un gráfico, la clase dice ser ggplot, >class(tmp) [1] "gg" "ggplot" > typeof(tmp) [1] "list" El 30 de enero de 2014, 13:00, daniel

Estimados compañeros Escribo para anunciar en la lista que el martes de la semana que viene (8 de Octubre) a la *13:00 h. *tendrá lugar en el *Instituto de Ciencias Matemáticas del CSIC-* *ICMAT (Universidad Autónoma, Cantoblanco,ver ubicacion <http://www.icmat.es/facilities/howtoarrive>) (Aula Gris 1), Madrid* un seminario sobre "Estimacion para areas pequeñas" ( "Small

Hola, Estoy tratando de leer una base de datos: tengo 39 columnas, de las cuales 38 son variables y una es el nombre de las empresas. Por lo que escribo el siguiente codigo: Variables <- read.csv("C:/Users/usuario/Documents/variables2.csv", header=TRUE, sep=";", comment.char="" , colClasses=c(Empresas="character",

si se tiene un data.table (DT), supongamos de 100 filas por 3 columnas de datos numericos, como puedo hacer para obtener el correspondiente valor de la columna 1 si busco, por ejemplo, el maximo de la columna 2 agrupado por la columna 3 ? para buscar el maximo de la columna 2 escribo. DT[ , max(c2), by=c3 ] muchas gracias, saludos, eric. -- Forest Engineer Master in Environmental and

Hola amigos, Supongamos que se quiere ejecutar un modelo no lineal con nls. Pensemos en el ejemplo de la ayuda: DNase1 <- subset(DNase, Run == 1) fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1) summary(fm1DNase1) Aquí se está modelando la densidad óptica de un ensayo relacionada de forma no lineal (logística) con (el logaritmo) de la concentración de una proteína.

Hola a todos, estoy procesando unos datos con las siguientes intrucciones que creo son del paquete stats: cloud(Imp89 ~ ImpTot89 * NDA | Edad, data = mtif, +       screen = list(x = -90, y = 45), distance = .4, zoom = .7) y obtengo el error: Error en eval(substitute(groups), data, environment(formula)) :  argumento ''envir'' númerico no tiene longitud uno La verdad es que no se

Hello, I have this code to plot a certain normal distribution and represent the pnorm value for a certain "x": x<-300 xx <- seq(2.5,7.5, by=0.1) yy <- dnorm(xx,5.01,0.77) d<-signif(pnorm(log(x), 5.01,0.77),4) xpts <- round(exp(0:8)) par(bg = "antiquewhite") plot(xx,yy, type="l", col="blue", lwd=2, xaxt="n",