On Tue, Oct 18, 2011 at 03:40:27PM +0200, Martin Batholdy wrote:
Ok, I think that would work – thanks!
However, in my case I read a data.frame via read.table().
So some of the columns get transformed to factors automatically – I don't
generate the factor-variables as in the example, so I
where the first column is the chromosome location and the second column is
some value. What I'd like to do is have a histogram created for each chr
location (i.e. a separate histogram for chr1, chr2, chr3, chr7, chr9, and
chr22). I am just having a hard time getting everything to work out
On Sun, Aug 21, 2011 at 09:18:25AM -0700, danielepippo wrote:
Sweave(example.Rtex) in R it seems working
[...]
* ...le/Desktop/dati/LaTeX1.Rtex*
Sounds like you are first running sweave on the file 'example.Rtex'
and later LaTex on 'LaTeX1.Rtex'
Two points:
1) why do these files have
On Thu, Aug 18, 2011 at 04:52:58PM +0700, Rut S wrote:
I tried to recode some complex multiple variables and run into a problem that
r can change only some column that I want to change.
I can reproduce the problem with this
idfortest -
guess is correct you should be able to fix your problem by
including quote= in both your read.table commands.
If it doesn't, also try comment.char= - another popular source of
problems.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
? Does your file ONLY contain
sequences - or are there also identifiers, annotations etc.?
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http://webclu.bio.wzw.tum.de
[, 'Sequence']
[1] NNATTAAAGGGC
str(dat)
'data.frame': 1 obs. of 1 variable:
$ Sequence: chr NNATTAAAGGGC
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http
.
http://lmgtfy.com/?q=r+volcano+plot
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
without resorting to lossy compression. You may have hit that
limit.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de
/expect4win/
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R
than you will
care for.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
which is the maximum large of digits that R has?, because SQL work
with 50 digits I think. and I need a software that work with a lot
of digits.
The .Machine() command will provide some insight into these matters.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte
On Fri, May 06, 2011 at 09:17:11AM -0400, David Winsemius wrote:
On May 6, 2011, at 4:03 AM, Philipp Pagel wrote:
The .Machine() command will provide some insight into these matters.
On my device (and I suspect on all versions of R) .Machine is a
built-in list and there is no .Machine
need to let us know where your
Hist() function came from.
hist() from package graphics does not have a scale parameter and
honours ylim without a problem.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum
message:
In open.connection(file, rt) : unsupported URL scheme
?url told me Note that ‘https://’ connections are not supported.
Case closed, problem unsolved...
Dirty workaround: use system() and wget or whatever command is available on
Windows for this.
cu
Philipp
--
Dr. Philipp Pagel
, on a single figure).
That's what I was about to suggest. Don't loop over the tasks - use
the power of lattice. I think this should be close to what you want:
bwplot(factor(maxthreads) ~ time | factor(tasks), x)
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
or fill='header' may exist. Did I just
not find it or does it not exist? And if it doesn't - does anyone
else think it would be a nice item for the wishlist?
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum
:1
Median :3.000 Median :6 NA's:1
Mean :2.667 Mean :6
3rd Qu.:3.500 3rd Qu.:7
Max. :4.000 Max. :8
NA's :1.000 NA's :2
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
but is slow because I
have used several loops to achieve this simple task.
Oh - and in case you ONLY wnt the number of NAs in each column this
should be pretty efficient:
lapply(foo, function(x){sum(is.na(x))})
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
On Mon, Apr 18, 2011 at 04:11:57PM +0530, Ramnath R wrote:
Hai
From which CRAN mirror can get the package ?LPP2005REC?
As the first hit of a google search for LPP2005REC told me it
is not a package but a dataset in package timeSeries.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für
.
see ?history for details.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
this can get you in trouble if y is a floating point
number because of the issues with identity of such numbers, as
discussed in ?all.equal and FAQ 7.31 Why doesn't R think these
numbers are equal?.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing
does not make
it less wrong: rm() will delete objects from the workspace but has
absolutely no effect on the history.
cu
Philipp
/...Tanvir Ahamed
━━━
From: Philipp Pagel p.pa...@wzw.tum.de
to find out what part of the program is actually eating up the most
time.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de
what you mean by integrating with these systems.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
$group2))
Perfect! Using the levels directly had not occured to me. Thanks!
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
the list of elements of
ineffFilesList?
as.data.frame(ineffFilesList)
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de
Of course I can now use the col argument to manually assign the colors
by group but is there a more elegant way that I missed?
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org
,]123
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
to the same file?
Sweave an odfWeave are very nice methods for generating reports with
both text, R code, Results from R and Graphics.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von
requires at least a filename which
seems to be missing from your comamnd above.
In addition, I recommend to post your question on the bioconductor
mailing list.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum
as a password protected csv file.
They are just text files, after all. So I guess you have something
else. How or what did the presumed pasword protection?
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
x description
1 0.4 my first measurement
2 1.6 Normal 5 object
3 0.4 Some measuremetn
4 0.7 A 4 long sample
R thinks that the description in row 2 ends in row 4 and you loose
data.
Try read.delim(..., quote=).
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für
is the way out.
And of course, I agree that jpeg should never be used for graphs.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
forget to close the device with dev.off() after the plot.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing list
- data.frame(x=c(0,0,0.9,0.73,0.78,1,0.3,0.32),
freq=c(14,3,12,15,2,15,2,8))
hist(with(dat, rep(x, times=freq)))
My example did not take special binning wishes into account but you
can easily customiye that with the breaks argument to hist.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für
won't have to do it again if I need to produce
another revision of the figure - all it takes is anoother run of my
script. And I can re-use good solutions in the future. Any manual
touch-ups have to be done manually every single time = not my idea of
efficiency.
cu
Philipp
--
Dr. Philipp
)
dat2 - as.data.frame(t(dat))
If the comlumns are of differnt types (e.g. some numeric, some
factors) I don't think you can do this at all, because columns of a
data.frame represent vectors, i.e. all value sin a column need to be
of the same type.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl
an alternative solution that doesn't require the
transposition.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
names
rownames(frame) = frame[,1]
# remove the old name column
frame - frame[,2:ncol(frame)]
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
tell us some more about your actual analysis,
more help can be provided.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org
is a named vector things are even easier:
names(v1) - v2
barplot(v1)
As I said, there are a bunch of other ways - e.g. using the lattice function
barchart() which works a bit differently.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
(equilogs = FALSE))
pln[,1] is just one column, so you are not plotting the values vs the
year but vs. their index. As xlim is set to the interval 1993-2011 you
simply don't see your data...
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
)
dev.off()
cat(\\only, slidenum, {\includegraphics{, plotfilename ,}}\n
,sep=)
@
If you want toget really fancy, you could wrap most of this in a
conveniance function...
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
Oops - have to comment my own answer:
echo=F=
For this to work it needs to be
echo=F, results=tex=
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354
Hi!
max.col does what you want. Example:
dat - data.frame(a=rnorm(20),b=rnorm(20),c=rnorm(20))
dat
a b c
1 1.17910304 -0.56951219 -0.2243664
2 -1.43840866 -0.99013855 -0.1613536
3 1.08515152 -0.77975274 0.3734530
4 -0.92154605 -0.20318367
not sure what else you could do if you
are not interested in a specific distance emasure etc.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
both matrices to files and compare them in an external
viewer (Excel, less, ...).
If I remember correctly, the object browser/data viewer of JGR allows
editing several objects at once.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
Assume that the 5th iteration (subject=5) leads to the error
message. How can I tell R to continue with the 6th iteration?
try or tryCatch are probably what you want.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org
/qbeta.c
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help
() to fit whatever function you find suitable
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
something that would like to be indexed starting with 1 -
just the other way round.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
some thing like this
jpeg('somefilename.jpg')
... plot... commands...
dev.off()
In addition to savePlot, which has already been recommended, you may
also want to look at dev.copy2eps and dev.copy2pdf.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
at the manual page confirms that Epi's version does not
have an origin option.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
...):
as.Date(36525, origin=1900-01-01')
[1] 2000-01-02
At first I thought the excess single quote maight be causingyour
problem, but it doesn't for me.
Maybe you need to upgrade R? Possibly it's an already fixed issue?
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte
not 100% sure, but if I understand your problem correctly,
loess() may be applicable.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
of levels of each factor in my data.frame seems to
make sense.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
par(mfrow=c(3,2))
The 6 graphs are coming out quite all right, but now I would like to
put a title on top of the page - i.e. something that is common for
all 6 graphs - how can I do that?
title(main=My title, outer=TRUE)
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für
there is no other choice. So I convert to a
e.g. 600dpi png. (Has this changed in recent versions of PowerPoint?)
But mind you: I don't do that in R, so I always have a vector format
master figure.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische
with.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing list
of essentually infinite
resolution for all practial purposes. So your problem is not with
the R-generated eps but somewhere downstream from that. Any
postprocessing, conversion or editing?
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
cases which, at least in my opinion, is the correct
value.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
like
write.table may be of interest to you. Load and save are meant for use
in R, only.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
, there is no magical link between the two. Of course, when you
load objects from a file this will overwrite any objects of the same
names (object names, not file names!) that happen to live in your
workspace before the load command.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für
it, but this does what
you asked for:
# plot data omitting the ylab
plot(1:10,1:10, ylab='')
# add the ylab myself using flushleft (adj=0.0)
mtext('foo', side=2, line=3, adj=0.0)
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing list
https
will most likely not fix this (unless a change in
the format of the .RData files had occurred - but to my knowledge no
such thing has happened, recently.)
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum
on the complexity of the task.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
(2,1,5))
Did I misinterpret the sentence from the help page or is this a bug?
Yes - I know that his works fine:
stripplot(yield~year|site, barley, layout=c(2,1))
Just curious...
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
to get but if you explicitly set the limits that's
what R is going to use.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
, where I would like to have them.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
is solved.
Nevertheless, take a look at the version in R-devel (pre-2.12.0)
which give you more options.
Thanks for that information - I will definitely have a look at that.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
in the manual page that explains this
behaviour or says how to change it. So I'm asking the community: Is
there a way to fix this or do I have to code around it?
This is in R 2.11.1 (Linux), by the way.
Thanks in advance
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
It works fine.
Could you explain to me why it did not worked for read.table?
Because of what Gavin already explaied in his reply: read.table
returns a data.frame and hist needs a vector.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische
in *.eps; do convert -density 600 $file `basename $file .eps`.png; done
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing list
https
and it is not a very good idea to
use function names as variable names.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
concept that I am
missing here...
I'm confused - what result where you expecting? p.adjust will need to
know the number of test you are trying to adjust for - either by
giving explicitly giving the number or by handing a vector of p-values
to the function.
cu
Philipp
--
Dr. Philipp
into a database and process it in
pieces. Or use sqldf or bigmemory or something like that.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http
.
[...]
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Maximus-von-Imhof-Forum 3
85354 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r
the example section is useful for a jumpstart.
Here is an example of computing row means:
apply(foo, 1, mean)
Instead of 'mean' you can insert whatever function you'd like to apply.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
feasible?
From what I gather fomr ?memory.limit it does not tell you how much
memory it currently available. So my guess is that you have som rather
large objects in your workspace already and thus there is not enough
space left for you vectors.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl
(determinant) :
no applicable method for determinant
det(as.matrix(foo))
[1] -0.1092489
So probably your problem is somewhere else.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Freising
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing list
https
(mrfow=c(2,1))
plot(d$a, d$b, col=red)
plot(d$a, d$c, col=red, pch=|)
As both plots use the same data for X you are set. If you need to
force two datasets with different x-ranges into the same range, you
can use the xlim parameter to define the desired range.
cu
Philipp
--
Dr. Philipp
the graph is a concern, you may also want to look at the
grid() command.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http://webclu.bio.wzw.tum.de/~pagel/
__
R-help@r-project.org mailing list
https
foo$B - as.numeric(sub(',', '', as.character(foo$B)))
foo
AB
1 1 12300.0
2 2 256001.0
3 3900.1
4 4 80.0
str(foo)
'data.frame': 4 obs. of 2 variables:
$ A: int 1 2 3 4
$ B: num 12300 256001 900 80
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte
in previous postings, as the data fram eused
there has 9 rows, not 14. Without the details (code) on what you did
its all guesswork. Perhaps you are mixing two data.frames of differnt
shape or ...
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität
On Fri, Sep 18, 2009 at 12:59:16PM +0200, Philipp Pagel wrote:
foo$B - as.numeric(sub(',', '', as.character(foo$B)))
Thinking about it some more, you should use gsub instead of sub here.
Otherwise only the first occurrence of the thousands separator will be
removed.
cu
Philipp
- sorry you are right. There is a Problem with inserting the
object. Try this instead:
dt$date - as.Date(dt$datetime, %d%b%Y)
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising
... are there good ways to generate tables in R?
Have a look at Sweave and xtable - I think that's what you want.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
http
On Thu, Sep 17, 2009 at 10:08:57AM +0200, Philipp Pagel wrote:
is it possible to convert latex code to pdf in R (like a
latex-program would do it)?
Is there a package that comes with this capabilities?
My problem is that I want to generate tables automatically -
and I can't use
the impression that your data file uses a mixture of
separation characters and fixed-width formatting. Maybe I
misinterpreted your example. Have a look at read.fwf() an if that does
not solve your problem maybe explain the Structure and expected result
a little further.
cu
Philipp
--
Dr. Philipp
hard to
parse the data - and it may even be impossible depending on what
characters are allowed in your free-text fields.
cu
Philipp
--
Dr. Philipp Pagel
Lehrstuhl für Genomorientierte Bioinformatik
Technische Universität München
Wissenschaftszentrum Weihenstephan
85350 Freising, Germany
1 - 100 of 212 matches
Mail list logo