similar to: Using RasterBricks

Hola Jorge, > My idea was to use 'stackApply': Para esto que necesitas es más sencillo usar calc. Si quieres usar stackApply, debes indicar rep(1, 3) como índices en lugar de sólo 1. Así, la función será (los ... son necesarios en stackApply): myFUN <- function(y, ...) { x <- seq_along(y) trapz(x, y) } Y aplicada con calc: FINAL <- calc(s, myFUN) o con

Gracias a todos, finalmente instalando primero el libxml2-dev y luego el dplR, funciono ----- Mensaje original ----- De: r-help-es-request en r-project.org Para: r-help-es en r-project.org Enviados: Miércoles, 17 de Julio 2013 12:59:47 Asunto: Resumen de R-help-es, Vol 53, Envío 26 Envíe los mensajes para la lista R-help-es a r-help-es en r-project.org Para subscribirse o anular su

I tried to update my package {pracma} on CRAN from 2.4.2 (2022-09-21) to version 2.4.4 (2023-11-08). This package reverse depends / imports / suggests on 350 packages on CRAN and 25 packages on Bioconductor. The only changes are small corrections on some help files, a new function for stereographic projection, and `gcd` and `Lcm` require integer inputs now (these functions are not used in the

Necesitaría ayuda con una cuestión. Pongamos por ejemplo estos valores aleatorios para generar un RasterBrick s de 3 capas. library(raster) r <- raster(ncol=10, nrow=10) r[]=1:ncell(r) s <- brick(r,r,r) s <- s * 1:3 Habría alguna forma sin hacer un bucle for, que me devolviese un vector con los valores medios de cada una de las 3 capas. mean_s<-NULL for (i in 1:3) {

  Hola:   Perdonad por la insistencia, pero me gustaría saber si cualquier proyecto en R debe afrontarse desde el uso de clases, o si sólo procede ante determinados casos. De ser así, ¿cuáles serían?.   Gracias.   Eva [[alternative HTML version deleted]]

Hello, I need to integrate the absolute difference between two lines measured on different points. # For example : x <- seq(0, 1, 1/100) f_x <- runif(101) + x y <- seq(0, 1, 1/23) f_y <- runif(24) + (1 - y) plot(x, f_x, type="l") lines(y, f_y) Then I would like to compute Integral( | f_x - f_y | )dx. (This is not the same as | Integral(f_x)dx - Integral(f_y)dx |.)

I am trying to rasterize a SpatialPolygon object by a RasterBrick object. The documentation of the raster::rasterize function explicitly says this is allowed. Here's what I am doing # load the raster package library("raster") # create a raster brick object using the example from the brick function documentation b <- brick(system.file("external/rlogo.grd",

Hola He hecho el típico mapa con gadm y spplot: library(sp) ES <- local(get(load(url("http://biogeo.ucdavis.edu/data/gadm2/R/ESP_adm2.RData ")))) acerco las canarias: A <- which(ES en data$NAME_1 == "Islas Canarias") L <- length(ES en polygons[[A]]@Polygons) L <- length(ES en polygons[[A]]@Polygons) for (i in 1:L){ ES en polygons[[A]]@Polygons[[i]]@coords

Hola a todos, a<-c(1,2,3,4,5) b<-c(1.3,2.4,5,6,8) etiquetas<-c("etiq1","etiq2","etiq3","etiq4","etiq5") c<-rbind(a,b) Me gustaría hacer un gráfico en el que se represente en el eje Y las series a y b que están unidas en c y en el eje X las etiquetas ¿cómo se podría hacer? Espero haberme explicado. Sería el equivalente a la figura de

Buenas tardes a todos, Estoy interesado en realizar varios gráficos y luego guardarlos individualmente como PDF. Hasta ahí todo bien. Sin embargo, el archivo resultante tiene aproximadamente 20MB. Podría alguien por favor indicarme cómo reducir su tamaño? A continuación un ejemplo de lo que intento hacer y mi sessionInfo(). Estoy usando Mac OS X. Muchas gracias, Jorge Ivan Velez # R code

Hola, Si no te entiendo mal, necesitas la conversión a matriz como paso intermedio para agregar por meses (o cualquier unidad temporal). Si es así, te aconsejo que uses el paquete zoo y el método aggregate para la clase zoo. Por ejemplo: library(zoo) dats <- data.frame(date=seq(as.Date('2012-01-01'), by='day', length=3*365), A=rnorm(3*365),

Hola, Yo lo resuelvo usando grid.rect. Tienes una posible solución aquí: https://github.com/oscarperpinan/spacetime-vis/blob/master/choropleth.R#L216 (que es una versión actualizada de este artículo https://procomun.wordpress.com/2012/02/18/maps_with_r_1/) Saludos. Oscar. ----------------------------------------------------------------- Oscar Perpiñán Lamigueiro Dpto. Ing. Eléctrica,

Hola, Estoy intentando dibujar una curva, pero que esté divivida por la mitad en dos colores: plot( seq(0, 100,by=1), seq(0,1,by=0.01), bty = "n",col = "white") x <- c(50: 100) y <- 20*x^-1 polygon( c(50,50, c(52:100),100 ) , c(1,y[1:50],1), col = "steelblue2", border = NA) x <- c( 0: 50 ) y <- 20*x^-1 polygon( c(c(21:50), 50 ) , c(y[21:50],1),

Hola ¿Alguno tiene un código para sumarizar medidas? my.df <- data.frame( "ID" = c( rep("A", 3 ), rep("B",2),"C", rep("D",3) ), "Obs.1" = rnorm( 9,0,1 ), "Obs.2" = rnorm( 9,0,3 ) ) Algo que quede mas bonito que esto: by( my.df, my.df$ID, mean ) Gracias un saludo! -- Patricia García González [[alternative HTML

Is there packaged code to convert geographical coordinates (e.g., longitude, latitude, elevation) to Cartesian coordinates in 3-space? I can see how to do this using 1. a spherical-to-Cartesian conversion like pracma::sph2cart(tpr) http://cran.r-project.org/web/packages/pracma/ 2. a geographical-to-spherical conversion. This seems to involve (in roughly increasing order of difficulty or

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120614/518ed68f/attachment.html>

Hola, ¿qué tal? Escribo en esta lista para recabar la ayuda de socios de la asociación y la comunidad de usuarios de R. Resulta que el nombre de la asociación que se eligió originariamente, "Comunidad de Usuarios de R", ha sido rechazada por el registro español de asociaciones a instancias del Ministerio de Sanidad, Servicios Sociales e Igualdad (¡nada menos!) por incumplir el artículo

Estimada comunidad, imagino que debe ser simple, pero llevo ya varias horas buscando como hacer para que el eje Y de un grafico hecho con lattice use solo enteros. Los valores que se grafican son enteros entre 0 y 3, pero el grafico incluye 0.5, 1.5 ... Alguna idea ? ... ya lei la ayuda de xyplot pero no veo donde podria encontrarse la opcion. Gracias, Slds, eric. -- Forest Engineer

Dear all, I am a new user to R and I am using pracma and nloptr libraries to minimize a numerical integration subject to a single constraint . The integrand itself is somehow a complicated function of x and y that is computed through several steps. i formulated the integrand in a separate function called f which is a function of x &y. I want to find the optimal value of x such that the

Hola buenos días: Hace algún tiempo mandé una duda que sobre fichero tipo shapefile que publica el Instituto Geográfico Nacional de España y que detallo en el word Duda_sobre_shp. El tema es que al final no obtuve respuesta (es cierto que lo mandé en Agosto y en esa fechas entiendo que sea normal, por tanto por recomendación de un conocido vuelvo de nuevo a remitir la cuestión) Como no podía