[R] Packages for Nonlinear ARIMA Estimation

2019-09-04 Thread Michael Howell
Hello, I am looking for an R package that uses nonlinear least squares to fit ARIMA models. Initially I was using tseries but according to the documentation for the Arima function in tseries: *The exact likelihood is computed via a state-space representation of the ARIMA process, and the innovatio

Re: [R] creating a line plot for time series data for several years

2019-09-04 Thread Rui Barradas
Hello, Thanks for the reproducible example. All you have to do is to use the aesthetic group = year in either ggplot() or geom_line(). This works: library(ggplot2) daT$date <- as.Date(daT$date) daT$DATE <- format(daT$date, format="%m/%d") Ab <- ggplot(daT, aes(x = as.factor(DATE),

Re: [R] [effects] allEffects does not accept integer value for xlevels

2019-09-04 Thread Fox, John
Dear Rolf, Thanks for trying to help. The bug wasn't in AnalyzeModel(). There was a bug in Effect.lm(), Effect.multinom(), and Effect.polr() in how xlevels=n (e.g., xlevels=4) was handled, now fixed in the development version of the effects package on R-Forge, from which it can be installed via

Re: [R] [effects] allEffects does not accept integer value for xlevels

2019-09-04 Thread Rolf Turner
I'm obviously not understanding something here, but it seems to me that the conjecture It appears to me that the cause is buried in effects:::Analyze.model in or close to the the lines if (is.numeric(xlevels) & length(xlevels) == 1L) { levs <- xlevels for (name in focal.predictors)

Re: [R] Change Y-axis labels

2019-09-04 Thread David Winsemius
On 9/4/19 7:35 AM, Medic wrote: The Y scale is divided (by default) as: 0.0 ... 0.2 ... 0.4 ... 0.6 ... 0.8 ...1.0 But I would like so: 0 ... 20 ... 40 ... 60 ... 80... 100 (with rotating axis labels) When I use par function (marked as comment here) it turns out correctly for ONLY ONE picture?!

[R] creating a line plot for time series data for several years

2019-09-04 Thread Marna Wagley
Hi R Users, I have been getting a trouble to create a time series plot (line) as I was trying to create a line graph for each year using month and date (in x axis). I would like to put only one x axis (month and date) for three years using the facet_grid, but it has not been creating a lineplot for

[R] Change Y-axis labels

2019-09-04 Thread Medic
The Y scale is divided (by default) as: 0.0 ... 0.2 ... 0.4 ... 0.6 ... 0.8 ...1.0 But I would like so: 0 ... 20 ... 40 ... 60 ... 80... 100 (with rotating axis labels) When I use par function (marked as comment here) it turns out correctly for ONLY ONE picture?! Help me, please. (This is the code

Re: [R] [effects] allEffects does not accept integer value for xlevels

2019-09-04 Thread Fox, John
Dear Gerrit, Yes, that appears to be a bug in Effect() -- too bad that it wasn't discovered earlier because a new version of the package was submitted yesterday, but thank you for the bug report. We'll fix the bug, but until then a work-around is to specify the number of levels for each numeri

[R] [effects] allEffects does not accept integer value for xlevels

2019-09-04 Thread Gerrit Eichner
Dear list, citing from allEffects' help page (of package effects 4.1-2): "If xlevels=n is an integer, then each numeric predictor is represented by n equally spaced values rounded to 'nice' numbers." However, adapting the first example from allEffects' help page throws an an error: mod

Re: [R] Durbin- Levinson algorithm

2019-09-04 Thread Stephen Ellison
Google "levionson-Durbin in R The first reference is the levinson function in the signal package. > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Faheem > Jan via R-help > Sent: 04 September 2019 11:29 > To: R-help Mailing List > Subject: [R] Durbin-

[R] Durbin- Levinson algorithm

2019-09-04 Thread Faheem Jan via R-help
hi, i have a problem that how i could use Durbin- levinson algorithm for prediction in case of multiple time series? how i could do this in R... [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and mor

[R] Functional final prediction error

2019-09-04 Thread Faheem Jan via R-help
Hi, my question is related functional data analysis, what is functional final prediction error how we can use to transform vector autoregressive model to its functional form...need help in this regard i will be thanks ful.. [[alternative HTML version deleted]] __

Re: [R] Efficient way to update a survival model

2019-09-04 Thread Frank S .
Chris, thank you so much for your answer!! Best, Frank S. De: Andrews, Chris Enviado: martes, 3 de septiembre de 2019 14:14 Para: Frank S. Cc: r-help@r-project.org Asunto: Re: [R] Efficient way to update a survival model library("survival") set.seed(1) v <- ru

Re: [R] Efficient way to update a survival model

2019-09-04 Thread Frank S .
Charles, thank you for your suggestion! Frank S. De: Berry, Charles Enviado: s�bado, 31 de agosto de 2019 19:21 Para: Frank S. Cc: Andrews, Chris ; r-help@r-project.org Asunto: Re: Efficient way to update a survival model The i^th model is included in the Cox[