If the daya set is not too big, I believe the easiest will be:
stepAIC
from the MASS package...
On Thu, Jul 2, 2009 at 11:21 PM, Angela Radulescu
angela.radule...@gmail.com wrote:
Hi,
I have performed a binary logistic regression using the glm ( family =
binomial(logit) )command on 6
2009/7/2 Peter Dalgaard p.dalga...@biostat.ku.dk:
Mark Wardle wrote:
[...]
Otherwise, it sounds like a recipe for obfuscation! What are you
trying to do really?
Not really obfuscation, there are legitimate reasons to want this sort of
behaviour. See demo(scoping) for a related example
Dear R-User,
Is it possible to simulate Impulsive differential equation using R?
If yes, please let me know how?
Thanking you,
Best Regards,
Mini Ghosh
See the Web#39;s breaking stories, chosen by people like you. Check out
Yahoo! Buzz. http://in.buzz.yahoo.com/
My point is your problem is not entirely related to this package in
particular but more should I trust additional R packages, ie should I
use package cron when R can handle date objects, should I use ggplot
when i could do it all in grid, ...
operators comes with unit tests that allow me to
On Thu, Jul 2, 2009 at 11:57 PM, seeliger.c...@epamail.epa.gov wrote:
Thanks for replying, Francois.
To directly answer your question, the difference between using base R
functions and a library comes down to code correctness and stability, as
well as future support. Portions of the
From yesterday I cannot connect anymore to wiki.r-project.org
Is the problem known?
Thanks!
mario
--
Ing. Mario Valle
Data Analysis and Visualization Group| http://www.cscs.ch/~mvalle
Swiss National Supercomputing Centre (CSCS) | Tel: +41 (91) 610.82.60
v.
Hi all,
Does anyone know if Lo's modified R/S statistic is implemented in R?
I found rsFit function from fArma package (LrdModelling), but I guess this
is to calculate the classical R/S statistic.
Thank you,
M.
--
View this message in context:
MUHC-Research wrote:
I'm currently having trouble with the implementation of a groupedData
object in the lme() function.
Executing the following function
applyScalingSimp - function(input.population)
{
## GA is a time value
varInOrder - c(GA,weight,grouping,sex)
Yes, it is known. It is a major power supply fealure during all night...
and UPS were exhausted. The site is now back again. Sorry for the
inconvenience.
Best,
Philippe Grosjean
Mario Valle wrote:
From yesterday I cannot connect anymore to wiki.r-project.org
Is the problem known?
Thanks!
My guess is that the original
poster is expecting a matrix
not a scalar. If so, then:
Reduce('+', a)
[,1] [,2] [,3]
[1,]300
[2,]030
[3,]003
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of The R Inferno and
Hi,
I'm stuck with a fairly basic re-ordering problem. I want to extract part of a
matrix as a vector, and reorder it to match with an erratic sequence of x,y
coordinates. Here's an example that shows what I want and how close I've got
(but for some reason my mind just can't get the final
Hi,
You could have a look at the sp-package, it provides R-classes for
spatial data. In addition, it also provides very nice plotting
facilities. To read your data into these sp classes you can use the
rgdal package. There is also an R mailing list specifically for
geographic data
Thanks Gabor for this reply. However can please clarify one more thing? How I
want to create a list,
wherein each member of that list is the monthly observations. For example,
1st member of list contains daily observation of 1st month, 2nd member
contains daily observation of 2nd month etc.
barbara.r...@uniroma1.it wrote:
Is there a command as mat.or.vec(nr,nc) to create an array that I
must calculate with more cicle?
Hi Barbara,
The approach that I usually use is something like this:
new.matrix-matrix(NA,nrow=3,ncol=3)
This creates a 3x3 matrix named new.matrix that is
Hi all
On page 39 of this paper [1] by Andrew Lo there is a very interesting
correlation network diagram (sorry I dont have a more convenient link to
the type of diagram I'm talking about). does anyone know of any package in
R that can generate these types of diagrams?
Cheers
-- Rory
[1]
Sorry - In my previous mail, I said page 39, when I actually meant page 36.
Cheers
-- Rory
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Is there an specific function
to eliminate blank space within
strings?
Thanks
--
Dr. Agustin Lobo
Institut de Ciencies de la Terra Jaume Almera (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: agustin.l...@ija.csic.es
http://www.ija.csic.es/gt/obster
Data from Fisher's paper: Confidence Limits for a Cross-Product Ratio.
y
col1 col2
[1,] 103
[2,]2 15
fisher.test(y)
Fisher's Exact Test for Count Data
data: y
p-value = 0.0005367
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence
On 07/03/2009 10:34 AM, Agustin Lobo wrote:
Is there an specific function
to eliminate blank space within
strings?
Thanks
gsub( [[:space:]]+, , string )
--
Romain Francois
Independent R Consultant
+33(0) 6 28 91 30 30
http://romainfrancois.blog.free.fr
|- http://tr.im/qJ8V : RGG#153: stars
How much blank space do you want to remove? all to one, all to none,
leading, trailing?
You can use regular expressions:
x -mytest string
# leading
sub(^ +, , x)
[1] mytest string
# trailing
sub( +$, '', x)
[1]mytest string
# multiple to one
Dear Mailinglist,
I just set up an R 2.9.1 environment with Rcmdr 1.4-6 on Ubuntu Jaunty
9.04. As I'm from Germany my $LANG variable is set to de_DE.UTF-8.
Now, when I open up Rcmdr and try to edit a new datamatrix there is no
edit window appearing:
Datenmatrix - edit(as.data.frame(NULL))
ERROR:
Agustin Lobo wrote:
Is there an specific function
to eliminate blank space within
strings?
No, but gsub could do that.
Duncan Murdoch
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Hi
I have a string 56 Fe [1] that I would like to use as a variable name by
using assign however I think the spaces and brackets might be causing R some
trouble.
How can I change the string so that it just becomes 56Fe1 and can be used as a
variable name.
Thank You
Kurt
On 07/03/2009 12:56 PM, Kurt Smith wrote:
Hi
I have a string 56 Fe [1] that I would like to use as a variable name by using
assign however I think the spaces and brackets might be causing R some trouble.
Not really;
R assign( 56 Fe [1], 10 )
R `56 Fe [1]`
[1] 10
How can I change the
Dear Mailinglist,
I just set up an R 2.9.1 environment with Rcmdr 1.4-6 on Ubuntu Jaunty
9.04. As I'm from Germany my $LANG variable is set to de_DE.UTF-8.
Now, when I open up Rcmdr and try to edit a new datamatrix there is no
edit window appearing:
Datenmatrix - edit(as.data.frame(NULL))
ERROR:
Some suggestions in increasing complexity:
1. Cut-and-paste from the R console. (Cut-and-paste has got to be the
best computer invention ever.)
2. If you want to generate the output automatically from an unattended
script, try help(cat) and help(file)
3. Try help(Sweave) for even more
On 07/03/2009 02:13 PM, Kurt Smith wrote:
Are you sure it would work?
It works when I physically enter 56 Fe [1] but fails when I try to
enter that anything other then directly.
instrument.input[1,6]
[1] 56 Fe [ 1 ]
43 Levels: 10840277.06 109014288.37 11055630.67 11456522.47 ... CPS
Thanks for the help, I really appreciate it.
However, I have one more small issue and I don't really want to post a second
message in such a short space of time, so I will put it up here if someone is
willing to help.
Basically I need to transfer two vectors into a data.frame/matrix (with a
Hi!
Is there some R version for flashdisk, pendrive etc ?
--
Regards,
Raphael Saldanha
saldanha.plan...@gmail.com
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Quicvk and dirty
Insall R to PC and copy the appropriate directory to flash disk.
Regards
Petr
r-help-boun...@r-project.org napsal dne 03.07.2009 14:34:59:
Hi!
Is there some R version for flashdisk, pendrive etc ?
--
Regards,
Raphael Saldanha
saldanha.plan...@gmail.com
R-helpers,
I thank Jonathan Greenberg and David Winsemius for their responses. I
will keep R64.app in mind but I found that by deleting some large
objects that I didn't need I was able to do my computations using R
2.9.1. (This is consistent with Winsemius's experiment on a 10GB
Hi
I have a dataset that looks like this (dput'd below):
head(x, 20)
time status
1 2009-07-02 10:32:37 1
2 2009-07-02 10:32:43 0
3 2009-07-02 10:32:43 1
4 2009-07-02 10:32:44 0
5 2009-07-02 10:32:44 1
6 2009-07-02 10:32:48 0
7 2009-07-02 10:32:48 1
8 2009-07-02 10:32:54 0
9 2009-07-02 10:33:04
Hi Sir,
Thanks for making package available to us. I am facing few problems if
you can give some hints:
Problem-1:
The model summary and residual deviance matched (in the mail below) but
I didn't understand why AIC is still different.
AIC(m1)
[1] 532965
AIC(m1big_longer)
[1] 101442.9
Try the network package from the statnet suite of packages. See
also http://www.jstatsoft.org/v24/i02 with page 27 having similar
graphs. The statnet tutorial in the same issue is recommended. Home
page is http://csde.washington.edu/statnet/.
Allan
On 03/07/09 11:43,
Maybe not address your question directly. There is a package for
correlation visualization in r-forge
(http://r-forge.r-project.org/projects/corrplot/).
2009/7/3 rory.wins...@gmail.com:
Hi all
On page 39 of this paper [1] by Andrew Lo there is a very interesting
correlation network diagram
Hello,
I'd like to say that it's clear when an independent variable can be ruled
out generally speaking; on the other hand in R's AIC with bbmle, if one
finds a better AIC value for a model without the given independent variable,
versus the same model with, can we say that the independent variable
Try this:
d - diff(DF$time)
d[rep(c(FALSE, TRUE), length = length(d))]
Time differences in secs
[1] 0 0 0 10 1 6 23 2 1
d[rep(c(TRUE, FALSE), length = length(d))]
Time differences in secs
[1] 6 1 4 6 0 2 4 8 1 0
On Fri, Jul 3, 2009 at 9:21 AM, rory.wins...@gmail.com wrote:
Hi
I
Hello,
On Fri, Jul 3, 2009 at 2:15 PM, Allan Engelhardtall...@cybaea.com wrote:
1. Cut-and-paste from the R console. (Cut-and-paste has got to be the
best computer invention ever.)
First option is probably interesting when using in Word a monospace
font for R results. For tables, you can also
We were not told the OS, but if this is Windows, see the rw-FAQ
'Can I run R from a CD or USB drive?'
If it is a Unix-alike: the scheme below will not work and you need to
edit R_HOME in the 'R' script (including the copy in R_HOME/bin) if
you relocate the R file tree. That can be a problem
How about export it to a HTML and then copy to MS word.
I have packaged some code to export lm object to HTML in AJS style.
You can have a look at ASAtable
(http://r-forge.r-project.org/R/?group_id=328).
Ronggui
2009/7/3 Suyan Tian st...@mail.rockefeller.edu:
Hi, everyone:
I forget on how to
On Jul 3, 2009, at 9:21 AM, rory.wins...@gmail.com wrote:
Hi
I have a dataset that looks like this (dput'd below):
head(x, 20)
time status
1 2009-07-02 10:32:37 1
2 2009-07-02 10:32:43 0
3 2009-07-02 10:32:43 1
4 2009-07-02 10:32:44 0
5 2009-07-02 10:32:44 1
6 2009-07-02 10:32:48 0
7
Better yet use two graphs on the same page. It is damned hard to read and
understand a two axes display.
--- On Thu, 7/2/09, Mark Knecht markkne...@gmail.com wrote:
From: Mark Knecht markkne...@gmail.com
Subject: Re: [R] Plot two graphs with different ranges in one
To: Rolf Turner
I noticed that Gabor's interpretation of your desires was the opposite
of what mine was, and also that you were thinking tapply might be a
useful strategy. Here is a tapply implementation that lets you choose:
x$wanted - c(0,-diff(x$status))
x$dif.tm - c(diff(x$time),NA)
sum(x[x$wanted==1,
One technique that I use a screen capture that you can paste into the
document. This works nicely since it will also capture the color/size
of the fonts. Works well for the graphics that you want to include.
If you have MS OneNote, you can use this. There is a freeware package
(MWSnap) that I
I often use Dotplot() from the Hmisc package to construct dotplots.
Typically, I plot multiple symbols on each line of a dotplot.
Sometimes these symbols overlap. In these cases, I would like to
eliminate the overlap by slightly adjusting the vertical placement
of the plotted symbols. Is there
Hi,
the igraph package will be able to do the job.
Yours sincerely
Sigbert Klinke
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
Hi.
I have the following two ecdf plots in one graph:
plot(
ecdf(),
do.points=FALSE,
verticals=TRUE,
main=paste(Ecdf of distances ,DIM,sep=),
col=red
);
lines(
ecdf(),
do.points=FALSE,
verticals=TRUE
);
How do I change the color of the resulting graph?
Hi,
Is there an specific function
to eliminate blank space within
strings?
Thanks
gsub( , , txt)
Sigbert
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Dear All,
I did a new installment of R on my office computer.
However, I encounter one technical problem with my R Help. When I write down
?, followed by the function I am interested in, the next message appears in
Help Menu:This program cannot display the webpage. Please, find attached
did you load the package with
library(Rdonlp2)
before running the example code?
On Fri, Jul 3, 2009 at 9:53 AM, Yoonus CAyoonus...@ifmr.ac.in wrote:
Dear Sir,
I installed successfully using the install.packages(Rdonlp2, repos=
http://R-Forge.R-project.org;). I tried on the example which is
On Thu, Jul 2, 2009 at 7:11 PM, Sean Carmody seancarm...@gmail.com wrote:
I have been plotting the same charts using png on a Windows machine
and on a Mac OSX and the quality of the resulting images, particularly
in relation to the fonts, look far superior in the plots produced on
the Mac. Is
Dear Sir,
I installed successfully using the install.packages(Rdonlp2, repos=
http://R-Forge.R-project.org;). I tried on the example which is defined on
the Rdonlp help file ,p - c(-10,10)par.l - c(0,0); par.u -
c(100,100)lin.u - 1; lin.l - 1A - t(c(1,1))fn - function(x){
x[1]^2+x[2]^2}ret -
Hallo,
How can I have the control of the width of each bin in the histogram?
I would like to plot my data modifying the number and the width of the bins.
Is there a good pdf or manual where i can learn to use the plot/hist in R
very well?
I find the help not always clear
Thanks a lot
Ale
--
Dear List,
My data consist of nine columns and about 50,000 rows. It looks like
this.
-9.0225 3.46464 2.80926 -0.3847 3.73735 1.1058 -2.989361.38901 -8.1846
-2.4315 -5.1189 1.8225 3.3798 1.7874 4.693 -3.9286 1.4266 5.7849
-3.4894 -4.0305 3.7879 3.5195 2.9186 2.8685
Check out the R-Graphviz package on Bio-Conductor.
On Fri, Jul 3, 2009 at 5:43 AM, rajesh jakshay.raj...@gmail.com wrote:
Hi,
I need to draw something like the attached Image.Where can I do this in R?
--
Rajesh.J
__
R-help@r-project.org mailing
aledanda wrote:
Hallo,
How can I have the control of the width of each bin in the histogram?
I would like to plot my data modifying the number and the width of the bins.
Is there a good pdf or manual where i can learn to use the plot/hist in R
very well?
See
?hist
and its argument breaks.
Please provide a sensible subject line as you are asked in the posting
guide which you obviously have not read before posting.
Please also specify your OS, the R version and in this case which help
system your are going to use. HTML? Plain text? Compiled html (if under
Windows)? PDF?
My
On 7/3/09, aledanda danda.ga...@gmail.com wrote:
Hallo,
How can I have the control of the width of each bin in the histogram?
I would like to plot my data modifying the number and the width of the bins.
Is there a good pdf or manual where i can learn to use the plot/hist in R
very
Lars Bergemann wrote:
Hi.
I have the following two ecdf plots in one graph:
plot(
ecdf(),
do.points=FALSE,
verticals=TRUE,
main=paste(Ecdf of distances ,DIM,sep=),
col=red
);
lines(
ecdf(),
do.points=FALSE,
verticals=TRUE
);
How do I change the color
On Jul 3, 2009, at 11:56 AM, Pooja Jain wrote:
Dear List,
My data consist of nine columns and about 50,000 rows. It looks like
this.
-9.0225 3.46464 2.80926 -0.3847 3.73735 1.1058 -2.98936
1.38901 -8.1846
-2.4315 -5.1189 1.8225 3.3798 1.7874 4.693 -3.9286 1.4266 5.7849
-3.4894
Hi David,
Sorry for not being absolutely clear.
x-axis will be the values ranging from (-50 to 50) and y-axis will be
the density distribution as you would get after plotting every column
fitted to the normal distribution. However, what I need here is to
reduce the nine column file to a
Reading in your data
==
xx - t(read.table(textConnection(x), header=FALSE, as.is=TRUE)); xx
mm - apply(xx, 1, mean)
mhi - apply(xx,1, max)
mlow - apply(xx, 1, min)
plot(mm, ylim=c(min(xx)-1, max(xx)+1))
points(mhi, col=red)
Hi,
I have two questions regarding the cloud function in the lattice
package:
1) Is there a way to not print the surrounding frame (i.e. the square
surrounding the entire plot)?
2) Is there a way to italicize the text displayed with the key argument?
Some sample code:
data(iris)
On Jul 3, 2009, at 12:52 PM, Pooja Jain wrote:
Hi David,
Sorry for not being absolutely clear.
x-axis will be the values ranging from (-50 to 50) and y-axis will
be the density distribution as you would get after plotting every
column fitted to the normal distribution. However, what I
On Jul 3, 2009, at 1:02 PM, Martin Eklund wrote:
data(iris)
cloud(Sepal.Length~Petal.Length*Petal.Width,data=iris,
groups=Species,screen=list(z=20,x=-70), col=c(red,blue,green),
pch=c(20,20,20),
#try adding
par.box=list(col=white),
perspective=FALSE,
Hi,
Two easy questions I'm sure.
1) As an example if I use the code
require(zoo)
then once it's loaded is there a command that lists everything that
zoo provides so that I can study the package?
Certainly help(zoo) gives me some clues about what zoo does but I'd
like a list. Maybe
R.version
_
platform i486-pc-linux-gnu
arch i486
os linux-gnu
system i486, linux-gnu
status
major 2
Dear R-Helpers,
I need to compute cross-correlation on two signals wich may contain
missing values.
One cannot pass Na.action=na.pass to the ccf() function.
So, I wrote a naive function of my own (see below).
Unsurprisingly, this function is not very fast.
Do you think that it is
library(help = zoo)
will list all the help files and, at the bottom, all the vignettes.
It does not list demos so to get those you have to issue a
second command:
demo(package = zoo)
On Fri, Jul 3, 2009 at 1:21 PM, Mark Knechtmarkkne...@gmail.com wrote:
Hi,
Two easy questions I'm sure.
1)
Does not work for me with rnorm(20) as an argument. (And supplying a
dummy value for DIM.) The error after fixing the missing DIM issue was:
Error in plot.stepfun(x, ..., ylab = ylab, verticals = verticals, pch
= pch) :
argument 4 matches multiple formal arguments
The help page for
On Fri, Jul 3, 2009 at 10:02 AM, Martin
Eklundmartin.ekl...@farmbio.uu.se wrote:
Hi,
I have two questions regarding the cloud function in the lattice
package:
1) Is there a way to not print the surrounding frame (i.e. the square
surrounding the entire plot)?
There is an example in the
Hi , I have created a parallel coordinate plot using rrgobi. I like to make
a powerpoint slide. I am not able to copy and paste. Is there any way to do
this? Thanks for any help.
Chetty
--
Professor of Family Medicine
Boston University
Tel: 617-414-6221, Fax:617-414-3345
emails:
On 7/3/2009 1:21 PM, Mark Knecht wrote:
Hi,
Two easy questions I'm sure.
1) As an example if I use the code
require(zoo)
then once it's loaded is there a command that lists everything that
zoo provides so that I can study the package?
ls(package:zoo) will list all the exported items in
On Jul 3, 2009, at 1:21 PM, Mark Knecht wrote:
Hi,
Two easy questions I'm sure.
snipped as already answered
2) Related to the above, how do I tell what packages are currently
loaded at any given time so that I don't waste time loading things
that are already loaded? search() tells me
To see all functions in a package:
library(zoo)
ls(package:zoo)
On Fri, Jul 3, 2009 at 2:21 PM, Mark Knecht markkne...@gmail.com wrote:
Hi,
Two easy questions I'm sure.
1) As an example if I use the code
require(zoo)
then once it's loaded is there a command that lists everything that
Thanks to everyone for your answers. They were great and I think I
have what I was looking for.
Cheers,
Mark
On 7/3/09, Duncan Murdoch murd...@stats.uwo.ca wrote:
On 7/3/2009 1:21 PM, Mark Knecht wrote:
Hi,
Two easy questions I'm sure.
1) As an example if I use the code
On Jul 3, 2009, at 1:35 PM, Veerappa Chetty wrote:
Hi , I have created a parallel coordinate plot using rrgobi. I like
to make
a powerpoint slide. I am not able to copy and paste. Is there any
way to do
this? Thanks for any help.
Save it (specifics may depend on OS) and then in PPT (this
Duncan Murdoch-2 wrote:
The Windows GUI was written before the Mac GUI, and was written using a
fairly strange toolkit, which are two reasons the Mac looks prettier.
But the Windows GUI has a few advantages over the Mac:
(snip)
You can copy from the console, prompts and all, and
I've never tried it but I would have thought the package DescribeDisplay can
help you,
http://www.ggobi.org/describe-display/
baptiste
2009/7/3 David Winsemius dwinsem...@comcast.net
On Jul 3, 2009, at 1:35 PM, Veerappa Chetty wrote:
Hi , I have created a parallel coordinate plot using
On Fri, Jul 3, 2009 at 2:14 PM, cls59ch...@sharpsteen.net wrote:
Duncan Murdoch-2 wrote:
The Windows GUI was written before the Mac GUI, and was written using a
fairly strange toolkit, which are two reasons the Mac looks prettier.
But the Windows GUI has a few advantages over the Mac:
On 7/3/2009 2:14 PM, cls59 wrote:
Duncan Murdoch-2 wrote:
The Windows GUI was written before the Mac GUI, and was written using a
fairly strange toolkit, which are two reasons the Mac looks prettier.
But the Windows GUI has a few advantages over the Mac:
(snip)
You can copy from the
Hi,
how can I get a more descriptive text
instead of the variable names in my XY-lattice plot,
according to the table below?
Variabletext
acet = Acetylaspartate Thalamus
chol = Choline Thalamus
acetp= Acetylaspartate parieoc
ino = Inositole Thalamus
I could not find a
On Fri, Jul 3, 2009 at 3:07 PM, Duncan Murdochmurd...@stats.uwo.ca wrote:
On 7/3/2009 2:14 PM, cls59 wrote:
Duncan Murdoch-2 wrote:
The Windows GUI was written before the Mac GUI, and was written using a
fairly strange toolkit, which are two reasons the Mac looks prettier. But
the Windows
Does this do what you want?
b.test - 3
f - function(a,b) {
+ a+b.test
+ }
f(10,20)
[1] 13
environment(f) - baseenv()
f(10,20)
Error in f(10, 20) : object 'b.test' not found
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
Duncan Murdoch-2 wrote:
cls59 wrote:
I recently noticed that the Mac GUI contains a way to do block
selection- so
you can select just the code you want to copy and exclude the prompts.
Just
hold down the option key and the text cursor will turn into a crosshair
which allows you to
As for the editor debate...
I'm not touching that one with a ten hundred foot pole :p
You don't have to: in emacs, C-x M-p will do it for you.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
?strip.custom
p -
xyplot(acet+chol+ino+acetp ~ zp,
group=grp,
data=data,
type=l,
scales=list(relation=free),
auto.key=list(title=
Neurotransmitters, border=TRUE))
update(p, strip=strip.custom(factor.levels=letters[1:4]))
HTH,
baptiste
2009/7/3
Hi all,
OK, I'm struggling a bit with getting good basic plots so I put
together this little example. It works - sort of - but it seems a bit
hokey to me. Can someone suggest how a pro would do it?
The goal is to take data from 4 columns of a data.frame (data read
in by read.csv) and
On Fri, Jul 3, 2009 at 1:02 PM, baptiste
auguiebaptiste.aug...@googlemail.com wrote:
?strip.custom
p -
xyplot(acet+chol+ino+acetp ~ zp,
group=grp,
data=data,
type=l,
scales=list(relation=free),
auto.key=list(title=
Neurotransmitters, border=TRUE))
update(p,
On Fri, Jul 3, 2009 at 1:38 PM, Dylan
Beaudettedylan.beaude...@gmail.com wrote:
Hi,
how about something like this:
a - 1:10
b - cumsum(a)
c - cumsum(b)
d - cumsum(c)
X- data.frame(a,b,c,d)
plot(b ~ a, data=X, type=l, col=blue, ylim=c(0,max(X)))
lines(c ~ a, data=X, col=green)
Hi,
how about something like this:
a - 1:10
b - cumsum(a)
c - cumsum(b)
d - cumsum(c)
X- data.frame(a,b,c,d)
plot(b ~ a, data=X, type=l, col=blue, ylim=c(0,max(X)))
lines(c ~ a, data=X, col=green)
lines(d ~ a, data=X, col=red)
legend('topleft', legend=c('a', 'b', 'c'), col=c('blue', 'green',
Thanks a lot
Baptiste and Deepayan!
That was very helpful-things can be so easy...
-Udo
Quoting Deepayan Sarkar deepayan.sar...@gmail.com:
On Fri, Jul 3, 2009 at 1:02 PM, baptiste
auguiebaptiste.aug...@googlemail.com wrote:
?strip.custom
p -
xyplot(acet+chol+ino+acetp ~ zp,
I have hundreds of megabytes of price data time series, and perl
scripts that extract it to tab delimited files (I have C++ programs
that must analyse this data too, so I get Perl to extract it rather
than have multiple connections to the DB).
I can read the data into an R object without any
Thanks for the reponse Winston. Your chart
(http://stdout.org/~winston/X/r-antialias/pch.html) illustrates my
problem perfectly. I did experiment with Cairo, with no success in
improving the display. In the blog post you linked to, the following
comment appeared:
If you don't have anti-aliasing
Hi Rory,
There are several. Have a look at the gR Task Views. There you will also
find a link to the statnet suite, where you will find links to a dedicated
set of jstatsoft articles.
Regards, Mark.
Rory Winston wrote:
Hi all
On page 39 of this paper [1] by Andrew Lo there is a very
On Fri, Jul 3, 2009 at 2:48 PM, Ted Byersr.ted.by...@gmail.com wrote:
I have hundreds of megabytes of price data time series, and perl
scripts that extract it to tab delimited files (I have C++ programs
that must analyse this data too, so I get Perl to extract it rather
than have multiple
alexander russell-2 wrote:
Hello,
I'd like to say that it's clear when an independent variable can be ruled
out generally speaking; on the other hand in R's AIC with bbmle, if one
finds a better AIC value for a model without the given independent
variable,
versus the same model with,
Dear experts,
I'm attempting to solve a constrained optimization problem using the Rdonlp2
package.
I created a Lagrange function (L=f(x)-lambda(g(x)-c)), where x is a vector
of 16 parameters. This is what I'm using as objective function in the code
below. In addition, I set bounds on these
JoK LoQ wrote:
Just a quickly beginner's question.
I wanna find the mean only from the values from a column related to
specific values from another one. Like, theres a 'region' column, i want
the mean of the value on 'profit' column only from south sells from
'region' column
Hi Mark
Thanks for replying.
Here is a short snippet that reproduces the problem:
library(PerformanceAnalytics)
thedata = read.csv(K:\\Work\\SignalTest\\BP.csv, sep = \t, header
= FALSE, na.strings=)
thedata
x = as.timeseries(thedata)
x
table.Drawdowns(thedata,top = 10)
1 - 100 of 120 matches
Mail list logo