Just to pick up
Third, the resulting line is not optimal for either predicting y for a
new x or x from a new y. It's hard to see why it is ever of much
interest.
It is not a regression (and hence the subject line was misleading), but it
does come up in errors-in-variables problems. Suppose
__
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, self-contained, reproducible code.
Hi everyone!
I have 1000 text files that contain data in a standard format separated by
semi colons. However, each file also has single header (first) and footer
(final) lines. The text files are all named with a sample number, which do
not follow a pattern. I need to delete the header and
Hallo,
my noisy time series represent a fading signal comprising of long
enough parts with a simple trend inside of each such a part.
Transition from one part into another is always a non-smooth
and very sharp/acute. In other words I have a piecewise
polynomial noisy curve asymptotically
Hi,all:
May be a pointless question
a - 1:10
b - matrix(1:8,nrow = 4)
c - letters[4:8]
……
ls()
[1] a b c
ls() can print the names of the objects in the memory ,
but I want to get
Dear all,
I tried out the example in the help document for mca (the multiple
correspondence analysis of the MASS package):
farms.mca - mca(farms, abbrev=TRUE)
farms.mca
plot(farms.mca)
But the graphic that I get seems unfeasible to me: I cannot recognize the
numbers (printed in black)
Sorry for inconvenience - but since I can't find any hints on this list at
stat.ethz.ch, I'll ask here: How do I unsubscribe from this list?
Thanks for answering - Richard
--
Richard Müller - Am Spring 9 - D-58802 Balve-Eisborn
www.oeko-sorpe.de
__
1. ? ls.str()
2. My favourite eval(parse(text= ... )) :)
sapply(ls(),FUN=function(x)eval(parse(text=x)))
--
View this message in context:
http://www.nabble.com/-R--How-to-print-the-objects-in-the-memory-tf3147472.html#a8725398
Sent from the R help mailing list archive at Nabble.com.
On Wed, 31 Jan 2007, Michael Reinecke wrote:
Dear all,
I tried out the example in the help document for mca (the multiple
correspondence analysis of the MASS package):
farms.mca - mca(farms, abbrev=TRUE)
farms.mca
plot(farms.mca)
But the graphic that I get seems unfeasible to me: I cannot
Dear List,
we're trying to install R on Solaris10 on a x86 (amd64).
During the installation we pass successfully the ./configure but we get
an error through the built-in function _isnan which we see existing in
/lib
When passing the command make we get :
gcc -std=gnu99 -I../../src/extra/zlib
On Wed, 31 Jan 2007, Vladimir Eremeev wrote:
[replying to something else, without any indication or context. That
'something' was about how to print all objects.]
1. ? ls.str()
2. My favourite eval(parse(text= ... )) :)
sapply(ls(),FUN=function(x)eval(parse(text=x)))
Note
Hmm,
I thought that learning to read was pretty wide-spread in your
country...
Richard == Richard Müller [EMAIL PROTECTED]
on Wed, 31 Jan 2007 10:01:34 +0100 writes:
Richard Sorry for inconvenience - but since I can't find
Richard any hints on this list at stat.ethz.ch, I'll ask
Dear List,
must be simple, but i got stuck on this...
I´ve been trying around for some time, but I could not find a
straigthforward way to do this:
How can one ask if a certain object is element of a certain group of
objects, e.g stored in a vector?
x - c(a, b, c)
# is there a simple
Peter Robinson wrote:
Dear List,
I have used the image() function to show a heat plot of a matrix of data
whose intensity is color-coded. I have two questions that I have not
been able to solve by using the help system or google.
1) How can one add a scale/legend that shows what
Christoph Heibl wrote:
Dear List,
must be simple, but i got stuck on this...
I´ve been trying around for some time, but I could not find a
straigthforward way to do this:
How can one ask if a certain object is element of a certain group of
objects, e.g stored in a vector?
x - c(a,
You need %in%
a %in% x
d %in% x
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics, methodology
This hasn't generated any feedback after a few days on R-devel, so I'm
forwarding it to R-help in case anyone here has any ideas...
Thanks,
Jon
-- Forwarded message --
From: Jon Clayden [EMAIL PROTECTED]
Date: 26-Jan-2007 11:25
Subject: readBin is much slower for raw input than
image.plot from fields plots an image and a legend (or scale)
Dr. med. Peter Robinson wrote:
I have used the image() function to show a heat plot of a matrix of data
whose intensity is color-coded.
1) How can one add a scale/legend that shows what numerical values a
given color
Chen, Xiao wrote:
Greetings -
I have a quick question that I hope someone will have a quick answer. I
have tried to use the R function system with the MS-DOS command type
to display the full content of a text file. But it always returns with a
message saying the text file is not found.
Hi,
I wrote some simple rpanel package script for visual spectral data comparison.
At this example i have a three samples and i want to zoom through x and y
axis to compare differences between samples. With my script below I can zoom
to some data region and add some other spetra to the plot,
Here is a workable solution:
df1 - data.frame(ar1)
df2 - data.frame(ar2)
cmn - intersect(names(df1),names(df2))
rbind(df1[,cmn],df2[,cmn])
Best
Bendix
__
Bendix Carstensen
Senior Statistician
Steno Diabetes Center
Niels Steensens Vej 2-4
DK-2820
This has come up several times, all with long obselete versions of gcc.
R does not use __builtin_isnan, but gcc does in its own (rather than
Solaris') header files. So it is an issue with the compiler, and note
that 3.4.2 is an old (Sept 2004) version of the compiler.
If you are happy to use
Hi,
I have two matrices A (m x 2) and B (n x 2), where m and n are large integers
(on the order of 10^4). I am looking for an efficient way to create another
matrix, W (m x n), which can be defined as follows:
for (i in 1:m){
for (j in 1:n) {
W[i,j] - g(A[i,], B[j,])
} }
I have tried the
Try 'file.show'
On 1/30/07, Chen, Xiao [EMAIL PROTECTED] wrote:
Greetings -
I have a quick question that I hope someone will have a quick answer. I
have tried to use the R function system with the MS-DOS command type
to display the full content of a text file. But it always returns with a
--- John Kane [EMAIL PROTECTED] wrote:
Date: Wed, 31 Jan 2007 08:32:03 -0500 (EST)
From: John Kane [EMAIL PROTECTED]
Subject: Re: [R] Simple Date problems with cbind
To: [EMAIL PROTECTED]
My thanks to Marc, Tony Plate, Phil Spector Mark
Leeds.
I totally missed the significance of
I think something like this should work:
outFile - file(output.txt, w)
for (i in ListOfFileNames){
x - readLines(i)
writeLines(x[2:(length(x)-1)], outFile)
}
close(outFile)
On 1/30/07, pif pif [EMAIL PROTECTED] wrote:
Hi everyone!
I have 1000 text files that contain data in a
Dear useRs,
an update of the ROCR package is available on CRAN.
ROCR helps in evaluating the performance of scoring classifiers using
ROC graphs, precision/recall plots, lift charts and many other
performance metrics.
For further information check http://rocr.bioinf.mpi-sb.mpg.de and
On Wed, 31 Jan 2007, Vladimir Eremeev wrote:
Chen, Xiao wrote:
Greetings -
I have a quick question that I hope someone will have a quick answer. I
have tried to use the R function system with the MS-DOS command type
to display the full content of a text file. But it always returns with a
Amy, et al,
I agree with you and the group that comparing test set classification
errors between the two methods is the way to go.
On interpretation, I find the partial dependence plots from
randomForest are useful - especially when talking to clients about
what the forest means. See slides 32
Hi usstata,
I think this will get you what you want:
mget(ls(),globalenv())
On 1/31/07, usstata [EMAIL PROTECTED] wrote:
Hi,all:
May be a pointless question
a - 1:10
b - matrix(1:8,nrow = 4)
c - letters[4:8]
……
ls()
[1] a b
I am using the rcor.test function available in the ltm package and
wanted to know if there was a way to use this function to compute the
correlation of two datasets that are not of equal length (similar to the
base cor function). I was thinking an apply function like by might work
but so far I
This problem also comes up in financial hedging problems,
but usually the 'errors' need not be of comparable size, so Errors in
Variables or Total Least Squares might be used.
David L. Reiner
Rho Trading Securities, LLC
Chicago IL 60605
312-362-4963
-Original Message-
From: [EMAIL
Hi Folks,
I just finished (the first draft of) a Wiki document to explains
options for producing graphics to incorporate in MS Word.
Thanks to all those who provided input, and who provided the wiki space.
Here is the direct link.
http://wiki.r-project.org/rwiki/doku.php?id=tips:using-
Hi,
Today I came upon unexpected R behaviour. I did some modelling and the
result was R object, about 28MB size (nested list, with matrixes as
list elements). When I was saving the session with save.image, the
resulting .RData file was 300MB. There were no other large objects:
To get the best results you need to transfer it using vector
graphics rather than bitmapped graphics:
http://www.stc-saz.org/resources/0203_graphics.pdf
There are a number of variations described here (see
entire thread). Its for UNIX and Windows but I think
it would likely work similarly on
Sir I am not finding the function to plot least square regression line on
type=o plot of two variables.
guid me in this regard.
--
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.
Email:
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[[alternative HTML
Sir I am new user of R. I am facing problems in solving matrices.
for example
a-array(c(3,4,1),dim=c(3,1))
b-a%o%a
b
, , 1, 1
[,1]
[1,]9
[2,] 12
[3,]3
, , 2, 1
[,1]
[1,] 12
[2,] 16
[3,]4
, , 3, 1
[,1]
[1,]3
[2,]4
[3,]1
solve(b)
Error in
Dear R-Help,
Thanks much.
I have received very good advice from a couple of experts. R-help is
just wonderful!
I have combined all the solutions that I got and they are shown below:
zz-file(d:/work/test/test.txt, w)
cat(this is a test\n, file=zz)
close(zz)
setwd(d:/work/test)
shell(type
On Wed, 2007-01-31 at 09:25 -0800, amna khan wrote:
Sir I am not finding the function to plot least square regression line on
type=o plot of two variables.
guid me in this regard.
Did you want something like this:
x - 1:50
y - rnorm(50)
plot(x, y, type = o)
abline(lm(y ~ x))
See ?abline
Hi,
I am pleased to announce the availability of the document, R for SAS
and SPSS Users, at
http://oit.utk.edu/scc/RforSASSPSSusers.doc
I've looked at the document and printed it. I think it will be very
useful to me, even if I will use it the reverse way : learn how to
use SAS from R...
As
Two comments:
1) ls() does not list all the objects: it has all.names argument.
2) save.image() does not just save the objects in the workspace, it also
saves any environments they may have. Having a function with a
large environment is the usual cause of a large saved image.
(And finally, a
Hi, all
I met a problem to query GSE3524, which cannot be open on my computer. I
hope some of you would be kind to give me some advice. Thanks!
The code is as follow:
##
library(GEOquery)
gsename=GSE3524
gse=getGEO(gsename)
##
The error information follows as
Previous subject:
bootstrap bca confidence intervals for large number of statistics in one model;
library(boot)
Jacob Wegelin asked for an easier way to do many bootstrap confidence
intervals for regression output.
The syntax would be easier with S+Resample, example below.
You create an ordinary
a-array(c(3,4,1),dim=c(3,1))
b-a%*%a
solve(b)
amna khan [EMAIL PROTECTED] a écrit : Sir I am new user of R. I am facing
problems in solving matrices.
for example
a-array(c(3,4,1),dim=c(3,1))
b-a%o%a
b
, , 1, 1
[,1]
[1,]9
[2,] 12
[3,]3
, , 2, 1
[,1]
[1,] 12
[2,] 16
Hello,
Is there a way to convert a character to a number with out getting a warning?
I have a vector that has both numbers and letters in it and I need to convert
it to only numbers. At the moment I'm using as.numeric but it is generating a
warning when it converts a letter. Is there another
Justin,
Did you actually run that code? :-)
a-array(c(3,4,1),dim=c(3,1))
b-a%*%a
Error in a %*% a : non-conformable arguments
I suspect that Amina wants:
a - array(c(3, 4, 1), dim=c(3, 1))
b - t(a) %*% a
solve(b)
[,1]
[1,] 0.03846154
or perhaps:
solve(crossprod(a))
suppressWarnings(a - as.numeric(c(1, 2, pi, a, 9, z)))
b
On Jan 31, 2007, at 2:35 PM, Konrad wrote:
Hello,
Is there a way to convert a character to a number with out getting
a warning? I have a vector that has both numbers and letters in it
and I need to convert it to only numbers. At
Konrad wrote:
Hello,
Is there a way to convert a character to a number with out getting a warning?
I have a vector that has both numbers and letters in it and I need to
convert it to only numbers. At the moment I'm using as.numeric but it is
generating a warning when it converts a
Hello all,
I have a population of 112 servers that are experiencing different levels of
packet loss. I don't want to poll all 112 of them (the analytical tools must
be manually run on each individually) so it seems best to sample among them;
then I plan on using R to run comparisons of the
On Wed, 2007-01-31 at 14:35 -0500, Konrad wrote:
Hello,
Is there a way to convert a character to a number with out getting a
warning? I have a vector that has both numbers and letters in it and
I need to convert it to only numbers. At the moment I'm using
as.numeric but it is generating a
On 31-Jan-07 Konrad wrote:
Hello,
Is there a way to convert a character to a number with out getting a
warning? I have a vector that has both numbers and letters in it and I
need to convert it to only numbers. At the moment I'm using as.numeric
but it is generating a warning when it
On 31-Jan-07 Benilton Carvalho wrote:
suppressWarnings(a - as.numeric(c(1, 2, pi, a, 9, z)))
Of course! Much better!
Ted.
E-Mail: (Ted Harding) [EMAIL PROTECTED]
Fax-to-email: +44 (0)870 094 0861
Date: 31-Jan-07
On Wed, 2007-01-31 at 20:02 +, [EMAIL PROTECTED] wrote:
On 31-Jan-07 Benilton Carvalho wrote:
suppressWarnings(a - as.numeric(c(1, 2, pi, a, 9, z)))
Of course! Much better!
Ted.
In the context of my prior reply:
# Bear in mind that the above vector is of class character, not mixed...
Hi all,
This is probably a blindingly obvious question: Why does it matter in
the uniroot function whether the f() values at the end points that you
supply are of the same sign?
For example:
f - function(x,y) {y-x^2+1}
#this gives a warning
uniroot(f,interval=c(-5,5),y=0)
Error in uniroot(f,
[EMAIL PROTECTED] wrote:
This is probably a blindingly obvious question:
Yes, it is.
Why does it matter in the uniroot function whether the f() values at
the end points that you supply are of the same sign?
Plot some graphs.
Think about the *name* of the function
The last two times I have originated message threads on R or
Bioconductor I have received the message included below from someone
named Patrick Connolly. Both times I was the originator of the message
thread and used what I thought was a unique subject line that explained
as best I could what my
On 1/31/2007 5:38 PM, Kimpel, Mark William wrote:
The last two times I have originated message threads on R or
Bioconductor I have received the message included below from someone
named Patrick Connolly. Both times I was the originator of the message
thread and used what I thought was a unique
Kimpel, Mark William wrote:
The last two times I have originated message threads on R or
Bioconductor I have received the message included below from someone
named Patrick Connolly. Both times I was the originator of the message
thread and used what I thought was a unique subject line that
Peter,
Thanks you for your explanation, I had taken Mr. Connolly's message to
me to imply that I was not changing the subject line. I use MS Outlook
2007 and, unless I am just not seeing it, Outlook does not normally
display the in reply to header, I was under the mistaken impression
that that
Julien Barnier wrote: ... I think it will be very useful to me, even if
I will use it the reverse way : learn how to use SAS from R...
I hadn't thought of using the document in reverse to learn SAS or SPSS
if you already know R. I'll have to reread it from that perspective
see if there are any
Hi Mark,
I'm sending this off-list because I don't want to unnecessarily fill
people's inbox if I'm wrong.
I think you are misunderstanding Patrick's problem. The point is kind of
subtle and requires a bit of knowledge of how e-mail works.
What he is complaining about is that people take a
Hi,
is there a 3D equivalent for hist2d {gplots}
that would give me the 3-dimensional matrix
which holds the counts of the number of observed
(x,y,z) triplets that fall into each bin ?
Many thanks in advance,
georgia
[[alternative HTML version deleted]]
See below for Bert Gunter's off list reply to me (which I do
appreciate). I'm putting it back on the list because it seems there is
still confusion regarding the difference between threading and sorting
by subject. I thought the example I will give below will serve as
instructional for other
Your final paragraph has the take-home message for everyone (not just MS
Outlook users): just create, from scratch, a new message when
initiating a new subject.
Viewing threads can be completely different to sorting based on the
subject line. Your initial post with the subject regexpr and
Hi list,
I am trying to print a rgl scene. I can do this with
rgl.snapshot(test.pnt,fmt=png), but
rgl.postscript(test.pdf,fmt=pdf) does crash R and returns to the
linux shell after extented time period.
I am running R 2.4.0 on i686 linux mandrake 10.2. Do I need any other R
external program to
You can see how it looks to most readers by viewing it on gmane:
http://thread.gmane.org/gmane.comp.lang.r.general/78065
Note that even though the subject has been changed its still
listed as a child of another message rather than the start
of a new thread.
On 1/31/07, Tony Plate [EMAIL
Tony,
I went to the MS link that you suggested (see below) and it indeed says
that The Arrange by Conversation arrangement shows your e-mail items
grouped by message subject or 'thread.' Instead of arranging by
subject, I arranged my view by conversation and got exactly the same
result that I had
I need to sum the columns of a sparse matrix according to a factor -
ie given a sparse matrix X and a factor fac of length ncol(X), sum
the elements by column factors and return the sparse matrix Y of size
nrow(X) by nlevels(f). The appended code does the job, but is
unacceptably
Hello everyone,
Here is the setup.
z is a 119 x 15 matrix, m_index is a 119 x 5 matrix
What I am trying to do is return the results from fitCopula by sequentially
binding all 15 columns of z to the first column of m_index,
(cbind(z[,1],m_index[,1]),(cbind(z[,2],m_index[,1]), etc.
On 1/31/2007 8:15 PM, [EMAIL PROTECTED] wrote:
Hi list,
I am trying to print a rgl scene. I can do this with
rgl.snapshot(test.pnt,fmt=png), but
rgl.postscript(test.pdf,fmt=pdf) does crash R and returns to the
linux shell after extented time period.
I am running R 2.4.0 on i686 linux
Hi all,
This information must be out there, but I can't seem to find it. What I
want to do is to store functions I've created (as .R files or in
whatever form) and then load them when I need them (or on startup) so
that I can access without cluttering my program with the function code.
This
This is really off-topic for both BioC and R-help, so I'll
keep it short.
From: Kimpel, Mark William
See below for Bert Gunter's off list reply to me (which I do
appreciate). I'm putting it back on the list because it seems
there is still confusion regarding the difference between
Dear all R users,
Suppose I have two data sets:
data1: 1,3,5,7,9
data2: 2,4,6,8
Now I want to combine these two data sets like that:
data = 1,2,3,4,5,6,7,8,9
Can anyone tell me how to do that?
Thanks and regards,
[[alternative HTML version deleted]]
The recommended approach is to make a package for your functions that will
include documentation, error checks etc.
Another way to accomplish what you want is to start a new R session and
'source' your .R files and then to save the workspace in a .RData file, e.g.
myFunctions.RData.
Finally
Hello sir:
I wanna get such kind of plot: a line whose start point is(1,10),end point
is(5,10)
In other words:
How can I draw a line if I only know the coordinate of the start point and end
point? Thanks! My best
__
R-help@stat.math.ethz.ch
Dear R-help,
I'm using prop.test() to compute a confidence interval for a proportion
under R version 2.4.1, as in:
prop.test(x = 340, n = 400)$conf
[1] 0.8103309 0.8827749
I have two questions:
1) from the source code my understanding is that the confidence
interval is computed according to
Hi
see ?segments
segments(1,10,5,10)
HTH
Petr
On 1 Feb 2007 at 14:21, XinMeng wrote:
From: XinMeng [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Date sent: Thu, 01 Feb 2007 14:21:34 +0800
Subject:[R] line plot
Send reply to:
Hi,
On 1/31/07, Professor Brian Ripley [EMAIL PROTECTED] wrote:
Two comments:
1) ls() does not list all the objects: it has all.names argument.
Yes, I tried it with all.names, but the effect was the same, I forgot
to mention it in a letter.
2) save.image() does not just save the objects in
78 matches
Mail list logo