similar to: Bug? Avoiding options in subdirectories

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

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

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: >

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í.

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

Estimad en s eRRer en s, Tengo un txt que quiero importar a R. Pero no tiene un formato adecuado para usar cosas normales, como por ejemplo read.csv() El formato es algo así: time 1 col1 col2 col3 col4 dato dato dato dato dato dato dato dato dato dato dato dato dato dato dato dato dato dato dato dato end time 2 col1 col2 col3 col4 dato dato dato dato dato dato dato dato dato dato dato dato dato

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, 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

Estimados herreros, Sabéis si hay algún sitio donde pueda hacer una pregunta sobre estadística, en concreto un diseño experimental. Alguna lista, foros, etc. Ya se que este no es el sitio adecuado para preguntas de estadística, pero por si alguien se pica os la pongo a continuación. Dos tratamientos C y H. Se aplican a una serie de especies (p.e. 10) Se observa la respuesta de las especies a

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, 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. -

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 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",

Hi! Is there any way of having a greyed ("ghosted") text (i.e, 2006) in the background of a plot? I'm making a dynamic plot and would like to show the year of each time step as a big greyed text in the background. (the idea comes from Hans Rosling video: http://video.google.com/videoplay?docid=4237353244338529080&sourceid=searchfeed ) Thanks Agus -- Dr. Agustin Lobo Institut

Hi! Given a vector (or a factor within a df),i.e. v1 <- c(1,1,1,2,3,4,1,10,3) and a dictionary cbind(c(1,2,3),c(1001,1002,1003)) is there a function (on the same line than recode() in car) to get v2 as c(1001,1001,1001,1002,1003,4,1001,10,1003) ? I'm using myself a function based on match() since long ago (I think that thanks to advice by Prof. B. Ripley), but would like to know if there

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

Hi! Is it possible to set a file to which both commands and output would get automatically saved? I've tried with sink(), but only get the output. I mean something like a combined history and sink, as you get with File/Save to File.. in the windows GUI. Tis is done with diary filename in Matlab, and you can state diary on and diary off to control what is being saved to the file. Thanks Agus