Search Archives

Hello, Does anybody know if (and how) it is possible to obtain standard errors of estimated variances from StructTS? (R 2.10.0). Thank you in advance, Giovanni

consult the help page for || ?"||" -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

I ran an instant experiment... > typeof(0) [1] "double" > typeof(-0) [1] "double" > identical(0, -0) [1] TRUE Best, Giovanni > By the way: > > Are there difference between -0 and 0?

Thank you to all who responded. In the end the simple solution suggested by Peter Dalgaard worked for us: install the package and zip the resulting directory. (So, incidentally, this still works with R 1.8.0) Thanks, Giovanni

> > xs <- lapply(xdata, function(x){(x - mean(x))/sqrt(var(x))}) Incidentally, the function 'scale' does just that. -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630...

> Hello > I am trying to generate data in say 2D "x,y" for a > circle. ^^^^^^^^ I think data are usually collected... GP -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479...

I am using Sweave to produce handouts for teaching. Is there a way of making Sweave keep the comments following the `#' in the code chuncks? Thanks, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas...

Is there a way of getting line numbers in Schunks? Ideally, I would like to have numbers printed every two or five lines. Thank you in advance, Giovanni -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences...

Hello, How can I see a function called "+.dlm"? > methods("+") [1] +.Date +.dlm* +.POSIXt Non-visible functions are asterisked > getAnywhere("+.dlm") Error in grep(pattern, x, ignore.case, extended, value, fixed, useBytes) : invalid regular expression '+\.dlm' Thanks in advance, Giovanni...

Goodmorning, and sorry for the off-topic question. Does anybody know if the data sets used in Harvey's book "Forecasting, structural time series models and the Kalman filter" are available online? Thanks in advance, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at...

I'm pretty sure I've seen some examples of a function printing zero entries in a matrix as dots, but I'm not able to find it now... Any suggestions...? Thanks in advance. (Of course, I might have dreamt...

> I thought all the values are false, all none of them, because there > aren't any that are true: > > any(logical(0)) > [1] FALSE > This is for the same reason why a product over an empty set of factors is...

> > You can live with it, but > > be aware that it is there. My suggestion is to start the optimization > > from several different initial values and compare maximized values of > > the likelihood. Simulated annealing may be used to better explore the...

Thanks to everybody who responded to my query. I got many useful suggestions about books and editors, plus notes and other material online. Summarizing, the books suggested were - Monahan, Numerical Methods of Statistics - Lange, Numerical analysis for statisticians In terms...

Hello, After reading the help for predict.lm and predict.glm, it is not clear to me what are the values returned by predict( ..., type="terms"). Anybody willing to enlighten me? The example provided by Peter Dalgaard in a recent...

Check the apply function: > ?apply Giovanni > i've a matrix say a(30,10,5). how can a function (for example histogram) > be used > to calculate along the first axis, or to create a 10x5 matrix of > histograms? -- __________________________________________________ [ ] [ Giovanni Petris...

Hello! I am wondering if there is a function in R that returns quantiles of a weighted sample, i.e., a set of values, each coming with a different weight. The function quantile() does that only for the case when...

I am creating lists of vectors withing a loop. I also would like to change the dim attribute to the vectors in order to make them matrices. I have tried the following, but it doesn't work... > sim <- c('simMeans...

Has anybody else experienced something like the example below? Any clues about where I could start looking? Thank you in advance, Giovanni > version _ platform sparc-sun-solaris2.7 arch sparc os solaris2.7 system sparc, solaris2.7 status major 1...

Hi all, Using Beamer, in order to highlight a piece of R code I do something like this - note the "\structure" and "\alert" commands: \begin{semiverbatim} > mleOut <- \structure{dlmMLE}(Nile, + parm = c(0.2, 120), # initial values for optimizer + lower...