Hi all,
I feel very comfortable with R in Windows, and am trying to switch over to
Linux. I have R 2.8.1 on Linux and it works great (so far), except for
instalilng/updating packages. For example, using Emacs with ESS as a text
editor, I entered:
> update.packages()
Here are some error message
Hi,
I got some problems while was trying to download data from Yahoo using
yahoo.get.hist.quote() function. My script is as follows :
app <- yahoo.get.hist.quote("aapl", start="02/07/09", end="02/07/06",
quote="close")
However I got following error :
trying URL
'http://chart.yahoo.com/table.cs
On Fri, 6 Feb 2009 19:12:09 -0800 (PST) Paul Heinrich Dietrich
wrote:
PHD> I feel very comfortable with R in Windows, and am trying to switch
PHD> over to Linux. I have R 2.8.1 on Linux and it works great (so
PHD> far), except for instalilng/updating packages. For example, using
PHD> Emacs with
On Fri, 6 Feb 2009, Paul Heinrich Dietrich wrote:
Hi all,
I feel very comfortable with R in Windows, and am trying to switch over to
Linux. I have R 2.8.1 on Linux and it works great (so far), except for
I suspect you didn't install it in any of the ways described in the
R-admin manual. An
Pele yahoo.com> writes:
> I have been trying to output all my results (text, plots, etc) into the same
> postscript file as
>
> one document, but have been unable to...Can anyone help me improve my code
> below so that I can
>
> accomplish this? Currently I have to output them separately the
Jason Rupert wrote:
...
Also, do you know how to make the radial rings "black" and a little
thicker? A co-worker requested that the be made a little darker and
easier to see.
Hi Jason,
You can make the radial grid black with grid.col="black", but
unfortunately there is no option for a th
Hi Pele,
>> I have been trying to output all my results (text, plots, etc) into the
>> same postscript file...
Surely I missed something? Is it not simply that you are turning off the
postscript device in the wrong place? At the moment you only seem to be
using the postscript device in the plott
On 06/02/2009 4:48 PM, Paul Warren Simonin wrote:
Hello,
I am attempting to create plots using two continuous variables and
it seems I should be able to use the simple "plot" command. However,
my x-axis values are dates, and I believe this could be the reason I
am receiving error message
On Sat, Feb 7, 2009 at 6:10 AM, Duncan Murdoch wrote:
> On 06/02/2009 4:48 PM, Paul Warren Simonin wrote:
>>
>> Hello,
>> I am attempting to create plots using two continuous variables and it
>> seems I should be able to use the simple "plot" command. However, my x-axis
>> values are dates, an
Hello R users,
Below is the code and output of what I am trying to do. My goal is to
insert/print all items in the chart function into a pdf document. Only the
acf and pacf charts gets printed.
Again, thanks in advance for any help I can get!
options (scipe
I don't see anywhere that you opened a pdf device. When I try :
pdf("test.pdf") and then run your code I get what looks like the
desired output sitting in my working directory:
test.pdf is attached below but my Mac email client sometimes does
unexpected (to me anyway) things to graphics att
I didn't look into the yahoo.get.hist.quote, but the quantmod package has a
working simple solution with getSymbols:
library(quantmod)
getSymbols("AAPL")
> head(AAPL)
AAPL.Open AAPL.High AAPL.Low AAPL.Close AAPL.Volume AAPL.Adjusted
2007-01-03 86.29 86.5881.90 83.80
Dear all, i have two problems
if i have a principal components analysis report as follows:
?princomp
problem one.
> pc.cr <- princomp(USArrests, cor = TRUE)
> pc.cr
i want to export a txt file (*.txt) to save "pc.cr" report
###
Call:
princomp(x = USArr
Hello R-users,
My first consult is that I would like to tackle the relationship between a
dependent variable and four environmental factors using a glm binomial
model. The problem is that I want to test the interaction between factors
too and I would like to learn some way to do it automatically.
Hi guRus,
I'm looking for advice on a good way to approach analysis of some
multi-level data I've obtained. I had humans classify words and
measured response time. Words were 10 positive words ("happy", "joy",
etc) and 10 negative words ("sad","grumpy", etc). Words were also
presented in either wh
Sorry. I now see that you did open a pdf device and attached results.
So your request is for results of summary and AIC to appear in the pdf
output? The plotrix library has functions capable of adding tables or
textboxes to plots. The Sweave and odfWeave packages provide a general
approach
Hi David, Pele:
David Winsemius wrote:
>
>>> I don't see anywhere that you opened a pdf device. When I try :
>
>>> pdf("test.pdf") and then run your code I get what looks like the
>>> desired output sitting in my working directory:
>
> Pele does open a PDF device (previously it was a postsc
Sorry, the message seems to have got botched. Here it is again:
Pele does open a PDF device (previously it was a postscript device). It
looks like what Pele is trying to do is "plot" the printed results of the
summary of a model + the AIC, together with the acf() and pacf() plots.
As Dieter ind
In R using:
formula = (ENF/TOT) ~ (VAR1 + VAR2 + VAR3 + VAR4)^2 would give you all
the main effects and 2 way interactions.
What do you mean by the R^2 for a binary model? It is not constructed
simply by minimizing the squared distances from the 0,1 outcomes.
There are analogues to R^2 (a
On Feb 7, 2009, at 9:07 AM, Nash wrote:
Dear all, i have two problems
if i have a principal components analysis report as follows:
?princomp
problem one.
pc.cr <- princomp(USArrests, cor = TRUE)
pc.cr
i want to export a txt file (*.txt) to save "pc.cr" report
Mike Lawrence thatmike.com> writes:
Thanks for the excellent reproducible sample set!
> I'm most interested in the interaction between color and type, but I
> know that there is likely an effect of word. Yet since word is not
> completely crossed with type, simply adding it to an aov() won't wor
Hi Experts,
I would like to present time series data in meaningful way in building
some graphics. I've tried:
(1) plot(ts(mbaye3))
and
(2) plot(ts(mbaye3), start=1990)
But I always get this error-message:
Fehler [error] in plotts(x = x, y = y, plot.type = plot.type,
Would it improve things if "type" were a continuous variable rather
than categorical? I chose words at the extreme ends of a valence
rating scale but I still have the raw valence ratings for each word.
On Sat, Feb 7, 2009 at 12:02 PM, Dieter Menne
wrote:
> Mike Lawrence thatmike.com> writes:
>
>
I am not an expert, but I will sometimes take a shot at answering a
question when the construction of an example it not too difficult. In
your case that does not occur. I would need to do multiple cut and
pastes and spend a bunch of time doing unnecessary work. I have zoo
and tseries packag
Mike Lawrence thatmike.com> writes:
>
>Would it improve things if "type" were a continuous variable rather
>than categorical? I chose words at the extreme ends of a valence
>rating scale but I still have the raw valence ratings for each word.
>
> >
> > With the interaction, the extreme would be
Hi David,
GOOD POINT! dput() seems to be very useful. I will use it in my future
post!! THX
dput(mbaye3)
structure(c(0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 3, 1, 0, 0,
0,
Each series should be a column, not a row. See ?t
On Sat, Feb 7, 2009 at 11:12 AM, Dominik Hattrup
wrote:
> Hi Experts,
>
> I would like to present time series data in meaningful way in building
> some graphics. I've tried:
>
> (1) plot(ts(mbaye3))
>
> and
>
> (2) plot(ts(mbaye3), start=
Thanks. It did give me the opportunity to discover one reason why
dput() may not be a totally general solution to the problem of
minimizing R-help responder time. That structure depends on an object
that is not included, namely the mbaye dataset that you constructed
mbaye3 from. I was able
And if I decided to ignore the "type" variable altogether and simply
use the continuous "valence" variable, this is what I'd use?
summary(lme(
fixed = rt~valence*color
, data = a
,random = ~1|id
))
I also have continuous luminance measurements of each color that vary
from
Mike Lawrence thatmike.com> writes:
>
> And if I decided to ignore the "type" variable altogether and simply
> use the continuous "valence" variable, this is what I'd use?
>
> summary(lme(
> fixed = rt~valence*color
> , data = a
> ,random = ~1|id
> ))
>
> I also have continuo
Hello all,
I have a *.csv file that looks like this (actual file is orders of magnitude
larger):
site taxa no.ind meadow LMA 2 meadow LCY 1 meadow MSA 2 forest LMA 1
forest LCY 1 forest MSA 1 forest MSX 1
I am interested in, but have failed to create, code that efficiently
converts it to a si
It's working now! Thank You Guys!!
t.mbaye3 <- t(mbaye3)
plot(ts(t.mbaye3))
Now I am searching for an alternative to ylim, which is not working. I
am not the first one :):
https://stat.ethz.ch/pipermail/r-help/2004-October/059376.html
https://stat.ethz.ch/pipermail/r-help/2006-J
Try
library(zoo)
plot(zoo(t.mbaye3), ylim = c(0, 20))
or to put them all on the same plot:
plot(zoo(t.mbaye3), screen = 1, col = 1:10)
On Sat, Feb 7, 2009 at 2:07 PM, Dominik Hattrup
wrote:
> It's working now! Thank You Guys!!
>
>t.mbaye3 <- t(mbaye3)
>plot(ts(t.mbaye3))
>
> No
Perfect! the zoo-library is great. Thank You Gabor!
Dom
Gabor Grothendieck schrieb:
> Try
>
> library(zoo)
> plot(zoo(t.mbaye3), ylim = c(0, 20))
>
> or to put them all on the same plot:
>
> plot(zoo(t.mbaye3), screen = 1, col = 1:10)
>
> On Sat, Feb 7, 2009 at 2:07 PM, Dominik Hattrup
> wro
Hello all,
I have a *.csv file that looks like this (actual file is orders of magnitude
larger):
Site taxa no.ind
forestLMA1
forestLCY1
forestSCO1
meadow
There are various ways, but for simple convenience that's exactly
the job that crosstab() in ecodist was written for.
Sarah
On Sat, Feb 7, 2009 at 2:47 PM, Drew Garey wrote:
> Hello all,
>
>
>
> I have a *.csv file that looks like this (actual file is orders of magnitude
> larger):
>
>
>
> Site
on 02/07/2009 01:47 PM Drew Garey wrote:
> Hello all,
>
>
>
> I have a *.csv file that looks like this (actual file is orders of magnitude
> larger):
>
>
>
> Site taxa no.ind
>
> forestLMA1
>
> forestLCY
Hi, I have a package for which I'm writing a vignette. The vignette
looks fine as indicated by R CMD check. However after installing the
package (in my personal R library location) doing
vignette('rcdk')
says
Warning message:
vignette 'rcdk' *not* found
But I can see the PDF file located u
Hi,
I have a question about ROCR package. I got the ROC curve plotted
without any problem following the manual. However, I don't know to
extract the values, e.g. y.values ( I think it is the area under the
curve auc measure). The return is an object of class "performance"
which have Slots and o
Hi,
Have you tried :
> vignette( package = "rcdk" )
Romain
--
Romain Francois
Independent R Consultant
+33(0) 6 28 91 30 30
http://romainfrancois.blog.free.fr
Rajarshi Guha wrote:
Hi, I have a package for which I'm writing a vignette. The vignette
looks fine as indicated by R CMD check. Ho
Hello,
I am plotting my data as side by side box plots.
Could somebody tell me how I can plot the points from which the box plots are
derived inside the box plots?
I imagine some use of the plot(x,y) function with x = the data and y = location
of box plot, but how do I distribute the points abo
On Feb 7, 2009, at 3:13 PM, Romain Francois wrote:
Hi,
Have you tried :
> vignette( package = "rcdk" )
Yes, I had tried that
> vignette(package='rcdk')
no vignettes found
---
Rajarshi Guha
GPG Fingerprint: D070 5427 CC5B 79
Dear George,
You could use jitter(x) in place of x to randomly perturb the horizontal
coordinates. See ?jitter for details.
I hope this helps,
John
--
John Fox, Professor
Department of Sociology
McMaster University
Hamilton, Ontario, Canada
web: socserv.mcmaster.ca/j
If you post the code that you used, prehapsy using the sample data in
the same package, you may get more useful replies:
"commented, minimal, self-contained, reproducible code."
Gee, ... where have I seen that before?
Also could see the thread;
https://stat.ethz.ch/pipermail/r-help/2009-F
Dear Waverley,
Try this:
unlist(slot(perf,"y.values"))
See ?slot for more details.
HTH,
Jorge
On Sat, Feb 7, 2009 at 3:17 PM, Waverley wrote:
> Hi,
>
> I have a question about ROCR package. I got the ROC curve plotted
> without any problem following the manual. However, I don't know to
> e
This is not really an answer to the original post, but an alternative
possibility that might benefit the R community as well as offering some
commercial revenue to someone with entrepreneurial instincts.
We have several web-based and similar interfaces to R. It should be
possible to have a server
Hi Waverley,
I forgot to tell you that "perf" is your performance object. Here is an
example from the ROCR package:
## computing a simple ROC curve (x-axis: fpr, y-axis: tpr)
library(ROCR)
data(ROCR.simple)
pred <- prediction( ROCR.simple$predictions, ROCR.simple$labels)
perf <- performance(pred,"t
I am new to 'R' and also new to the concept of a 'Hessian' with non-linear
optimization. I would like to avoid going through all of the reference articles
given with ?optim as access to a library is not handy. Would someone be able to
elighten me on what is in the Hessian matrix if 'hessian = TR
Drew, would look at the reshape package
On Sat, Feb 7, 2009 at 2:47 PM, Drew Garey wrote:
> Hello all,
>
>
>
> I have a *.csv file that looks like this (actual file is orders of
> magnitude
> larger):
>
>
>
> Site taxa no.ind
>
> forestLMA
Paul,
On Fri, Feb 6, 2009 at 3:25 PM, Paul Warren Simonin
wrote:
> Hello,
> I am writing regarding log transformation of data in a single matrix
> column, and subsequent use of these data in a glm model fit. I have a data
> matrix in which I am using the log function to transform the values. Th
Quick newbie question please:
I am trying to turn a dataframe into xts with the function;
As.xts
But it returns the error;
Error in as.POSIXlt.character(x, tz, ...) :
character string is not in a standard unambiguous format
could someone give me some pointers please
the
Hi,
can someone help me;
I don't understand why this code doesn't do what it's supposed to do;
x <- c("F+", "F+")
x <- sub("F+", "F", x)
x
[1] "F+" "F+"
(I want "F" "F")
when I try this;
x <- c("F+", "F+")
x <- sub("+", "", x)
x
I get an error message ("invalid regular expression")
Dear Jörg,
Try this:
> gsub('\\+',"",x)
[1] "F" "F"
Take a look at "Basic Regular Expressions" in ?regex.
HTH,
Jorge
On Sat, Feb 7, 2009 at 8:55 PM, Jörg Groß wrote:
> Hi,
>
>
> can someone help me;
> I don't understand why this code doesn't do what it's supposed to do;
>
> x <- c("F+", "F
Hi Group,
Here is some data.
p <- runif(1000) # sample data
groups <- rep(c(1,2),each=500) #conditioning variable
mydata <- cbind(p,groups)
n <- length(p)
u <- (1:n)/(n + 1) # uniform distribution reference for qqplot
logp <- -log(p,base=10)
logu <- -log(u,base=10)
qqplot(logp,logu)
How can I ma
This might seem like a strange question but is there any way to compress an
R object (such as a matrix) and know its resulting size in bytes ?
Clearly, I could implement this in the following way (if x is my matrix):
zz <- gzfile(fname,"w");
write.table(x,zz);
close(zz);
fil
Dear fellow R users,
I read through the "Writing R Extensions" document and am able to now create
my own packages/libraries which so far are just well documented collections
of my own R functions. I use package.skeleton() and the tools package to
build these packages.
However, it is not clear to me
Hi,
This is a resend because the attachment to the previous message was removed.
Sorry.
I wrote a bash script that can either perform unattended upgrades of R as a
cron job or by hand in a terminal. It has run in several environments and
fits my purposes but I'd like to submit it to the commun
What do you want the compressed R object to be? (It is not an R
object.)
Omegahat package Rcompression may help you, but it returns a raw
vector (and that has overheads such as the header: you could use its
length if appropriate).
On Sat, 7 Feb 2009, Markus Loecher wrote:
This might seem
On 2/7/09, Juliet Hannah wrote:
> Hi Group,
>
> Here is some data.
>
> p <- runif(1000) # sample data
> groups <- rep(c(1,2),each=500) #conditioning variable
> mydata <- cbind(p,groups)
> n <- length(p)
> u <- (1:n)/(n + 1) # uniform distribution reference for qqplot
> logp <- -log(p,base=1
59 matches
Mail list logo
