search for: predecir

Hola a todos, 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

...,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 y time, predecir la probabilidad de germinación. Para eso uso: predict(m1.pile,newdata=data.frame(tem=15,pot=-0.3,time=3),type="response") Con esto me da valores de probabilidad de germinación lógicos y razonables. Por razones ajenas a mi voluntad, necesito poder hacer esto mismo usando los coeficientes...

...rmula que debería usar es esta? > 0,0006077 + 0,7043138*temp + 0,9962766*pot + 0,5060756*time + > 0,4923288*temp^2 + 0,4997649*temp^2*pot > Donde los números son los valores de los coeficientes y temp, pot y time > son los valores de mis variables, los que cambio en la formula para > predecir la germinación en distintas condiciones. > Al hacer el predict me da valores entre 0 y 1 que son probabilidades, pero > al aplicar esta formula me da valores mayores de 1, no tiene sentido. > Preguntado de otra forma ¿Si en R pone I(tem^2):pot eso equivale a > 0,4997649*temp^2*pot? Sien...

Hola a todos, feliz día del trabajador. Y bueno, les consulto: ¿alguno observa una diferencia en predecir con R 3.0? Copio parte del código, viejo y guardado, donde al ejecutarlo me sale el siguiente error: > predict(modelo, predecir, level=0) Error in X %*% fixef(object) : non-conformable arguments Originalmente, el resumen del modelo era con unos datos, donde había por ejemplo tres meses, ene...

..."jun","jul", "ago", "sep", "oct","nov", "dic") ) if(!is.null(contr)) cyc <- C(cyc, contr) model.matrix( ~ cyc ) } > modelo1 <- function(xx, yy) { # xx es serie explicativa (economia) # yy es serie para predecir --- gasolina.eq name.x <- deparse(substitute(xx)) name.y <- deparse(substitute(yy)) xx <- ts(rep(xx,rep(3, length(xx))), frequency=12, start= c(start(xx)[1], (start(xx)[2]-1)*3+1) ) common <- ts.intersect(xx, yy) colnames(common) <- c(name.x, name.y) xx <- common...

Hola Comunidad, Estoy intentado usar una application web (http://rsscalculatorforpets.royalcanin-world.com/) con R. La idea es tener los datos en R y importalos a la applicacion web y bajar los resultados otra vez a R. Alguna idea de donde puedo encontrar informacion relacionada con este tema? Muchas gracias, Marcal Plans Pujolras UPC-ESAB. Campus del Baix Llobregat . Edifici D4. C. Esteve

...1] stats graphics grDevices utils datasets methods base Lo que yo prentendo es clasificar unos datos que están en ficheros de texto, el tamaño de cada uno de estos ficheros va desde los 10mb a 90mb. En un fichero de texto (input) le indico los ficheros que tiene que ir cogiendo para predecir y en el archivo output los nombres de los ficheros de salida con las clasificaciones. Como os comentaba antes, con el clasificador randomForest no tengo problema pero si con support vector. Cuando trabajo con los datos de entrenamiento y test no hay ningún problema. El problema surge al intentar cl...

Buenos días, Soy Cristian Rodelo y quería haceros una pregunta. Tengo un conjunto de datos en los que a un grupo de sujetos se les realizaron medidas de tensión arterial a lo largo de 24 horas. En promedio 4 determinaciones horarias. Es decir que cada sujeto se le realizaron 4 determinaciones horarias durante 24 horas. Bien, quiero determinar la trayectoria longitudinal de las mediciones de

...La red neuronal no produce 0,0,1 como resultado, este podría ser 0, 0,00009, 0,999998. Yo puedo decir a R que como 0,00009 y 0,999998 están próximos a 0 y 1 estos vales 0 y 1, trasladando el resultado a 0,0,1 que significa una categoría (palabra de objeto de la casa). Hasta ahí todo correcto, puedo predecir la categoría. Pero ¿Qué pasa si estas categorías son la presencia de objetos de la casa para clasificar la casa? Quiero decir, techo, pieza es casa común. Otro es techo, pieza, pileta es casa grande. Pero otro usuario solo ingresa pieza y pileta (supone que hay techo) y también es casa grande. En e...

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]))

Estimado Oliver Nuñez Envío un ejemplo reproducible. Javier Marcuzzi # de donde tomo datos, y tiene el modelo (en el pdf) library(MCMCglmm) # librería con las funciónes que voy a usar library(nlme) datos0<-ChickWeight # creo algunos datos que agrego a los origonales Factor<-as.numeric(datos0$Chick) Factor[Factor > 0 & Factor <= 10] <- 'A' Factor[Factor > 10

Es justo ese ejemplo el que estoy mirando, pero no sale la prediccion He probado a cambiar la funcion generadora, haciendo que devuelva como lista solo los input, pero sigue devolviendo error: Error in py_call_impl(callable, dots$args, dots$keywords) : ValueError: Error when checking model : the list of Numpy arrays that you are passing to your model is not the size the model expected.

...mo resultado, este podría ser 0, > 0,00009, 0,999998. > > Yo puedo decir a R que como 0,00009 y 0,999998 están próximos a 0 y 1 > estos vales 0 y 1, trasladando el resultado a 0,0,1 que significa una > categoría (palabra de objeto de la casa). > > Hasta ahí todo correcto, puedo predecir la categoría. > > Pero ¿Qué pasa si estas categorías son la presencia de objetos de la casa > para clasificar la casa? > > Quiero decir, techo, pieza es casa común. > > Otro es techo, pieza, pileta es casa grande. > > Pero otro usuario solo ingresa pieza y pileta (supone...

Hola a todos, en estos momentos me encuentro inmerso en laresolución del siguiente problema. Resulta que dispongo de una variable categóricadenominada Severity y que consta de tres categorías (low, medium, high). Además,tengo otra variable que se denomina ZX y que puede tomar cualquier valorcomprendido entre 0 y 10. Quiero determinar los dos puntos de corte óptimos demanera que me dividan a la

Hola a todos, Estoy iniciandome en R y la verdad es que aun estoy muy muy verde.... Estoy intentando clasificar unos datos con support vector machines, pero me da fallo al usar la funcion predict. El código que uso es el siguiente: calibrate<-read.table("calibration.txt", header=TRUE) calibrate$calibration<-as.factor(calibrate$calibration)

Hola compañeros de la lista. Estoy aprendiendo a usar el paquete "rsm" para superficies de respuesta. Siguiendo este ejemplo todo va bien, hasta que trato de obtener el valor máximo predicho para la variable de respuesta en los valores de x1 y x2 estimados. ------------------------------------------------------------- library("rsm") ChemReact CR <- coded.data(ChemReact,

Muchas gracias, Javier, por tu explicación. Los dotos que manejo son todos datos linguisticos individualizados según variables geográficas y sociologicas; o sea que no hay que preocuparse demasiado por la peligrosidad de mis aseveraciones. El ejemplo lo hice a lo bruto para que se entendiera cómo se estructuran los datos, pensando que así no generaría distración sobre el objeto de mi problema. Lo