similar to: suma del resultado de multiplicar fila x columna

Les consulto por lo siguiente: Se me ocurrió usar datos de MCMCglmm en otros procedimientos de R, comparar los resultados y el tiempo como facilidad de uso, para lo cual use LME4, los resultados fueron parecidos, entendiendo que internamente los algoritmos tienen diferencias, es lógico que los resultados no coincidan exactamente, pero aproximadamente tienen que ser compatibles, en este caso

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

Hello all, I'm new in R, and I have a data-frame like this (dput information below): Specie Fooditem Occurrence Volume 1 Schizodon vegetal 1 0.05 2 Schizodon sediment 1 0.60 3 Schizodon vegetal 1 0.15 4 Schizodon alga 1 0.05 5 Schizodon sediment 1 0.90 6 Schizodon

Estimado Carlos Ortega Muchas gracias, al ver su correo recordé por donde estaba eso muy fácil, copio y pego la solución, era tan fácil que no lo recordaba. > table(d) , , Etablecimiento = Est1 Factor Dieta A B 1 220 0 2 0 120 3 0 0 4 0 0 , , Etablecimiento = Est2 Factor Dieta A B 1 0 0 2 0 0 3 0 120 4 0 118

In the text file pressione2008.csv I have the following "Data","MAX","MIN","Note" "07-01-2008 08:00:00", 135, 90, "Eccessi feste, inizio dieta" "07-01-2008 18:00:00", 135, 85, "" "08-01-2008 08:00:00", 125, 75, "" which is a collection of blood pressure data at different time of the day. I would

Estimados Tengo un error que me desconcierta, es un código que simplifiqué de otro trabajo donde no hay problemas, sin embargo me da un error. Una diferencia es que en mi otro trabajo uso spline y ahora polinomio, este es de segundo grado y se encuentra tanto en efectos fijos como aleatorios, el modelo es correcto, corre con MCMCglmm pero no con nlme. grid <- expand.grid(Tiempo=4:6,

Buenos días, *Muchas gracias, todas las aportaciones han sido bien útiles.* Las he tenido en cuenta y he pasado los datos con el R, siguiendo el siguiente comando: *modeloRTUN2<-lm(AVE.~ Tariff + d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8 + d9+ d10 + d11+ d12+ d13+ d14+ d15+ d16+ d17+ d18+ d19+ d20 +d21 + Tariff*d1 + Tariff*d2 + Tariff*d10)* *summary(modeloRTUN2)* Siendo: AVE. = Variable

Hi R users I want to use the lme package for a multilevel analysis on the following example: > math<-c(2, 3,2, 5, 6 ,7 , 7) > sex<-c(1, 2, 1, 2, 2, 2, 1) > school_A<-c(1,1,1,2,2,2,2) > school_B<-c(10,10,10,20,20,20,20) > mydata<-data.frame(math, sex, school_A, school_B) > mydata School_A and school_B are two different school characteristics, math is an

Estimados Estoy realizando algo tan sencillo que no se que forma sería la más apropiada, la situación es la siguiente, hay un data.frame con unos 500 registros, con algunas variables de las cuáles me interesa marcar solo tres, o dicho de otra forma, me interesan tres columnas, de las cuáles cada una tiene entre dos y cuatro factores. Nada complicado, me interesa colocar la combinación de los

Dear R users, I am having problems trying to fit quasipoisson and negative binomials glm. My data set contains abundance (counts) of a species under different management regimens. First, I tried to fit a poisson glm: > summary(model.p<-glm(abund~mgmtcat,poisson)) Call: glm(formula = abund ~ mgmtcat, family = poisson) . . . (Dispersion parameter

Hola a todos, queria saber si existe un comando que me permita sumar filas agrupandolas o algo asi y luego sacarle el promedio a eso. Tengo una matriz de datos con varias escuelas, cada escuela tiene distintos salones de clases, algunas 2 otras tres etc, y cada salon de clases tiene un numero de alumnos, lo que quiero es calcular el promedio de alumnos por escuela, y el promedio de salones por

Buen día, Tengo un data frame con los siguientes datos > Sexo<-c(1,1,1,2,2,2,2,1) > edad<-c(10,12,15,12,10,15,12,15) > frec<-c(150,125,147,123,156,174,185,195) > df<-data.frame(Sexo,edad,frec) Quisiera saber como elaborar una tabla en la cual seleccione el sexo y la edad y R realice la suma de la frecuencia debido a que la lista original incluye miles de

Hi all, When upgrading to 2.11.1 recently I noticed different results being produced by my code. After much digging I have finally narrowed it to a call to nlm(). This can be replicated by: FixedRemovals<-1836180125888 AbStageInitial<-2223033830403 Rates<- 0.3102445 nlm(function(rootM,Abund,Loss,OtherM) {(Loss-(rootM/(rootM+OtherM)* (1-exp(-(rootM+OtherM)))*

Dear all, I have run the confint() function according to below and I get the following error message: > confint(stepAIC.glm.spe.var.konn2.abund, level=0.95) Waiting for profiling to be done... Error: missing value where logical needed In addition: Warning message: NaNs produced in: sqrt((fm$deviance - OriginalDeviance)/DispersionParameter) or > confint(stepAIC.glm.spe.var.konn2.abund,

Hello all, I'm a beginner in R working on a script that will produce a set of models (linear, polynomial and logistic) for each location in a dataset. However, the self-starting logistic model often fails - if this happens I would like to just skip to the next iteration of the loop using tryCatch. I've looked at a few examples and read the help file, but didn't understand tryCatch

No sé si olvidé aclarar, pero quiero crear una lista para cada columna y así almacenar esos valores. On Wed, 18 Jan 2023 at 14:13, David Camilo Gomez Medina < dcgomezme en unal.edu.co> wrote: > Hola Carlos. > > Por ejemplo la segunda columna (16040050), quiero sumar todos los datos de > esa columna correspondientes al año 2000 (quiero relacionarlos con la > primera columna

Mandi! Fyodor Kravchenko via samba In chel di` si favelave... > Thank you, nesting is set to 1. Is there anything else we could check? AFAIK, no. Next step is use a privileged container, but really in this vase better using a VM. -- Berlusconi: "Da oggi sono a dieta" Il Paese lo ? gi? da 4 anni (Il Ruggito del Coniglio)

Hello everybody, In am trying to analyse a BACI experiment and I really want to do it with R (which I find really exciting). So, before moving on I though it would be a good idea to repeat some known experiments which are quite similar to my own. I tried to reproduce 2 published examples but without much success. The first one in particular is a published dataset analysed with SAS by

It is a trivial problem, but in the book I couln`t figure out how to put different colours at different regression lines plot(bif,abund,type="n", xlab= "number_bifurcations", ylab="abundances") sbif<-split(bif,stage) sabund<-split(abund,stage) points(sabund[[2]],sbif[[2]],pch=16, col="red") for(i in 1:2) abline(lm(sabund[[i]]~sbif[[i]])) Thanks in

Hola Carlos. Por ejemplo la segunda columna (16040050), quiero sumar todos los datos de esa columna correspondientes al año 2000 (quiero relacionarlos con la primera columna donde está la fecha) y ese valor almacenarlo en una lista y así sucesivamente con los demás años y con las demás columnas. Tengo pensado una lista así: est_16040050 [[2000]] [1] 2.3 [[2001]] [1] 1.7 [[2002]] [1] 4.8