on wiki for easy access
of some simple terms.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of Ista Zahn
Sent: Tuesday, February 02, 2010 12:03 PM
To: Peng Yu
Cc: r-h...@stat.math.ethz.ch
Subject: Re: [R] population variance
On Thu, Feb 4, 2010 at 4:34 PM, Henrik Bengtsson h...@stat.berkeley.edu wrote:
On Thu, Feb 4, 2010 at 9:18 PM, Liaw, Andy andy_l...@merck.com wrote:
From: Peng Yu
On Wed, Feb 3, 2010 at 10:01 AM, Peng Yu pengyu...@gmail.com wrote:
Some examples in the help page are too long to be copied
On Tue, Feb 2, 2010 at 11:22 PM, David Winsemius dwinsem...@comcast.net wrote:
On Feb 3, 2010, at 12:20 AM, Peng Yu wrote:
On Tue, Feb 2, 2010 at 11:14 PM, David Winsemius dwinsem...@comcast.net
wrote:
On Feb 3, 2010, at 12:11 AM, Peng Yu wrote:
On Tue, Feb 2, 2010 at 11:04 PM, David
On Wed, Feb 3, 2010 at 2:12 AM, Emmanuel Charpentier
charp...@bacbuc.dyndns.org wrote:
Le mercredi 03 février 2010 à 00:01 -0500, David Winsemius a écrit :
On Feb 2, 2010, at 11:38 PM, Peng Yu wrote:
?contrast in the contrast package gives me the following description.
However, I have
On Wed, Feb 3, 2010 at 8:41 AM, David Winsemius dwinsem...@comcast.net wrote:
On Feb 3, 2010, at 9:28 AM, Peng Yu wrote:
On Tue, Feb 2, 2010 at 11:22 PM, David Winsemius dwinsem...@comcast.net
wrote:
On Feb 3, 2010, at 12:20 AM, Peng Yu wrote:
On Tue, Feb 2, 2010 at 11:14 PM, David
On Wed, Feb 3, 2010 at 10:01 AM, Peng Yu pengyu...@gmail.com wrote:
Some examples in the help page are too long to be copied from screen.
Could somebody let me know some easy way on how to extract the example
to a file so that I can play with them?
I forget to mention. I use a terminal version
Some examples in the help page are too long to be copied from screen.
Could somebody let me know some easy way on how to extract the example
to a file so that I can play with them?
__
R-help@r-project.org mailing list
Suppose that I have the following list of lists of frames 'root'
(let's call it a 'tree' of frames). I want to flatten it to be a list
of frames. However, if I unlist(root), it will flatten the frames as
well. Is there a simply way to flatten the tree to certain depth?
On Wed, Feb 3, 2010 at 12:29 PM, Steve Lianoglou
mailinglist.honey...@gmail.com wrote:
Hi,
On Wed, Feb 3, 2010 at 1:19 PM, Peng Yu pengyu...@gmail.com wrote:
Suppose that I have the following list of lists of frames 'root'
(let's call it a 'tree' of frames). I want to flatten it to be a list
Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Peng Yu
Sent: Wednesday, February 03, 2010 10:20 AM
To: r-h...@stat.math.ethz.ch
Subject: [R] How to flatten a tree (based on list) to a certain depth?
Suppose that I have the following list
-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Peng Yu
Sent: Wednesday, February 03, 2010 11:40 AM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] How to flatten a tree (based on list) to a certain depth?
On Wed, Feb 3, 2010 at 12:26 PM, Bert Gunter gunter.ber...@gene.com wrote
Hi Chad,
There are not many significant transcript clusters after p.adjust.
sum(adjusted_pvalue[,'C-D'].05)
[1] 0
sum(adjusted_pvalue[,'A-D'].05)
[1] 0
sum(adjusted_pvalue[,'D-E'].05)
[1] 2
sum(adjusted_pvalue[,'A-B'].05)
[1] 8
However, there are may TCs that are significant before
On Wed, Feb 3, 2010 at 3:48 PM, Max Kuhn mxk...@gmail.com wrote:
Wait, what were we talking about? Right...it is a typo. It should be
c'beta in the numerator.
Peng: As the package maintainer, you really should send me a quick
email about it instead of posting to the list.
Let's not waste
pe...@morgan:~/projects/APP_Hippo_Exon_Jankowsky/analysis/pvalue_contrast_table.RData
__
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
On Wed, Feb 3, 2010 at 4:53 PM, David Scott d.sc...@auckland.ac.nz wrote:
Peng Yu wrote:
On Wed, Feb 3, 2010 at 3:48 PM, Max Kuhn mxk...@gmail.com wrote:
Wait, what were we talking about? Right...it is a typo. It should be
c'beta in the numerator.
Peng: As the package maintainer, you
bcp()'s help says it implements the paper by Barry and Hartigan
(1993). But it has 4 other citations which are later than 1993. Could
somebody who have experience on this function let me know why these 4
later citations matters to this particular function? Or these
citations are just for the topic
On Mon, Oct 19, 2009 at 12:53 PM, Kingsford Jones
kingsfordjo...@gmail.com wrote:
sum((x-mean(x))^2)/(n)
[1] 0.4894708
((n-1)/n) * var(x)
[1] 0.4894708
But this is not a built-in function in R to do so, right?
hth,
Kingsford
On Mon, Oct 19, 2009 at 9:30 AM, Peng Yu pengyu...@gmail.com
?contrast in the contrast package gives me the following description.
However, I have no idea what Type II and III contrasts are. Could
somebody explain it to me? And what does 'type' mean here?
*‘type’*: set ‘type=average’ to average the individual contrasts
(e.g., to obtain a Type II or
It seems that Eq (2) in the vignettes for the 'contrast' packages is
not correct. That is, the numerator on the right hand side should be
$c' \beta$ rather than $c' \lambda$, right? If I'm correct, could
somebody notice the author to fix it?
__
library(contrast)
contrast
function (fit, ...)
UseMethod(contrast)
environment: namespace:contrast
I guess the above function is S3. Could somebody let me know how to
show the function body in an R session?
__
R-help@r-project.org mailing list
contrast.glscontrast.lm
[5] contrast.lme
methods(contrast.lm)
no methods were found
Warning message:
In methods(contrast.lm) :
function 'contrast.lm' appears not to be generic
contrast.lm
function (fit, ...)
contrastCalc(fit, ...)
environment: namespace:contrast
On 2/2/2010 8:48 PM, Peng Yu
I read the vignette of contrast package. I don't think that I
understand how to use it.
I made the following simpler example (contrast between '3' and '4').
Could somebody let me know what is the correct way to compute the
contrast in the following example?
library(contrast)
a=3
n=4
A =
On Tue, Feb 2, 2010 at 11:04 PM, David Winsemius dwinsem...@comcast.net wrote:
On Feb 2, 2010, at 11:42 PM, Peng Yu wrote:
It seems that Eq (2) in the vignettes for the 'contrast' packages is
not correct. That is, the numerator on the right hand side should be
$c' \beta$ rather than $c
On Tue, Feb 2, 2010 at 11:14 PM, David Winsemius dwinsem...@comcast.net wrote:
On Feb 3, 2010, at 12:11 AM, Peng Yu wrote:
On Tue, Feb 2, 2010 at 11:04 PM, David Winsemius dwinsem...@comcast.net
wrote:
On Feb 2, 2010, at 11:42 PM, Peng Yu wrote:
It seems that Eq (2) in the vignettes
Could somebody recommend some good nonparametric statistical inference
textbooks for a beginner? And what are pros and cons of each book?
Nonparametric statistical methods by Hollander seems to be more
difficult for a beginner, but is great as a reference, right? Are
there any books that are
I'm wondering if there is a way to make blocks of code independent
from each other. Please see the following example for what I mean.
x=1
## is there a way to make the following assignment not affect the above x?
## in C++, I can use {} to make the effect local. Is there an
equivalent construct
By default, my R functions run silently if NA is encountered? I would
like them to give me some message when NA is encountered by default.
Note that I could check each intermediate variable by something like
is.na(). But this is not manageable for big programs, so I'd like a
way to check NA
On Fri, Jan 22, 2010 at 3:33 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
On Thu, Jan 21, 2010 at 1:06 PM, Peng Yu pengyu...@gmail.com wrote:
On Thu, Jan 21, 2010 at 2:41 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
On Thu, Jan 21, 2010 at 12:29 PM, Peng Yu pengyu...@gmail.com wrote
On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris gpet...@uark.edu wrote:
I have found a good refernce to be S Programming by Venables and
Ripley.
I'll take a look at this book.
Since S and R are not completely the same, there are delicate
differences between S and R, which an S book may be
I found this paper on ANOVA on unequal error variance. Has this be
incorporated to any R package? Is there any textbook that discuss the
problem of ANOVA on unequal error variance in general?
http://www.jstor.org/stable/2532947?cookieSet=1
__
start with, if I want to
understand the method?
Do you have any simple explanation that may help me understand what is
the difference between the method in 'Exact Statistical Methods for
Data Analysis' and the method in gls()?
HTH,
Dennis
On Thu, Jan 21, 2010 at 12:03 PM, Peng Yu pengyu
On Thu, Jan 21, 2010 at 2:41 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
On Thu, Jan 21, 2010 at 12:29 PM, Peng Yu pengyu...@gmail.com wrote:
On Thu, Jan 21, 2010 at 2:16 PM, Dennis Murphy djmu...@gmail.com wrote:
Hi:
This paper was a prelude to his first book 'Exact Statistical
Please see the following example. I can not write '' to a csv file
successfully. Could somebody let me if it is possible to write '' to
a csv file with the default setting of write.csv?
my_home$ Rscript main_quote.R
x=rbind(
+ \A\
+ , \B\
+ )
x
[,1]
[1,] \A\
[2,] \B\
write.csv(x,
.
I have looking for a way so that the resulted csv file can be read by excel.
http://www.creativyst.com/Doc/Articles/CSV/CSV01.htm#FileFormat
On Thu, Jan 21, 2010 at 6:29 PM, Peng Yu pengyu...@gmail.com wrote:
Please see the following example. I can not write '' to a csv file
successfully
I got the following warnings when I install R.oo. Are these warnings
normal? Should I reinstall the package as mentioned in the warnings?
How to reinstall? The sessionInfo() is at the end.
install.packages(R.oo, dependencies=T)
Warning in install.packages(R.oo, dependencies = T) :
argument
Apparently, Bengtsson.pdf is not current any more. For example, the
installation method is outdated. R.oo.pdf is a reference manual rather
than a tutorial. I'm wondering if there is a better and more current
R.oo tutorial so that I can quickly get started on R.oo.
R-lang.pdf has the following description in Section 3.1.1.
Any number typed directly at the prompt is a constant and is evaluated.
1
[1] 1
Perhaps unexpectedly, the number returned from the expression 1 is a
numeric. In most
cases, the difference between an integer and a numeric value will be
According to R-lang.pdf (Section 2):
Function mode gives information about the mode of an object in the
sense of Becker,
Chambers Wilks (1988), and is more compatible with other
implementations of the S
language. Finally, the function storage.mode returns the storage mode
of its argument
in the
I don't find a tutorial on S3. Bengtsson.pdf cites MASS (1999
edition). However, I don't think that MASS (2002 edition) clearly
explain what S3 is and help a user who knew very little about S3 to
quickly understand it. Could somebody let me know if there are some
better learning materials to help
On Wed, Jan 20, 2010 at 5:04 PM, Peng Yu pengyu...@gmail.com wrote:
I don't find a tutorial on S3. Bengtsson.pdf cites MASS (1999
edition). However, I don't think that MASS (2002 edition) clearly
explain what S3 is and help a user who knew very little about S3 to
quickly understand it. Could
On Wed, Jan 20, 2010 at 5:30 PM, Steve Lianoglou
mailinglist.honey...@gmail.com wrote:
Hi,
On Wed, Jan 20, 2010 at 6:05 PM, Peng Yu pengyu...@gmail.com wrote:
On Wed, Jan 20, 2010 at 5:04 PM, Peng Yu pengyu...@gmail.com wrote:
I don't find a tutorial on S3. Bengtsson.pdf cites MASS (1999
I don't find where in the R document the discussion of nested
namespace is. If there is nested namespace supported in R, could
somebody let me know whether the document is?
__
R-help@r-project.org mailing list
On Thu, Jan 7, 2010 at 11:10 AM, Peng Yu pengyu...@gmail.com wrote:
I don't find where in the R document the discussion of nested
namespace is. If there is nested namespace supported in R, could
somebody let me know whether the document is?
Could somebody let me know if there is nested name
read.table terminates the program if the input file is empty. Is there
way to let the program continue and return me a NULL instead of
terminating the program?
$ Rscript read_empty.R
read.table(empty_data.txt)
Error in read.table(empty_data.txt) : no lines available in input
Execution halted
$
.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
/Henrik
On Fri, Jan 1, 2010 at 12:41 PM, Peng Yu pengyu...@gmail.com wrote:
read.table terminates the program if the input file is empty. Is there
way to let the program continue and return me
))
print(g)
V1 V2 V3 V4
1 name title1 title2 title3
2 row1 10.1 20.1 30.1
3 row2 12.1 22.1 32.1
On Fri, Jan 1, 2010 at 12:41 PM, Peng Yu pengyu...@gmail.com wrote:
read.table terminates the program if the input file is empty. Is there
way to let the program
PM, Ted Harding ted.hard...@manchester.ac.uk
wrote:
On 01-Jan-10 20:41:52, Peng Yu wrote:
read.table terminates the program if the input file is empty. Is there
way to let the program continue and return me a NULL instead of
terminating the program?
$ Rscript read_empty.R
read.table
On Fri, Jan 1, 2010 at 3:19 PM, Peng Yu pengyu...@gmail.com wrote:
On Fri, Jan 1, 2010 at 2:54 PM, Dylan Beaudette
dylan.beaude...@gmail.com wrote:
?try
This works. Thank you!
f=try(read.table(empty_data.txt))
Error in read.table(empty_data.txt) : no lines available in input
One more
I don't see where describes the implementation of '[]'.
For example, if x is a matrix or a data.frame, how the lookup of
'colname1' is x[, 'colname1'] executed. Does R perform a lookup in the
a hash of the colnames? Is the reference O(1) or O(n), where n is the
second dim of x?
R internal resolve
which columns to use if I specify the column names rather than the
indexes. If R does it by search, the access time should be O(log(n))
where n is the number of columns. If R does it by hash, the access
time should be O(1).
On Thu, Dec 31, 2009 at 5:15 PM, Peng Yu pengyu
wrote:
Dear Peng Yu
how big is the RAM of your computer? You could try with closing all other
applications before running this script. You could try on a server with more
RAM.
I tried downloading the file whose URL who give below, but gave up after
some failed rounds with the extraordinarily
, Benilton Carvalho bcarv...@jhsph.edu wrote:
pm(data)
b
On Dec 26, 2009, at 2:21 PM, Peng Yu wrote:
I use the following code to do RMA. I'm wondering how get the probe
level values before the summary to the probeset level values.
library(oligo)
data-read.celfiles(list.celfiles(data_dir
I only want to load a limited number of rows by dbReadTable(). I don't
see an option in the help. Is there an option to do so?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Suppose that I 'library()' a package, in the 'NAMESPACE' file of the
package, it export some functions to the global namespace. Is there a
way to 'library()' the package without export the functions in
'NAMESPACE' to the global namespace?
__
On Fri, Dec 18, 2009 at 3:33 AM, Gustaf Rydevik
gustaf.ryde...@gmail.com wrote:
On Thu, Dec 17, 2009 at 4:33 PM, Peng Yu pengyu...@gmail.com wrote:
On Thu, Dec 17, 2009 at 5:33 AM, Gustaf Rydevik
gustaf.ryde...@gmail.com wrote:
On Wed, Dec 16, 2009 at 10:13 PM, Peng Yu pengyu...@gmail.com
I don't find a function to print a string to file. Would somebody let
me know what function I should use?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On Thu, Dec 17, 2009 at 5:33 AM, Gustaf Rydevik
gustaf.ryde...@gmail.com wrote:
On Wed, Dec 16, 2009 at 10:13 PM, Peng Yu pengyu...@gmail.com wrote:
Currently, I load the RData file then ls() and str(). But loading the file
takes too long if the file is big. Most of the time, I only interested
Currently, I load the RData file then ls() and str(). But loading the file
takes too long if the file is big. Most of the time, I only interested what
the variables are in the the file and the attributes of the variables (like
if it is a data.frame, matrix, what are the colnames/rownames, etc.)
On Tue, Dec 15, 2009 at 10:32 PM, hadley wickham h.wick...@gmail.com wrote:
I don't understand what these addresses mean. Would you please help me
understand it?
Did you try reading the documentation?
When an object is traced any copying of the object by the C
function ‘duplicate’
2009/12/15 tlum...@u.washington.edu:
On Wed, 16 Dec 2009, Peng Yu wrote:
On Tue, Dec 15, 2009 at 10:32 PM, hadley wickham h.wick...@gmail.com
wrote:
I don't understand what these addresses mean. Would you please help me
understand it?
Did you try reading the documentation?
When
Please see the following example. If I change the length of a vector
twice below, I don't understand why 'tracemem' is shown in the first
case but not the second case. Could somebody help explain it to me?
x=1
tracemem(x)
[1] 0x1bad8c8
length(x)=100# The 'tracemem' string is shown as expected.
On Tue, Dec 15, 2009 at 2:34 PM, Peng Yu pengyu...@gmail.com wrote:
2009/12/15 tlum...@u.washington.edu:
On Wed, 16 Dec 2009, Peng Yu wrote:
On Tue, Dec 15, 2009 at 10:32 PM, hadley wickham h.wick...@gmail.com
wrote:
I don't understand what these addresses mean. Would you please help me
I'm wondering if lazy copy is available in R or not. For example, in
the following code, I'm wondering if the memory for y is allocated in
the 2nd line or the 3rd line. Is there a documentation for this?
x=1:1
y=x
y[[10]]=5
__
R-help@r-project.org
PM, Benilton Carvalho bcarv...@jhsph.edu wrote:
use tracemem() to figure out... and read its documentation in detail.
b
On Dec 15, 2009, at 1:03 AM, Peng Yu wrote:
I'm wondering if lazy copy is available in R or not. For example, in
the following code, I'm wondering if the memory for y
The following examples are confusing to me. It is OK, to assigned NULL
to one element in a list. The result is still a list. However, a list
of NULL's are reduced to NULL. I don't understand how this conversion
occurs. Could somebody let me know what is going on?
X=matrix(1:8, nr=4)
apply(X,1,
http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-set-components-of-a-list-to-NULL_003f
The explanation on this FAQ entry is not clear. It says '... similarly
for named components...'. What I understood was x[i]-list(NULL) is
the same as x$a_name-list(NULL). But, they are not. As the example
On Fri, Dec 11, 2009 at 11:01 AM, William Dunlap wdun...@tibco.com wrote:
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu
Sent: Friday, December 11, 2009 8:44 AM
To: r-h...@stat.math.ethz.ch
Subject: [R] Why a list
On Fri, Dec 11, 2009 at 11:27 AM, Steve Lianoglou
mailinglist.honey...@gmail.com wrote:
Hi,
On Dec 11, 2009, at 12:07 PM, Peng Yu wrote:
http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-set-components-of-a-list-to-NULL_003f
The explanation on this FAQ entry is not clear. It says
On Fri, Dec 11, 2009 at 11:43 AM, William Dunlap wdun...@tibco.com wrote:
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Peng Yu
Sent: Friday, December 11, 2009 9:18 AM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] Why a list
On Fri, Dec 11, 2009 at 11:51 AM, Steve Lianoglou
mailinglist.honey...@gmail.com wrote:
On Dec 11, 2009, at 12:36 PM, Peng Yu wrote:
[snip]
What seems confusing to me is:
even 'x[i]-list(NULL)' and 'x[[i]]-list(NULL)' are different, why
x[i]-NULL and x[[i]]-NULL are the same?
Shouldn't
On Fri, Dec 11, 2009 at 12:05 PM, hadley wickham h.wick...@gmail.com wrote:
A very common situation is that the users don't know all the possible
return types of 'some_third_party_function()'. If the users don't know
all the return types, he/she can not make sure the return type of
function(x)
mode and a non strict mode. If it is always better to return
the same type, which is the case in other strong typed language, then
it might better that there is such a mode in R to enforce the same
return type.
At 12:24 PM -0600 12/11/09, Peng Yu wrote:
On Fri, Dec 11, 2009 at 12:05 PM, hadley
On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
ch...@sharpsteen.net wrote:
On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu pengyu...@gmail.com wrote:
How do you figure out all the possibilities?
Well, the Value section of the third party function's help page should
outline the return types
On Fri, Dec 11, 2009 at 3:49 PM, Steve Lianoglou
mailinglist.honey...@gmail.com wrote:
On Dec 11, 2009, at 4:45 PM, Peng Yu wrote:
On Fri, Dec 11, 2009 at 2:37 PM, Charlie Sharpsteen
ch...@sharpsteen.net wrote:
On Fri, Dec 11, 2009 at 10:24 AM, Peng Yu pengyu...@gmail.com wrote:
How do you
I'm wondering where I can find the detailed descriptions on R memory
management. Understanding this could help me understand the runtime of
R program. For example, depending on how memory is allocated (either
allocate a chuck of memory that is more than necessary for the current
use, or allocate
:
Related...
Rule of thumb:
Pre-allocate your object of the *correct* data type, if you know the
final dimensions.
/Henrik
On Thu, Dec 10, 2009 at 8:26 AM, Peng Yu pengyu...@gmail.com wrote:
I'm wondering where I can find the detailed descriptions on R memory
management. Understanding this could
On Thu, Dec 10, 2009 at 3:46 AM, Karl Ove Hufthammer k...@huftis.org wrote:
On Wed, 9 Dec 2009 19:20:47 -0600 Peng Yu pengyu...@gmail.com wrote:
Is there a way to figure out which of these variants is actually
dispatched to when I call split? I know that if the argument is of the
type
In findInterval's help page, it says 'v[i[j]] = x[j] v[i[j] + 1]'.
I'm wondering if there is a variant of it such that 'v[i[j]] x[j] =
v[i[j] + 1]'.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
If I use system.time() to measure the runtime of an expression, I will
not get the result. Is there a way to measure the runtime and get the
result as well?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
a variation of this; preallocate and then extend when you read a
limit
x - numeric(2)
for (i in 1:100){
if (i length(x)){
# double the length (or whatever you want)
length(x) - length(x) * 2
}
x[i] - i
}
On Thu, Dec 10, 2009 at 11:30 AM, Peng Yu pengyu...@gmail.com
On Thu, Dec 10, 2009 at 7:03 PM, Dirk Eddelbuettel e...@debian.org wrote:
On 10 December 2009 at 18:12, Peng Yu wrote:
| If I use system.time() to measure the runtime of an expression, I will
| not get the result. Is there a way to measure the runtime and get the
| result as well?
Use
Is there a way to profile an R program similar to valgrind
(valgrind.org), in the sense that I can easily see which function is
the bottleneck of an R program?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
The following code returns a list with the 2nd element as NULL. I'm
wondering what the best way to get rid of NULL element in an
'apply()'s result.
lapply(1:3, function(x) {
+ if(x==2) {
+ return(NULL)
+ } else {
+ return(x)
+ }
+ }
+ )
[[1]]
[1] 1
On Thu, Dec 10, 2009 at 10:00 PM, Marc Schwartz marc_schwa...@me.com wrote:
On Dec 10, 2009, at 9:44 PM, Peng Yu wrote:
The following code returns a list with the 2nd element as NULL. I'm
wondering what the best way to get rid of NULL element in an
'apply()'s result.
lapply(1:3, function(x
I have following the message dim(refdata) and dimnames(refdata) no
longer allow parameter ref=TRUE, use dim(derefdata(refdata)),
dimnames(derefdata(refdata)) instead when I loaded data.table. Is it
from the package ref? Could it be fixed? Or there is something wrong
with my installation?
On Tue, Dec 8, 2009 at 11:06 PM, David Winsemius dwinsem...@comcast.net wrote:
On Dec 9, 2009, at 12:00 AM, Peng Yu wrote:
On Tue, Dec 8, 2009 at 10:37 PM, David Winsemius dwinsem...@comcast.net
wrote:
On Dec 8, 2009, at 11:28 PM, Peng Yu wrote:
I have the following code, which tests
On Tue, Dec 8, 2009 at 11:05 PM, David Winsemius dwinsem...@comcast.net wrote:
On Dec 8, 2009, at 11:37 PM, Peng Yu wrote:
I want to split a matrix, where both 'u' and 'w' are results of
possible ways. However, whenever 'n' changes, the argument passed to
mapply() has to change
On Wed, Dec 9, 2009 at 11:20 AM, Charles C. Berry cbe...@tajo.ucsd.edu wrote:
On Wed, 9 Dec 2009, Peng Yu wrote:
On Tue, Dec 8, 2009 at 11:06 PM, David Winsemius dwinsem...@comcast.net
wrote:
On Dec 9, 2009, at 12:00 AM, Peng Yu wrote:
On Tue, Dec 8, 2009 at 10:37 PM, David Winsemius
On Tue, Dec 8, 2009 at 10:37 PM, David Winsemius dwinsem...@comcast.net wrote:
On Dec 8, 2009, at 11:28 PM, Peng Yu wrote:
I have the following code, which tests the split on a data.frame and
the split on each column (as vector) separately. The runtimes are of
10 time difference. When m
I see 'library(stats)' at the beginning of
R-2.10.0/src/library/stats/tests/nls.R.
I'm wondering if I am developing my own package 'mypackage' whether I
should put 'library(mypackage)' in a .R file in mypackage/tests/? If I
do, then it seems awkward to me, because to use 'library(mypackage)',
I
There are a number of functions that are dispatched to from split().
methods('split')
[1] split.data.frame split.Date split.defaultsplit.POSIXct
Is there a way to figure out which of these variants is actually
dispatched to when I call split? I know that if the argument is of the
type
I make a version for matrix. Because, it would be more efficient to
split each column of a matrix than to convert a matrix to a data.frame
then call split() on the data.frame. Note that the version for a
matrix and a data.frame is slightly different. Would somebody add this
in R as well?
Sorry. I sent this to r-help by mistake. Could somebody help delete it
from the archive?
On Wed, Dec 9, 2009 at 9:29 PM, Peng Yu pengyu...@gmail.com wrote:
I make a version for matrix. Because, it would be more efficient to
split each column of a matrix than to convert a matrix to a data.frame
()
print(n)
Uwe Ligges
Peng Yu wrote:
I know that R is a dynamic programming language. But I'm wondering if
there is a way to make the assignment in a for-loop not affect
variables outside the loop.
n=10
for(i in 1:n){
+ n=3
+ print(n)
+ }
[1] 3
[1] 3
[1] 3
[1] 3
[1
I did a search on www.rseek.org to look for the function to test if a
vector is ordered or not. But I don't find it. Could somebody let me
know what function I should use?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
I'm doing some data manipulations. I thought originally that I should
use data.frame, as the elements in a data.frame can have multiple
types but the elements in a matrix has to be the same, which all the
data have to convert to strings if there is a single column that is
string.
However, when I
something productive by now.
See further comments interspersed below:
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of Peng Yu
Sent: Friday, December 04, 2009 2:06 PM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] grep() exclude
On Fri, Dec 4, 2009 at 6:12 PM, Gavin Simpson gavin.simp...@ucl.ac.uk wrote:
On Fri, 2009-12-04 at 15:18 -0600, Peng Yu wrote:
On Fri, Dec 4, 2009 at 3:06 PM, Peng Yu pengyu...@gmail.com wrote:
On Fri, Dec 4, 2009 at 2:35 PM, Greg Snow greg.s...@imail.org wrote:
The invert argument seems
I have the following code, which tests the split on a data.frame and
the split on each column (as vector) separately. The runtimes are of
10 time difference. When m and k increase, the difference become even
bigger.
I'm wondering why the performance on data.frame is so bad. Is it a bug
in R? Can
I want to split a matrix, where both 'u' and 'w' are results of
possible ways. However, whenever 'n' changes, the argument passed to
mapply() has to change. Is there a way to pass elements of a list as
multiple arguments?
m=10
n=2
k=3
set.seed(0)
x=replicate(n,rnorm(m))
f=sample(1:k, size=m,
1 - 100 of 351 matches
Mail list logo