similar to: Unable to create/index a zoo irregular timeseries

I have the following code #################################################################### library(zoo) miedate <- yearmon((2006)+seq(0,23)/12) tab <- zoo(cbind(A = c(79.47, 89.13, 84.86, 75.68, 72.82, 78.87, 93.46, 78.18, 82.46, 77.25, 80.95, 84.39, 81.7, 74.76, 65.29, 60.3, 66.59, 73.19, 92.39, 65.76, 77.45, 74.22, 101.36, 100.01), B = c(77.95, 76.73, 51.2, 51.86, 51.29, 49.45,

Buenas tardes, tengo 2 dfs: Dieta de (108x11) y Abund de (591x108). Necesito multiplicar cada columna de la 1ª (108 elementos) por cada fila de la 2ª (108 elementos) y crear una nueva df con las sumas de esas multiplicaciones. He hecho esto, pero no sale y creo que está lejos de estar bien: Res <- matrix(nrow=nrow(Abund),ncol=ncol(Dieta)) Res <- as.data.frame(Res) for(i in

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

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,

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

ng, I have a matrix (x) with binary content. Each row of the matrix holds exactly one 1, and the rest of the row is zeros. The thing is that I need to 'collapse' the matrix to one column where each row holds the original column index of the 1's (y). Sometimes, the matrix is quite large, so I have a perfomance problem. x <- matrix(c(1,0,0, 0,0,1, 0,1,0, 0,0,1, 0,1,0,

hi, I have some session data in a dataframe, where each session is recorded with a start and a stop date. Like this: session_start session_stop =================== 2009-01-03 2009-01-04 2009-01-01 2009-01-05 2009-01-02 2009-01-09 A session is at least one day long. Now I want a dataframe with 'active sessions' per date. Like this: date active_sessions ============= 2009-01-01 1

I'm trying to create a new column in my data.frame where subjects are categorized depending on values on four other columns. In any other case I would just nest a few ifelse statements, however, in this case i have 4*6*2*3=144 combinations and i get weird 'context overflow' errors. So I wonder if there is a more efficient way of doing this. For illustrational purposes, let's say

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

hi everybody I have to do a lot of Anova with R and I would like to have another type of coefficients coding.. I explain. by default if I have 2 temperatures for an experience. 100°C or 130°C and I want to see the temperature effect on the presure I want to estimate the coefficient of each temperature. I will obtain ,with the anova, juste one coefficients for example +3,56 (for 100°C), and the

Hi R folks, In my previous post I forgot to mention that I was new to R. I was really grateful for your quick help. I have two further questions: 1) In the graph of a regression line I would like to show one specific residual yi obs - yi pred (let's take the person whose residual is 76). How do I add a bracket to this vertical distance and name it? I'am getting stuck after the

Hi everyone, I'm new in this list and I hope I'll write my problem in a right way. So, this is my problem: Every an undefined time (It could be 5 min or 20 min) my pc restarts the samba service and It writes this line on the screen: "Reloading /etc/samba/smb.conf smdb only" I have Ubuntu Server 9.10 x32. I am looking for fix this problem for a long time but I don't know

Hi, A small (but annoying) problem with RMySQL library. When a connection is closed, it echoes 'TRUE' to the console. Like this: R> mysqlCloseConnection(con) [1] TRUE The real problem comes when I use RMySQL with Sweave, since the TRUE echo gets into my final Sweave document. Even when I explicitly state <<include=FALSE, echo=FALSE>>= in the header of the R code chunk.

Hi, I've set up a very simple R repository. Just a single source library. Everything works fine. I can install the package on my client using: install.packages(repos='http://www.myServer.se/myRepo/', pkgs='myLib', dep=TRUE) However, I want to protect the repo, so I use a .htaccess, placed directly under 'myRepo' on the server. I use 'Authentication Basic' and

Hello there, Could anybody please help on how to correctly use ggplot when printing irregural time series, by irregural here I mean for example the absence of some dates in the whole timespan of a dataset. To be more precise in the following example I generated some random data which spans the whole November up to now and dropped weekend days but for some reasons ggplot continues to plot the

Hello Marco, On 05.02.2024 16:44, Marco Gaiarin wrote: > Mandi! Fyodor Kravchenko via samba > In chel di` si favelave... > >> Have to add about the environment - this is an unprivileged TurnKey >> Fileserver Linux container run under Proxmox. The extensive googling for >> the problem suggests Samba will not work in such environment because of >> ACL and such,

Greetings -- I've got some sensor data of the form t1_1, t1_2 t2_1, t2_2 ... tN_1,tN_2 -- time intervals measuring starts and stops of sensor activity. I'd like to see whether there's any regularity in it. Seems natural to consider these data timeseries -- except most of the timeseries packages and models assume regular ones, with a fixed frequency. I wonder what's a

I saw an email about 8.0.1 being not quite out, and planning for 9.0.0. When is 8.0.1 supposed to come out? I installed llvm at develop on x86 on March 30, and its clang says it is version 9.0.0 (trunk 357348). I tried to do another spack install llvm at develop, but it said that I already had it installed. Does that mean there were no changes since March 30, or does spack not look for any