Hello,
I have a set of data frames, generated by an SQL query that I am working
with. Because of the way the query was written, zero values for the
dependent variable (V2 in the example) are not recorded. Up until now
this has not been a problem.
I would like to be able to fill all absent data
Hello,
I am plotting a set of data where subsets are of different length and
not all sites are noted in the data sets. To indicate the extent of
sites, I am marking this with a line on the xyplot being produced:
require(lattice)
lengths-as.integer(runif(10)*100)
results - matrix(ncol=3)
for (i
Thank you, that works perfectly and is more elegant - I'm sorry not to
have included complete code/data.
I have another question that I have been having problems with recently,
described by this code:
reads-as.data.frame(cbind(c(1:25),rnorm(25),rnorm(25),rnorm(25),rnorm(25)))
Hi,
I am trying to plot a set of means+/-SD calculated by an external
program (an RDBMS) since the data set is too large to happily fit in R
(740M x 100 values - which are summarised to 100 means/SD by the DB). I
want to have a mean with whiskers at +/-1SD.
Can anyone suggest a way to do this?
Thank you for that, you have helped amazingly.Are you also able to point
me to a barplot equivalent?
Dan
On Thu, 2009-11-19 at 18:53 -0500, David Winsemius wrote:
install.packages(plotrix)
library(plotrix)
?plotCI
__
R-help@r-project.org
chromosomes
id refseq namelength
1 0 NC_01.9 Homo sapiens chromosome 1 247249719
2 1 NC_02.10 Homo sapiens chromosome 2 242951149
3 2 NC_03.10 Homo sapiens chromosome 3 199501827
4 3 NC_04.10 Homo sapiens chromosome 4 191273063
5 4
mistake.
On Thu, 2009-11-19 at 22:23 -0500, David Winsemius wrote:
On Nov 19, 2009, at 7:38 PM, Dan Kortschak wrote:
Thank you for that, you have helped amazingly.Are you also able to
point
me to a barplot equivalent?
It is a great puzzle to me that people do not work through
probably not the nicest way to do it, but it works.
hits-cbind(hits,chromosomes$name[hits$Chromosome])
names(hits)[5]-'name'
hits$name-factor(hits$name,ordered(chromosomes$name[1:dim(chromosomes)[1]]))
xyplot(weighted~bin|name,data=hits)
On Fri, 2009-11-20 at 12:16 +1030, Dan Kortschak wrote
Hi,
Am am plotting aggregated frequency data (extracted from an RDBMS)
relating to DNA sequence features for each of the human chromosomes as
described in the table chromosomes below (the frequency data is in a
table 'hits' that has a value (or not) for each of a set of bins across
each
Hi Deepayan,
I'm building the pipeline to generate the data, so that can easiy be
inserted into the code upstream.
Thanks for your help.
Dan
On Sun, 2009-10-04 at 12:16 -0700, Deepayan Sarkar wrote:
On Thu, Oct 1, 2009 at 11:30 PM, Dan Kortschak
dan.kortsc...@adelaide.edu.au wrote:
Thanks
-29 at 23:17 -0700, Deepayan Sarkar wrote:
On Tue, Sep 29, 2009 at 11:09 PM, Dan Kortschak
dan.kortsc...@adelaide.edu.au wrote:
Hi Deepayan,
Thanks for that, I had a think (a few hour too late) and came to the
same conclusion. I had wanted to have vertical bars for each of the
data
at 22:59 -0700, Deepayan Sarkar wrote:
On Wed, Sep 16, 2009 at 2:02 PM, Dan Kortschak
dan.kortsc...@adelaide.edu.au wrote:
Hi, I trying to produce a bar chart describing hits to specific bins by
chromosome for a large data set (I am asking here because
experimentation with options is precluded
).
My question is how I can get a sane number of labels along the x axis
(about 5-10 would be good). I suspect that 'scales' might be the way to
go, but I can't see an obvious option to do this.
thanks
--
Dan Kortschak dan.kortsc...@adelaide.edu.au
Hello R users,
I have a data set which is a set of lengths and types of objects. I want
to calculate the mean length for each type of object as opposed to the
mean of all the objects in the set.
This is in order to make a comparison between the lengths of each type
of objects and the number of
Perfect. In conjunction with Jorge's contrib that works a treat. Thanks.
Dan
On Tue, 2009-03-24 at 19:00 -0400, David Winsemius wrote:
?try
?tryCatch
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
I'm not really looking for a needle in a haystack, there are a small
number of the 60 tests (about 20) that are likely to concord with other
experiments I have, and in a particular pattern. Since I already have
the data in tables for graphic depiction, I was hoping to have a
reasonably easy way to
Hi R users,
I have a very large data set that has two conditioning variables for the
test I want to perform.
A toy set can be simulated:
type-sample(1:3,100,replace=TRUE)
class-sample(1:20,100,replace=TRUE)
value-rnorm(100)
data-cbind(type,class,value)
(though type and class are alphanum)
I
Hi Jorge,
That is exactly what I wanted - I should have given a reasonable
number of observations (my set has *almost* all paired observations,
so it will still break with that approach unless I manicure the data
set). Is there a way to fail nicely on a single one of the tests
without
That is a valid point, the number of samples I expect to be different
is actually quite small, but it is supportable (or otherwise) by
other experimental data.
Unfortunately the question I really want answered is pretty much
covered by doing this.
thanks
Dan
On 25/03/2009, at 10:25 AM,
Hi,
I'm making a multipanel lattice densityplot figure with 2 factors (3 and
20 classes in each factor) with the following statement (the
type=percent is there to prevent plotting the actual points which
detract from the figure - is there another way of doing this?):
densityplot(~End-Begin |
Hi,
This is another question relating to my 2 factor figure.
densityplot(~End-Begin | Type * Chromosome, data=Mon, layout=c(5,12),
xlab=Element Length,type=percent, col=grey60,
strip=strip.custom(style=3, bg=grey90, par.strip.text=list(cex=0.5)))
I would like to flip the plot so those at the
Perfect! Thanks.
Dan
On Mon, 2009-03-23 at 08:35 +0100, Dimitris Rizopoulos wrote:
try setting the 'as.table' argument to TRUE, e.g.,
densityplot(~End-Begin | Type * Chromosome, data=Mon, layout=c(5,12),
xlab=Element Length, type=percent, col=grey60,
strip=strip.custom(style=3, bg=grey90,
That's part way to my desired solution, but not completely there.
Here's an example:
chromosome-rep(c(A,X,Y),time=20)
type-rep(c(1:20),times=3)
length-rnorm(60)
densityplot(~length | type * chromosome, layout=c(5,12))
What I would like to see is the chromosome strip (A, X, Y) once on the
left
23 matches
Mail list logo