I found the problem. Operator error :(
Thank you.
-Original Message-
From: Kevin Burton [mailto:rkevinbur...@charter.net]
Sent: Tuesday, January 31, 2012 11:53 AM
To: 'Jim Lemon'
Cc: 'r-help@r-project.org'
Subject: RE: [R] Different type of legend?
Sorry this is the
Sorry this is the commands I am using:
colors <- grey(0:10/10)");
leg.txt <- c("-1", "0", "1")
mo <- map('county', projection='polyconic')");
color.legend(mo$range[1]-1,mo$range[2]+1,mo$range[3]-3,mo$range[3]-1," +
"leg.txt[c(1,length(leg.txt)/2
There must be a trick because I am not seeing a legend:
> library(maps)
> library(plotrix)
> colors <- grey(0:10/10)
> mo <- map('county', projection='polyconic')
> leg.txt <- c("-1", "0", "1")
> color.legend(mo$range[1]-1,mo$range[2]-3,mo$range[2]+1,mo$range[2]-2,
+ leg.
12 1:12 PM
To: David Winsemius
Cc: Kevin Burton; r-help@r-project.org
Subject: Re: [R] Interactive map graphics.
Hi Kevin,
I haven't watched through the video yet, but I guess It's not ggplot2, just
like David pointed out, it's actually using a different graphic engine,
imp
stall these packages?
From: Tengfei Yin [mailto:yinteng...@gmail.com]
Sent: Tuesday, January 24, 2012 1:12 PM
To: David Winsemius
Cc: Kevin Burton; r-help@r-project.org
Subject: Re: [R] Interactive map graphics.
Hi Kevin,
I haven't watched through the video yet, but I guess It'
I just watched
http://www.youtube.com/watch?v=iSXNfZESR5I
and there is a section where Hadley Wickham showed demonstrated a county map
of the US that was interactive. This is exactly what I would like. I have
downloaded ggplot2 but still failed to find out how to even read in and plot
some
I have some data where the frequency is heavily weighted on the lower end.
So I have lots of low values with very few higher values. I would like to
find breakpoints that cover the data with as much detail as possible. I find
that if I use hist() to automatically find the breaks for me it finds bre
Thanks for the tip. I will try it out. I am more after a legend of colors than
readable labels. In fact I would prefer getting rid of the text labels
altogether.
On Jan 22, 2012, at 12:46 AM, Jim Lemon wrote:
> On 01/22/2012 06:56 AM, Kevin Burton wrote:
>> Of course by entering &
colors are generated by:
colors <- rainbow(100)
How would I make such a legend?
Thanks for the suggestions. They are tricks definitely worth knowing.
On Jan 21, 2012, at 10:16 AM, Uwe Ligges
wrote:
>
>
> On 21.01.2012 14:47, Kevin Burton wrote:
>> I can put a legend o
:
> Try something like this:
>
> legend('bottom', leg.txt, horiz = TRUE, cex=.75)
>
> A workable example is requested and would have been helpful.
>
>
> ----- Original Message -
> From: Kevin Burton
> To: r-help@r-project.org
> Cc:
> Sent: S
-
From: Ray Brownrigg [mailto:ray.brownr...@ecs.vuw.ac.nz]
Sent: Thursday, January 19, 2012 1:55 PM
To: r-help@r-project.org
Cc: Kevin Burton
Subject: Re: [R] Alaska and Hawaii map data?
On Fri, 20 Jan 2012, Kevin Burton wrote:
> I can plot each county of the contiguous 48 states or all of
I can put a legend on a plot with something like:
legend('bottom', leg.txt, horiz = TRUE, fill = colors)
But what if the arrays leg.txt and colors are too big? I would still like to
provide a legend but to save space I would like to just show small boxes
with the color filled in so it will
ase that listed
FIPS codes for the counties in Alaska and Hawaii. Similar to 'county.fips'.
Thank you.
Kevin Burton
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/list
iginal Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: Sunday, November 27, 2011 4:24 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Missing data?
On Sun, Nov 27, 2011 at 4:08 PM, Kevin Burton
wrote:
> I admit it isnt reality but I was hoping through
I was just trying to be complete. Why is the frequency argument and
attribute available?
-Original Message-
From: R. Michael Weylandt [mailto:michael.weyla...@gmail.com]
Sent: Saturday, November 26, 2011 2:40 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Missing data
11 at 6:50 PM, Kevin Burton
> wrote:
>> Void of any other suggestions this approach makes sense but for my case I
>> think I need to use zoo objects rather than xts. If I sequence the data
>> generally I don't know if there will be 365 days in the year or 366. So I
>> h
Seems to work fine. Thank you.
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: Saturday, November 26, 2011 10:11 AM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Time series merge?
On Sat, Nov 26, 2011 at 10:55 AM, Kevin Burton
wrote:
>
I have two time series
a <- ts(1:10, start=c(1,6), end=c(2,5), frequency=10)
b <- ts(1:5, start=c(2,1), end=c(2,5), frequency=10)
Obviously 'b' is a subset of 'a'. I want a single index value indicating
where that start of 'b' lines up with the start of 'a'. So in this simple
example I wou
rable work to solve the original
problem. If I create a time series and fill in the values that are
appropriate I still could have NA in the series it seems to.weekly has a
problem with NA in the time series:
> ds <- xts(rep(NA,365), seq(from=as.Date("2011-01-01"),
to=as.Date(
e first column should be consistent.
Notice that 10-09 and 10-10 are properly considered different weeks because
the 9th is a Sunday and the 10th is a Monday (the beginning of the week).
-Original Message-
From: R. Michael Weylandt
[mailto:michael.weyla...@gmail.com]
Sent: Tuesday, Novembe
--
From: R. Michael Weylandt
[mailto:michael.weyla...@gmail.com]
Sent: Tuesday, November 22, 2011 3:10 PM
To: Kevin Burton
Cc:
Subject: Re: [R] Missing data?
Couldn't you use seq.Date() to set up the time index and then just fill as
appropriate?
Alternatively, to.weekly if you are starting with
I was wondering what the best approach is for missing data in a time series.
I give an example using xts but I would like to know what seems to be the
"best" method. Say I have
library(xts)
xts.ts <- xts(1:4,as.Date(c("1970-01-01", "1970-1-3", "1980-10-10",
"2007-8-19")), frequency=52)
I w
Say I have the following data:
> s <- list()
> s[["A"]] <- list(name="first", series=ts(rnorm(50), frequency=10,
start=c(2000,1)), category="top")
> s[["B"]] <- list(name="second", series=ts(rnorm(60), frequency=10,
start=c(2000,2)), category="next")
If I use unlist since this is a list of lists
What is wrong with the following?
x <- 1:2
if(x[1] > 0)
{
if(x[2] > 0)
{
print("1 & 2 > 0")
}
else
{
print("1 > 0")
}
}
els
I hadn't considered altering the font. Thank you I will try that.
-Original Message-
From: Sarah Goslee [mailto:sarah.gos...@gmail.com]
Sent: Tuesday, November 15, 2011 1:53 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Plot alignment with mtext
Hi Kevin,
On Tue, N
I would like the text plotted with 'mtext' to be alighned like it is for
printing on the console. Here is what I have:
> print(emt)
ME RMSE MAE
MPE MAPE MASE
original -1.034568e+07 1.097695e+08 2.433160e+07 -31.30554 3
Thank you. I mainly didn't know about the vector/matrix printing rules.
Kevin
-Original Message-
From: William Dunlap [mailto:wdun...@tibco.com]
Sent: Tuesday, November 15, 2011 10:43 AM
To: Kevin Burton; r-help@r-project.org
Subject: RE: [R] Controlling the precision of the d
Has anyone come across the right combinations to print a limited number of
digits? My trial and error approach is taking too much time. Here is what I
have tried:
> op <- options()
> a <- c(1e-10,1,2,3,.5,.25)
> names(a) <- c("A", "B", "C", "D", "E", "F")
> # default
> a
A B
I can get an array of strings for the data that I want using 'paste()' as
follows:
paste('ma', 1:am$arma[2], '=', coef(am)[1:am$arma[2] + am$arma[1]], sep='')
This results in a vector of strings like:
[1] "ma1=1.17760133668255" "ma2=0.649795570407939" "ma3=0.329456750858276"
What I
I am sorry to ask this group but the maintainer of this package did not
leave an email address.
Has anyone used or is using the 'rugarch' package with time-series data
(ts)? I try to fit a GARCH model to my data using the following:
> gf <- ugarchfit(data=l[["MEN"]]$series, spec=spec)
I am also using statConn so I will let you know if I hear anything new.
-Original Message-
From: Cem Girit [mailto:gi...@comcast.net]
Sent: Monday, November 14, 2011 8:52 AM
To: 'Kevin Burton'
Cc: r-help@r-project.org
Subject: RE: [R] Upgrade R?
Hello Kevin,
Thank y
It seems that there is a bug in the forecast::tslm function. I have forwarded
what I think is the bug (the call to get() should supply the argument
'envir=parent.frame()').
Thank you.
On Nov 11, 2011, at 12:11 PM, Joshua Wiley wrote:
> Seems like it could work---can you save your script as a
ve a
version
> of 'R' from my machine? This is a Window PC running Windows 7.
>
>
>
> Thank you.
>
>
>
> Kevin
>
>
>
> From: Kevin Burton [mailto:rkevinbur...@charter.net]
> Sent: Monday, November 07, 2011 2:23 PM
> To: 'r-
ully remove a
version
> of 'R' from my machine? This is a Window PC running Windows 7.
>
>
>
> Thank you.
>
>
>
> Kevin
>
>
>
> From: Kevin Burton [mailto:rkevinbur...@charter.net]
> Sent: Monday, November 07, 2011 2:23 PM
> To: 'r-help@r-pr
I can get multiple plots on a page like:
op <- par(mfcol = c(3, 1))
What I was wondering is if there is a way to have a title for the whole
page? I can specify the title for each individual plot like:
plot(xxx, main=".")
But I would like a 'title' for the group of plots. Is this pos
: Thursday, November 10, 2011 9:17 AM
To: Kevin Burton
Cc: r-help@r-project.org; r-help-boun...@r-project.org
Subject: Re: [R] Upgrade R?
Why don't you just download the latest release from CRAN - R. It is the
recommended approach to installing R.
Steve Friedman Ph. D.
Ecologist / Spatial Statis
Will traceback() work in the error routine specified in tryCatch?
error <- function(e)
{
traceback()
}
tryCatch(..., error=error)
-Original Message-
From: William Dunlap [mailto:wdun...@tibco.com]
Sent: Wednesday, November 09, 2011 5:09 PM
To: Thomas Lumley; rkevinbur...@charter.net
Since apparently there is no one familiar with this error message let me
rephrase the question. Is there a 'manual' process to fully remove a version
of 'R' from my machine? This is a Window PC running Windows 7.
Thank you.
Kevin
From: Kevin Burton [mailto:rkev
]
Sent: Tuesday, November 08, 2011 6:20 PM
To: Kevin Burton
Cc: R. Michael Weylandt; r-help@r-project.org
Subject: Re: [R] window?
The ets() function in the forecast package requires either a numeric vector
or a Time-Series object (produced from ts()). The frequency argument in ts()
refers to the
a...@gmail.com]
Sent: Tuesday, November 08, 2011 2:28 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] window?
I'm not entirely sure that your request makes sense: what do you expect the
frequency to be? It makes sense to me as is...Might your troubles be because
53 is prime?
Mor
This doesn't seem to work:
> d <- rnorm(2*53)
> ds <- ts(d, frequency=53, start=c(2000,10))
> dswin <- window(ds, start=c(2001,1), end=c(2001,10), frequency=1)
> dswin
Time Series:
Start = 2001
End = 2001
Frequency = 1
[1] 1.779409
> dswin <- window(ds, start=c(2001,1), end=c(200
Can someone enlighten me on why the following doesn't work?
setwd('C:/Temp/R')
d <- rep(1:53,2)
(s <- ts(d, frequency=53, start=c(2000,10)))
n <- length(s)
k <- n%/%3
for(i in (n-k):n)
{
st <- c(start(s)[1] + (start(s)[2] + i)%/%frequency(s), (start(s)[2] +
i) %% frequency(s
I am trying to upgrade to R 2.14 from R 2.13.1 I have compied all the
libraries from the 'library' directory in my existing installation (2.13.1)
to the installed R 2.14. Now I want to uninstall the old installation (R
2.13.1) and I get the error:
Internal Error: Cannot find utCompiledCode reco
I started to check what I thought I knew with autocovariance and it doesnt
jive with the the calculations given by R. I was wondering if there is
some scaling or something that I am not aware of.
Take the example
Ø d <- 1:10
Ø (a <- acf(d, type="covariance", demean=FALSE, plot=FALSE))
I was able to solve it by supplying the 'mar' argument. Even when I spell it
right this argument did not solve it.
Thanks again.
Kevin
-Original Message-
From: Eik Vettorazzi [mailto:e.vettora...@uke.de]
Sent: Monday, October 24, 2011 8:46 AM
To: Kevin Burton
Cc: r-help@r-p
This is very basic but I have not been able to find an answer. Basically I
want to find the length of a string.
length("Text")
returns 1 so I know that is not right.
Thank you.
Kevin
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
I get the following error trying to install the package:
> install.packages("grid")
--- Please select a CRAN mirror for use in this session ---
Warning message:
In getDependencies(pkgs, dependencies, available, lib) :
package grid is not available (for R version 2.13.2)
>
-Original Messag
Thank you. This works pretty well. I am having some trouble with the text on
the left-hand side getting cut off. I have tried haling="center" with not
luck. How can I make the text seem wider than it really is to avoid the
truncation. It is only a few characters but still it is annoying. Thank you
48 matches
Mail list logo
