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
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
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’
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
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
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.
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
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
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
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
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
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(
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
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
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
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
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
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,
> 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
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 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
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
> 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
23 matches
Mail list logo