Dear all R users,
Can anyone please tell me how to add grid lines in any plot in R including
in Histogram, QQ plot etc?
Thanks and regards,
Arun
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
On 8/21/06, Martin Maechler [EMAIL PROTECTED] wrote:
Hi Göran,
GB == Göran Broström [EMAIL PROTECTED]
on Mon, 21 Aug 2006 11:12:49 +0200 writes:
GB A new version, 0.65-1, of glmmML is now on CRAN. It is a major rewrite
GB of the inner structures, so frequent updates (bug
Please see the posting guide: low-level programming questions to R-devel.
On Mon, 21 Aug 2006, Sender wrote:
Hello:
I was hoping to get some advice about how to return a tree (basically a
linked list -- with each node containing a parent, left, and right node
pointers) from a C routine
Arun Kumar Saha wrote:
Dear all R users,
Can anyone please tell me how to add grid lines in any plot in R including
in Histogram, QQ plot etc?
Have you ever typed ?grid before posting?
Uwe Ligges
Thanks and regards,
Arun
[[alternative HTML version deleted]]
Perhaps ?grid will help you.
--
Bjørn-Helge Mevik
__
R-help@stat.math.ethz.ch 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 provide commented, minimal,
Hi all,
I've added a new package 'weaver' to the BioC repository:
http://www.bioconductor.org/packages/1.9/bioc/html/weaver.html
The weaver package provides extensions to the Sweave utilities
included in R's base package. The focus of the extensions is on
caching computationally expensive
I'd like to pick every imbricated five character long subsets from a
vector. I guess there is some efficient way to do this without loops...
Here is a for-loop-version and a model for output:
VECTOR=c(1,4,2,6,5,0,11,10,4,3,6,8,6);
ADDRESSES=c();
for(i in 1:(length(VECTOR)-4)){
Hi
Thanks again. I hope not to waste to much of your time.
I delete some lines of your answer
Each time myfun is run a new environment is created to hold
its local variables. The parent of that environment is e in
this example by construction. So e and the environment that
is temporarily
Finney's method for finding the confidence interval for a ratio of
parameters is quite simple and is probably described in his book
'Probit analysis'. It is also known as Fieller's method so an
RSiteSearch on 'Fieller' might show something useful.
On Mon, Aug 21, 2006 at 08:46:24AM -0700, carlos
embed(VECTOR, 5)[, 5:1]
gives the subsets, so something like
apply(embed(VECTOR, 5)[, 5:1], 1, paste, collapse=)
does the job.
The following is a bit more efficient
ind - 1:(length(VECTOR)-4)
do.call(paste, c(lapply(0:4, function(j) VECTOR[ind+j]), sep=))
but by looking at
VECTOR=c(1,4,2,6,5,0,11,10,4,3,6,8,6)
x - lapply(seq(length(VECTOR)-4),function(z)paste(VECTOR[z:(z+4)],
collapse=''))
unlist(x)
[1] 14265 42650 265011 6501110 5011104 0111043 1110436
104368 43686
On 8/22/06, kone [EMAIL PROTECTED] wrote:
I'd like to pick every imbricated five
Dear Robert,
Thanks for your time.
I have downloaded Rgraphviz (windows binary) from www.bioconductor.org
and put inside R2.3.0 library then i installed from the local zip
its says package 'graph' couldnot be loaded.
Am i doing the installation correctly? Still the new user.
Can you guide me
Hi
Can I get R to handle really big numbers?I am not interested
in more than (say) 10 sig figs, but I would like to deal with numbers
up to, say, 10^1.
If
a - 10^1
b - pi* a
I would like a+b to return 3.1415926e1.
Toy example, illustrating why I can't deal with log(a) and
You also need to install the 'graph' package, i think it is also available
from bioconductor.org. Other packages might be needed as well
Gabor
On Tue, Aug 22, 2006 at 06:42:31PM +0800, j.joshua thomas wrote:
Dear Robert,
Thanks for your time.
I have downloaded Rgraphviz (windows binary)
Thanks!
I have used tons of for- and while-loops (I'm ashamed to reveal these scripts,
but I'm primarily a musician;-)
http://users.utu.fi/attenka/SetTheoryScripts.r), taken some or more cup of
cocoa and mostly been happy ;-) Now I got so many new ways to do these things,
that it takes a
Hi,
I'm having trouble applying the matrix function. I'd like to be able to
create a matrix of vectors filled in by rows, which are not all the same
length, and so I need it to fill in NAs where applicable.
It's easiest to explain with a simple example:
Suppose vec = c(3,4,5). How can I
--- MARK LEEDS [EMAIL PROTECTED] wrote:
these people/experts provide all these packages and
documentation as a FAVOR
and for the fact that they enjoy spreading
knowledge/statistical computing
abilities etc. It's not their job so I think
criticism of the docs and the
fact that they use
hello
i drew a boxplot with:
boxplot(voxit$AGE ~ voxit$A02X)
and the boxes are from left to right:
Ja Leer Nein wn k.A.
levels(voxit$A02X)
[1] Ja Leer Nein wn k.A. 98
there are no entries with 98
but i want:
Ja Nein Leer wn k.A.
how can i change the order of the boxes
bye
thomas
If I start a computation in R, how can I interrupt it?
I' using R 2.1.0.
Thanks for your help,
Lothar
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On 8/22/06, Sergio Martino [EMAIL PROTECTED] wrote:
Hi
Thanks again. I hope not to waste to much of your time.
I delete some lines of your answer
Each time myfun is run a new environment is created to hold
its local variables. The parent of that environment is e in
this example by
The loop method took 195 secs. Just assigning to an answer of the correct
length reduced this to 5 secs. e.g. use
ADDRESSES - character(length(VECTOR)-4)
Moral: don't grow vectors repeatedly.
Other languages (eg. Java) grow the size of the vector independently
of the number of
On 8/22/06, Lothar Schmid [EMAIL PROTECTED] wrote:
If I start a computation in R, how can I interrupt it?
I' using R 2.1.0.
Hi Lothar,
What do you mean by interrupt - cancel completely and
return to the prompt? Or pause, then resume? Your operating
system would be helpful too.
If the former,
Lothar,
Which system do you use? Windows, Linux, Mac?
Stefan
Lothar Schmid schrieb:
If I start a computation in R, how can I interrupt it?
I' using R 2.1.0.
Thanks for your help,
Lothar
__
R-help@stat.math.ethz.ch mailing list
On Tue, 2006-08-22 at 14:27 +0200, Thomas Kuster wrote:
hello
i drew a boxplot with:
boxplot(voxit$AGE ~ voxit$A02X)
and the boxes are from left to right:
Ja Leer Nein wn k.A.
levels(voxit$A02X)
[1] Ja Leer Nein wn k.A. 98
there are no entries with 98
but i want:
Ja Nein Leer
Here are two solutions. seq(length = ...) instead of
just seq(...) is so that v can possibly contain zeros.
# data
v - 3:5
# solution 1 - rbind/lapply
f - function(n) {
s = seq(length = n)
replace(rep(NA, max(v)), s, s)
}
do.call(rbind, lapply(v, f))
# solution 2 - loop
mat -
Hi
f - function(a,n){(1:a)[1:n]}
t(sapply(c(2,3,4,4,4,5,6),f,n=5))
[,1] [,2] [,3] [,4] [,5]
[1,]12 NA NA NA
[2,]123 NA NA
[3,]1234 NA
[4,]1234 NA
[5,]1234 NA
[6,]12345
[7,]1
Gabor makes a good point about seq() vs a:b [a common gotcha
for me].
I'll revise my original function to:
f - function(a,n){(seq(length=a))[1:n]}
t(sapply(c(2,3,4,4,4,5,6,0),f,n=5))
[,1] [,2] [,3] [,4] [,5]
[1,]12 NA NA NA
[2,]123 NA NA
[3,]12
Gabor == Gabor Grothendieck [EMAIL PROTECTED]
on Mon, 21 Aug 2006 21:03:49 -0400 writes:
Gabor It is worthwhile to note that what is being
Gabor illustrated here is aggregating a numeric matrix by a
Gabor factor using the aggregate.default method and, of
Gabor course, a
We had set R2WinBugs to use the burnin as the adaptive phase. I think
this was changed very slightly in the latest version of R2WinBUGS so
that the adaptive phase would equal burnin minus 1. This was to allow
DIC to be calculated. So I think it should work OK now. If you send an
example I
I'm using Linux.
And I'd like just to cancel a running computation, not the entire R prompt.
Lothar
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Thanks for the tip!
unfortunately though conf limits calculated with Fieller and delta methods do
not seem to be in agreement with (and seem to be worse than..) my SPSS
results.. Am i doing something wrong?
thanks a lot in advance for your help!!
An RSiteSearch on 'Fieller' gave me
Dear All,
I'm looking for code that does circular-linear correlations as proposed
by Mardia (1976).
Thank you,
Agnes
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Then Ctrl-C should do the trick.
Sarah
On 8/22/06, Lothar Schmid [EMAIL PROTECTED] wrote:
I'm using Linux.
And I'd like just to cancel a running computation, not the entire R prompt.
Lothar
--
Sarah Goslee
http://www.stringpage.com
__
Hi Folks,
I've encountered something I hadn't been consciously
aware of previously, and I'm wondering what the
explanation might be.
In (on another list) using R to demonstrate the difference
between different contrasts in 'lm' I set up an example
where Y is sampled from three different normal
The 'gmp' package may be of use here, but I'm not sure.
-roger
Robin Hankin wrote:
Hi
Can I get R to handle really big numbers?I am not interested
in more than (say) 10 sig figs, but I would like to deal with numbers
up to, say, 10^1.
If
a - 10^1
b - pi* a
I would
Hello,
Version 2.2 of The R Guide is available for download in
the Contributed Documents section on CRAN. The R Guide
is written for the beginning R user. I use the guide in my
undergraduate probability and math stat sequence, but anyone
with a basic understanding of statistics (who wants to
Its a very nice document. Here is the link:
http://cran.r-project.org/doc/contrib/Owen-TheRGuide.pdf
--
View this message in context:
http://www.nabble.com/new-version-of-%22The-R-Guide%22-available-on-CRAN-tf2146496.html#a5926389
Sent from the R help forum at Nabble.com.
there is
no factor in the dataset but why there is not one and
why a call to another dataset is totally opaque.
The reason is purely historical. The state dataset is about
10 years older than the data.frame concept. At the time the
state.* variables were constructed it was not possible to
Here is another variation using Robin's idea of t(sapply(...))
v - c(3:5, 0)
t(sapply(lapply(v, function(n) seq(length = n)), length-, max(v) ))
# which can be shortened even further for the case where
# there are no zeros in v
v - 3:5
t(sapply(lapply(v, seq), length-, max(v) ))
On 8/22/06,
On Mon, 21 Aug 2006, Damien Moore wrote:
For very large regression problems there is the biglm package (put you
data into a database, read in 500,000 rows at a time, and keep updating
the fit).
thanks. I took a look at biglm and it seems pretty easy to use and,
looking at the source,
Dear marc: thank you for your tip.
cat(x.new, \n)
3\',5\'-cyclic-nucleotide phosphodiesterase activity
here cat is printing on screen.
how can I direct the output to an object.
I cannot do:
y - cat(x.new, \n)
is there any other way.
thanks
srini
Try this:
x -
The software for my recent book
Statistical Analysis and Data Display
Richard M. Heiberger and Burt Holland
http://springeronline.com/0-387-40270-5
is now available as an R package from the book's website.
The package has been submitted to CRAN.
Description: Support software for Statistical
Hello!
Does anybody know how to compare pre-post-effect sizes of different
variables from the same sample with statistical tests? I have the same
problem with areas under the curve (AUC) from ROC-Analysis.
Any recommendations on references or programms are welcome!
Thank you,
Will
Hi again,
Sorry for the noise, but I need to make a correction:
Seth Falcon [EMAIL PROTECTED] writes:
To install:
source(http://bioconductor.org/biocLite.R;)
biocLite(weaver)
At present, the above install sequence will _only_ work if you are
using a development version of R.
If you are
On Tue, 22 Aug 2006, hadley wickham wrote:
The loop method took 195 secs. Just assigning to an answer of the correct
length reduced this to 5 secs. e.g. use
ADDRESSES - character(length(VECTOR)-4)
Moral: don't grow vectors repeatedly.
Other languages (eg. Java) grow the size
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce
our September courses:
(1) Regression Modeling Strategies in R/Splus --- by Prof Frank Harrell
http://xlsolutions-corp.com/Rstats2.htm
*** Washington DC, September 28-29, 2006
Dear R-helpers,
For those of you who wish to use the useful and interesting HH
package on the Mac: I successfully downloaded the package from http://
astro.temple.edu/~rmh/HH/HH_1.4.tar.gz
Safari routinely upacks archives, so after downloading I had to (in
the Terminal)
gzip
Hello!
I am running R-2.3.1-i386-1 on Slackware Linux 10.2. I am a former matlab user,
moving to R. In matlab, via the cftool, I performed nonlinear curve fitting
using the method nonlinear least squares with the Trust-Region algorithm
and not using robust fitting. Is it possible to perform
On Tue, 2006-08-22 at 07:32 -0700, Srinivas Iyyer wrote:
Dear marc: thank you for your tip.
cat(x.new, \n)
3\',5\'-cyclic-nucleotide phosphodiesterase activity
here cat is printing on screen.
how can I direct the output to an object.
I cannot do:
y - cat(x.new, \n)
is
Me and some colleagues are planning to write a textbook together
(Statistics using R) where the target audience for the book is
psychologists and students of psychology.
We thought that it might be a good idea to use a Wiki when writing the
text. Is that a good idea? Does anybody have any
On Tue, 22 Aug 2006, [EMAIL PROTECTED] wrote:
Hi Folks,
I've encountered something I hadn't been consciously
aware of previously, and I'm wondering what the
explanation might be.
Try
contr.helmert(letters[1:3])
[,1] [,2]
a -1 -1
b1 -1
c02
On Mon, 21 Aug 2006, Damien Moore wrote:
For very large regression problems there is the biglm package (put you
data into a database, read in 500,000 rows at a time, and keep updating
the fit).
thanks. I took a look at biglm and it seems pretty easy to use and,
looking at the source,
j.joshua thomas [EMAIL PROTECTED] writes:
Dear Robert,
Thanks for your time.
I have downloaded Rgraphviz (windows binary) from www.bioconductor.org
and put inside R2.3.0 library then i installed from the local zip
its says package 'graph' couldnot be loaded.
Am i doing the installation
Like this:
do.call( paste, c( list(sep=), lapply(1:5,function(x)
VECTOR[x:(length(VECTOR)-5+x)]) ))
[1] 14265 42650 265011 6501110 5011104 0111043 1110436
104368 43686
HTH,
Chuck
On Tue, 22 Aug 2006, kone wrote:
I'd like to pick every imbricated five character long subsets from a
I have a dataframe with the following structure
iddate value
-
122/08/2006 48
124/08/2006 50
128/08/2006 150
130/08/2006 100
101/09/2006 30
211/08/2006 30
222/08/2006 100
228/08/2006 11
2
Is there a way (simple or not) to get the marginal prediction from lme
(in nlme) and/or lmer (in lme4)?
Rick B.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Here is a solution that uses gsub with a negative lookahead perl-style
regexp to do it:
VECTOR - c(1,4,2,6,5,0,11,10,4,3,6,8,6)
e -
([[:digit:]]+),(?=([[:digit:]]+),([[:digit:]]+),([[:digit:]]+),([[:digit:]]+))
out - gsub(e, \\1\\2\\3\\4\\5 , paste(VECTOR, collapse = ,), perl = TRUE)
Dear Sir / Madam
Hi, I have written some code For R that uses for loops to do
2-dimensional grid searches for maximum likelhood combined with
iterated GLS estimation. As can be expected, depending on the szie of
the grid, estimation can take quite some time. However, I have noticed
that the same
Hi,
thanks for this. I'll keep it in mind next time in teaching/referring
someone to R. BTW, before the R-core guys get you ;) Just replace
all places where you use library to refer to a package (see all
comments on the definition of these on r-help/r-devel), e.g.
Page 17: FYI, .GlobalEnv is
Hi there,
I am using SAS Proc NLMIXED to maximize a likelihood with
multivariate normal random effects. An example is the two part random
effects model for repeated measures semi-continous data with a
cluster at 0. I use the model y ~ general(loglike) statement in
Proc NLMIXED, so I can
--- Richard M. Heiberger [EMAIL PROTECTED] wrote:
there is
no factor in the dataset but why there is not one
and
why a call to another dataset is totally opaque.
The reason is purely historical. The state dataset
is about
10 years older than the data.frame concept. At the
time
Dear R users,
The R Project website doesn't seem to have any links devoted to R training.
Are there any R trainers out there?
Thank you for your help!
Jeffrey V. Ludlow
-
[[alternative HTML version deleted]]
Rick,
if by marginal prediction, you mean the prediction without random
effects, then use the level argument. See ?predict.lme or ?fitted.lme
If not then I don't know :)
Cheers
Andrew
On Tue, Aug 22, 2006 at 02:27:06PM -0400, Rick Bilonick wrote:
Is there a way (simple or not) to get the
Dear all,
I would like to know how to run an analysis of covariance in R. For
example, I have a data frame (data) consisting of two second-degree
categorical variables (diagnosis and gender), one continous
independent variable (age) and one continous dependent variable
(response).
I ran a
Dear all,
I got a wierd problem when using lm.ridge() in MASS library. When my X
matrix has few columns, there is no problem. But when my X matrix gets
larger (over 1000 columns), I got the following error:
Error in Xs$v %*% a : non-conformable arguments
In addition: Warning messages:
1: longer
I have always find heaveyweight solutions while use R. If I were rich, I would
have made a great donation to R fundation ...
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
On Wed, 2006-08-23 at 06:43 +1000, Andrew Robinson wrote:
Rick,
if by marginal prediction, you mean the prediction without random
effects, then use the level argument. See ?predict.lme or ?fitted.lme
If not then I don't know :)
Cheers
Andrew
Thanks. I'm familiar with level in
Greetings -
Is there a way to automatically perform what I believe is called rule
extraction (by Quinlan and the machine learning community at least) for
the leaves of trees generated by rpart? I can use path.rpart() to
automatically extract the paths to the leaves, but these can be
On Tue, 22 Aug 2006, Rick Bilonick wrote:
On Wed, 2006-08-23 at 06:43 +1000, Andrew Robinson wrote:
Rick,
if by marginal prediction, you mean the prediction without random
effects, then use the level argument. See ?predict.lme or ?fitted.lme
If not then I don't know :)
Cheers
Andrew
Try this:
# data
DF - structure(list(id = c(1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 4,
4, 4, 4, 4, 4, 4), date = structure(c(8, 9, 10, 11, 3, 7, 8,
10, 4, 1, 2, 3, 8, 9, 10, 11, 3, 5, 6), .Label = c(01/07/2006,
01/08/2006, 01/09/2006, 02/09/2006, 03/09/2006, 06/09/2006,
11/08/2006, 22/08/2006,
Simply using Rtangle(file.Rnw) gives the following error message in R-2.3.1.
It used to work fine with the earlier versions of R.
Rtangle(qtlbim.Rnw)
Error in Rtangle(qtlbim.Rnw) : unused argument(s) ( ...)
RtangleSetup(file=qtlbim.Rnw,output=qtlbim.R)
Writing to file qtlbim.R
Error in
Check out with XLSolutions Corp
[EMAIL PROTECTED]
--- jeffrey ludlow [EMAIL PROTECTED] wrote:
Dear R users,
The R Project website doesn't seem to have any
links devoted to R training. Are there any R
trainers out there?
Thank you for your help!
Jeffrey V. Ludlow
Dear friends,
When i clean my dataset , i met a difficulty
suppose my data set is :
* data-data.frame(x=c(1:5,1,2,3))
data
x
1 1
2 2
3 3
4 4
5 5*
6 1
7 2
8 3
Now i need to add the data which are less than 3.5 at the bottom, not
including the top data, so the results should be :
x
1 1
2 2
3
Hi list,
I've searched in R-help and found some related discussions but still could
not understand this type of error. My own function is pretty complex, so I
would not put it here, but the basic algorithm is like this:
myfun-function(k){
mydata-...#by someway I create a data frame
Mike Wolfgang asks:
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Mike Wolfgang
Sent: Wednesday, 23 August 2006 1:31 PM
To: R-help list
Subject: [R] glm inside one self-defined function
Hi list,
I've searched in R-help and found some related discussions but still
could
Hi All,
I have also similar feelings for R. I really thank each and every one in
the R community for joining together to Create R and Spread R.
thanks
Sayonara With Smile With Warm Regards :-)
G a u r a v Y a d a v
Senior Executive Officer,
Economic Research Surveillance
aov.out - aov(response~diagnosis*gender,data)
Just add it where you think it belongs in the
sequential sum of squares
To adjust the factors for the covariate use
aov.out - aov(response ~ age + diagnosis*gender, data)
To adjust the covariate for the factors
aov.out - aov(response ~
77 matches
Mail list logo