similar to: Leer un txt a trozos

Estimad en s errer en s He hecho este modelo glm m1.pile<-glm(ger~tem+pot+time+I(tem^2)+I(tem^2):pot ,family="binomial" ,data=long.PILE ) Que nos da la probabilidad de germinación de una semilla en función de tem (Temperatura), pot (Humedad del suelo) y time (Tiempo que la semilla pasa en esas condiciones). Ahora quiero, para diferentes tem, pot

Un ejemplo con un modelo más simple: He especificado este modelo: >formula(m2.pile) ger ~ tem + pot + time Si hago predict me da: >predict(m2.pile,newdata=data.frame(tem=25,pot=0,time=3),type="response") 0.08243262 Extraigo los coeficientes: > coef(m2.pile) (Intercept) tem pot time -1.89521331 -0.02303313 4.74499714 0.02043222 Ahora calculo la

Desde el terminal de Ubuntu, he conseguido descargar los datos usando esto: DESCARGAR LISTA DE ESTACIONES wget --method GET --header 'cache-control: no-cache' --no-check-certificate --output-document - https://opendata.aemet.es/opendata/api/valores/climatologicos/inventarioestaciones/todasestaciones/?api_key=*vuestra api key* DESCARGAR DATOS HISTÓRICOS Diciembre 2017, ALICANTE wget

Hola, No tengo mucha idea de esto y hace tiempo que no lo uso. Unas ideas: - ¿Puedes hacer un summary(prueba) para ver si los datos que tienes dentro del raster son los que te está representado? Eso te permitirá saber si el problema esta en raster() o en plot() - Puede que al hacer plot() te haga categorías en lugar de asignar un color para cada valor, igual el parámetro col de plot() te ayuda. -

Muchas gracias, Probaré eso también y ya os cuento. Jaume. El lun., 17 feb. 2020 a las 22:10, Javier Marcuzzi (< javier.ruben.marcuzzi en gmail.com>) escribió: > Estimado Jaume Tormo > > En lo personal yo utilizo un enfoque como el que comenta Carlos Ortega, se > me ocurre que posiblemente funcione si a su código le coloca algo de > formato, me refiero a esta forma: >

Carlos, muchas gracias, voy a probarlo. Pero me sigue intrigando por que no puedo ponerlo como elemento de un vector... Misterios del R. SI lo averiguo os lo digo. Jaume. El sáb., 15 feb. 2020 a las 19:08, Carlos Ortega (<cof en qualityexcellence.es>) escribió: > Hola, > > Una alternativa que te puede ayudar es enfocar el problema de esta otra > forma. > Puedes ir guardando

Hola eRReros, Necesito extraer datos de un dataframe según uno de los factores. La cosa sería algo así: df.nuevo <- iris[iris$Species=="setosa"|iris$Species=="virginica",] Que me crea un df nuevo solo con las filas en que Species es setosa o virginica. El problema es que mi factor de interés tiene 22 niveles y a veces he de seleccionar 13, 10, 8 niveles o números así.

