Try source(beams.txt) which reads and executes the commands in
beams.txt.
It might make more sense in future to store and work with dataframes,
especially with big datasets.
Regards, Adai
On Tue, 2005-03-22 at 17:00 +0800, chin wei wrote:
hi, all.
I am a beginner. Could you tell me how to
See Andy Liaw's and my suggestion to this post
http://files.protsuggest.org/biocond/html/7818.html
On Fri, 2005-03-18 at 11:40 +, Pingping Zheng wrote:
Suppose I have a nx2 matrix of data, X, the following code generate
density estimation for each column and plot them
denlist - apply(X,
If you know the exact formulae for the distribution, replace it with 'f'
function below. You may want to use the log=x in the plot.
f - function(x) 1 - exp( -x/20 );
plot( f, xlim=c(0,100), ylim=c(0.5, 1) )
Otherwise generate sufficient realisations from it and fit a line as
below
x -
Yes, I was thinking of the trivial problem of 2 matrices. Nice to know I
am not the only one who made the same error.
Thanks to Thomas Lumley, John Fox, Dimitris Rizopoulos for pointing this
out and many others for providing the correct solution.
Regards, Adai
On Wed, 2005-03-16 at 11:07
As promised here is my reply. I am assuming that your problem is
changing the order of the xyplots by idno.
The plotting order in lattice is determined by the levels of md$idno. By
default it was reading the levels from top to bottom of the dataset as
unique(md$idno) would do.
library(nlme);
We need more info to help you. What do you mean by crash, does it
generate an error or something ?
I presume you have just installed R. How did you install it ? If you did
from source, did you do a make check ? I had this problem once and it
was traced to a broken gcc version
Well, one way you can try is to define the different styles you want in
your $HOME/.Rprofile file (see ?Startup). For example
---
library(graphics)
op - par(no.readonly = TRUE) # store original par
par0 -
mylist - list( matrix(1:6, nc=3), matrix(7:12, nc=3) )
do.call(+, mylist)
[,1] [,2] [,3]
[1,]8 12 16
[2,] 10 14 18
Regards, Adai
On Wed, 2005-03-16 at 18:21 +0200, Vicky Landsman wrote:
Dear all,
I think that my question is very simple but I failed to solve it.
I have a
You will need to _apply_ the t-test row by row.
apply( genes, 1, function(x) t.test( x[1:2], x[3:4] )$p.value )
apply() is a C optimised version of for. Running the above code on a
dataset with 56000 rows and 4 columns took about 63 seconds on my 1.6
GHz Pentium machine with 512 Mb RAM. See
If you use par(new=TRUE), you are overlaying one graph on top of the
other. In which case, make sure your xlim and ylim are correctly set.
Another way is to split the plotting window. For example
par(mfrow=c(2,3))
for(i in 1:6) hist( rnorm(100), main=paste(Histogram, i))
See help(par) for
You will need to capture the value of ss at the end of each 'i' as such
z4 -function(w){
output - numeric(w)
for (i in 1:w){
set.seed(i+6) # this is redundant line
ss-0
for (j in 1:5){
set.seed(j+1+(i-1)*6)
r-rnorm(1)
ss-ss+r
}
output[i] -
Not answering your question here directly but here is an alternative.
Interactive search and replace TRUE for FALSE with M-% after
highlighting the region that you are interested in. And then repeat for
changing FALSE to TRUE.
You can bind some function keys in your configuration file to
interface of Stata for windows does not exist.
Shige
On Fri, 04 Mar 2005 00:19:23 +, Adaikalavan Ramasamy
[EMAIL PROTECTED] wrote:
There is an ESS mailing list that might be more appropriate.
I use emacs rather than Xemacs, but I think you need to point the path
to ess in init.el
See help(set.seed).
set.seed(1)
rnorm(5)
[1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078
rnorm(5)
[1] -0.8204684 0.4874291 0.7383247 0.5757814 -0.3053884
set.seed(1)
rnorm(5)
[1] -0.6264538 0.1836433 -0.8356286 1.5952808 0.3295078
On Fri, 2005-03-04 at 13:07 +0200, Clark
Use the collapse argument in paste.
paste( c(a, b, c), collapse= )
[1] abc
On Fri, 2005-03-04 at 03:32 -0800, Matthieu Cornec wrote:
Hello,
I would like to convert c(a,b,c) into abc.
Anyone could help?
Thanks,
Matthieu Cornec
__
Are you using Windows operating system ? If so, then you will need to
download the executable not the source codes from
http://www.cran.r-project.org/bin/windows/base/rw2001.exe
Regards, Adai
On Fri, 2005-03-04 at 12:54 +0100, Ramseyer Amandine wrote:
Hello,
I'm a student in biology,
How will you deal with multiple word searches such as
help.search(eps dev)
One way to implement would be ??eps dev but this looks awkward to me.
Regards, Adai
On Thu, 2005-03-03 at 10:21 +, Yan Wong wrote:
On 3 Mar 2005, at 10:08, Duncan Murdoch wrote:
That's not a bad
Try scatterplot3d() in the scatterplot3d package.
Alternatively try searching http://maths.newcastle.edu.au/~rking/R/
Regards, Adai
On Thu, 2005-03-03 at 17:45 +, Stephan Freyberger wrote:
Hi,
Trials have generated a vast number of points, which I would like
to have plotted in 3D.
Why not simply read it as an csv file, then transpose it. If you also
store it as a data frame, you can use attach() or detach() the object to
the search path whenever you want to access the variables directly.
df - read.csv( tmp.txt, header=FALSE, row.names=1 )
df - data.frame( t( df ) )
df
One solution is to cut() 'x' according to the breaks defined by 'y'.
Using cut with labels=FALSE is really fast. See a simulation below.
However the accuracy depends on the number of ties you have in your
empirical distribution. I have tried to simulate with the round()
function below.
#
There is an ESS mailing list that might be more appropriate.
I use emacs rather than Xemacs, but I think you need to point the path
to ess in init.el file which is located on the home directory. E.g. :
(setq ess-icon-directory C:/Programme/xemacs-packages/etc/ess)
(require 'ess-site)
Just out of curiosity, what is the difference between the terms for
package and library ? Why are we loading a package with the library()
command ?
If this is a case of RTFM, I would be happy to do so if pointed in the
right direction. I have searched the FAQ and mail archives and only came
up
You have to decided whether you want one single plot or 23 plots (one
for each chromosome). Since I suspect that, say, copy_num at the end of
Chromosome 1 and the copy_num at the start of Chromosome 2 are not
related (like time series data) and would suggest 23 plots instead.
# simulate data #
Schwartz wrote:
On Wed, 2005-03-02 at 13:42 +, Adaikalavan Ramasamy wrote:
Just out of curiosity, what is the difference between the terms for
package and library ? Why are we loading a package with the library()
command ?
If this is a case of RTFM, I would be happy to do so if pointed
Here is one way
x - c( VGT1_CONTR_B020020301.H0V0.MIR,
VGT1_VGT2_CONTR_B020020611.H0V0.MIR,
VGT1_VGT2_CONTR_B020030401.H0V0.MIR,
VGT1_VGT2_CONTR_B020030711.H0V0.MIR,
VGT1_VGT2_CONTR_B020031211.H0V0.MIR )
( tmp1 - sapply( strsplit(x, split=_), tail, n=1 ) )
[1]
Might be slightly more interesting. If we want to generate values which
are completely missing at random, then we can just simply sample all
available index of a 2-d array.
# simulate data #
set.seed(1) # for reproducibility
m - matrix( rnorm(12), nr=4, nc=3 )
m
[,1] [,2]
I am assuming you are talking about a Windows machine. If so, reading
http://cran.r-project.org/bin/windows/base/rw-FAQ.html#Installation-and-
Usage
will point you the following windows installer
http://cran.r-project.org/bin/windows/base/rw2001.exe
However, you may need to ask your IT person
You might want to read (or re-read) the posting guide about giving a
simple example. See comments below.
On Mon, 2005-02-28 at 23:03 +0800, Feng Chen wrote:
Dear all,
I have something about function outer() that I can't understand. Just see the
following example. The two NaNs are due to
Reading the posting guide http://www.R-project.org/posting-guide.html
will tell you to specify an informative subject line in your postings.
You can change working directory by using setwd(), see help(setwd).
For example setwd(c:\My Documents) might work. I am not sure if you
need to truncate the
Searching for graph publication on
http://maths.newcastle.edu.au/~rking/R/ gave me the following hit :
http://tolstoy.newcastle.edu.au/R/help/04/03/0202.html which suggests
postscript().
Have you tried printing other documents in black and white on the same
printer or tried different printers for
Sorry, all backslashes have to be doubled in R as mentioned in FAQ 2.14
http://cran.r-project.org/bin/windows/base/rw-FAQ.html#R-can_0027t-find-
my-file
Regards, Adai
On Tue, 2005-03-01 at 04:33 +, Adaikalavan Ramasamy wrote:
Reading the posting guide http://www.R-project.org/posting
Forgive me for I do not fully comprehend the idea of classes and methods
but I was wondering if someone could help explain why the function args
() behaves the way it does.
Why does args(cut) show the simplified version instead of the more
complete one as in help(cut). This is true for few other
, 25 Feb 2005, Adaikalavan Ramasamy wrote:
Forgive me for I do not fully comprehend the idea of classes and methods
but I was wondering if someone could help explain why the function args
() behaves the way it does.
Why does args(cut) show the simplified version instead of the more
://www.med.kuleuven.ac.be/biostat/
http://www.student.kuleuven.ac.be/~m0390867/dimitris.htm
- Original Message -
From: Adaikalavan Ramasamy [EMAIL PROTECTED]
To: R-help r-help@stat.math.ethz.ch
Sent: Friday, February 25, 2005 1:14 PM
Subject: [R] display full form in args
I am confused. Are you saying that your two data frames are of different
dimensions ?
In any case what I think what you are looking for is which.
# generate the conditioning matrix
a - matrix( sample(0:1, 9, replace=TRUE), nc=3 )
a
[,1] [,2] [,3]
[1,]111
[2,]100
say 10,
of element in test1$'cm' , if not then 0
Luis
Adaikalavan Ramasamy [EMAIL PROTECTED] 23/02/2005 15:07:55
I am confused. Are you saying that your two data frames are of
different
dimensions ?
In any case what I think what you are looking for is which.
# generate
See the details and references section of help(wilcox.test). A small
section from the details might be relevant :
By default (if 'exact' is not specified), an exact p-value is
computed if the samples contain less than 50 finite values and
there are no ties. Otherwise, a normal
Two other possibilities come to my mind :
1) If memory serves, when you un-subscribe using the web interface, you
will be sent an email asking for confirmation. Till you click on the
link in that email, you will not be un-subscribed. Did you you receive
any such mail ? If not, you might want to
AFAIK, it does not matter as long as you got the correct paths set in
your init.el file.
Please see the excellent installation instructions from John Fox's
http://socserv.mcmaster.ca/jfox/Books/Companion/ESS/
Regards, Adai
On Wed, 2005-02-23 at 17:16 -0600, Laura Holt wrote:
Dear R
Which plotting function are you using ?
I think most of plotting can accept xaxt=n which is the command to
supress the x-axis. If this works, at least you do not have to redefine
the function. Examples
plot(1:10, xaxt=n)
hist( rnorm(100), xaxt=n )
boxplot( rnorm(10), rnorm(10), rnorm(10),
Here is an another way
count - is.na(x) + is.na(y)
which( count == 1, arr.ind=TRUE )
'count' gives you the number of missing values at for each row and
column. Then you can find out how many occurances of both missing, none
missing and one missing.
On Mon, 2005-02-21 at 15:48 +0100,
The sigToEnv error has been reported several times on the BioConductor
mailing list in this month. See :
https://stat.ethz.ch/pipermail/bioconductor/2005-February/007632.html
https://stat.ethz.ch/pipermail/bioconductor/2005-February/007674.html
You can read in the data using read.delim() or read.table(). For
illustration let us generate some artificial data and suppose that you
are interested in equal sized breaks of 5 (you can define your own break
points instead).
x - rchisq(50, df=10, ncp=5)
brk - seq(0,
1) Why do you think there is a bug in rnorm. What is your expected
numbers for cases 4 or cases -4 ? If you can show your calculations
then maybe we can see if there is an error in how you calculated it or
with R.
2) Please give a reproducible example. It might be low or high in one
instance
Here is a generalisation of the function that others have suggested to
take take more than 2 vectors.
my.barplot - function(...){
my.list - list(...)
lev - sort( unique( unlist(my.list) ) )
tmp - t(sapply( my.list, function(v) table(factor(v, levels=lev))) )
barplot(tmp, beside=T)
}
w
Please try to use something other than l2 for showing an example
because it is looks awfully similar to the number 12.
Here is another of doing this but this will only work if you have equal
lengths in each of your elements.
# simulate data
mylist - lapply( 1:5, rnorm, n=2 )
names(mylist) -
Assume that you have stored the lm object as 'fit' and the summary as
fit.summ as such
x - rnorm(100)
y - rnorm(100)
fit - lm( y ~ x )
fit.summ - summary( fit )
fit.summ$coefficients and fit.summ$adj.r.squared gives you the
coefficients and adjusted R-square.
names( fit.summ ) or str(
Either set the 'main' or 'xlab' in the hist(). See help(par) for more
information on graphical arguments or help(hist).
mat - matrix( rnorm(1000), nc=5 )
colnames(mat) - LETTERS[1:ncol(mat)]
for( i in 1:ncol(mat) ){
hist( mat[ ,i],
main=paste( Histogram of data from column ,
yes but add 'replace=TRUE' into that statement to sample with
replacement.
On Wed, 2005-02-16 at 14:00 -0500, Liaw, Andy wrote:
Wouldn't sample(length(V), prob=V) do?
Andy
From: Ben Hyde
I have a vector V. sum(V) = 100, i.e. it's percentages.
length(V) is
large. I wish to
On Sun, 2005-02-13 at 20:53 -0800, Cuichang Zhao wrote:
Hello,
right now, i have a program to collect data into a table. right now, my table
is
table1 - data.frame(trial = NA, x = NA, y = NA)
One often uses the term 'table' as in tabulate or cross-tabulate
discrete values.
for each time
See help(cut) for general or help(cut.POSIXt) for time related cut.
On Mon, 2005-02-14 at 13:08 -0500, Omar Lakkis wrote:
I am collecting one price for market daily.
date1 price1
date2 price2
date3 price3
I have a roll schedule where if
date1 is between d1 and d2 then market is
This can be simplified slightly by
unique( t( sapply( permn( c(1,1,0) ), c ) ) )
Here is another possibility :
a - expand.grid( 0:1, 0:1, 0:1 )
a[ which( rowSums(a) == 2 ), ]
which gives
Var1 Var2 Var3
4110
6101
7011
Regards, Adai
On Sun,
not knoe why it is doing that
Is this something that was meant and not documented or something?
Jean
On Fri, 11 Feb 2005, Adaikalavan Ramasamy wrote:
See below.
On Thu, 2005-02-10 at 14:53 -0800, T. Murlidharan Nair wrote:
Just to explain my previous mail, here
How about data.lm$model ?
On Wed, 2005-02-09 at 16:32 +0100, [EMAIL PROTECTED] wrote:
Petr,
It works perfectly! But I still have a question;
I have fit the following data;
x,y,z
1,10,11
2,11,15
3,12,21
4,13,29
5,14,39
6,15,51
7,16,65
8,17,81
9,18,99
10,19,119
dat.lm -
Please use a sensible subject line. In short, you need help().
Most R functions have a section called Examples in the help file. E.g.
help(plot) or even help(help). Granted this will give you only
snippets of a whole analysis but it is helpful nonetheless.
Since you have already identified the
Looking at the output from this and the previous response using barplot
worries me for two reasons :
a) The bars in front may obscure the bars at the back. The order of
plotting becomes important here.
b) IMHO, it is distracting and does not convey information well unless
one is willing to risk
When you start you R session, see if it says
[Previously saved workspace restored]
immediately before the first command prompt.
I think the only way a workspace can be restored automatically is if you
used save.image() at some point, in which case it would save all objects
as a hidden file
Please give a simple example of the input data and output that you
desire. It is difficult to understand from you partial codes what you
mean. For example what is Y ?
Are you trying to find add values from pairs of rows ? If so, please see
my posting pairwise difference operator where I wanted to
Please learn to wrap your emails at about 72 characters.
See below for other comments.
On Thu, 2004-12-16 at 08:29, vasilis pappas wrote:
Hello R users,
I have three questions and I would be grateful if someone could give me an
answer to each of these.
1) I have constracted a
://www.R-project.org/posting-guide.html
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/
Cancer Research UK Tel : 01865 226 677
Old Road Campus, Headington, Oxford Fax : 01865 226 962
! http://www.R-project.org/posting-guide.html
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/
Cancer Research UK Tel : 01865 226 677
Old Road Campus, Headington, Oxford Fax : 01865 226 962
://www.R-project.org/posting-guide.html
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/
Cancer Research UK Tel : 01865 226 677
Old Road Campus, Headington, Oxford Fax : 01865 226 962
__
[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
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine http
]]
__
[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
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine
://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/
Cancer Research UK Tel : 01865
You have not called legend() in your codes below, so we do not know what
your problem is. See other comments below.
On Mon, 2004-11-15 at 01:08, Sean David Richards wrote:
R : Version 1.9.1
Hi,
Am having trouble adding a legend to scatterplot. R code is shown below.
I have tried various
What do you mean by power curves ? Is it the power of a study as the
effect size varies or power output of a machine with some other
parameter ?
I usually generate a sequence of numbers (for the x-axis) that spans the
range of interest and calculate its output.
# Example 1
f - function(x)
Sorry typo. The last line should read
legend(1500, 9000, legend=paste(Data from, sfiles), pch=1:n, col=1:n )
^^^
On Mon, 2004-11-15 at 11:39, Adaikalavan Ramasamy wrote:
You have not called legend() in your codes below, so we do
Try with less ambitious numbers such as my.choose1(60,20). I think it
works fine.
I think the problem is that gamma(6001) and prod(1:6000) are so large
that it gives Inf as the answer. Hence the numerator and denominator
approaches Inf and division of two Inf gives NaN.
You could use the natural
get(myvar1) + get(myvar2)
[1] 295 411 230
On Thu, 2004-11-04 at 13:58, Oskar Villani wrote:
Hello list,
I'd like to use a variable (or a column of a data frame) by using its name as a
string. E.g.:
Data2003 - c(150,200,120)
Data2004 - c(145,211,110)
myvar1 - Data2003
myvar2 -
x - 1.0001
all.equal(x, 1)
[1] TRUE
x == 1
[1] FALSE
Reading help(all.equal) will tell you that the tolerance level by
default is around 1.490116e-08.
On Fri, 2004-11-05 at 13:48, Kay Pilz wrote:
Hello All.
I am running R 2.0.0 with a Win XP operating system. The same problem occured
Check http://cran.r-project.org/src/contrib/PACKAGES.html.
AFAIK, there is tree, rpart and knnTree.
The simplest wat to install, say the package tree, is to type in
install.packages(tree) in the R command line.
On Fri, 2004-11-05 at 15:07, Xin Qi wrote:
Dear all R users and helpers:
I
typed warning() (in most cases the group
I am testing is relatively large). But all the data
are from the same data set.
So what is the distinction between these two?
Many thanks,
Fang
--- Adaikalavan Ramasamy [EMAIL PROTECTED]
wrote:
Reading the posting guide will tell you (among
Sorry if this is a question more on regular expressions. I am dealing
with several files which have been badly named. For example the files
are given either the extensions txt, TXT or Txt. I wish to select all
those files ending with 'txt' ignoring case.
Here is how I would do it in bash (Redhat
:
On Thu, 4 Nov 2004, Adaikalavan Ramasamy wrote:
Sorry if this is a question more on regular expressions. I am dealing
with several files which have been badly named. For example the files
are given either the extensions txt, TXT or Txt. I wish to select all
those files ending with 'txt
Reading the posting guide will tell you (among other things) to :
a) use a sensible subject line
b) give a reproducible example when possible
See other comments below.
On Wed, 2004-11-03 at 19:04, fang lai wrote:
Hi there,
When I ran the wilcox.exact test, it always shows
number of items
/. You system
administrator should be able to help.
Regards, Adai
On Mon, 2004-11-01 at 12:28, [EMAIL PROTECTED] wrote:
Hi, Adaikalavan Ramasamy:
Thanks your message.
I installed R in our SunRay. It seems it is only commend line, could I run
RGui in SunRay? Thanks.
Have a good day
Francisco, a more reproducible example would have helped but see if the
following helps in your understanding
# Create dataset
df - data.frame( type=1:5, species=LETTERS[1:5] )
df
type species
11 A
22 B
33 C
44 D
55 E
df[ ,2] == E
I cannot make sense of your question. There is no reproducible example
to help either. Please read the posting guide.
However, I think you might find using arr.ind=TRUE option in which()
when dealing with matrices to be useful. See help(which). Example :
mat - matrix( 1:9, nc=3 ) + 100
mat[3,1]
See comments below.
On Thu, 2004-10-28 at 18:49, S Peri wrote:
Dear Group,
I am using DEAL package for modeling signal
transduction nets. This process is deal slow on a
SunFire server with 8 gigs ram.
we have a grid that can process much faster that one
individual server.
I presume R when executed from its bin directory works fine.
Are you actually using bash shell. Try 'echo $SHELL'.
If yes, then try 'which R' and see if it pointing to the right path.
On Mon, 2004-10-25 at 15:55, Alicia Amadoz wrote:
Hello,
I'm new to the installation R on Linux and I've
You will need to get R working first before installing BioConductor
packages.
Please provide a more useful output (like the error message) if you want
useful help.
Have you tried reading the R manual or searching the mail archives [2].
[1] http://cran.r-project.org/doc/manuals/R-admin.pdf
[2]
Perhaps save.image(file=lala.rda, compress=T) followed by
load(lala.rda) to re-load it.
On Fri, 2004-10-22 at 18:05, Andreas Buness wrote:
Hello,
I like to save the complete state of R, i.e. including
all environments, objects/workspaces, loaded packages etc..
This wish has arisen since
1. This question is more appropriate for the BioConductor mailing list.
2. Any Affymetrix pre-processing can be split into background
correction, normalisation and summary stages. The combination of rma
background, quantile normalisation and median polish summary gives rise
to RMA. Therefore
Have you checked Section 5.14 (page 64) of R-extension titled Using
these functions in your own C code.
Section 5.7 tells you about some of the available subroutines.
On Fri, 2004-10-15 at 01:31, doktora v wrote:
Hey everyone,
I have been looking for a while for ways to integrate R's
This is not answering your question directly.
I usually use the BATCH command for running R non-interactively. You can
also use commandArgs() to get any arguments from the command line. For
more information, see help(BATCH) or help(commandArgs).
On Thu, 2004-10-07 at 16:14, Kunal Shetty wrote:
r[ order(o) ] will give you the answer. More generally,
x - rnorm(100)
identical(x, x[ order(x) ][ order(order(x)) ])
[1] TRUE
On Mon, 2004-10-04 at 15:21, Wolfram Fischer wrote:
I have:
d - sample(10:100, 9)
o - order(d)
r - d[o]
How I can get d (in the original order), knowing
Alternatively you can use the combinations() function in gregmisc.
library(gregmisc)
n - 3
sapply(1:n, function(x) apply(combinations(n, x, LETTERS[1:n]), 1,
paste, collapse=) )
Here is the same code written with a for loop
n - 3
out - NULL
for(i in 1:n){
tmp - combinations(n, i,
1) Read the posting guide http://www.R-project.org/posting-guide.html
which tells you to
a) use a meaningful subject line
b) give a reproducible example. I cannot find plotData and there is a
syntax error (comma after the xyplot line).
2) The following example works fine for me
for(i in
m - matrix(1:9, nc=3)
a - list(m, m+10, m+100, m+1000)
sapply(a, function(mat) mat[1,1])
[1]1 11 101 1001
On Thu, 2004-09-16 at 12:33, Perez Martin, Agustin wrote:
DeaR useRs:
I have a list with 500 elements, in each other there are data.frames and I
want to take the first row and
456274 12.2 984024 26.3
Vcells 10123014 77.3 10538396 80.5
On Mon, 2004-09-13 at 18:47, Prof Brian Ripley wrote:
On Mon, 13 Sep 2004, Adaikalavan Ramasamy wrote:
I am comparing two different algorithms in terms of speed and memory
usage. I can calculate the processing time
See ?density or try help.search(kernel)
On Mon, 2004-09-13 at 16:05, Alexandre Bournery wrote:
Hello,
Does anyone could tell me if R sofware can be used for smoothing, using the
kernel distribution ?
I will appreciate
Alex
[[alternative HTML version deleted]]
I am comparing two different algorithms in terms of speed and memory
usage. I can calculate the processing time with proc.time() as follows
but am not sure how to calculate the memory usage.
ptm - proc.time()
x - rnorm(100)
proc.time() - ptm
I would like to be within R itself since
Just finished updating and installing new packages from CRAN and
BioConductor (including annotation data) and am happy to say that my R
has just exceeded the 1 GB mark.
On Fri, 2004-09-10 at 10:11, Jonathan Baron wrote:
On 09/10/04 03:54, Adaikalavan Ramasamy wrote:
There is another issue
Try setting something like cex=0.5 in the plot. It works for at least
hclust. Here is a quick example using hclust to show how you can get
away with set the colours.
data(USArrests)
hc - hclust(dist(USArrests), ave)
plot(hc, hang=-1, cex=0.5)
labels - rownames(USArrests)[ hc$order]; n -
Out of curiosity, why is srt not respected in mtext ? See examples below
plot(1:10, xaxt=n);
par(srt=45)
mtext( paste(Point, LETTERS[1:10]), side=1, at=1:10 )
plot(1:10, xaxt=n);
mtext( paste(Point, LETTERS[1:10]), side=1, at=1:10, srt=45 )
It would be nice to have 'srt' as an argument in
))
plot(1:10, main=paste(x, eval(expression(Delta)), values))
plot(1:10, main=paste(x, expression(Delta, values )))
plot(1:10, main=paste(x, expression(paste(Delta, values
Many thanks.
Regards,
--
Adaikalavan Ramasamy[EMAIL PROTECTED]
Centre for Statistics in Medicine
), list(x = x)))
-roger
Wolski wrote:
Hi!
plot(1:10, main=expression(paste( x, , Delta, values )))
/E
*** REPLY SEPARATOR ***
On 9/6/2004 at 3:50 PM Adaikalavan Ramasamy wrote:
I have been struggling with this problem for a while and I hope
Using append=TRUE in write.table only appends rows (see example below).
You might be interested in a recent tread about appending to save()
http://tolstoy.newcastle.edu.au/R/help/04/07/1467.html
write.table( matrix(1:9, 3), file=aaa.txt, sep=\t )
write.table( matrix(101:109, 3), file=aaa.txt,
Look into the code of power.t.test in the stats package. For example,
the sample size for two-sample t-test, two-tail testing and strict
interpretation of tail probability can be found by solving the following
equation iteratively :
\begin{equation}
1 - \beta = \Pr ( t_{v,ncp}^{*} t_{v,
201 - 300 of 412 matches
Mail list logo