search for: introduzco

...decir que llevo muchas horas buscando información y peleándome con los manuales, tanto en inglés como en español, pero no sé como resolver mi problema que es el siguiente: ¿Cómo puedo hacer para crear un programa con el que, al darle una función y un punto, me diga que valor tiene? Es decir, si yo introduzco x^2 y el punto 2, que me diga que vale 4 (2^2=4). La idea sería poder poner cualquier tipo de función, pero si me dáis una solución que sólo sirva para introducir funciones polinómicas también me sirve. Lo que yo he hecho es lo siguiente, en el que la x seria la función a introducir y a sería el pu...

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])) #Valores observados realmente obs<-table(IB$nijt)[1:14] #Tabla comparativa rbind( pred, obs) Y obtengo la siguiente tabla: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 pre...

...uasipoisson(link = "log"),data=as.data.frame(cbind(DYO,Conjunto))) Os escribo el codigo del bucle: for(i in (1:nrow(VarCombS))){ Formula <- as.formula(paste("DYO",sep=" ~ ",paste(as.matrix(VarCombS)[i,1:NVar],sep="",collapse=" + "))) ### Cómo introduzco el offset?? GAM <- gam(Formula,family=quasipoisson(link = "log"),data=as.data.frame(cbind(DYO,Conjunto))) Summ <- summary(GAM) Var.p.values <- Summ$s.table[,"p-value"] MaxSignif <-c(MaxSignif,max(Var.p.values)) Significance <- all(Var.p.values<0.1) #...

...mp$difID[i] <- 0 df.tmp$difYE[i] <- 0 } else { df.tmp$difID[i] <- df.tmp$ID[i] - df.tmp$ID[i-1] df.tmp$difYE[i] <- df.tmp$YEAR[i] - df.tmp$YEAR[i-1] } } } df.tmp #------- Segundo bucle para introducir filas en los saltos de años # Introduzco filas cuando el salto de años sea mayor que 2... df.new <- 0 for(i in 1:nrow(df.tmp)) { #Copio la fila tal cual cuando la diferencia en años es 0 o menor que dos. if(df.tmp$difYE[i] < 2) { df.new <- rbind(df.new, c(df.tmp$ID[i] , df.tmp$YEAR[i] ,df.tmp$cusum[i])) } else {...

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

...gt; 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])) >> #Valores observados realmente >> obs<-table(IB$nijt)[1:14] >> #Tabla comparativa >> rbind( pred, obs) >> >> Y obtengo...

...ID[i] <- df.tmp$ID[i] - df.tmp$ID[i-1] > df.tmp$difYE[i] <- df.tmp$YEAR[i] - df.tmp$YEAR[i-1] > } > } > } > #df.tmp (Lo comento porque no es el resultado finalmente buscado) > > #------- Segundo bucle para introducir filas en los saltos de años > # Introduzco filas cuando el salto de años sea mayor que 2... > df.new <- 0 > for(i in 1:nrow(df.tmp)) { > > #Copio la fila tal cual cuando la diferencia en años es 0 o menor que dos. > if(df.tmp$difYE[i] < 2) { > > df.new <- rbind(df.new, c(df.tmp$ID[i] , df.tmp$YEAR[i...

Hola a todos: Les pregunto, desconociendo si es posible de forma fácil. Utilizando la librería lme4, introduzco el modelo, se ejecuta, y obtengo resultado por ejemplo con el comando ranef(modelo). De esta forma entiendo que obtengo la soluciones. Pero me gustaría conocer el PEV (prediction error variance) junto a las soluciónes. Hay algún comando adecuado para este caso (tipo ranef). Desde ya muchas...

...uasipoisson(link = "log"),data=as.data.frame(cbind(DYO,Conjunto))) Os escribo el codigo del bucle: for(i in (1:nrow(VarCombS))){ Formula <- as.formula(paste("DYO",sep=" ~ ",paste(as.matrix(VarCombS)[i,1:NVar],sep="",collapse=" + "))) ### Cómo introduzco el offset?? GAM <- gam(Formula,family=quasipoisson(link = "log"),data=as.data.frame(cbind(DYO,Conjunto))) Summ <- summary(GAM) Var.p.values <- Summ$s.table[,"p-value"] MaxSignif <-c(MaxSignif,max(Var.p.values)) Significance <- all(Var.p.values<0.1) #...

> Hola a tod@as a ver si me podéis ayudar: > > tengo una base de datos en excel con 120 variables y 58 observaciones para > cada una de ellas y no encuentro la manera de que reconozca la base de > manera correcta, es decir que reconozca que son 120 variables y que se trata > de series temporales. He usado el paquete ts, zoo pero tampoco así consigo > indexarlas de manera

Hola Alfonso, buenos días: Lo he probado pero no me va, eso sí ha cambiado el error, el código utilizado es: library(ROAuth)library(twitteR) download.file(url="http://curl.haxx.se/ca/cacert.pem", destfile="cacert.pem") reqURL <- "https://api.twitter.com/oauth/request_token"accessURL <- "http://api.twitter.com/oauth/access_token"authURL <-