Hola eRRer en s, He creado un modelo mixto (adjunto) con lme4 He dibujado la predicción usando: df <- ggpredict(m.glob.pre.anu , terms = c("pre", "area_m2", "sum.vs.win") ) plot(df) Luego he arreglado el gráfico usando: df <- ggpredict(m.glob.pre.anu , terms = c("pre", "area_m2",

Hola eRReros. Os lo explico con un ejemplo: # Cargamos los paquetes y el ejemplo install.packages("SPEI") library("SPEI") data(wichita) # los primeros 12 datos head(wichita,12) # mi subset de los primeros 12 datos meu<-wichita[c(1:12),] meu # como veis los valores de TMED son iguales en ambos dataframes. # ahora viene el problema # calculamos

Hola, Estoy encontrando un problema al intentar poner un objeto Dates en un vector o dataframe. Mi ejemplo # preliminares install.packages( lubridate ); library( lubridate ) v <- c(0, 0, 0) original<-c(2019,308,1700, 25) # c(año, día del año, hora, temperatura) esto sale así de un sensor de temperatura # convertimos los datos originales en algo que entienda R a <-

Hola, en la misma definici?n de la funci?n: # P15. Precipitation Seasonality(Coefficient of Variation) # the "1 +" is to avoid strange CVs for areas where mean rainfaill is < 1) p[,15] <- apply(prec+1, 1, cv) Un saludo, Jorge On Martes, 19 de Junio de 2018 13:07:27 Marcelino de la Cruz Rot escribi?: > Hola Jaume: > > Si miras el c?digo de biovars() ver?s que la

Hola eRReros, Estoy trabajando con datos geográficos unos en formato raster y otro en formato shape. Lo que quiero es hacer las estadísticas (sobretodo la media) de los píxeles del raster que caen dentro de algunos de los polígonos del shape. Lo he intentado con sextante pero se flipa con los valores perdidos, nuestros queridos NA. Y he pensado: estando el paquete raster en R seguro que se puede

Estimad en s, Junto con saludar, quería saber si alguien me podría ayudar con este gráfico, tengo un análisis de datos del año 2016 de 3 variables con dos dimensiones (muestra objetivo y muestra lograda), y quería agregar la comparación con el año 2015, me gustaría agregar un gráfico de linea y puntos para el año 2015, pero no he podido lograrlo , la idea es generar un gráfico de esta forma.

Hola a todos, Estoy trabajando con modelos lineales generalizados (GLM), en particular con las funciones glm y anova.glm de R. Tengo una pregunta que es más bien técnica en el sentido estadístico (y un poco Off Topic y probablemente un tanto naive). No tengo claro si es correcto decir análisis de varianza o de devianza al utilizar GLMs, veo que en el help de anova.glm dice análisis de

Estimados erreros, Estoy intentando entender como calcula el paquete dismo ( https://cran.r-project.org/web/packages/dismo/index.html) un coeficiente de variación. Os pongo un ejemplo: tmin <- c(10,12,14,16,18,20,22,21,19,17,15,12) # temperatura mínima media mensual de un año tmax <- tmin + 5 # temperatura máxima media mensual de un año prec <- c(0,2,10,30,80,160,80,20,40,60,20,0)

Hola a todos. Quería saber si conocéis alguna librería que haga esto: SPSS permite enlazarse con PowerPoint de manera que si quiero modificar un gráfico del ppt, con solo pinchar en el gráfico me aparece los datos del spss y al modificar estos cambia la gráfica directamente sin tener que hacer la gráfica de nuevo. Hay algo parecido en R???? Saludos José Miguel Contreras [[alternative

Estimados Estoy realizando algunos cambios, migrando desde windows 7 64 bit a linux SUSE 12.02. No encontré inconvenientes hasta que: Error en plot.new() : márgenes de figura muy grandes Este error se da en Rstudio, ¿alguna idea? En windows sabía que usaba la consola en varias ventanas, pero en linux ¿...?. ¿Hay alguna configuración que desconozco? Desde ya muchas gracias. Javier Marcuzzi

Hola a todos. Un problema dos pasos. Quiero desde RStudio trabajar con dos versiones de R (pongamos R2.13 y R2.15) desde Linux, Ubuntu 10.04 (lucid). Ya tengo instalada la R.2.15.1. Por defecto trabajo con ella. Pasos que creo he de dar: 1º. Instalar dos versiones de R en la misma maquina. ¿Como se hace?. Me estoy metiendo con el PATH y tal... hay algún protocolo clarito por ahí? 2º.

Hola a todos. En linux uso rkward o gedit con el plugin Rgedit, pero me gustaría saber qué editores se pueden usar en windows. He probado Tinn-R , ¿qué usáis vosotros? Gracias

Hola a todos: Estoy intentado realizar un modelo senoidal de unos datos de radiación solar con el fin de afrontar el relleno de la serie y aplicar la prueba de Thom para verificar su homogeneidad [0]. De momento me encuentro con los siguientes problemas: 1- ¿Existe la prueba de Thom en R? ¿O debo crearme mi propia función? 2- Para la realización del modelo senoidal estoy siguiendo los pasos