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 wi
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
Hello,
I am trying to use lattice barchart so that it will draw a barplot where
the bars do no reach the frame of the chart, but only extend to 0 as the
minimum. With barplot, I can set this by using the offset value, but I
cannot see how to do this with barchart.
Using the example in the barchar
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)))
names(reads)<-c('id
It's 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:
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 t
> 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
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
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?
t
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 chromosome
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
> wrote:
> > Thanks for
e, 2009-09-29 at 23:17 -0700, Deepayan Sarkar wrote:
> On Tue, Sep 29, 2009 at 11:09 PM, Dan Kortschak
> 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 ea
, 2009-09-29 at 22:59 -0700, Deepayan Sarkar wrote:
> On Wed, Sep 16, 2009 at 2:02 PM, Dan Kortschak
> 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
> > experimentatio
density of the labels).
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
_
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 th
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
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 t
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 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 the
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)
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 th
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,
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
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 | Ty
24 matches
Mail list logo