ermore National Laboratory
> 7000 East Ave., L-627
> Livermore, CA 94550
> 925-423-1062
>
>
>
>
>
> On 8/17/12 9:28 AM, "Scott Duke-Sylvester"
> wrote:
>
>>I have a quick question: It appears that in rgdal v0.7-12 (R version
>>2.15.1, OSX
1] "ID" "ANID" "F_AREA" "Avg_z" "Manning"
[6] "IniWL" "IniSal" "DispCoeff" "HydRad" "veg1_DW"
[11] "veg2_IW""veg3_SW""veg4_DC&q
reg Snow
Sent: Saturday, January 28, 2012 8:12 PM
To: Duke; r-help@r-project.org
Subject: Re: [R] percentage from density()
If you use logspline estimation (logspline package) instead of kernel density
estimation then this is
simple as there are cumulative area functions for logspline fits.
If you ne
Hi folks,
I know that density function will give a estimated density for a give
dataset. Now from that I want to have a percentage estimation for a
certain range. For examle:
> y = density(c(-20,rep(0,98),20))
> plot(y, xlim=c(-4,4))
Now if I want to know the percentage of data lying in (-20
Thanks for your response, Marc. HG and LG are high-grade/low-grade tumors.
The data has not been collected yet, but will be soon. It's all archived
data that will be pulled from computer records. The IRB wants some mention
of power or sample size, but doing it for this scenario has been a bit o
useR's,
I am trying to do a power calculation for a survival analysis using a
logrank test and I need some help properly doing this in R. Here is the
information that I know:
- I have 2 groups, namely HG and LG
- Retrospective analysis with subjects gathered from archival data over 20
years. No ne
On 9/15/11 4:40 PM, Duke wrote:
Hi folks,
Please let me know what I am doing wrong. I want to have a legend with
symbols that are filled with same color as the drawn line, but I
failed to do that:
plot(1:100, 1:100, pch=21, bg="red")
legend("bottomright", "test&q
Hi folks,
Please let me know what I am doing wrong. I want to have a legend with
symbols that are filled with same color as the drawn line, but I failed
to do that:
plot(1:100, 1:100, pch=21, bg="red")
legend("bottomright", "test", bty='n', pch=21, bg="red", col="red")
It looks to me that th
On 2/24/11 1:21 PM, David Winsemius wrote:
On Feb 24, 2011, at 1:03 PM, Duke wrote:
Hi folks,
I am wondering if the following is possible: I want to control the
arguments sent to a function by string variables. For example,
instead of
> heatmap.2( A, col=greenred(75) )
I would want
Hi folks,
I am wondering if the following is possible: I want to control the
arguments sent to a function by string variables. For example, instead of
> heatmap.2( A, col=greenred(75) )
I would want to have something like:
> heatmap.2 ( paste(A, "col=greenred(75)", sep=",") )
Is this possible
Hi all,
I am new to heatmap function in stats package, so any advice will be
greatly appreciated. So I have a dataset that I want to draw a heatmap,
for example:
> dput(a)
structure(c(-8.256023, -14.764429, -10.145779, -5.123104, -11.974729,
-16.329029, -16.329029, -16.329029, -2.764114, -4.9
On 1/12/11 6:44 PM, Duke wrote:
Thanks so much for your suggestion Martin. I had Bioconductor
installed but I honestly do not know all its applications. Anyway, I
am testing GenomicRanges with my data now. I will report back when I
get the result.
I got the results. My code took ~ 580
On 1/12/11 6:12 PM, Martin Morgan wrote:
The Bioconductor project has many tools for dealing with
sequence-related data. With the data
k <- read.table(textConnection(
"chr132375463237547rs523104280+
chr132375493237550rs520975820+
chr24513326451332
Hi folks,
I am working on a project that requires subsetting of a found file based
on some known file. The known file contains several lines like below:
chr132375463237547rs523104280+
chr132375493237550rs520975820+
chr245133264513327rs2976928
most subsetting/extractions
should be in less than a second unless there is something funny in
your data. So provide the 'str' so we can see.
On Thu, Dec 30, 2010 at 11:28 AM, Duke wrote:
Hi Jim,
Is this really a problem for me to use [1] instead of [[1]]? Will this make
it run slow
Actually there are different ways of doing subsetting:
[1]
[[1]]
[,1]
$V1
Please let me know which one is the fastest (and most used) one. Thanks.
D.
On 12/30/10 11:28 AM, Duke wrote:
Hi Jim,
Is this really a problem for me to use [1] instead of [[1]]? Will this
make it run slower? Also
.
$ runif.n..1: num 0.691 0.783 0.357 0.619 0.717 ...
$ runif.n..2: num 0.705 0.584 0.279 0.789 0.192 ...
$ runif.n..3: num 0.873 0.678 0.455 0.843 0.383 ...
On Thu, Dec 30, 2010 at 10:23 AM, Duke wrote:
Hi all,
First I dont have much experience with R so b
Hi all,
First I dont have much experience with R so be gentle. OK, I am dealing
with a dataset (~ tens of thousand lines, each line ~ 10 columns of
data). I have to create some subset of this data based on some certain
conditions (for example, same first column with another dataset etc...).
H
Hi Laurent,
On 12/7/10 5:01 PM, Laurent Gatto wrote:
Dear Duke,
I'm not sure if limma can generate colourful Venn diagrams. As an
alternative, I can however recommend Vennerable [1].
Thanks. I also found out that there are so many packages drawing venn
diagram now. I am playing
Hi all,
I am trying to plot a venn diagram using vennDiagram in limma package,
but I have no idea how to have it colorful such as the color one created
by the following website: http://www.pangloss.com/seidel/Protocols/venn.cgi.
Any one has any advice?
Thanks,
D.
__
ike matlab, R or python.
"Duke" wrote:
On 9/16/10 5:00 PM, Rolf Turner wrote:
On 17/09/2010, at 8:51 AM, Duke wrote:
Hi Duncan,
On 9/16/10 3:47 PM, Duncan Murdoch wrote:
On 16/09/2010 3:40 PM, Duke wrote:
Hi all,
I am writing a function (fun.R), but I dont know how
On 9/16/10 5:00 PM, Rolf Turner wrote:
On 17/09/2010, at 8:51 AM, Duke wrote:
Hi Duncan,
On 9/16/10 3:47 PM, Duncan Murdoch wrote:
On 16/09/2010 3:40 PM, Duke wrote:
Hi all,
I am writing a function (fun.R), but I dont know how to code the
function so that the Help Text will be shown
Hi Duncan,
On 9/16/10 3:47 PM, Duncan Murdoch wrote:
On 16/09/2010 3:40 PM, Duke wrote:
Hi all,
I am writing a function (fun.R), but I dont know how to code the
function so that the Help Text will be shown up when one types ?fun (of
course, after he loads it up). Anyone has any advice
Hi all,
I am writing a function (fun.R), but I dont know how to code the
function so that the Help Text will be shown up when one types ?fun (of
course, after he loads it up). Anyone has any advice for me how to do that?
Thanks,
D.
__
R-help@r-pr
On 9/10/10 4:24 PM, Gabor Grothendieck wrote:
On Fri, Sep 10, 2010 at 4:20 PM, Duke wrote:
On 9/10/10 2:49 PM, Gabor Grothendieck wrote:
On Fri, Sep 10, 2010 at 1:24 PM, Dukewrote:
Hi all,
I have to filter a tab-delimited text file like below:
"GeneNames""val
On 9/10/10 2:49 PM, Gabor Grothendieck wrote:
On Fri, Sep 10, 2010 at 1:24 PM, Duke wrote:
Hi all,
I have to filter a tab-delimited text file like below:
"GeneNames""value1""value2""log2(Fold_change)"
"log2(Fold_change) n
Hi Phil,
On 9/10/10 1:45 PM, Phil Spector wrote:
Duke -
One possibility is to check the help files for the functions
involved to see if there are options to control this behaviour.
For example, the check.names= argument to read.table, or the quote=
argument to write.table. How about
Yes
Hi all,
I have to filter a tab-delimited text file like below:
"GeneNames""value1""value2""log2(Fold_change)"
"log2(Fold_change) normalized""Signature(abs(log2(Fold_change)
normalized) > 4)"
ENSG0209350435-3.81131293562629-4.14357714689656TRUE
ENSG0
On 9/6/10 12:18 PM, Stephen Liu wrote:
Hi folks,
Debian 504 64-bit
I found following document;
http://www.biostat.wisc.edu/~kbroman/Rintro/
Whether it is the right document for installing Emacs+ESS and R so that R can
run on Emacs?
A lot on the web. Google gave me for example:
http://www.u
On 4/28/10 10:16 AM, peter dalgaard wrote:
On Apr 28, 2010, at 4:03 PM, Duke wrote:
I just did it, and it said it ignored the arg:
$ R R_PROFILE_USER=try.R
ARGUMENT 'R_PROFILE_USER=try.R' __ignored__
...
What am I doing wrong?
Using Windows syntax
On 4/28/10 9:53 AM, Duncan Murdoch wrote:
On 28/04/2010 9:38 AM, Duke wrote:
On 4/27/10 8:04 PM, Duncan Murdoch wrote:
> On 27/04/2010 6:15 PM, Duke wrote:
>> Hi all,
>>
>> I am testing R by trying R on command line (ie without opening R
>> console) on a Mac Pro wit
On 4/27/10 8:04 PM, Duncan Murdoch wrote:
On 27/04/2010 6:15 PM, Duke wrote:
Hi all,
I am testing R by trying R on command line (ie without opening R
console) on a Mac Pro with 10.6. So I created a file plot.R:
x11(width=5,height=3.5)
plot(1:10,1:10)
So I have two questions below (please
On 4/27/10 6:29 PM, mtmor...@fhcrc.org wrote:
Welcome to R!
But actually, I should be saying 'welcome to Bioconductor!' and you
should be on the Bioconductor (or Bioc-sig-seq) mailing list
http://bioconductor.org/docs/mailList.html
because ShortRead is a Bioconductor package.
BWA produces
Hi all,
I am testing R by trying R on command line (ie without opening R
console) on a Mac Pro with 10.6. So I created a file plot.R:
x11(width=5,height=3.5)
plot(1:10,1:10)
So I have two questions below (please forgive me if the questions have
been asked already)
1. What is the different
Dear folks,
Please welcome a newbie both to R and the mailing list :). I am
currently working on a sequencing project, and heard about R as well as
some of its packages for next gen sequencing, and decided to give it a
try. Starting with ShortRead, I found a document
(http://www.bioconductor.
de or outside the polygons defined in the shapefile.
Many thanks,
scott.
--
Scott M. Duke-Sylvester
Assistant Professor
Department of Biology
Office : 300 E. St. Mary Blvd
Billeaud Hall, Room 141
Lafayette, LA 70504
Mailing address : UL Lafayette
Depar
36 matches
Mail list logo
