similar to: Problem with R 1.9.1

Hi, I make this model: > moths.m6 <- glm(nind~metros+especie+habitat+metros:habitat+habitat:especie,family=poisson) > anova(moths.m6,test="F") Analysis of Deviance Table Model: poisson, link: log Response: nind Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F) NULL 81 488.32

Buen día, Tengo un data frame con los siguientes datos > Sexo<-c(1,1,1,2,2,2,2,1) > edad<-c(10,12,15,12,10,15,12,15) > frec<-c(150,125,147,123,156,174,185,195) > df<-data.frame(Sexo,edad,frec) Quisiera saber como elaborar una tabla en la cual seleccione el sexo y la edad y R realice la suma de la frecuencia debido a que la lista original incluye miles de

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

Hi, I am trying to run analyzes incorporating sample weight, strata and cluster (three-stage sample) with PNS data (national health survey) and is giving error. I describe below the commands used. I could not make the code reproducible properly. Thanks, ################################################# library(survey)####change to 0 and 1 variable

Hi, I make an analysis and depending of the order of the variables, the significance change, look. m1 <- glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot) > anova(m1,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: (infec/ntot) Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F)

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

Hello, I?m starting with my PhD and I have to stop because i got a little knowledge in R and statistics. I?ve got a model of this kind: binary response variable: prevalence of infection (0/1) 3 categorical independent variables: sex, month and name of the area I was trying with a full model like this, before the simplification model<-aov(prevalencia~sex*month*area) but the Fligner test

Un saludo a todos. Es posible que con el asunto no quede claro lo que quiero. Para eso añado una captura. Dispongo de 3 variables (dos discretas y una continua) y me interesaría sacar la media tanto del grupo como de los subgrupos. Esta creo que es la segunda vez que escribo preguntando algo parecido. En la primera me hablasteis de tapply, pero con este comando sólo sé obtener la media de una de

Hi, I have this data set: obitoss = c( 5.8,17.4,5.9,17.6,5.8,17.5,4.7,15.8, 3.8,13.4,3.8,13.5,3.7,13.4,3.4,13.6, 4.4,17.3,4.3,17.4,4.2,17.5,4.3,17.0, 4.4,13.6,5.1,14.6,5.7,13.5,3.6,13.3, 6.5,19.6,6.4,19.4,6.3,19.5,6.0,19.7) (dados = data.frame( regiao = factor(rep(c('Norte', 'Nordeste', 'Sudeste', 'Sul', 'Centro-Oeste'), each=8)), ano =

Hola. Estoy intentando crear una función pero no logro que termine de cerrar y mi manejo no permite que el google me ayude... n.reg.dep <- function(x, y) { etiqueta <- str_replace(nombres[Variable == y, Descripcion], "[ ]", "") tabla <- datos[, .N, by = .(etiqueta = A3a)] # tabla <- tabla[, Porc := round(N/sum(N)*100,1)] # % tabla[, PorcAc :=

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

Hi everyone, Does anyone know if there is a package to compute raking weights using R? What I need is to create a variable with weights base in some demographic variables (e.g. sexo, age group, area) using the raking procedure. Thank you in advance! -- Sebasti?n Daza

Hola, Quiero especificar una ecuación con varias variables explicativas de una manera eficiente sin necesidad de escribir todas y cada una. Tengo un conjunto de variables (junto con otras) dentro de una base de datos que se llaman pot23 pot311 pot312 pot 316 pot317........... pot80. No necesariamente están secuenciadas. Quisiera saber cómo indicar que incluya todas las variables de pot23 a pot80

Hola a todos, Por fin pude terminar de preparar la descarga y procesamiento de los datos por edad y sexo. Ya están en el repositorio: https://github.com/rubenfcasal/COVID-19 Ya puestos también incluí la descarga e importación de los datos acumulados del ISCIII y actualicé las tablas en https://rubenfcasal.github.io/COVID-19/COVID-19-tablas.html empleando el paquete DT. Ahora tiene tres

Hola a todos, Soy nueva usando R y tengo unos problema en procesado de datos que quiero hacer en R. Tengo un fichero de datos en formato texto separado con ; (punto y coma). Este fichero contiene las columnas: Poblacion;Sexo;Nota; EdadP, EdadM Tengo 2 problemas que no se como tratar. 1.- Problema La columna *Poblacion* debo tenerla siempre con la primera letra en mayúsculas y tengo

Estimado Mauricio Monsalvo R es un lenguaje medio complicado, no es orientado a objetos, aunque hay formas para un trabajo con objetos, por otro lado se puede definir una función o emplear 5 paquetes para lo mismo en dos líneas de código. Usted tiene 1500 líneas de código, no es fácil de comprenderlas, tampoco hay recetas mágicas. Yo ordenaría como si fuesen cajas, una entrada, un proceso

I am using pyramid.plot() from the plotrix package. I have something like this ############################################ xy.pop<-dados$masfr xx.pop<-dados$femfr #agelabels<-dados$femlab xycol<-color.gradient(c(0,0,0.5,1),c(0,0,0.5,1),c(1,1,0.5,1),11) xxcol<-color.gradient(c(1,1,0.5,1),c(0.5,0.5,0.5,1),c(0.5,0.5,0.5,1),11) xylab<-dados$maslab xxlab<-dados$femlab

Buenas tardes a todos los participantes del foro. Me dirijo a vosotros porque estoy atascado con una duda de programación respecto al data frame: > dd # Data frame de 5 variables, leído de un archivo txt id sexo nacim origen final 1 1 0 02/09/1955 01/04/1985 01/02/2014 2 2 1 29/10/1951 15/08/1996 01/05/2009 3 3 0 30/10/1942 02/08/2000 01/02/2014 4 4 1

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 amigos, Estoy familiarizándome con los modelos lineales generalizados en R. Estoy interesado en realizar un análisis lig linear y me gustaría saber cuáles son o como extraer los valores correspondientes al chi cuadrado en el análisis para cada grupo y para las interacciones. Desde ya muchas gracias y disculpas si la pregunta es muy básica, adjunto los comandos que estoy utilizando. Si