search for: fijas

Hola todos Cuando realizo las predicciones del modelo multinomial con el paquete nnet, estas cambian cada vez que lo ejecuto ... saben por qué pasa esto ?? Gracias por la ayuda. [[alternative HTML version deleted]]

...en principio, algo así: lmer.efi<-lmer(mean_efi_week ~plot*year+(1|transect),data = efi_sem, method = "REML") lme.efi<-lme(efi ~year* plot,random=~1|year/week | plot/transect,data = efi_sem, method = "REML") Mi variable respuesta será el número de grillos, como variables fijas he puesto parcela y año (solo hay dos años de muestreo), como variable aleatoria he puesto el transecto (cada transecto se corresponde con una única parcela). El problema es que no sé cómo definir el muestreo semanal como medida repetida dentro de cada año, ¿podríais ayudarme? Muchas Gracias de an...

Necesito realizar un an?lisis del tama?o de muestra para un experimento con las siguientes caracter?sticas: 1) Una variable independiente intrasujeto es el tipo de queso con dos niveles, Queso A (1) y Queso B (2). 2) La otra variable independiente intrasujeto es la ronda de degustaci?n con tres niveles: ronda de degustaci?n a ciegas (1), ronda de degustaci?n de expectativas (2) y ronda de

...Randomforest. Encuentro que: Partial dependence plot gives a graphical depiction of the marginal effect of a variable on the class probability (classification) or response (regression) que nos indica como varía la VR en función de la variable considerada, manteniendo el resto de variables fijas. No encuentro lo que es esa VR por ningún sitio (varianza?), ni la explicación de qué son los dos ejes. Gracias, Manuel -- Dr Manuel Mendoza Department of Biogeography and Global Change National Museum of Natural History (MNCN) Spanish Scientific Council (CSIC) C/ Serrano 115bis, 28006 MADRI...

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

...3. Un saludo, Guillermo > Hola, ¿qué tal? > > Sigue http://www.datanalytics.com/2011/08/04/desarrollo-de-paquetes-con-r-iv-funciones-genericas/ > a rajatabla y lo tendrás. > > Un saludo, > > Carlos J. Gil Bellosta > http://www.datanalytics.com > > P.D.: Si te fijas bien, no estás siguiendo esa guía "a rajatabla". > > El día 23 de julio de 2015, 16:26, <Guillermo.Vinue en uv.es> escribió: > > Hola, > > > > Estoy tratando de crear un método S3 llamado "anthr" dentro del paquete > > que estoy desarrolland...

Otra manera sería usar la función split. Que devuelve una lista de data.frame. Para expotar luego cada elemento puedes usar lapply, ejemplo: data(iris) lista<-split(iris,iris$Species) lapply(names(lista),function(x) write.table(x = lista[[x]],file = paste0(x,".txt")) ) Saludos El 11 de septiembre de 2014, 14:44, Isidro Hidalgo <ihidalgo en jccm.es> escribió: > ¿En qué

Buena noche a todos, Nuevamente requiero de su ayuda con algo puntual: 1. Tengo en el vector Muestras: rm(list = ls()) set.seed(2085089)# Fija Datos data<-matrix(-1,30,ncol=10) > MuestraS [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1 0 0 0 1 1 0 1 1 1 Ahora, dentro de un bucle for voy a generar muestras aleatorias, pero requiero que

...datos. En primer lugar, trato de clasificar las observaciones en las categorías de Y1 mediante un glm y obtengo estimadores de los coeficientes que son muy similares a los valores simulados, como era esperable. Sin embargo, cuando trato de utilizar un glmer con un intercepto aleatorio y pendientes fijas para tratar de distinguir entre categorías, obtengo estimadores muy alejados de los valores simulados. He leído mucha de la bibliografía específica (por ejemplo, los capítulos relacionados con los glmer de los libros de Bates o Goldstein) y he buscado por internet pero no logro encontrar el error...

Hola a todos, Tengo un archivo que va desde el año 2002 hasta el 2010; lo que quiero hacer es separar ese archivo por años.Crear un archivo para cada año, alguien sabe como podria hacerlo? Muchas gracias, un saludo [[alternative HTML version deleted]]

