similar to: Interpretación de salida de un GLM

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

Dear Consultant I've done linear regression successfully on R a few times before. But this time it keeps telling me:- "Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases" The model is:- fm1 <- lm(TS.CM ~ AGE + SEX + HFE.Y.01 + TFC2B.01 + HFE.Y.01*TFC2B.01, data = IRONresults, subset = DIAG2.1D == 0) summary (fm1) TS.CM is a

Hola a todos. Me encuentro realizando un gráfico en R y en los ejes quisiera que los valores aparezcan con separador de miles. Es decir... 1000000 1,000,000 Gracias... -- Paula Altamar Estratificación-DIG-DANE Cel: 3005137541 [[alternative HTML version deleted]]

Buenas. Abajo pongo la salida de un modelo de cox , dónde he estratificado por una variable de país (Countryb) y por otra (Q6). Además hay interacción entre la variable mobilityPDurG2 (es una variable 0,1, y 0 es la categoría de referencia) país. La categoría de referencia para país es "united kingdom". Mi duda surge si quiero calcular el hazard ratio para los que tienen un 1

Full_Name: Harald Weedon-Fekj?r Version: 2.0 OS: Windows Submission from: (NULL) (80.202.87.143) Dear R core Thanks again for a superb program. This is really one of the rare instances I find a bug. The bug is easy to bypass, but can be quite confusing; "Rcmd check" seems to fail without any good error report if the package contains a certain amount of R code and is run from a

Buenas, ¿Alguien ha hecho algun curso sobre bigdata que merezca la pena? La verdad es quequiero adquirir competencias sobre el manejo de grandes volumenes de información, hadoop, spark, etc y no se muy bien que estudios cursar. He estado mirando, pero veo que las escuelas de negocio lo orientan a big data + business analytics y no se si pueden pasar muy por encim adel tema o estar centrados en la

Hi all, After running this code (attaches is the input file): dat <- read.table("gene_prob.txt", sep = "\t") n <- length(dat$V1) print(n) print(dat$V1) I get this print out. ...... [8541] LOC552889 GPR15 SLC2A11 GRIP2 SGEF [8546] PIK3IP1 RPS27 AQP7 8548 Levels: 3.8-1 A2M A4GALT A4GNT AAAS AAK1 AAMP AANAT AARSD1 AASS ...

Full_Name: Harald Fekjær Version: 1.22 OS: Windows Submission from: (NULL) (158.36.132.220) Dear R developers Actually, this is not a bug, but a small remark about a way R is working that makes trouble in porting the code to S. If I run the following i R: --- midl <- 4 attr(midl,"Object created") <- date() dump("midl","midl.R") --- I get a file like: ---

When I try to adjust a mixed model with random effects I can make this order without problem > lm.FA<-lme(absFA~trait*condition,random=~1|individual) But if I try to fit a model in which the response (absFA) is not the same in all individuals at different levels of "trait" factor , but varies randomly from each. That is, this order >

Buenas Tengo un dataset como el que viene: datos1 datos2 grupo 0 122 A 1 128 A 2 129 A 3 241 A 0 201 A 1 185 A 2 174 A 3 125 A ... .... .... ... .... .... ...

Gracias Beatriz, efectivamente, lo que indicas en tu ejemplo es lo que obtengo al final de mi proceso. En todo caso pruebo tu opción con mis datos, si es como imagino seguro que es más rápida de montar y más elegante que tratar el resultado de un sink() (recuerdo que en su momento lo intenté con fotify pero no supe bien como atacarlo, no conocía el enlace que me mandas, por tanto lo pruebo de

Me funciono perfecto, muchas gracias Carlos ... si no es mucha la molestia, alguien podria explicarme la diferencia tecnica entre ambos signos, y porque == produce el resultado que produce ? en que casos deberia usar == y en que casos %in% ... cual es el significado de cada uno ? Disculpen que no me conforme solo con obtener el resultado adecuado, es solo que no me suena logico que == no

Estimados, tengo el siguiente problema: estoy trabajando con un data.frame que tiene 892 filas y 9 columnas, el cual adjunto para que puedan probar. Cuando consulto por las filas que cumplen una cierta condicion, por ejemplo, ... graph[graph$NCar==c("160s"),] obtengo el siguiente conjunto de datos: sol con dia NCar NIns isom area rep tipo 1 con 0.001 1 160s 0 s

Tengo un vector de fechas en formato R: > fechas <- c(?2016-07-15?, ?2016-10-12?, ?2017-02-11?) quiero que las fechas de 2017 (la tercera del vector del ejemplo) pasen a ser NA > fechas2 <- ifelse(fechas > ?2016-12-31?, NA, fechas) El problema que me encuentro es que el resultado (fechas2) es un vector numérico y no un vector de fechas. No encuentro la solución. Se agradecen

Muchas fracias por la rapida contastacion Javier. Si, ya me he dado cuenta que a veces es mas limpiar y formatear que calcular. Lo que yo necesito hacer no es crear un data.frame, es crear una lista por cada grupo/condicion de muestras. Y queria hacerlo en un for loop, para que no tuviese que indicar cuantas listas necesito si no que cogiese esa informacion del usuario (y) Creo que no me

Hola, que tal? Sin saber nada mas de los datos, por lo que dices, mi aprioximacion seria usar la funcion split, e.g.: splittedData <- split(myData, "variablegrupo/condicion") Esto te crearia una lista con tantos objectos como condiciones tengas. Saludos On 6/27/24 10:15, Leticia G Leon wrote: > Muchas fracias por la rapida contastacion Javier. > > Si, ya me he dado

Hola Neo, esto es lo mejor que me ha salido. A ver si por lo menos te ayuda. Un saludo colores<-c("black", "green", "red", "steelblue", "purple") #creamos el vector de 5 colores que usaremos para unificar luego la leyenda y el grafico xyplot(V5 ~ dia | sol, groups=con, layout=c(1, 3), type= "l", pch=1,

I''ve been stumbling on the tut at oriontransfer (http://wiki.oriontransfer.co.nz/main/show/SortableTree) I''ve got it working in a one dimensional list and I''ve got the drag and drop working on children, but as can''t get it to save the tree. Here''s my view: <h2>Menu Prototype 2</h2> <%= flash[:notice] %> <ul

Hi, I have encountered this error when attempting a One-way Repeated-measure ANOVA with my data. I have read the "Anova in car: SSPE apparently deficient rank" thread by I'm not sure the within-subject interaction has more degrees of freedom than subjects in my case. I have prepared the following testing script: rm(list = ls())