similar to: replication time series using permutations of the each "y" values

Hola a todos, Estoy trabajando con modelos lineales generalizados (GLM), en particular con las funciones glm y anova.glm de R. Tengo una pregunta que es más bien técnica en el sentido estadístico (y un poco Off Topic y probablemente un tanto naive). No tengo claro si es correcto decir análisis de varianza o de devianza al utilizar GLMs, veo que en el help de anova.glm dice análisis de

Hello everybody, I recently wrote a "program" that to generate AR1 time series, here the code: #By Jomopo. Junio-2007, Leioa, Vizcaya #This program to create the AR1 syntetic series (one by one) #Where the mean is zero, but the variance of the serie AR1 and #the coef. of AR1 are be changed. If var serie AR1 = 1 then is standarized! #Final version for AR1 time series program #Mon Jul

Hello everybody, I want to save a graphic using the "pdf environment", everything in order, but I can not to save the graph in "horizontal" way. I have read the manual (to plot, par), the FAQ, google, but nothing I found that can help me. Anyone can say me a hint? Thanks a lot, -- Josue [[alternative HTML version deleted]]

Hi everybody, I am reading the Lomb paper (Lomb, 1976) and I found an interesting equation, and I wish to resolve it using R. I am wondering if anybody has a hint. The equation is: Yis[i] = a*cos((2*pi/T)*(Times[i] - Tau)) + b*sin((2*pi/T)*(Times[i] - Tau)) ... (1) Where T and Tau are constants. I know the "Times" and "Tis" values (in fact these values come from a Time

Hi, I have recently used the CTS package in order to use the Lomb-Scargle periodogram (spec.ls) function. I have noticed an issue that I hoped you may be able to explain. If a regularly spaced time series has two points removed, one at either side of a single data point (thus making an irregularly spaced time series), a spectrum with a very large peak at the highest frequencies is produced. An

Hi all, I am trying to do wavelets and I got an error message saying "The length of data is not a power of 2" Is there a way of handing that? or should the data length be exactly the power of 2? I am using R version 2.9.2 (2009-08-24) The is library(wavethresh). wds <- wd(ds$v,filter.number=1) Thanks

Hola, quisiera que me dijeran como puedo mandarle unos datos al R, de una aplicación en java. "La educación es lo que queda cuando uno se olvida de todo lo que aprendió en la escuela." [[alternative HTML version deleted]]

Dear all, Before I file this as a bug, I wanted to check if I didn't miss something. The help page of lag() says that the function returns a time series object. It actually does return something that looks like a ts object (the attribute tsp is set). But when using a vector, the class "ts" is not added to the result: > avec <- 1:10 > lag(avec) [1] 1 2 3 4 5 6 7 8

Hola, estoy haciendo un estudio de la desaceleración de la mortalidad y tengo dudas con la función nls(), adjunto un archivo donde he redactado mi inquietud. saludos cordiales Boris Polanco ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20141015/b9761c95/attachment-0001.html> ------------ próxima

Dear List, today I turn to you with a next problem. I'm trying to compare species richness between various datasets (locations) using species accumulation curves (Chapter 4, page 54 in Tree diversity analysis<http://www.worldagroforestry.org/treesandmarkets/tree_diversity_analysis.asp>by Kindt & Coe). To accomplish this I'm using package BiodiversityR. My data is comprised of

Hello, I am trying to place a legend in the title area of a plot. I can not seem to find a way to get legend() to do this. Nor can I use the title() function. Does anyone know a way to do this? I would greatly appreciate any advice. Thanks Josue Samayoa Center for Biomolecular Science & Engineering School of Engineering Division of Physical & Biological Sciences University of

Hello, I have two series, one with stream stage measurements every 5 minutes, and the other with barometric pressure measurements every hour. I want to subtract each barometric pressure measurement from the 12 stage measurements closest in time to it (6 stage measurements on either side of the hour). I want to do something like the following, but I don't know the syntax. "If the

Dear Ladies and Gentlemen how are you? I would like to wish everyone and all their families, may God continue to bless and always illuminating your steps, that this Christmas the joy of our Lord Jesus invade your home and participate in their family on this Christmas. May God bless you always! Merry Christmas and a Happy and Prosperous 2014. With Best Regards! Josue -------------- next part

I've noticed that I get a warning message every time a package masks some functions from Hmisc. The warning message says : Warning message: In identical(get(., i), get(., lib.pos)) : ignoring non-pairlist attributes This happens with eg: library(plyr) library(xtable) I think I've seen this passing by before, but I'm not sure any more. Just thought I'd mention it. Cheers Joris

Hi, I would like to have an index for a column in a matrix encoded in a cell of the same matrix. For example: x = matrix(c(11,12,13,1, 21,22,23,3, 31,32,33,2),byrow=T,ncol=4) In this case, column 4 is the index. I then access the column specified in the index by: > for (i in 1:3) print(x[i,x[i,4]]) [1] 11 [1] 23 [1] 32 > > for (i in 1:3) {x[i,x[i,4]] <- x[i,x[i,4]] + 5} > x

?typeof? is your friend here: > typeof(`[`) [1] "special" > typeof(mc[[1]]) [1] "symbol" > typeof(mc2[[1]]) [1] "special" so mc[[1]] is a symbol, and thus not a primitive. - Lukas > On 28 Mar 2017, at 14:46, Michael Lawrence <lawrence.michael at gene.com> wrote: > > There is a difference between the symbol and the function (primitive >

Dear All, I am currently using the MI package (Su, Gelman, Hill and Yajima) to make multiple Imputations of my dataset with missing values. After fitting a model, I can use display(model) to visualize the pooled estimates as well as estimates of each imputed dataset. I can visualize these also by typing print(model). However I would like to know how I can extract estimates of single imputed

On 31/01/2018 6:33 AM, Joris Meys wrote: > 3. given your criticism, I'd like your opinion on where I can improve > the documentation of https://github.com/CenterForStatistics-UGent/pim. > I'm currently busy updating the help files for a next release on CRAN, > so your input is more than welcome. After this invitation I sent some private comments to Joris. I would say his

Dear gurus, I was utterly surprised to learn that one of my examples illustrating the need of match.fun() doesn't give me the expected result. center <- function(x,FUN) FUN(x) center(1:10, mean) mean <- 4 center(1:10, mean) Used to give me the error message "could not find function FUN". Now it just works, even though I didn't expect it to. I believe this is at least

I would like to prepare the data for barplot. But I only have the data frame now. x1=rnorm(10,mean=2) x2=rnorm(20,mean=-1) x3=rnorm(15,mean=3) data=data.frame(x1,x2,x3) If there a way to put data within a specific range? The expected result is as follows: range x1 x2 x3 -10-0 2 5 1 (# points in this