search for: ajustado

...dos, He realizado un dredge (para obtener todos los modelos GAM posibles a parir de un full model), luego he seleccionado un confidence set (los modelos que no se diferencian en 2 en AIC) y he hecho un model averaging con ese confidence set. Ahora me gustaría aplicar ese modelo "average" ajustado sobre otro set de datos pero no se como especificar en R que use el mismo modelo average en otro set de datos. tanto el dredge como el model average lo he realizado con el paquete MuMIn Alguien me puede ayudar? Joan -- *Joan Giménez Verdugo* *PhD Student* *Severo Ochoa* Estación Biológica de...

Una consulta desde un Uruguay muy lluvioso, como se realiza una prueba Chi-cuadrado ajustando por otra variable. Muchas gracias! Saludos! Sebastián. [[alternative HTML version deleted]]

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 ... la pregunta es: hay alguna forma de obtener el gráfico de superficie de respuesta con las unidades originales habiendo ajustado la superficie de...

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 pregunta ahora, es como extraigo los residuos desde inf.ninf.Aov...

...ia). El problema que tengo es que los gamlss no toman el efecto aleatorio como random(Familia/Genero) o (1|Familia/Genero), sino como: m2<-gamlss(Teleosteos ~ cs(LT_max)+ random(Familia,lambda=T)+random(Familia:Genero,lambda=T),family=BEINF, data = datos) El modelo con el código en amarillo es ajustado con gamlss sin problemas. ¿Es correcto la forma randon(A)+random(A:B) para considerar un efecto aleatorio anidado como (1|A/B)? Muchas gracias de antemano y saludos cordiales a toda la comunidad, Santiago. -- Santiago A. Barbini *Laboratorio de Ictiología* Instituto de Investigaciones Marinas y...

Hola Eric, A la que me refiero es a la que se utiliza para probar la independencia de dos variables entre sí, mediante la presentación de los datos en tablas de contingencia. Pero no se si se puede realiza esa prueba Chi-cuadrado ajustando por otra variable. Y no me refiero a si quiero probar la independencia de X1 y X2 y tengo además la variable *sexo*, realizar la prueba para entre X1 y X2

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

...; encuentra almacenado en el objeto fit. Si lo que necesitas para la > variables A y B depende de la informacion contenida en ese objeto, la > parte 2.1 y 2.2 que mencionas pueden resolverse facilmente creando una > funcion que tome como argumento el objeto que contiene el modelo > ajustado y el nombre de la variable que necesitas calcular (A o B en > este caso). Esta funcion debe tener la posibilidad de no solo > calcular lo que necesitas para cada variable, sino tambien la opcion > de poder exportar el resultado obtenido, ojala con el nombre de la > varible de entra...

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,

...del mismo obtengo por ejemplo la significación (valor p) de la variable independiente, así como el valor del AIC para comparar mi modelo con otro similar, pero no sé cómo obtener un valor del grado de ajuste de la recta de regresión resultante a los datos observados, algo así como el R al cuadrado ajustado que me da el modelo "lm". ¿Alguien sabe cómo obtenerlo? Muchas gracias de antemano, Salu2, Javier

Hola a tod en s, Tengo unos datos y he ajustado un modelo de regresión polinómica de segundo orden. Quiero representar la recta de regresión en el scaterplot y he visto que sería con la función lines() pero me no sale una linea, sino varias. yo únicamente lo que quiero representar es la curva que representa la ecuación de regresión. ¿Alquien...

...ae el libro con el paquete lmer. En este libro se indica que en la comprobación de los supuestos del modelo hay que estudiar los gráficos de residuales tanto del nivel 1 como de los superiores. En el primer conjunto de datos que se analiza hay dos niveles. Obtener el los residuales y los valores ajustados por al nivel inferior no ha sido problema con las funciones residuals y fitted. Para el segundo nivel, la función ranef me extrae los residuales. El problema lo tengo con los valores predichos, y no puedo realizar el gráfico de los residuales frente a los valores predichos. Mi pregunta es: ¿c...

