> Ok. I made a copy of the arima.r function called earima.r to put in some
> print
> statements. Fair enough.
>
> Now when I run earima, the .Call statements call find the C subroutines.
>
> I know that this should be a really simple fix
Not quite simple, as it's a nasty-looking namespace p
Duncan Murdoch said:
> environment(arima) <- environment(stats::arima)
Doh!, in a word.
That is a really useful trick; same effect as copying and fooling about with
body() to get an object with the right environment, but I'd completely
overlooked the fact that you can assign an environment from
On 10/6/2014 6:36 AM, S Ellison wrote:
Duncan Murdoch said:
environment(arima) <- environment(stats::arima)
Doh!, in a word.
That is a really useful trick; same effect as copying and fooling about with
body() to get an object with the right environment, but I'd completely
overlooked the fact
Yes, you should train 5 different models, or find an outcome that can
combine them together, since caret only accepts a list or vector as
outcome.
On Sun, Oct 5, 2014 at 1:51 PM, Lorenzo Isella
wrote:
> Thanks a lot.
> At this point then I wonder: seen that my response consists of 5
> outcomes f
> On Oct 5, 2014, at 4:51 PM, Lorenzo Isella wrote:
>
> Thanks a lot.
> At this point then I wonder: seen that my response consists of 5
> outcomes for each set of features, should I then train 5 different
> models (one for each of them)?
> Cheers
caret can only model one outcome at a time so
Hi David ,
coef alone won't do as it still contains the description. As per the
example in ?coef shows ,
> x <- 1:5; coef(lm(c(1:3, 7, 6) ~ x))(Intercept) x
-0.7 1.5
Unless I use cat again ,
> x <- 1:5; cat(coef(lm(c(1:3, 7, 6) ~ x))[1])-0.7
Thanks!
Billy
On Mon,
My wife is trying to run some bootstrap model validation. She is using
code to the effect of
library(rms)
residual <- sasxport.get("~/R_boot/residual.xpt")
c1.odat <- somers2(residual$risk.5, residual$caco)["C"]# c - statistic
for original risk estimate in full sample #
fitboot <- fu
I tried to convert a monthly dataset to quarterly by averaging
dt_mz<-zoo(dt_m, dt_m$date_m)
dt_mz2q<-aggregate(dt_mz, as.yearqtr, FUN=mean)
However, the resulting zoo object are just NAs
How can I get a quarterly average dataset?
Thanks!
> dput(dt_mz)
structure(c("Jan 1981", "Feb 1981", "Mar 1
Yes, I have, but the scripts would need some tweaking.
Hadley
On Mon, Oct 6, 2014 at 12:28 PM, Greg Snow <538...@gmail.com> wrote:
> Hadley, have you tried producing the book in other electronic formats
> (other than pdf)? such as epub? I tried and ended up with a file that
> worked, but all the
On Oct 5, 2014, at 6:50 PM, billy am wrote:
> Hi David ,
>
> coef alone won't do as it still contains the description.
It's not a "description", ; in R parlance it's a "name" although that parlance
can be confused since there is also a language class called "name". In this
case it's an attri
On Oct 6, 2014, at 2:50 AM, jpm miao wrote:
> I tried to convert a monthly dataset to quarterly by averaging
>
> dt_mz<-zoo(dt_m, dt_m$date_m)
> dt_mz2q<-aggregate(dt_mz, as.yearqtr, FUN=mean)
>
> However, the resulting zoo object are just NAs
Because you put those month values inside the core
Hello!
I have about the arimax function, please:
If you see the example page, you see the following:
# Exhibit 11.6
air.m1=arimax(log(airmiles),order=c(0,1,1),seasonal=list(order=c(0,1,1),
period=12),xtransf=data.frame(I911=1*(seq(airmiles)==69),
I911=1*(seq(airmiles)==69)),
transfer=list(c(0,0),c(
Hadley, have you tried producing the book in other electronic formats
(other than pdf)? such as epub? I tried and ended up with a file that
worked, but all the example code was missing (which defeats the
convenience of having it on an ebook reader), I did not check if
everything else was there or
Hi all -- I just encountered a behavior that I believe has changed
from previous versions, though I haven't chased back the last version
that behaves as my existing code expects quite yet.
Perhaps this is a bug, though perhaps I'm missing a subtle detail
somewhere in the documentation...
Here's so
New to R. When trying to install packages always get the following message:
> install.packages("HSAUR2")
Installing package(s) into ‘C:/Users/sbueno/Documents/R/win-library/2.14’
(as ‘lib’ is unspecified)
--- Please select a CRAN mirror for use in this session ---
Warning: unable to access index f
Another (easier) kludge is augmenting a call to each par("plt" = ...) call:
par("plt" = some_plt_coordinates); box(lty = 0)
The box(lty = 0) addition makes downstream calls work as expected, but
yeah... this is a kuldge.
-m
On Mon, Oct 6, 2014 at 12:00 PM, Murat Tasan wrote:
> Hi all -- I just
On 06/10/2014 3:17 PM, Santiago Bueno wrote:
New to R. When trying to install packages always get the following message:
> install.packages("HSAUR2")
Installing package(s) into ‘C:/Users/sbueno/Documents/R/win-library/2.14’
(as ‘lib’ is unspecified)
--- Please select a CRAN mirror for use in thi
For me, it output html from rstudio complete with code etc., with no changes
except a tinker to a call to setwd(), downloading some libraries, and no other
intervention. I haven't explored the rmarkdown capabilities before, and I was
just amazed.
> On Oct 6, 2014, at 1:28 PM, Greg Snow <538.
Thanks Pikal -
I'd like to actually get the names of the variables for those indices both
above and below .3/-.3 with the intention of generating a grid of
scatterplots (which I will use ggplot2 for).
I played around with your code but can't seem to get the indice names.
Best,
Edward
On Fri, O
I believe that what is happening is that the clipping region is being
reset when you call box, but not when you call rect. If you insert
the command "par(xpd=NA)" (or TRUE instead of NA) after the plot.new
and use the rect commands then you can see both rectangles (because
this turns the clipping
On Mon, Oct 6, 2014 at 12:29 PM, Duncan Murdoch
wrote:
> On 06/10/2014 3:17 PM, Santiago Bueno wrote:
>>
>> New to R. When trying to install packages always get the following
>> message:
>>
>> > install.packages("HSAUR2")
>> Installing package(s) into ‘C:/Users/sbueno/Documents/R/win-library/2.14’
In package TSA?
You may need to do some studying for yourself, this is complicated stuff.
As I read the help page, the intention is that the transfer= bit is to allow a
_covariate_ to affect the process in an ARMA-like fashion. So c(1,0) would be
AR(1)-like which if I remember correctly corres
This is great...thanks so much!
On Mon, Oct 6, 2014 at 5:56 PM, peter dalgaard wrote:
> In package TSA?
>
> You may need to do some studying for yourself, this is complicated stuff.
>
> As I read the help page, the intention is that the transfer= bit is to
> allow a _covariate_ to affect the pr
23 matches
Mail list logo
