//people.ee.ethz.ch/~oetiker/lshort/lshort.pdf
It really was a not too short intro. I'll have a look at it.
Yes definitly not too short. But it states in LaTeX in 133 min ...
Seems to be for Linux only. My server is Windows, even if I have the
rest of the components.
Hm at the projects'
Dear list,
Previously I posted these question to R-SIG-GUI. Perhaps here is a better place.
1. Is there some option for mantaining the tktoplevel window always on top?
2. Is there some option to eliminate the border icons maximize, minimize and
close of a tktoplevel window?
3. Is possible to
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
Right. the grid call shouldn't be necessary, axis.line controls the
panel borders. And tck can be a vector, to get rid of the ugly bumps
on top:
xyplot(x ~ x | g, type = l, lwd = lwd,
scales = list(tck = c(-1, 0)),
Hi all,
We've just had a paper accepted for publication in Nature. We used R for
95% of our analyses (one of my co-authors sneaked in some GenStat when I
wasn't looking.). The preprint is available from the Nature web site, in
the open peer-review trial section. I searched Nature for previous
Dear List,
after updating the exactRanksumTests package I receive a warning that
the package is not developed any further and that one should consider
the coin package.
I don't find the signed rank test in the coin package, only the Wilcoxon
Mann Whitney U-Test. I only found a signed rank test
Piet Bell skreiv:
The publisher of our paper has requested:
1. all tick marks should point inwards instead of outwards.
Point him to William S. Cleveland’s excellent book /The Elements of Graphing
Data/, where Cleveland strongly recommends that tick marks should point
*outwards* ‘because
Dear all,
I would like to thank everybody who replied for their useful
suggestions. Maybe, I am going through the book statistics to teach
(fresh) myself.
Wish you have a nice weekend.
Regards, Muhammad Subianto
On this day 24/08/2006 18:59, Muhammad Subianto wrote:
Dear all,
I apologize if
Dear Tom,
To add a few things to explore:
- I'd definitely go with LaTeX. Depending on how much formatting control you
want, though, and if your coworkers are reluctant to jump into LaTeX, you
might start with reStructuredText (http://docutils.sourceforge.net/rst.html)
or text2tags
On Fri, 25 Aug 2006, Stefan Grosse wrote:
Dear List,
Stefan,
after updating the exactRanksumTests package I receive a warning that
the package is not developed any further and that one should consider
the coin package.
I don't find the signed rank test in the coin package, only the
Simon Blomberg wrote:
Hi all,
We've just had a paper accepted for publication in Nature. We used R for
95% of our analyses (one of my co-authors sneaked in some GenStat when I
wasn't looking.). The preprint is available from the Nature web site, in
the open peer-review trial section. I
Andrej Kastrin [EMAIL PROTECTED] writes:
Simon Blomberg wrote:
Hi all,
We've just had a paper accepted for publication in Nature. We used R for
95% of our analyses (one of my co-authors sneaked in some GenStat when I
wasn't looking.). The preprint is available from the Nature web
Thomas Lumley wrote:
On Thu, 24 Aug 2006, Jason Liao wrote:
I recently coded a recursion algorithm in R and ir ran a few days
without returning any result. So I decided to try a simple case of
computing binomial coefficient using recusrive relationship
choose(n,k) = choose(n-1,
--- Gabor Grothendieck [EMAIL PROTECTED]
wrote:
Use cbind to create a two column matrix, mat,
and multiply that by the appropriate inflation
factors.
Then use rowsum to sum the rows according to the
id grouping factor.
inf.fac - list(year1 = 1, year2 = 5, year3 = 10)
mat - cbind(s1 =
Hi everyone,
I have what may appear to be a newbie question, but I have looked
everywhere I can think to look and I cannot find an answer. On page 35
of An Introduction to R the following command appears:
plot(ecdf(eruptions), do.points=FALSE, verticals=TRUE). What is the
do.points argument?
Take a look at ?plot.stepfun.
'ecdf()' returns an object of class ecdf inheriting from class stepfun and
'plot.ecdf()' calls 'plot.stepfun'.
-roger
Catherine Carter wrote:
Hi everyone,
I have what may appear to be a newbie question, but I have looked
everywhere I can think to look and I
Thank you, Thomas and Joerg! Joerg's example is extremely useful. In
fact I had wanted to compare R with other interpreting language myself.
I heard that Ruby is not fast among scripting languages. I thus believe
that there is room to improve R's recursive function.
Jason
--- Joerg van den
On Fri, 2006-08-25 at 09:08 -0400, Catherine Carter wrote:
Hi everyone,
I have what may appear to be a newbie question, but I have looked
everywhere I can think to look and I cannot find an answer. On page 35
of An Introduction to R the following command appears:
plot(ecdf(eruptions),
Thank you! That is exactly what I needed.
Roger D. Peng wrote:
Take a look at ?plot.stepfun.
'ecdf()' returns an object of class ecdf inheriting from class
stepfun and 'plot.ecdf()' calls 'plot.stepfun'.
-roger
Catherine Carter wrote:
Hi everyone,
I have what may appear to be a
Hi R-users,
I am trying to calculate the correlation between 3 vectors of numbers.
Is there any other solution then passing by kendall.w function which
need ranks in input?
Thanks,
Bianca Dimulescu
__
R-help@stat.math.ethz.ch mailing list
i tried fib(30):
R: 8.99s
Python (interpreted): 0.96s
Python (interpreted but using psyco library): 0.03s
C++: 0.015s
cheers
Damien
--- On Fri 08/25, Joerg van den Hoff [EMAIL PROTECTED] wrote:
From: Joerg van den Hoff [mailto: [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Cc: [EMAIL
Dear Rense,
(This question was posted a few days ago when I wasn't reading my email.)
So-called effect decompositions are simple functions of the structural
coefficients of the model, which in a model fit by sem() are contained in
the $A component of the returned object. (See ?sem.) One
On Fri, 25 Aug 2006, Joerg van den Hoff wrote:
maybe someone's interested:
I made the same observation of seemingly very slow recursion recently: just
for fun I used the (in)famously inefficient
fib - function(n = 1) {
if (n 2)
fn - 1
else
fn - fib(n - 1) + fib(n - 2)
Dear R-helpers,
using snow on socket connections, I ran into the following error
cl - makeSOCKcluster(hosts)
Error in socketConnection(port = port, server = TRUE,
blocking = TRUE : all connections are in use
with showConnections(all=T) showing 50 open connections.
As
biglm fits linear and generalized linear models to large data sets, using
bounded memory.
What's New: generalized linear models.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
Hi All,
I want to plot y~ x under the condition of variable a and b. Followed is the
dataset:
plotid lndenlnvol source
369 9037.0 10.419002 -4.101039226 S
370 9037.0 9.840548 -2.432385723 S
371 9037.0 8.973351 -1.374842169 S
372 9037.0 8.242756
It's the |source in your formula that tells lattice to separate them.
If you drop that, you'll get all points without S and P distinguished at
all. If you add a groups argument, you should get them presented with
different colors/symbols/etc. depending on your trellis settings (warning:
untested
I apologize if this is redundant. I've been Googling, searching the
archive and reading the help all morning and I am not getting closer
to my goal.
I have a series of data( xi, yi). It is not evenly sampled and it is
messy (meaning that there is a lot of scatter in the data). I want to
On Fri, 25 Aug 2006, Marc Kirchner wrote:
As - for administrative reasons - I would prefer to use snow's
SOCK capabilities (instead of MPI and the like), I wonder if there is a
way to increase the number of simultaneous open connections allowed in
an R session (~100 would be sufficient).
On Fri, 25 Aug 2006, Marc Kirchner wrote:
Dear R-helpers,
using snow on socket connections, I ran into the following error
cl - makeSOCKcluster(hosts)
Error in socketConnection(port = port, server = TRUE,
blocking = TRUE : all connections are in use
with
Dear List,
I try to make a xyplot with different colors and symbols, I came this far:
library(DAAG)
xyplot(csoa~it|sex*agegp,data=tinting,groups=target,pch=list(1,2),auto.key=list(space
= right))
this produces a plot with different colors and symbols but unfortunately
the legend does only
On Fri, 25 Aug 2006, [EMAIL PROTECTED] wrote:
Hi,
I have a casting/formatting question on hist.POSIXt:
The histogram plot from POSIXct works perfect (with help of Prof. Ripley
-thanks!).
When processing the hist(plot=FALSE) output and then plotting the
results over the x-axis (bins)
yep, it works!
impressingly elegant. R seems worth to dig into the details. I'll try harder.
Thank You very much
Peter
-Ursprüngliche Nachricht-
Von: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Gesendet: Freitag, 25. August 2006 17:53
An: Marx Peter
Cc: r-help@stat.math.ethz.ch
On Fri, 25 Aug 2006 17:51:24 +0200,
Stefan Grosse [EMAIL PROTECTED] wrote:
Dear List, I try to make a xyplot with different colors and symbols, I
came this far:
library(DAAG)
xyplot(csoa~it|sex*agegp,data=tinting,groups=target,pch=list(1,2),auto.key=list(space
= right))
this produces a
Adam D. I. Kramer wrote:
Hello,
I have an HLM-like situation: a data frame with rows representing
individuals, and columns representing some statistics for the individual
(e.g., a mean and variance for that individual, plus the number of
observations for that person).
I'm
Thanks Matthew. This does help. But from the discription of the function and
examples, we are not able to use groups argument and the trellis settings you
mentioned.
Jen.
- Original Message -
From: Wiener, Matthew, [EMAIL PROTECTED]
Sent: 2006-08-25, 09:46:44
To: jennystadt;
Hi
and if speed is an issue and object is numeric matrix something like
function(x) which(colSums(abs(diff(x)))==0)
is a little bit quicker
Cheers
Petr
On 25 Aug 2006 at 13:39, [EMAIL PROTECTED] wrote:
Date sent: Fri, 25 Aug 2006 13:39:48 +1000
From: [EMAIL
I didn't see a response to this, and am not an expert, but wanted to
point you in the right direction.
You should ask this question on the bioconductor mailing list.
The site is very much alive at this end (though I'm very close to it,
physically), so you'll need to provide some more information
R trellis.par.set(superpose.symbol=list(pch=c(1, 2)))
Warning message:
Note: The default device has been opened to honour attempt to modify trellis
settings in: trellis.par.set(superpose.symbol = list(pch = c(1, 2)))
R
On 8/25/06, Sebastian P. Luque [EMAIL PROTECTED] wrote:
On Fri, 25 Aug 2006 17:51:24 +0200,
Stefan Grosse [EMAIL PROTECTED] wrote:
Dear List, I try to make a xyplot with different colors and symbols, I
came this far:
library(DAAG)
Hello all,
I am using the function lm to do my weighted least
square regression.
model-lm(Y~X1+X2, weight=w)
What I am confused is the r.squared.
It does not seem that the r.squared for the weighted
case is an ordinary 1-RSS/TSS.
What is that precisely?
Is the r.squared measure comparable to
Michael Koppelman wrote:
I apologize if this is redundant. I've been Googling, searching the
archive and reading the help all morning and I am not getting closer
to my goal.
I have a series of data( xi, yi). It is not evenly sampled and it is
messy (meaning that there is a lot of scatter
On Fri, 25 Aug 2006, Charles wrote:
Hello all,
I am using the function lm to do my weighted least
square regression.
model-lm(Y~X1+X2, weight=w)
What I am confused is the r.squared.
What r.squared? There is no r.squared in that object, but it is
calculated by the summary method.
It
Hi,
Here is a vector and the result from the embed-command:
VECTOR=c(0,3,6,3,11,2,4,3,7,6,4,5,10,2,3,5,8)
embed(VECTOR, dimension=5)
[,1] [,2] [,3] [,4] [,5]
[1,] 113630
[2,]2 11363
[3,]42 1136
[4,]342 113
[5,]
--- Prof Brian Ripley [EMAIL PROTECTED] wrote:
savePlot is just an internal version of dev.copy, part of the support for
the menus on the windows() graphics device.
It is described in `An Introduction to R' (the most basic R manual).
the most basic R manual doesn't quite answer my
Hi,
I'm making an interface, where a Tcl/Tk window have few listbox widgets.
I need to select separate parameters from separate listboxes.
It is clear how to get cursor selection value, once you know which listbox
widget you clicked.
The problem is I can't figure out which one tcltk command to
Hello, Ted,
Thank you for your help!
So I can not contrast the mean yields between sections 1-8 and 9-11
under Trt but I can contrast mean yields for sections 1-3 and 6-11
because there exists significant interaction between two factors
(Trt:section4, Trt:section5). Could I use the commands
Dear List,
The following code produces a list, which is what I what:
set.seed(123)
tmpf - function() {
x - rpois(rpois(1,4),2)
}
n - 3
m - replicate(n,tmpf())
m
[[1]]
[1] 3 2 4
[[2]]
[1] 0 2 4 2 2 5 2
[[3]]
[1] 2 0 4 1 0
Now I need something that would to extract iteratively (or as many
I would like a tick mark for each month; for example,
xyplot(runif(365)~I(as.Date(1999-01-01) + 1:365),
scales=list(x=list(format=%b %Y,tick.number=12)))
I know I could make x numeric and use 'at' and 'labels', but I was
wondering if there is a more direct route I'm missing. (In
Dear list,
I know it must be obvious and I did my homework. (In fact I've
RSiteSearched with keyword remove AND library but got timed
out.(why?))
How do I unload a library? I don't mean getting ride of it permanently
but just to unload it for the time being.
A related problem : I have some
Thank you. Yes, I do feel that I am under-qualified to even ask
questions of y'all. Plus I'm an astronomer, which doesn't help! ;)
I'll try again.
I have two columns of data, the first column (x) is a distance (or
length or separation) and the second column (y) is a flux (or number
of
On 8/24/2006 4:34 PM, Gaspard Lequeux wrote:
Hej,
On Thu, 24 Aug 2006, Duncan Murdoch wrote:
On 8/24/2006 11:19 AM, Gaspard Lequeux wrote:
On Wed, 23 Aug 2006, Duncan Murdoch wrote:
On 8/23/2006 5:15 PM, Gaspard Lequeux wrote:
When exporting a image from rgl, the following error is
II am translating some gauss code into R, and gauss has a matrix
product function called the horizontal direct product (*~), which is
some sort of variant on the Kronecker product.
For example if x is 2x2 and y is 2x2
the horizontal direct product, z, of x and y is defined (in the Gauss
I don't know if mount -l shows something interesting below.
We're wondering if there is an issue with the gcc compiler, which has been
updated since we last updated R. Do you see anything that could be wrong with
the info below?
We plan on updating to R 2.3.1, hoping that will clear things
see ?detach
Horace Tso [EMAIL PROTECTED] wrote:
Dear list,
I know it must be obvious and I did my homework. (In fact I've
RSiteSearched with keyword remove AND library but got timed
out.(why?))
How do I unload a library? I don't mean getting ride of it permanently
but just to unload
Dear Janet,
Performing a traceback after the error gives a hint:
tmp.pcc-polychor(tmp.mat, ML=T, std.err=T)
traceback()
8: stop(at least one element of , sQuote(lower), is larger than ,
sQuote(upper))
7: checkmvArgs(lower = lower, upper = upper, mean = mean, corr = corr,
sigma
try detach(package:zoo)
Sachin
Horace Tso [EMAIL PROTECTED] wrote:
Sachin,
I did try that, ex
detach(zoo)
Error in detach(zoo) : invalid name
detach(zoo)
Error in detach(zoo) : invalid name
But zoo has been loaded,
sessionInfo()
Version 2.3.1 (2006-06-01)
i386-pc-mingw32
Sachin,
I did try that, ex
detach(zoo)
Error in detach(zoo) : invalid name
detach(zoo)
Error in detach(zoo) : invalid name
But zoo has been loaded,
sessionInfo()
Version 2.3.1 (2006-06-01)
i386-pc-mingw32
attached base packages:
[1] methods datasets stats tcltk utils
On Fri, 25 Aug 2006, Horace Tso wrote:
Dear list,
I know it must be obvious and I did my homework. (In fact I've
RSiteSearched with keyword remove AND library but got timed
out.(why?))
Probably because the site is offline.
How do I unload a library? I don't mean getting ride of it
Aah, that works. The missing package:...
H.
Sachin J [EMAIL PROTECTED] 8/25/2006 1:16 PM
try detach(package:zoo)
Sachin
Horace Tso [EMAIL PROTECTED] wrote:
Sachin,
I did try that, ex
detach(zoo)
Error in detach(zoo) : invalid name
detach(zoo)
Error in detach(zoo) : invalid name
maybe something like:
%*~% - function(x, y){
n - nrow(x)
out - matrix(0, n, ncol(x) * ncol(y))
for(i in 1:n)
out[i, ] - c(y[i, ] %o% x[i, ])
out
}
x - matrix(1:4, 2, 2, TRUE)
y - matrix(5:8, 2, 2, TRUE)
x %*~% y
I hope it helps.
Best,
Dimitris
Dimitris
On Fri, 25 Aug 2006, Prof Brian Ripley wrote:
On Fri, 25 Aug 2006, Marc Kirchner wrote:
Dear R-helpers,
using snow on socket connections, I ran into the following error
cl - makeSOCKcluster(hosts)
Error in socketConnection(port = port, server = TRUE,
blocking = TRUE : all
Hi,
There are two ways that I know of:
detach(package:zoo)
or use detach() with the number of that package's
place in the search list.
library(zoo)
search()
[1] .GlobalEnvpackage:zoo package:methods
[4] package:graphics package:grDevices package:utils
[7] package:datasets
hi michael : ( i stupidly didn't send the initial email to the whole list so
they will have to read from the bottom ), it's clearer now but fitting a
gaussian still may not be the right thing to do in this case. someone else
can answer better but it sounds to me like you want to do some kind
Try this:
xyplot(runif(365)~I(as.Date(1999-01-01) + 1:365),
scales=list(x=list(at=seq(as.Date('1999-1-1'), length=12, by='1
month'), labels=NULL)))
On 8/25/06, Benjamin Tyner [EMAIL PROTECTED] wrote:
I would like a tick mark for each month; for example,
try this:
set.seed(123)
tmpf - function() {
+ x - rpois(rpois(1,4),2)
+ }
n - 3
m - replicate(n,tmpf())
m
[[1]]
[1] 3 2 4
[[2]]
[1] 0 2 4 2 2 5 2
[[3]]
[1] 2 0 4 1 0
lapply(m, function(x)x[x2])
[[1]]
[1] 3 4
[[2]]
[1] 4 5
[[3]]
[1] 4
On 8/25/06, xpRt.wannabe [EMAIL PROTECTED] wrote:
event.nab.2 is 0/1 and I dichotomized va to get va.2 to see if I could
get geeglm to work. glm has no problem with the data but geeglm chokes.
Each subject (patient.id) has at most 2 observations and more than 3/4
of the subjects have 2 observations. I have even worse problems trying
to use
On Fri, 25 Aug 2006, Luke Keele wrote:
II am translating some gauss code into R, and gauss has a matrix
product function called the horizontal direct product (*~), which is
some sort of variant on the Kronecker product.
For example if x is 2x2 and y is 2x2
the horizontal direct product, z,
AAh. Then my hypothesis has been rejected. Oh well!
Cheers,
Simon.
Simon,
Congratulations!
It used to be that
R Ihaka and R Gentleman
R: A Language for Data Analysis and Graphics
Journal of Computational and Graphical Statistics, 1996
was used to cite R.
I see a
Thank you for your answers yesterday. I now have another question!
Suppose that instead of creating a formula for a regression model I
simply wanted to add the variables. I believe I cannot use the
as.formula anymore. Again I tried to use expression to no avail. I get
an expression but I can't
Try this and note we must pay attention to which environment
the expression is evaluated in.
fit.sum - function(x, env = parent.frame())
eval(parse(text = paste(x, collapse = +)), env = env)
# test
x1 - x2 - 1
fit.sum(c(x1,x2))
On 8/25/06, Maria Montez [EMAIL PROTECTED] wrote:
Thank you
You likely want the answer that Sarah has already given but in
addition you might also want to look at the thread below,
the point being that detaching a package still leaves portions:
https://www.stat.math.ethz.ch/pipermail/r-help/2006-July/109056.html
On 8/25/06, Sarah Goslee [EMAIL
Dear all,
I have the following problem:
- I have written a program in R which runs out of physical memory
on my MacBook Pro with 1.5 GB RAM
- The memory usage is much higher then I would expect from the
actual data in my global environment
- When I save the workspace, quit R, restart
Jim and Patrick,
Both of you made the same suggestion, which works great!
A follow-up question: Suppose I change the condition 'x2' in 'lapply'
to 'x4', as follows:
set.seed(123)
tmpf - function() {
x - rpois(rpois(1,4),2)
}
n - 3
m - replicate(n,tmpf())
m
sub.m - lapply(m, function(x)x[x4]) #
Yet another question:
Let's say I do the following:
set.seed(123)
tmpf - function(){
x - rpois(rpois(1,4),2)
}
n - 3
m - replicate(n, tmpf())
m
sub.m - lapply(m, function(x)x[x2])
'sub.m' gives me:
[[1]]
[1] 3 4
[[2]]
[1] 4 5
[[3]]
[1] 4
The question is: What do I need to do such that I
74 matches
Mail list logo