similar to: Evaluar función?

Hola, quiero hacer un estudio de no inferioridad para ver si tres tratamientos quirúrgicos son equivalentes en el tratamiento de una patología. Me gustaría saber qué librerías se pueden utilizar y ver algún ejemplo. Gracias J Martín

Compañeros, Una de primero de funciones pero que no soy capaz de resolver. Necesito evaluar como texto el parámetro de mi función. Ejemplo: datos.A<-rnorm(10,3) datos.B<-rnorm(100,30) datos.C<-rnorm(1000,300) datos.D<-rnorm(10000,3000) describe <- function(letra){ summary(datos.letra) } describe(B) Quiero que haga summary(datos.B) y no soy capaz. Nunca me había visto en una

Buenos días, Me gustaría aplicar una regresión múltiple a los datos con los que trabajo pero no se como introducir los datos en R. He probado introducir el siguiente comando: modeloM<-lm(AVE.~ d1 + as.factor(d1T)*Tariff) summary(modeloM) Pero me da el siguiente error: > summary(modeloM) Error in if (attr(z$terms, "intercept")) sum((f - mean(f))^2) else sum(f^2) : argument is

Hola a todos, Me gustaria saber como puedo evaluar / obtener el resultado guardado en un objeto cuando este se ejecuta. Quiero decir, el resultado de la siguiente expresion: > paste("res_UP_sep_dec_","16_17_18",sep="") [1] "res_UP_sep_dec_16_17_18" es un objeto (una lista) cuyo nombre es precisamente res_UP_sep_dec_16_17_18 >

Buenos días a todxs, Estoy comparando la predicción de los valores (0, 1, 2, 3,.....hasta 13) frente a los observados. Con la idea de comparar el modelo Zero inflated y el Binomial negativo y ver cual presenta mas distancia frente a las predicciones observadas. Para ello introduzco los códigos en la consola: #Modelo ZIM pred<-round(colSums(predict(zeroinfl, type="prob") [,1:14]))

Hola a todos. ¿Cómo puedo ajustar un modelo de regresión logística en el que la varianza de los errores no son iguales, sino que esta puede ser modelada por otras variables ? ¿Utilizando modelos mixtos o voy desencaminado? Gracias..

Hola, he estado buscando por todos lados, como puedo saber que traigo en un una variable que es un list, y pues nomas no encuentro nada. alguien sabe alguna función como dim, en caso de los data.frame que nos da las dimensiones, o length en un vector que nos da cuantos elementos tiene algo como names para saber el nombre de los objetos dentro del list, en esencia me gustaría saber si alguien

Buenos días, Estoy realizando un desarrollo sobre un problema de programación lineal. Tengo una función en la que añado al modelo todas las restricciones y quiero controlar que esa parte no da error (por longitudes de vectores, por valores inexistentes, etc). He visto que existe, como en otros lenguajes, un tryCatch, pero no sé bien como utilizarlo. Estoy realizando pruebas a pequeña escala pero

Amigos de la erre. He creado mi primer bucle con for para entrenar unos modelos con GAM. La respuesta es quasipoisson porque estoy trabajando con densidades de peces. Sin embargo, tengo un problema, no se muy bien como añadir el offset a la formula siguiente cuando creo el bucle. GAM.A1 <-gam ((DYO)~s(DMA,k=4)+ s(WOD,k=4)+s(CIN,k=4)+s(DRA,k=4)+s(DBR,k=4)

> > Estimados muy buenas quería hacerles unas consulta: Estoy trabajando en mi tesis sobre mejoramiento animal y mi objetivo es evaluar la habilidad predictiva de modelos estadísticos mediante validación cruzada. Pero antes la intención es dividir mi base de datos en 3 partes y quisiera que todos los efectos incluidos en el estudio y cada uno de sus niveles, estén lo más equitativamente

Hola, Hola a todos, Os escribo porque no consigo finalizar el script necesario para realizar lo que a continuación planteo. Partiendo de un data frame (2 millones de casos), tengo: > datos2 ID FECHA YEAR CANTIDAD 1 100 2005-08-02 2005 1 2 100 2005-10-19 2005 2 3 100 2007-02-09 2007 1 4 100 2007-10-25 2007 1 5 100 2007-10-29 2007 1 6 120 2006-05-11

Buenas. Alguien sabe como representar juntos , de forma fácil, los boxplot de todas las columnas de un data.frame. A mí lo único que se me ha ocurrido es unir las columnas en un vector numérico y crear un factor que me diga a que columna pertenece cada uno . ¿hay alguna forma más fácil? indice <- rep(1:ncol(datos),each=nrow(datos)) indice <- factor(indice) levels(indice) <-

Hola!! Estoy intentando evaluar mi modelo de series temporales (uso auto.arima). Para ello he implemetado el método "rolling window" que se basa en ir añadiendo progresivamente datos al conjunto de train para testar el modelo. Por ejemplo: - Train: 1 año, test: día 1 (24 observaciones, una por hora) --> evalúo ese día (RMSE por ejemplo) - Train: 1 año + 1 día, test: día 2 -->

Hola, Vaya, en el código que he enviado, cusum no se incrementaba.. Y has indicado que se introduce un año más, con el mismo ID que el anterior y con la misma cantidad acumulada. Pero si el siguiente año es del mismo ID, acumula el valor de la cantidad que hemos introducido en esa fila... Con el siguiente código se resuelve este error y además ya está preparado para contemplar cualquier tipo de

Hola! Necesito evaluar la normalidad de una variable para la que tengo más de 5000 observaciones. Shapiro-Wilks no funciona para muestras mayores 5000 observaciones. AAlshap <- lapply(AAdf, shapiro.test) Error in FUN(X[[i]], ...) : sample size must be between 3 and 5000 Alguna sugerencia? Gracias. [[alternative HTML version deleted]]

Hola Comunidad, Tenia al siguiente duda ¿hay alguna forma de que R me de la Esperanza y Varianza de una distribución según sus parámetros? es decir, si tengo una weibull de parámetros "X" y "Y" para la esperanza hay una formula, pero me gustaría saber si hay una función en R que le diga algo como Esperanza(weibull, shape = "X", scale = "Y") y ya devuelva

Buenas a todos, acá estoy yo de nuevo con problemas de loops. Tengo el siguiente problema: un vector de datos (y) y una serie de efectos. El loop lo que intenta es evaluar el R² de modelos incrementando por vuelta una variable efecto. Seria algo así: for(i in 1:x) { modelo=lm(y~efectos[,1:i]) ... codigo para guardar R² y otros por cada vuelta... } apply() puede ir columna a columna pero no

Mira el conjunto de mandelbrot en r cómo mola https://www.google.es/amp/s/www.r-bloggers.com/the-mandelbrot-set-in-r/amp/ El 3 nov. 2017 4:20 PM, "eric" <ericconchamunoz en gmail.com> escribió: > hola luis, podrias explicarte un poco mas ? que quieres decir con "obtener > fractales de una imagen" ? > > Puedes construir una imagen con fractales, pero

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, LLevo buscando desde hace tiempo como hacer el Goodness of fit test en R. Es decir, me explico, intento hacer una cosa parecida que se hace en Minitab, por ejemplo, yo tengo un conjunto de datos, y lo que quiero es sabes que tipo de distibución es, en minitab se hace un histograma para ver si se ajusta bien o no a la campana de Gauss, luego vemos si aproximar la distribución de la muestra