John wrote:
Thank you, Uwe. I've found a way to do the job by
reading the FAQ 7.21 although it is not giving a
precise explanation to a novice or casual user at
first reading. For example, if you type the first two
But the corresponding help files do so, for sure, and the FAQ 7.21
points you to
Peter thanks for the response.
So the results from hist(mydata) and truehist(mydata, h = .5) are the
same. OK, but the sum of densities or intensities, for case that I gave,
don't sum to 1 but to 2. Look bellow. I have an example where these
density values are also up to 4 and sum to 5 (I have
UweL == Uwe Ligges [EMAIL PROTECTED]
on Sat, 27 Nov 2004 17:09:26 +0100 writes:
UweL Andrew Kniss wrote:
I am attempting to write an R function to aid in time series diagnostics.
The tsdiag() works well, but I would prefer to get a plot with ACF, PACF,
and Ljung-Box
Gregor GORJANC wrote:
Peter thanks for the response.
So the results from hist(mydata) and truehist(mydata, h = .5) are the
same. OK, but the sum of densities or intensities, for case that I gave,
don't sum to 1 but to 2. Look bellow. I have an example where these
density values are also up to 4
JohnF == John Fox [EMAIL PROTECTED]
on Sat, 27 Nov 2004 23:49:08 -0500 writes:
JohnF Dear Seth,
JohnF You don't say which variable is the explanatory
JohnF variable and which is the response, but assuming that
JohnF prob is to be regressed on effect, you can fit
JohnF
Peter,
Yes c(0,23) works.
Many thanks!
Matthew
-Original Message-
From: Peter Dalgaard [mailto:[EMAIL PROTECTED]
Sent: 26 November 2004 16:43
To: Peter Dalgaard
Cc: Matthew Dowle; '[EMAIL PROTECTED]'
Subject: Re: [R] Tcl error - brace in argument?
Peter Dalgaard [EMAIL PROTECTED]
Dear all,
Has anybody ever written some plot / hist functions that would return
PlotML code? [http://ptolemy.eecs.berkeley.edu/]
Regards,
Gregoire
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Dear Uwe,
I must say that I had thanked you for referring me to
the specific and exact FAQ 7.21 and I had solved my
simple problem from it. I alreadly had looked at some
of basic materials like 'An Introduction to R', 'R for
Beginners', 'R Data Import/Export as well as the
FAQ(that is, I know how
Hi,
I am building a 64-bit R 2.0.1 on Solaris 9. The compiler is Sun Studio
8. Make was successful but I have a core dump during a make check. By
the way, this problem also happens on my Solaris 8 machine though I did
not get a core dump. I do not have 64-bit versions of the readline,
tcl/tk,
Gregoire Thomas [EMAIL PROTECTED] writes:
Dear all,
Has anybody ever written some plot / hist functions that would return
PlotML code? [http://ptolemy.eecs.berkeley.edu/]
Not that I know of, and looking at the design of PlotML, it doesn't
look like a nice fit as an R device driver. PlotML
The following works fine with the x11 device, though it
may well be that an initial plot is overwritten. With a pdf
or postscript device, I get two plots, the first of which
still has the red border from having the focus, while the
second is the plot that I want.
library(lattice);
hi all
i'm new to the area of neural networks. i've been reading some
references and seem to understand some of the learning algorithms. i am
very familiar with regression and would just like to see how neural nets
handle this problem so i've been using the nnet package.
i simply want to use a 3
Yes. See ?abline or ?lines
-Original Message-
From: [EMAIL PROTECTED] on behalf of [EMAIL PROTECTED]
Sent: Sun 11/28/2004 8:35 PM
To: [EMAIL PROTECTED]
Cc:
Subject: [R]: Adding a line in the graph of 'plot()'
Dear All,
I wish to inform, especially Italian speaking R-users,
that on CRAN web site is now available a contribute
(in Italian language) about using R in ts analysis.
Any comments would be appreciated.
Best regards,
Vito
=
Diventare costruttori di soluzioni
Became solutions'
Dear all,
I am having trouble plotting a PCA result. The plot doesn't appear!!!
R goes through without any errors but doesn't make a plot appear!!
Could it be wrong window parameters? In this case how do I change them?
I am under red hat 9 with the latest version of R!
Thanks.
Hi
I have a package which was built using R 1.9.1 and everything worked
fine. I recently upgraded to R 2.0.1 and tried to re-install my package
- and I got:
Error in library(mypackage) : 'mypackage' is not a valid package --
installed 2.0.0?
So I tried rebuilding it using my new version of R:
Hi !
I am wondering how to plot data (e.g. f(x,y) ) in a coordinate system spanned
by two non-orthogonal basis vectors (e.g. hexagonal symmetry). The data is
given on an equally spaced grid in theses coords and i would like to do a
contour plot (e.g. with filled.contour).
Thanks for your
Please point me to the documentation explaining why some of the functions
returned
by calling methods() are marked as non-visible and whether there is indeed
no way of
viewing the R code of such functions
thanks
Steve
_
Stay in
non-visible functions are hidden in a namespace. You can view the
code by using getS3method().
-roger
steve houghton wrote:
Please point me to the documentation explaining why some of the
functions returned
by calling methods() are marked as non-visible and whether there is
indeed no way of
I am trying to use the Hmisc function transace to transform predictors
test-cbind(flowstress,pressres,alloy)
xtrans-transace(x,binary=pressres',monotonic='flowstress',
categorical='alloy')
and I am getting the following message¨
Error in ace(x[, -i], x[, i], monotone = im, categorical = ic)
I'd like to make just a couple of points:
R-help is considered by quite a few people to be high-traffic. As such,
many have low appetite for very basic questions. (I wouldn't call them
silly.) In many cases such questions are answered by pointing to a
particular function help page or manual
You mean something like ?methods, which says:
Value:
An object of class 'MethodsFunction', a character vector of
function names with an 'info' attribute. There is a 'print'
method which marks with an asterisk any methods which are not
visible: such functions can be examined
lots of good points from Andy and Uwe deleted
and perhaps the most important reason for the particular socratic form
of teaching on this list are the number of to-be, current, and former
faculty members who feel compelled to teach general solutions to the
problems (reading the FAQ is a rather
Roger == Roger D Peng [EMAIL PROTECTED]
on Mon, 29 Nov 2004 08:48:04 -0500 writes:
Roger non-visible functions are hidden in a namespace.
yes.
Roger You can view the code by using getS3method().
not always {namely when the hidden function is not an S3 method}
getAnywhere() is
Hello to everyone,
I have 2 problems and would be very pleased if anyone can help me:
1) When I use the package systemfit for SUR regressions, I get two
different variance-covariance matrices when I firstly do the SUR
regression (The covariance matrix of the residuals used for
estimation) and
On Monday 29 November 2004 16:42, [EMAIL PROTECTED] wrote:
Hello to everyone,
I have 2 problems and would be very pleased if anyone can help me:
1) When I use the package systemfit for SUR regressions, I get two
different variance-covariance matrices when I firstly do the SUR
regression
Hello!
I would like to know how do I citate R?
I have used it during my Master thesis but I dont know how to citate
during the text and on the references.
Ive looked for it on the web page but only found how to citate the FAQ.
Thank you in advance.
Tatiana Fernandes
Universidade Estadual
Hi William,
the 1st example given in ?screeplot works for me (R 2.0.0 on SuSE Linux 9.0):
(pc.cr - princomp(USArrests, cor = TRUE)) # inappropriate
screeplot(pc.cr)
plot appears
To help you we need more details.
Does normal plotting work? e.g.: plot(rnorm(20),rnorm(20))
Can you plot to a
See the function ?citation under 2.0.0
R citation()
To cite R in publications use:
R Development Core Team (2004). R: A language and environment for
statistical computing. R Foundation for Statistical Computing,
Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org.
HTH,
Did you see item 2.8 in the R FAQ?
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Citing-R
Tatiana Fernandes wrote:
Hello!
I would like to know how do I citate R?
I have used it during my Master thesis but I dont know how to citate
during the text and on the references.
Ive looked for it on the
Tatiana Fernandes wrote:
Hello!
I would like to know how do I citate R?
I have used it during my Master thesis but I dont know how to citate
during the text and on the references.
Ive looked for it on the web page but only found how to citate the FAQ.
Please the FAQ more carefully and find how
citation()
Hello!
I would like to know how do I citate R?
I have used it during my Master thesis but I dont know how to citate
during the text and on the references.
Ive looked for it on the web page but only found how to citate the FAQ.
Thank you in advance.
Tatiana Fernandes
michael watson (IAH-C) wrote:
Hi
I have a package which was built using R 1.9.1 and everything worked
fine. I recently upgraded to R 2.0.1 and tried to re-install my package
- and I got:
Error in library(mypackage) : 'mypackage' is not a valid package --
installed 2.0.0?
So I tried rebuilding it
Tatiana Fernandes [EMAIL PROTECTED] writes:
Hello!
I would like to know how do I citate R?
I have used it during my Master thesis but I dont know how to citate
during the text and on the references.
Ive looked for it on the web page but only found how to citate the FAQ.
R is a
From: Chuck Cleland
Did you see item 2.8 in the R FAQ?
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Citing-R
Or the start-up message, which has, in part:
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or
HI
the scan function when only one item is read says:
scan()
1: 3.3
2:
Read 1 items
[1] 3.3
I hope my english is not playing a trick on me, but 1 items sounds
very strange
this makes me feel very anal, and it's really not important and I
apologize for it but here it goes anyway
On Thu, 25 Nov 2004, Peter Dalgaard wrote:
I haven't got all that much experience with it, but obviously, the
various algorithms for constrained optimization (box- or otherwise) at
least allow you to find a proper maximum likelihood estimator.
It's harder than it looks (well, my experience is
Has anybody ever written some plot / hist functions that would return
PlotML code? [http://ptolemy.eecs.berkeley.edu/]
http://ptolemy.eecs.berkeley.edu/%5D
Not that I know of, and looking at the design of PlotML, it doesn't
look like a nice fit as an R device driver. PlotML works at the level
Thomas Lumley [EMAIL PROTECTED] writes:
Also, the algorithm in glm.fit, while not perfect, is a little smarter
than a simple IRLS. It uses step-halving to back away from the edge,
and when the parameter space is convex it has a reasonable chance of
creeping along the boundary to the true MLE.
On Sun, 28 Nov 2004, HOME - Didier Ledoux wrote:
Je suis débutant en R. Je voudrais faire un tableau de statistiques
descriptives ou les moyennes seraient calculées en fonction de deux critères:
sexe, région
Cela donnerait ceci:
SEXEREGIONMOYENNE
femmeAmoy1
femmeB
How can I get a single backslash in a character string?
My goal is to escape dots in a string that will be used as a regular
expression. I thought I could do it this way:
gsub(., \\., x)
Unfortunately, \\ does not represent a literal backslash as I
expected, but rather a pair of backslashes:
On Mon, 29 Nov 2004, steve houghton wrote:
Please point me to the documentation explaining why some of the functions
returned
by calling methods() are marked as non-visible and whether there is indeed
no way of
viewing the R code of such functions
Luke Tierney's article in Volume 3 No 1 of the
Hi there.
I just upgraded to 2.01 on Mac OS 10.3.6. I used to use the command
(on R 1.9.x):
world(ylim=c(-30,30), xlim = c(0,360), shift=TRUE, add=TRUE)
to draw a world outline over my image plots. My data uses longitude
from (0, 360) so I need to use the shift function. After I upgraded,
Dan Lipsitt [EMAIL PROTECTED] writes:
How can I get a single backslash in a character string?
My goal is to escape dots in a string that will be used as a regular
expression. I thought I could do it this way:
gsub(., \\., x)
Unfortunately, \\ does not represent a literal backslash as I
[EMAIL PROTECTED] wrote:
I am trying to use the Hmisc function transace to transform predictors
test-cbind(flowstress,pressres,alloy)
xtrans-transace(x,binary=pressres',monotonic='flowstress', categorical='alloy')
and I am getting the following message¨
Error in ace(x[, -i], x[, i], monotone = im,
Hello,
I have a df, pp, with five variables:
nobs(pp)
q10_1 q10_2 q10_3 q10_4 actcode
16201620162016201620
I want to create a loop to run four xtabs (the first four variables above by
the fifth) and then store the results in a matrix. Below I make my intent
Dear group,
I am using justRMA in bioconductor to get expression
values. Now I computed fold changes. the fold
changes i get are log2 values.
How can I convert these to log10 values to see them as
actual fold change values.
thank you.
__
[EMAIL
Ah, I see. Thanks. ?print.default and ?cat do not mention this.
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Multiply by log(2)/log(10); e.g.,
log2(1024) * log(2)/log(10)
[1] 3.0103
log(1024, 10)
[1] 3.0103
Andy
From: Srinivas Iyyer
Dear group,
I am using justRMA in bioconductor to get expression
values. Now I computed fold changes. the fold
changes i get are log2 values.
How can I
Srinivas Iyyer [EMAIL PROTECTED] writes:
Dear group,
I am using justRMA in bioconductor to get expression
values. Now I computed fold changes. the fold
changes i get are log2 values.
How can I convert these to log10 values to see them as
actual fold change values.
Divide by
Does this do what you want?
foo.df - data.frame(x = rnorm(12), y = runif(12), z = factor(rep(1:3,4)))
bar.mat - matrix(NA, nrow = ncol(foo.df)-1, ncol = nlevels(foo.df$z))
for(i in 1:(ncol(foo.df)-1))
{
bar.mat[i,] - xtabs(foo.df[,i] ~ foo.df$z)
}
bar.mat
There's probably a slicker way with
Srinivas Iyyer wrote:
Dear group,
I am using justRMA in bioconductor to get expression
values. Now I computed fold changes. the fold
changes i get are log2 values.
How can I convert these to log10 values to see them as
actual fold change values.
First, a bit of etiquette. It is considered
Hi
Deepayan Sarkar wrote:
On Monday 29 November 2004 04:46, John Maindonald wrote:
The following works fine with the x11 device, though it
may well be that an initial plot is overwritten. With a pdf
or postscript device, I get two plots, the first of which
still has the red border from having the
Hi
I am trying to tune an svm by doing the following:
tune(svm, similarity ~., data = training, degree = 2^(1:2), gamma =
2^(-1:1), coef0 = 2^(-1:1), cost = 2^(2:4), type = polynomial)
but I am getting
Error in svm.default(x, y, scale = scale, ...) :
wrong type specification!
I have it working now, I think. Since it's going into a regular
expression, I have to escape each of the escape characters, resulting
in four backslashes altogether:
sub([.], x, a.b)
[1] axb
sub([.], \., a.b)
[1] a.b
sub([.], \\., a.b)
[1] a.b
sub([.], \\\., a.b)
[1] a.b
sub([.], ., a.b)
A.J. Rossini wrote:
and perhaps the most important reason for the particular socratic form
of teaching on this list...
Golly, anyone who read Plato's Dialogues would realize that the Socratic
method involves patiently leading the questioner stepwise through the
solution, not simply writing
I just upgraded to 2.01 on Mac OS 10.3.6. I used to use the command
(on R 1.9.x):
world(ylim=c(-30,30), xlim = c(0,360), shift=TRUE, add=TRUE)
to draw a world outline over my image plots. My data uses longitude
from (0, 360) so I need to use the shift function. After I upgraded, I
Hi
I am trying to tune an svm by doing the following:
tune(svm, similarity ~., data = training, degree = 2^(1:2), gamma =
2^(-1:1), coef0 = 2^(-1:1), cost = 2^(2:4), type = polynomial)
but I am getting
Error in svm.default(x, y, scale = scale, ...) :
wrong type specification!
Sending a help request once is enough! Even if you don't have an answer
after 1 hour.
I am trying to tune an svm by doing the following:
tune(svm, similarity ~., data = training, degree = 2^(1:2), gamma =
2^(-1:1), coef0 = 2^(-1:1), cost = 2^(2:4), type = polynomial)
I think you want to
Dan Lipsitt danlipsitt at gmail.com writes:
:
: I have it working now, I think. Since it's going into a regular
: expression, I have to escape each of the escape characters, resulting
: in four backslashes altogether:
:
: sub([.], x, a.b)
: [1] axb
: sub([.], \., a.b)
: [1] a.b
: sub([.],
Dear R-users,
I state that this is for beginners, so you may ignore
this in order not to be irritated.
By the way, patience is another important thing,
together with kindness, we should keep in mind when
we teach students and our own children as Jim Lemon
pointed out well in the context of the
You may be missing something. After you create all those objects,
you'll want to use them. Use get():
for (i in 1:10) ... get(paste(object,i,sep=)) ...
It took me about a week to find out how to do this. I waited for a
few days, but before I got to ask this basic/rtfm question, someone
else -
Dear members,
I am wondering if there are any functions to perform the Kernel Fisher
Discriminant method, especially for multi-class problems, in R. I would
appreciate any kind of information on this.
Thanks for your time.
Seungho Huh, Ph.D.
Research Statistician
RTI International
It was enough for me to use the 'assign' function
alone. But I'll remember the 'get' function for future
reference. Thanks a lot for the note.
John
--- bogdan romocea [EMAIL PROTECTED] wrote:
You may be missing something. After you create all
those objects,
you'll want to use them. Use
Hi All:
This may turn out to be very simply, but I can't seem to add the name of
the school to a chart. The loop I created is below that subsets a
dataframe and creates a chart for each school based on certain
variables. As it stands now, they title includes the school's ID number.
Instead, I
Doran, Harold [EMAIL PROTECTED] writes:
Hi All:
This may turn out to be very simply, but I can't seem to add the name of
the school to a chart. The loop I created is below that subsets a
dataframe and creates a chart for each school based on certain
variables. As it stands now, they title
Your statement seems innocent enough on the face of it, but there are two
facets that I think are worthy of note.
The first is that of time, and more specifically who's time. As a user of other
lists I can say that this is the best list in terms of getting the answer to my
problem, albeit
Dear Michael,
I had some time this evening, so I programmed the two-step procedure
described in the article that I mentioned. This isn't the ML estimate of the
correlation, but apparently it performs reasonably well and it is quite
fast. I checked the function on some examples and it appears to
On 29 Nov 2004 18:07:40 +0100, Peter Dalgaard
[EMAIL PROTECTED] wrote:
I had a situation where there was like
a (virtual) maximum outside the boundary, and the algorithm would
basically stay on the path to that peak, banging its little head
into the same point of the wall repeatedly, so to
Hi,
Suppose I have a 4 by 4 matrix as following:
10 6 3 1
8 4 3 2
6 2 4 3
9 3 4 2
The x axis is the column index of the matrix, the y axis is the row index of
the matrix and the z axis is the value in the corresponding position of the
matrix.
I tried to use contour and
On Monday 29 November 2004 21:07, Peter Yang wrote:
Hi,
Suppose I have a 4 by 4 matrix as following:
10 6 3 1
8 4 3 2
6 2 4 3
9 3 4 2
The x axis is the column index of the matrix, the y axis is the row
index of the matrix and the z axis is the value in the corresponding
position
I've just upgraded to 2.0.1 and was taken by surprise by the changes
in graphical parameter handling for lattice plots. I'd previously
been using 1.9.1. The old settings seem to have been replaced by a
daunting number of new options. I've poked around a bit and have not
seen any discussion of
Hi, Dear all R users:
Does someone know whether R can calculate the Receiver Operating
Characteristic (ROC) Curves? I didn't find it from the packages.
Thanks a lot.
--- Xin
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On Monday 29 November 2004 22:39, [EMAIL PROTECTED] wrote:
I've just upgraded to 2.0.1 and was taken by surprise by the changes
in graphical parameter handling for lattice plots. I'd previously
been using 1.9.1. The old settings seem to have been replaced by a
daunting number of new options.
package ROC from bioconductor, eg:
http://www.bioconductor.org/repository/release1.5/package/Win32/
Cheers!
Christoph
Xin Qi wrote:
Hi, Dear all R users:
Does someone know whether R can calculate the Receiver Operating
Characteristic (ROC) Curves? I didn't find it from the packages.
Thanks a
There is a package, Lexis, not officil though, which contains
a function ROC (and some other stuff for epidemiology).
You can find it in:
http://biostat.ku.dk/~bxc/SPE/library/
Bendix Carstensen
--
Bendix Carstensen
Senior Statistician
Steno Diabetes Center
Niels Steensens
76 matches
Mail list logo