similar to: critical value for HEGY test in uroot

Dear R-list, I try to use StructTS for some univariate time series, with a type="BSM". I have two questions: 1-in the "fitted" value of the object, it seems to me that for getting back fitted values, one have to add level, slope and seasonal (or "level" contains already the two other ones) 2-i am really circumspect about the results of "KalmanForecast"

Dear R-list, i have a problem, in the framework of simulations, i want to vectorize for earning time: a variable, say X, has values on intervals and an other variable, say Y, has values on other intervals. For example Inf Sup X 0 2 1 2 4 2 4 6 3 and Inf Sup Y 1 3 1 3 5 2 5 7 3 i want to create a matrix like this Inf Sup

Dear R-list, I try to adjust GLMM on incidence cancer data. Without random effect, in GLM the command is, for example with sex effect, glm(Observed~sex+offset(log(Expected)),family=poisson) because the observed are Poisson distribued with parameter Expected*incidence rate. But know I want to introduce random effect (for example spatial effect) and it seems to me that the "offset" does

Dear R-list, I wish a very happy new year and send you a little question: I have different variables which names are m1, m2, m4, .., m10, ... and want to obtain for example mean of each of them without typing each mean(m1), mean(m2), .... What is the solution for decomposing names in mXX? With my best regards, EA. ============================================ Erik-André SAULEAU SEAIM Hôpital du

Hello R-List I use DClusters package (I work in a cancer regestry). I have 2 questions about it: 1-how is it possible to get back the bootstrap pValue? I mean the pValue of the calculated statistic with respect of the distribution of this statistic under the null hypothesis. 2-how is it possible to test an overdispersion in the poisson model? for choosing a best model I need this mesure of

Hi, i want to create an object which contains different classes: for example i have some time series and test if ARIMA models are best than HoltWinters models: for each of my time serie i want to collect in an unique object which model was the best: for some it will be an HoltWinters class and for some other an Arima class. is there any solution?? thanks in advance, with my best, erik sauleau

Hi, I'd like to use the HEGY test from the uroot package (s. attachment) and get the following error message: error in dimnames(Mypi)[[2]] <- paste("Ypi", 1:s, sep = "") : Object 'Mypi' not found For the air passenger example on http://127.0.0.1:11997/library/uroot/html/HEGY.test.html it works, but for my time series it doesn't (giving names to the

Dear all, I want to use GraspeR (on R 1.6.0) and have to install tcltk library. When I try to do this, I get back this error message: Error in firstlib(which.lib.loc, package) : TCL_LIBRARY is not set Error in library(pkg, character.only = TRUE) : .First.lib failed But I have a tcl/tk package (Active TCL 8.4.1.0) correctly install on my PC (Win95) and in the autoexec.bat the

Dear R-list, I think it's really a newbie question but ... I try to model age-period-cohort models with polynoms in each effect. I have for each level of age and period some cases -k- and persons-years -py-. Models are A_xP_y: log(k_age,period/py_age,period)=f(polynom(x),age)+f(polynom(y),period). For exemple I try to adjust an A2P3 :

Dear R-listers, I have to explain "y", a continous variable, with two variables "t" and "p" which are "factor variables": t is [1:12] and p is [1:18] (a value for y for each combination t*p). OK, then I try to adjust models like: y=a_p + b*t (I mean as many intercepts as levels of p) or y=a_p + b_p*t (as many intercepts and coefficients for t as levels of

Dear All, I want to apply HEGY test for bi-monthly data.Is it possible if I'll do that by editing the code of uroot package. Could you please help me to solve this problem. Thank You [[alternative HTML version deleted]]

Hi I am still having problems with using the stl function, when I read the csv file into R into a file called tkr and use dim(tkr) the result is 132 1 which is fine. When coerce it into a trime series using ts either: tstkr <- ts(t(tkr), deltat=1/12) or tstkr <- ts(c(tkr), deltat=1/12) and use the stl function I get the following error: Error in

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20070220/cb89f312/attachment.pl

There seems to be a bug with var() when the argument is a vector with exactly three values of 1e30 (or close to this value). This does not happen with twice, four (or more) times this value, or another value. > var(rep(1e30, 3)) [1] 2.971056e+28 > var(rep(1.2e30, 3)) [1] 2.971056e+28 > var(rep(0.9e30, 3)) [1] 2.971056e+28 > var(rep(0.8e30, 3)) [1] 0 > var(rep(1e29, 3)) [1] 0 >

Dear friends, I installed the package “uroot” then I wrote library (uroot) and I entered a series entitled extp, in spite of that the problem persite. In short,here are what I wrote and the message that I obtained: > library(uroot) > extp=c(1,3,10,14,12,5,8,12,13,15,9,8,7,10,9,7,3,10,3,11,12,7,4,9,8,10,20,17) > CH.test (wts=extp, frec=c(1,1), f0=1, DetTr=FALSE, ltrunc=NULL) Erreur dans

Hi When I'm looding the library "uroot" something happens with the basic R code that handles numeric classes. Its looks like a bug but I cant isolate the problem. So hopfully one of you smart guys out there can give me some help. This is how the problem looks. Take for example the vector TEST <- c("1","2") which I want to transform it to a numeric vector. I

Hello friends. I'm wondering what happened to package uroot. I worked quite well with older versions of R, however with 2.12.1 version it simply doesn't work. The worst thing is that I look for it in the contributed packages and it simply doesn't appear. I want to develop ADF test with seasonal (centered) dummies, is there any other possible package that contains this test? Thanks

Hi R users! I've come across using kpss tests for time series analysis and i have a question that troubles me since i don't have much experience with time series and the mathematical part underlining it. x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279, 304, 275, 276, 313, 292, 302, 322, 281, 298, 305, 295, 286, 327, 286, 270, 289, 293, 287,

Hi, I got command to execute R code from Mark. but it was showing error.. Error in library(hopach) : there is no package called 'hopach' Does it mean this library hopach is not installed. But I have installed complete Bioconductor package. I am confused. can you help me in this thanks Regards, -- Prateek Singh Final year Bioinformatics(BTech) Vellore Institute Of Technology

Full_Name: M. Edward Borasky Version: 1.2.2 OS: both Windows 2000 and Red Hat Linux 7.1 Submission from: (NULL) (208.130.245.188) According to the R language manual, "%" is the modulus operator. However, on both Windows and Linux, I get a syntax error when I type records <- 100000 records % 1000 I tried wrapping the operands with "as.integer ()" but that didn't help