search for: aleatorio

...tware 23(7)). El modelo que intento correr es: m1<-gamlss(Teleosteos ~ cs(LT_max)+random(Familia/Genero,lambda=T),family=BEINF, data = datos) donde LT_max es la talla máxima de especies de peces y la variable respuesta es la proporción del consumo de peces teleósteos de cada especie. El efecto aleatorio es anidado de géneros y familias taxonómicas de especies (género dentro de familia). El problema que tengo es que los gamlss no toman el efecto aleatorio como random(Familia/Genero) o (1|Familia/Genero), sino como: m2<-gamlss(Teleosteos ~ cs(LT_max)+ random(Familia,lambda=T)+random(Familia:Gene...

Queridos R-users: Tengo una duda que hace mucho tiempo que estoy intentando resolver, os explico a modo de ejemplo: Tengo estos efectos: Año(5 niveles),Localidad (10 niveles) y genotipo (3 niveles), año y localidad son aleatorios y genotipo es fijo (los he escogido yo). Me gustaría hacer obtener una tabla parecida a la Tabla Anova donde aparezca cada factor y sus interacciones y finalmente hacer un comparación de medias del efecto fijo y del efecto Localidad. NO ser si es posible hacer tal cosa con los efectos al...

Reciban un saludo cordial. Por favor, ¿en cuál biblioteca de R existe una función para generar vectores aleatorios correlacionados, por ejemplo (X,Y), con X~Weibull y Y~Lognormal y corr(X,Y) = 0.7? Gracias desde ya por su atención/ayuda. César Escalante C. [[alternative HTML version deleted]]

Hola a todos, Tengo algunas dudas sobre como introducir en un modelo lineal factores con efectos fijos y aleatorios. Mi diseño es el siguiente: Factores Line: 40 líneas de trigo (Fijo) Rep: 3 Bloques (Aleatorio) Year: 2 Años (Aleatorio) Variable dependiente alpha.ug.mg Nota: Adjunto tabla de datos Dentro del diseño, 'Rep' y 'Year' se considera de efectos aleatorios y 'Line' e...

Estimados quería realizarles la siguiente consulta: 1) ¿es posible dentro de un modelo con enfoque bayesiano, anidar un efecto con distribución multivariada ("aleatorio" en términos frecuentistas) sobre un factor de distribución no multivariada (fijo en términos frecuentistas) en R? . Trabajo sobre datos de competencia en enduro y dispongo de una variable Año (33 niveles) y una variable Reglamento (3 niveles) y la distribución de datos es la siguiente:...

Hola de nuevo, ya empiezo a ser pesado ¿no? bueno, no importa porque aprendemos todos. Eso, al menos, me parece. Hoy estuve estudiando en R el tema de la aleatoriedad. Veo que hay múltiples posibilidades pero me están chocando mucho. Encuentro que el generador de números pseudo aleatorios es más pseudo de lo que debería. Me explico, quiero generar 0 y 1 aleatorios. Estoy trabajando con una martingala. Uso x <- runif(1,0,1) y después según cual sea la probabilidad de ganar elijo el punto de división del intervalo apropiado. Bien, los modelos que extraigo no me parecen coheren...

...on Producción neta, Peso de 100 semillas y Proteína. Tercero: Observo los mismos genotipos (G) en varios años (A) en las misma localidades (L). La especie es Phaseolus vulgaris L., en forma de judía seca, por lo que cada año se planta y se recolecta. Con el lme no ser introducir dos efectos aleatorios, ni su interacción; y en el caso de que no exista significación en ninguna interacción con G, como podría hacer la comparación múltiple de medias. Muchas gracias de nuevo y gracias por el comentario Olivier Nuñez. Marçal Plans -- UPC-ESAB. Campus del Baix Llobregat . Edifici D4. C. Esteve...

...rix and marginal distributions Saludos, Carlos Ortega www.qualityexcellence.es El 18 de diciembre de 2013, 12:42, Cesar Escalante <c.escalante.c@gmail.com>escribió: > Reciban un saludo cordial. > > Por favor, ¿en cuál biblioteca de R existe una función para generar > vectores aleatorios correlacionados, por ejemplo (X,Y), con X~Weibull y > Y~Lognormal y corr(X,Y) = 0.7? > > Gracias desde ya por su atención/ayuda. > > César Escalante C. > > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-e...

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 ejemplo ver una pdDiag, pero cuando abro una...

Saludo cordial para cada uno. Les pido ayuda para generar números aleatorios de una mixtura k-puntos. Sabemos que la función de distribución F es una mixtura k-puntos si es de la forma F(x) = p_1 F_1(x) + p_2 F_2(x) + … + p_k F_k(x), donde F_j es una función de distribución de probabilidad, p_j > 0 y suma(p_j) = 1, para j = 1, 2, …, k. En mi caso particular F es la...

...if(random_matrix[x_list[random_sq]-1, y_list[random_sq]]== 0 || random_matrix[x_list[random_sq] + 1, y_list[random_sq]] == 0 || random_matrix[x_list[random_sq], y_list[random_sq] + 1] == 0 || random_matrix[x_list[random_sq], y_list[random_sq] - 1] == 0) { aleatorio <- runif(1) if (aleatorio < 1) { if (random_matrix[x_list[random_sq],y_list[random_sq]-1]== 0 ) { random_matrix[x_list[random_sq], y_list[random_sq] - 1] <- 1 flag <- 1 range_counter <- range_counter + 1 x_list[range_counter] <...

...2: In pf(q, df1, df2, lower.tail, log.p) : NaNs produced 3: In pf(q, df1, df2, lower.tail, log.p) : NaNs produced Y hasta aquí llegue. Veo que algo está mal con la distribución de los grados de libertad, ¿Qué es lo que tendría que odificar para resolver el problema? ¿Están mal elegidos los efectos aleatorios? ¿Tendría que usar otra formulación del modelo u otro modelo? Muchas gracias a todos Saludos Santiago -- Santiago. L. Poggio Ing. Agr., Dr. CONICET - Cátedra de Producción Vegetal Departamento de Producción Vegetal Facultad de Agronomía Universidad de Buenos Aires Av. San Martín 4453 (C1417DSE...

Hola a todos. Estoy utilizando la función lmer del paquete lme4 para ajustar un modelo mixto. Tengo varias variables en mi data.frame, unas son a nivel individual y otras a nivel de comarcas. Listo algunas. ingre_6 : Ingresos (nivel individual) iscs_a : un indicador sintético resumen de otras variables, calculado mediante componentes principales. sau_com : superficie agraria útil de cada una

Hola comunidad, Les consulto mi duda. Quisiera generar N números aleatorios tal que su suma esté predeterminada a un valor k. Sé que con la distribución Dirichlet, por ejemplo, se podrían obtener números aleatorios tal que su suma es la unidad. Por ejemplo MCMCpack::rdirichlet(1,rep(10,7)) genera 7 números tal que su suma es 1 (el parámetro rep(10,7) es cualquier cosa)....

...cuado he simulado unos 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...

...ción con lo siguiente: Tengo un vector dicotomico (Binario) con la siguiente información que me surgio de algunos procesos anteriores: > MuestraS [1] 1 1 1 1 0 1 1 0 1 1 Ahora necesito hallar todas las posibles combinaciones que se puedan hacer con estos elementos para luego hacer un muestro aleatorio simple con reemplazo y seleccionar algunas de las posibles muestras que se obtengan. Agradezco me puedan ayudar con este asunto Saludos, DC. [[alternative HTML version deleted]]

Hola, Esta pregunta la hice en stackoverflow <http://stackoverflow.com/questions/31137940/randomly-assign-teachers-to-classrooms-imposing-restrictions/31143808#31143808>pero nadie pudo contestarla. 1. Quiero generar N escuelas, con G grados y C divisiones. 2. Quiero asignar cada uno de T maestros a 2 divisiones en un grado y escuela Si tengo C=4 divisiones, puedo lograr lo que quiero con

Hola a todos He estado trabajando con el paquete osmar y funciona muy bien.Puede alguien decirme donde encontrar informacion o algun script para generar una lista con N pares de puntos (longitud/latitud) aleatorios o siguiendo una funcion determinada dentro del marco de coordenadas definido por corner_bbox(left, bottom, right, top) para despues plotearlos. saludos eloy ortiz -- Este mensaje le ha llegado mediante el servicio de correo electronico que ofrece Infomed para respaldar el cumplimiento de la...

Hola quisiera saber de que manera puedo generar numeros aleatorios enteros entre 0 y 4, para completar una matriz de un automata celular Muchas gracias Juan Nicolás Corral GómezEstudiante Carrera de EcologíaPontificia Universidad Javeriana [[alternative HTML version deleted]]

Hola, Soy nuevo manejando R y no tengo mucha experiencia. Estoy intentando modelar una función que me relacione el nº de cebas (nº de presas que los padres traen a los pollos) con el tamaño de parche de un bosque (factor categórico; 2 niveles= grande y pequeño). Al ser un conteo (nº de cebas) he pensado utilizar familia= poisson con link= logarítmico. He construido un GLMM con: Nº de cebas