...deviance: 268.31 on 270 degrees of freedom* *AIC: 762.1* Como puede apreciarse, la Null deviance cambia. Entonces, cuando quiero construir la table de análisis de la deviance para reportar los resultados, no sé contra qué comparar los modelos no-nulos: a) contra la Null deviance del modelo Nulo ajustado explícitamente, o b) contra cada Null deviance que arroja el programa al ajustar cada modelo no-nulo. Con la opción (a) la residual deviance del modelo no-nulo es mayor que la del modelo Nulo (263.5 vs. 268.31, lo que constituye un resultado al menos contraintuitivo). Con la opción (b) obviamente l...

Buen dia!! Pues me encuentro trabajando con un conjunto de datos simulados que se ajustan a un modelo Ar(1) y pues queria saber si existe un comando para realizarle una prueba de homocedasticidad, pues la prueba que hay en el R es la bartlett.test pero pues no estoy muy seguro para usarla. Gracias [[alternative HTML version deleted]]

..., ...) { if (!class(x) == "NDfit") { stop("Object must belong to class ND") } print.default(x$estimate, print.gap = 2, quote = FALSE) invisible(x) } El problema es que al hacerlo así, cuando hago coef.ND(ajuste) siendo "ajuste" un modelo ajustado clase NDfit, me sale bien, pero si hago coef(ajuste) me devuelve NULL. Lo mismo me pasa con print.ND, en cuyo caso me devuelve el resultado de la función print.default, o con logLik, donde me devuelve Error en UseMethod("logLik") : no applicable method for 'logLik' applied...

hola soy nuevo usuario de R y necesito crear un objeto p tal que, p=lambda*integral[f´´(x)^2 dx ] donde "lambda" es uno de los parametros que resultan de la funcion "smooth.spline()" y la integral es sobre la derivada 2 de esa misma funcion... dos cuestiones: 1) como extraigo lambda de los resultados de smooth.spline() para usarlo como objeto cuando lo requiera y 2) como

...at=flowed; charset="iso-8859-1"; >        reply-type=original > > Hola Julio: >               en este momento no tengo tiempo de mirar el paquete pero a > penas pueda lo haré, porque yo también me he encontrado con la necesidad de > hacer estimaciones a partir de modelos ajustados con lme y no me resultaron > fáciles. >              ¿le diste una leída a los mensajes de la lista [R-sig-ME]? > uno de los últimos que leí  donde discutían este tema fue " lme and > prediction intervals", también en esa seguidilla de e-mails  recomiendan > leer "Exe...

Estimado Gracias por su ayuda José 2023-05-06 19:47 GMT-04:00, Carlos Ortega <cof en qualityexcellence.es>: > Hola, > > Si añades lo que destaco, defines el layout y queda más ajustado en la > pantalla. > > #-------------------------- > library(readxl) > Sistolica <- read_excel("./Sistolica1.xlsx") > library(meta) > m <- metacont( > Ne, Me, Se, Nc, Mc, Sc, > data = Sistolica, > studlab = paste...

Hola. Hay alguna función que haga un Friedman test (digamos 4 tratamientos o tiempos relacionados/dependientes) y que después haga una comparación de un tratamiento contra el resto, digamos el primero, como un contratase simple, o un Dunnett? o simplemente ¿como hago un Dunnett para unos tratamientos relacionados? -- Antonio M [[alternative HTML version deleted]]

Hola, Construyo un gráfico utilizando layout y manejando convenientemente los mar y oma para luego ubicar los plot y las líneas de texto. El caso es que si ejecuto ese mismo código en un netbook (indicando windows(width=8.27, height=11.69) igual que en el pc), me da el típico error de plot.new : márgenes de figura muy grandes. ¿Cómo puedo hacer para solucionarlo?. ¿O no tiene solución?. Gracias.