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
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
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
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
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
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
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,
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")
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
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
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
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
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
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
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
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 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
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
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
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
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, 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
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
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
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
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
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
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
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
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
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
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
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=
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,
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
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
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:
>
>
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,
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
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
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
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
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. 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 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
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.
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
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
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
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
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
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
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
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
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
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
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
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
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
59 matches
Mail list logo
