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,
+
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, +
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 commands I am using
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
?
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's not ggplot2, just
like 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,
implemented in packages qtbase/qtpaint
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 j...@bitwrit.com.au wrote:
On 01/22/2012 06:56 AM, Kevin Burton wrote:
Of course
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
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
[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 them
using variations
...@yahoo.ca wrote:
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 rkevinbur...@charter.net
To: r-help@r-project.org
Cc:
Sent: Saturday, January
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 lig...@statistik.tu-dortmund.de
wrote:
On 21.01.2012 14:47, Kevin Burton wrote:
I can put a legend on a plot
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/listinfo/r-help
PLEASE do read
, 2011 at 6:50 PM, Kevin Burton rkevinbur...@charter.net
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
have
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
...@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 rkevinbur...@charter.net
wrote:
I admit it isnt reality but I was hoping through judicious use of these
functions I could
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
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 rkevinbur
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 would like
. Michael Weylandt michael.weyla...@gmail.com
[mailto:michael.weyla...@gmail.com]
Sent: Tuesday, November 22, 2011 3:10 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Missing data?
Couldn't you use seq.Date() to set up the time index and then just fill as
appropriate?
Alternatively
is a Sunday and the 10th is a Monday (the beginning of the week).
-Original Message-
From: R. Michael Weylandt michael.weyla...@gmail.com
[mailto:michael.weyla...@gmail.com]
Sent: Tuesday, November 22, 2011 3:10 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Missing data
, weeks, name = name, ...) :
missing values removed from data
-Original Message-
From: R. Michael Weylandt michael.weyla...@gmail.com
[mailto:michael.weyla...@gmail.com]
Sent: Tuesday, November 22, 2011 3:10 PM
To: Kevin Burton
Cc: r-help@r-project.org
Subject: Re: [R] Missing data
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 I don't end up
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 C D
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 digits
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
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, Nov 15
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)
}
}
else
{
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 you. 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)
and
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
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 jwiley.ps...@gmail.com wrote:
Seems like it could work---can you
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:rkevinbur...@charter.net
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
: 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 Statistical
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
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-project.org'
Subject: Upgrade R?
I am trying to upgrade to R 2.14 from R 2.13.1 I have compied all
? 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-project.org'
Subject: Upgrade R?
I am trying to upgrade to R 2.14 from R 2.13.1 I have compied all the
libraries from
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))
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(2001,10))
, 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?
More generally, most people don't like
]
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
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
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))
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.
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
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
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-project.org
48 matches
Mail list logo