Hola, Estoy tratando de crear un método S3 llamado "anthr" dentro del paquete que estoy desarrollando, cuyo argumento principal es "res" que básicamente es una lista con un solo componente. Pero si el segundo argumento llamado "oneSize" es FALSE, "res" es una lista de listas. Lo que he escrito hasta el momento es lo siguiente: anthr <- function(res,

...l dependence plot gives a graphical depiction of the marginal effect >> of a variable on the class probability (classification) or response >> (regression) >> >> que nos indica como varía la VR en función de la variable considerada, >> manteniendo el resto de variables fijas. >> >> No encuentro lo que es esa VR por ningún sitio (varianza?), ni la >> explicación de qué son los dos ejes. >> >> Gracias, >> Manuel >> >> >> -- >> Dr Manuel Mendoza >> Department of Biogeography and Global Change >> Natio...

Hola Marta, Si el código exacto que utilizaste es el que has compartido, hay un error... Te falta cerrar un corchete. Y he incluido un detalle para que los ficheros tengan como nombre los años... ?#--------------------------? ??filename<-"test1.csv" DBa<-read.csv(filename,sep="*;*" ,header=TRUE) anios<-unique(DBa$Year) for(i in 1:length(anios)) {

...e guardan los datos en esas hojas excel, donde los usuarios, que no son estadísticos, que no conocen R y no lo van a conocer completan esos excel. Es más cómodo hacer que guarden los datos en esos ficheros fijos respetando ciertas reglas en cuanto a formatos, estando dichos ficheros en ubicaciones fijas y llamar los datos desde R directamente a los archivos en sí, que tener que ir Excel por excel abriéndolos. Otras veces, las AAPP cuelgan información en excel que te puedes descargar con R, a algún sitio. Cuando son muchos ficheros, es más cómodo el referenciarlos y buscar los datos (si responden a...

Carlos, mi intencion es entender el codigo: https://github.com/woobe/blenditbayes/blob/master/2013-06-street-level-crime/analysis.R Para luego hacer mi trabajo, algo similar a esto, pero adaptado a mi región. Si te fijas en las lineas 141-161, esta esto: ## Convert the base map into a ggplot object ## All added Cartesian coordinates to enable more geom options later on map.ggmap <- ggmap(map.base, extend = "panel") + coord_cartesian() ## Display a message cat("Creating a ggplot object...

Buen día, no sé si estoy utilizando bien la lista, es la primera vez. Si lo hago mal me corrigen por favor. Sobre tu comentario Pedro, muchas gracias. Lo qeu entiendo con tu sugerencia de set.seed es qeu de esa forma fijas los resultados, pero no estoy seguro si otra agrupación funcione mejor. Es decir me interesa un método de agrupación que genere la "mejor" agrupación y como los resultados de kmeans cambian, no sé cual agrupación elegir. Utilicé otros métodos de agrupación como mclust y pvclust ( http://w...

Alguien tiene experiencia en escribir una pdMat. Para aquellos que no lo recuerden son las matrices de covarianzas de los efectos aleatorios que ajusta la función lme de la librería nlme Estas matrices tiene especial importancia en aplicaciones de genética de poblaciones y en particular en mapeo de asociación. Pinheiro y Bates dicen que el usuario puede crear sus propias pdMat y sugiere como

Buenos días; Estoy intentando aprender algo de sobre los modelos lineales mixtos. Estoy siguiendo el libro "Multilevlel analysis: techniques and applications" de Joop Hox. Estoy intentado reproducir las tablas y gráficos que trae 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

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150128/9e23caec/attachment.html>

Para mi lo más práctico y rápido es prescindir de estos paquetes, seleccionar los datos que me interesan del excel, copiarlos y en R poner: datos<-read.table("clipboard",header=T,dec=",",sep="\t",na.strings="") Y con eso me pega los datos como un data.frame en excel. Un saludo Jesús [[alternative HTML version deleted]]