similar to: Model averaging en R

Muchas gracias Carlos, la lógica es perfecta pero no se como identificar con código las lineas en blanco entre el bloque 2 y el bloque 3. Para de esta forma quedarme solo con el bloque 2. Tienes alguna idea? Muchas gracias. Joan 2016-07-28 17:00 GMT+02:00 Carlos Ortega <cof en qualityexcellence.es>: > Hola, > > Se me ocurre esta solución en pseudo-código...: > > >

Hola a todos, tengo 170 .csv donde tengo que eliminar las primeras 20 lineas (primer bloque) y luego todo un último bloque de datos (tercer bloque) que está separado por dos filas sin datos del segundo bloque (que es el que me interesa). El tercer bloque empieza en cada .csv en una linea diferente por lo tanto no se si puedo automatizar en R quedarme tan solo con la información a partir de la

Hola, tengo una matriz de especies donde cada celda tiene datos con caracteres (son parentesis). Cómo la puedo reordenar considerando tanto filas como columnas? Ejemplo: ,Specie 1, Specie2, Specie3 Specie1, NA, 3(1-4), 8(6-9) Specie2, 5(2-6),NA, 5(4-6) Specie3, 2(1-3), 10(5-15), NA Quiero: ,Specie 2, Specie3, Specie1 Specie2, NA, 5(4-6) , 5(2-6) Specie3, 10(5-15) ,NA,

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

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

Estimados usuaRios: deseo hacer una consulta, debo realizar una prueba no paramétrica a un conjunto de datos que posee una variable de respuesta con tres factores que poseen interacción, o al menos debo probar ello. Es decir, realizar una prueba Scheirer-Ray-Hare o de Contrastes de rangos alineados. Estuve trabajando el paquete ARTool aunque no me queda claro

Estimados usuarios de R: escribo porque tengo un problema con el paquete rgl que no puedo resolver desde hace un tiempo mi sessionInfo es. sessionInfo() R version 2.15.3 (2013-03-01) Platform: i686-pc-linux-gnu (32-bit) el problema específico es el siguiente: library(rgl) > demo(rgl) demo(rgl) ---- ~~~ Type <Return> to start : >

Hi, I have recently installed R-1.2.2 for windows (16MB RAM, P-166) and I getting the following message after processing my data (6 variables and 1200 observations): >Error: cannot allocate vector of size 4 Kb >In addition: Warning message: >Reached total allocation of 15Mb: see help(memory.size) Then, the program close. With the last version, 1.1.1 (I think) I didn't have this

Hola a todos: Alguién me podría explicar cuales son las diferencias entre las diferentes escalas de viridis en ggplot2, es decir en qué se diferencian scale-colour_viridis_d de scale_colour_viridis_c, scale_colour_viridis_b. Y cual es la diferencia entre scale-fiil-viridis_d de scale-fill_viridis_c y de scale_fill_viridis-b. He estado buscando la diferencia en internet y no he encontrado nada al

Hi All, How do I can create a new vector, i.e. 'seasonal' or 'trend' from the resultant seasonal or trend component of the Time.Series object produced by STL, and how I could superpose in the same graphic, i.e. original data and trend or seasonality? Thanks in advance! Antonio Antonio Rodr?guez Verdugo CICEM Agua del Pino Huelva Oceanography and Coastal Resources, PhD Program,

Estimada Comunidad, hice un Anova (SS III) usando la biblioteca CAR y necesito obtener los residuos del ajuste ... el argumento de Anova fue un modelo ajustado con glm(), como indica el siguiente codigo: inf.ninf.glm <- glm(CH ~ pa.+lam+ps.+age+area+expo+inf+P+ppacum1mes, data=hum[hum$CH < 45,]) inf.ninf.Aov <- Anova(inf.ninf.aov, type=3, error.estimate="pearson") mi

estimados, ajusté una superficie de respuesta con el paquete "rsm" de R y hago un grafico 2D de la superficie de respuesta con contour(), pero como el argumento de esta función es la superficie de respuesta ajustada y ésta se ha ajustado codificando las variables, de acuerdo a la recomendación de la documentación del paquete rsm, obtengo un gráfico con los ejes en unidades codificadas

Hi, How do I can make a lagged scatterplot of two variables: Yt (nao) versus Xt-h (mei) if they have the following structure: >series mei nao Jan 1950 -1.036 0.55 Feb 1950 -1.133 3.31 Mar 1950 -1.259 0.81 Apr 1950 -1.027 1.60 May 1950 -1.399 -1.73 Jun 1950 -1.366 1.26 Jul 1950 -1.300 -0.87 . . . I've tried with lag.plot but I don't understanf how to use it Thanks in

Hi, I have a data frame with 13 time series (same length). Two of them have several NA's values, so when I've tried to run: >acf(flow, na.action=na.omit) I get >Error in na.omit.ts(as.ts(x)) : time series contains internal NAs I've read through help files of acf, na.action, na.omit, and I don't understand how to manage NA's values. If I list the variable (eg: rain,

Hola a tod en s, tengo una duda que la comunidad R me puede ayudar. Estoy trabajando con gamlss, porque tengo una variable respuesta con valores entre 0 y 1 e incluidos estos. La distribución que utilizo com gamlss para este caso es "beta inflated" (Stasinopulos and Rigby 2007. Journal of Statistical Software 23(7)). El modelo que intento correr es: m1<-gamlss(Teleosteos ~

Buenas erreros!! Una cuestión de las que tengo ciertas dudas es saber en que se diferencian Spark y H2o, si son competencia, si valen para lo mismo o no.... Según lo poco que se, Spark es una manera de agilizar el Map-Reduce, y con la libreria MLlib, puedes hacer datamining de grandes datasheets, y si lo conectas con R o con Python, puedes usar ese lenguaje. H2O es una herramienta que nos

Tengo un modelo de regresión de Cox y quiero obtener el plot ajustado por una covariable (sexo) con la función ?ggadjustedcurves?, pero me da el siguiente error: > cox2 <- coxph(os ~ imc_25 + sexo.1, data = datos) > cox2 Call: coxph(formula = os ~ imc_25 + sexo.1, data = datos) coef exp(coef) se(coef) z p imc_25 -0.621 0.537 0.299 -2.08 0.038 sexo.1M 0.714

R múltiple Estoy pensando en un problema que tendré que solucionar pero aún no comencé a escribirlo, por lo tanto no hay código en R como para compartir, sin embargo no tengo idea de cómo realizarlo desde R. El planteo es el siguiente: Los datos son en una cantidad necesaria para que el procesamiento estadístico demore (minutos, horas). Supongamos dos variables (serían más), la A y la B,

Estimadas y estimados miembros de la comunidad. Junto con saludar, les escribo este mensaje porque estoy interesado en desarrollar un taller de estadística descriptiva con R para estudiantes de secundaria. El objetivo es que a través de un conjunto de datos puedan calcular medidas de tendencia central, dispersión, calcular correlaciones y hacer gráficos. Actualmente, he realizado un taller

Please help, I'm having this problem at my job and it really needs to get solved. I'm trying to setup samba 3.0 to be a member fileserver in my Windows 2000 active directory domain. I followed the instruction on this website: http://www.linux-sxs.org/networking/nt4dom_samba.html#win_sysreq net RPC join -W <domain> -U <domain user> works wbinfo -t works wbinfo -G works