Dear Francisco,
thanks for your solution. It turns out that it's best for me to use
\setkeys{Gin}{width=0.15\textwidth}
directly before I call the plot - that seems to work just fine.
Andrwe
On Thu, Sep 08, 2005 at 05:44:59AM +, Francisco J. Zagmutt wrote:
> Others may propose more elegan
Francisco J. Zagmutt wrote:
> I don't have much experience in the subject but it seems that library(akima)
> should be useful for your problem. Try library(help="akima") to see a list
> of the functions available in the library.
>
> I hope this helps
>
> Francisco
Yes, function aspline() of pa
Others may propose more elegant solutions but, in windows one quick an dirty
option would be to change the argument 'pin' and 'fin' within par to get an
image of exactly 1 inch (2.54 cm) i.e.
y <- c(40, 46, 39, 44, 23, 36, 70, 39, 30, 73, 53, 74)
x <- c(6, 4, 3, 6, 1, 5, 6, 2, 1, 8, 4, 6)
par(pi
I don't have much experience in the subject but it seems that library(akima)
should be useful for your problem. Try library(help="akima") to see a list
of the functions available in the library.
I hope this helps
Francisco
>From: Gabor Grothendieck <[EMAIL PROTECTED]>
>Reply-To: [EMAIL PROTEC
On 7 Sep 2005, [EMAIL PROTECTED] wrote:
> # yum list installed mysql
> Installed Packages
> mysql.i3864.1.11-2 installed
> mysql.x86_64 4.1.11-2 installed
I would have thought that you need to have a mysql-dev.x86_64 rpm
package installed in order to
I have a batch of data in each line of data contains three values,
calcium score, age, and sex. I would like to predict calcium scores as a
function of age and sex, i.e. calcium=f(age,sex). Unfortunately the
calcium scorers have a very "ugly distribution". There are multiple
zeros, and multiple val
Hi everyone,
I was wondering if anyone has any code they could share for creating
thumbnail plots in Sweave. For example, I'd like a plot like the
following:
y <- c(40, 46, 39, 44, 23, 36, 70, 39, 30, 73, 53, 74)
x <- c(6, 4, 3, 6, 1, 5, 6, 2, 1, 8, 4, 6)
opar <- par(mar=c(3,3,0,0))
plot(x, y, x
Dear listers,
I have a piece of code which performs an ANOVA type of analysis on 2D GC
data. The code is shown below:
# ANOVA 2D GC analysis
# maxc <- number of samples
# nreps <- number of samples
maxc <- 2
nreps <- 4
sscl <- NULL
cmean <- NULL
#
# Initial stat. variable
#
d
On 9/7/05, David James <[EMAIL PROTECTED]> wrote:
> The purpose of this email is to ask for pre-built procedures or
> techniques for smoothing and interpolating missing time series data.
>
> I've made some headway on my problem in my spare time. I started
> with an irregular time series with lots
Salang Pan wrote:
> hi，
>
> Is it possible to draw a string text in a rectangle according the width of
> this rectangle? that is, the fontsize of this string text can be adjusted
> according the width of the rectangle.
> How to set the cex parameter in text function?
>
> text (x, y = NU
Hello,
I cant seen to find the equivalent function in version 2.1.1. This always
served as a nice preliminary tool for looking visually at data.
Thanks,
Matt
Matthew MacManes
PhD Student
Museum of Vertebrate Zoology
Department of Integrative Biology
UC-Berkeley
Berkeley, CA. 94720
Office
Dear useRs,
I'm having a hard time installing RMySQL on a FC4 x86_64 box (R 2.1.0
and MySQL 4.1.11-2 installed through yum). After an initial
configuration error ("could not find the MySQL installation include
and/or library directories") I managed to install RMySQL with
# export PKG_LIBS="-L -
Actually, I think period = 365 * 24 = 8760 is really what I need.
That crashes arima as well.
> I am using period = 365, which makes sense me to me. Is this wrong?
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing li
I cannot state this with the certainty that others might, but the Rd format is
a text format. If you want to produce something else then you need to choose an
alternative method. For instance, 1.4 of "Writing R Extensions" notes that
"Documents in 'inst/doc' can be in arbitrary format, however w
The following command crashes my Mac OS X version of R:
(I'm running R on a PowerMac G5, with 1 GB of RAM and dual processors.)
> arima.0 <- arima(w3.ts,order=c(1,0,0),seasonal=list(order=c
(1,0,0),period=365))
-David
Here is some background:
w3.ts is hourly temperature data with about 20%
The purpose of this email is to ask for pre-built procedures or
techniques for smoothing and interpolating missing time series data.
I've made some headway on my problem in my spare time. I started
with an irregular time series with lots of missing data. It even had
duplicated data. Thank
Hi,
In manual for optim() function, there are three control parameters for
Nelder-Mead algorithm: alpha (reflection),beta(contraction) and
gamma(expansion), but in the original paper, there is another parameter
delta which controls shrinkage, how can I set this shrink parameter? Or is
beta act
I looked through the "Writing R Extensions" pdf, and I don't see how graphics
can be input in help files. For example, if I had a .eps plot that I wanted
to include in a help file, what would the syntax be to include it in an R
help file?
If there is graphics support in help files, which format ar
On Wed, 7 Sep 2005, Clark Allan wrote:
>
> $MLE$message
> [1] "ERROR: ABNORMAL_TERMINATION_IN_LNSRCH"
>
>
> WHAT DOES THIS ERROR MESSAGE MEAN???
>
Looking at the code in optim() a little, it looks as though this error
comes when the optimiser tries to do a line search in the steepest descent
di
At 13:15 5/9/2005, you wrote:
>Hi,
>
> I'm using mda library to solve a discriminant
>analysis. I get results, but the thing is that I want
>to use Fisher's method to obtain the classification
>functions and I'm lost in what I should do: libraries
>to use, ... Can anybody give me a clue??
Hi Ca
Sebastien Durand wrote:
> Dear all,
>
> I am running
> R : Copyright 2005, The R Foundation for Statistical Computing
> Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0
> Under Mac os X, a french version!
>
> I am preparing a package and I got the following issue
>
> I am trying to read dates tha
Hi
Tim Churches wrote:
> Paul Murrell wrote:
>
>> Hi
>>
>> Deepayan Sarkar wrote:
>> > On 9/7/05, Tim Churches <[EMAIL PROTECTED]> wrote:
>> >
>> >> Version 2.1.1 Platforms: all
>> >>
>> >> What is the trellis parameter (or is there a trellis parameter) to
>> >> set the leading (the gap between
Luis Pineda wrote:
> I sent this email before, but I got a r-help-bounce message and I don know
> if it got to the m-list. Sorry if you had already seen it.
>
> I'm using the areg.boot function to do an ace regression. So far I've been
> able to do some simple running tests to fit a model with s
Paul Murrell wrote:
> Hi
>
> Deepayan Sarkar wrote:
> > On 9/7/05, Tim Churches <[EMAIL PROTECTED]> wrote:
> >
> >> Version 2.1.1 Platforms: all
> >>
> >> What is the trellis parameter (or is there a trellis parameter) to
> >> set the leading (the gap between lines) when long axis values
> >> labe
[EMAIL PROTECTED] writes:
> Has anyone written a matrix editor or data.entry() replacement using the Tk
> table widget? I've been playing around with the examples at
>
> http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/tktable.html
>
> and making some progress, but I'd rather not spend much
Hi
Deepayan Sarkar wrote:
> On 9/7/05, Tim Churches <[EMAIL PROTECTED]> wrote:
>
>> Version 2.1.1 Platforms: all
>>
>> What is the trellis parameter (or is there a trellis parameter) to
>> set the leading (the gap between lines) when long axis values
>> labels or panel header labels wrap o
Dear all,
I am running
R : Copyright 2005, The R Foundation for Statistical Computing
Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0
Under Mac os X, a french version!
I am preparing a package and I got the following issue
I am trying to read dates that are written in
english and have them reco
On 9/7/05, Tim Churches <[EMAIL PROTECTED]> wrote:
> Version 2.1.1
> Platforms: all
>
> What is the trellis parameter (or is there a trellis parameter) to set the
> leading (the gap between lines) when long axis values labels or panel header
> labels wrap over more than one line? By default, the
Has anyone written a matrix editor or data.entry() replacement using the Tk
table widget? I've been playing around with the examples at
http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/tktable.html
and making some progress, but I'd rather not spend much time on this if
someone else has alrea
I sent this email before, but I got a r-help-bounce message and I don know
if it got to the m-list. Sorry if you had already seen it.
I'm using the areg.boot function to do an ace regression. So far I've been
able to do some simple running tests to fit a model with some input data,
predict the
Revised versions of maptools, a package for reading geographical data from
shapefiles, and sp, a package with classes and methods for spatial data
handling, have been released on CRAN. They are maptools release 0.5-1 and
sp release 0.8-1.
The maptools package now depends on sp (>= 0.8), so that us
On Wed, 7 Sep 2005, Douglas Bates wrote:
> On 9/7/05, Doran, Harold <[EMAIL PROTECTED]> wrote:
> > You could use <>= and the code inside the chunk will not
> > be evaluated. I suppose two other options could be to comment out the
> > bad code inside the code chunk or to use verbatim to make it loo
On 9/7/05, John Wilkinson (pipex) <[EMAIL PROTECTED]> wrote:
> Ronaldo,
>
> Further to my previous posting on your Glycogen nested aov model.
>
> Having read Douglas Bates' response and Reflected on his lmer analysis
> output of your aov nested model example as given.The Glycogen treatment has
Last week I posted a question concerning the mCall function, which is
used to create self-starting functions and is described in the book by
Pinheiro, J.C. and Bates, D.M. (Mixed-effects models in S and S-PLUS).
On page 345 one finds the following call:
xy<-sortedXyData(mCall[["x"]], LHS,data)
On 9/7/05, Doran, Harold <[EMAIL PROTECTED]> wrote:
> You could use <>= and the code inside the chunk will not be
> evaluated. I suppose two other options could be to comment out the bad code
> inside the code chunk or to use verbatim to make it look like a code chunk in
> your output.
>
>
>
You could use <>= and the code inside the chunk will not be
evaluated. I suppose two other options could be to comment out the bad code
inside the code chunk or to use verbatim to make it look like a code chunk in
your output.
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAI
FYI
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce 2-day "R/S-plus Fundamentals and Programming
Techniques" in New York City.
www.xlsolutions-corp.com/Rfund.htm
New York City September 22nd-23rd, 2005
Reserve your seat now at the early bi
Peter Dalgaard <[EMAIL PROTECTED]> writes:
> "Francisco J. Zagmutt" <[EMAIL PROTECTED]> writes:
>
> > Check some of the threads at RSiteSearch("Hotelling")
>
> Or use anova(lm(X~g), test="Hotelling"), where X is the matrix of
> responses and g is the grouping factor.
Oops, sorry. That is in r-d
"Francisco J. Zagmutt" <[EMAIL PROTECTED]> writes:
> Check some of the threads at RSiteSearch("Hotelling")
Or use anova(lm(X~g), test="Hotelling"), where X is the matrix of
responses and g is the grouping factor.
> Cheers
>
> Francisco
>
> >From: Bill Donner <[EMAIL PROTECTED]>
> >To: R-help@
Check some of the threads at RSiteSearch("Hotelling")
Cheers
Francisco
>From: Bill Donner <[EMAIL PROTECTED]>
>To: R-help@stat.math.ethz.ch
>Subject: [R] Hotelling Test
>Date: Wed, 7 Sep 2005 06:48:06 -0700 (PDT)
>
>Hello R-users,
>
>I've been looking for a function performing one and two sample
On 9/7/05, Thomas Petzoldt <[EMAIL PROTECTED]> wrote:
> Omar Lakkis schrieb:
> > Using system() is theer a way to make the R interpreter not wait for
> > the command to finish?
>
> system("cmd", wait=FALSE)
>
> see ?system in online help.
>
Coincidentially I recently posted on r-devel a note p
Hi all,
Is there an option in Sweave to avoid it from stopping on a code chunk with an
error? (I purposefully want to include code with an error in class notes.)
I suspect the answer is "no" and that I will be pointed to options("error").
That'd be fine, but which error parameter will just "do
Omar Lakkis schrieb:
> Using system() is theer a way to make the R interpreter not wait for
> the command to finish?
system("cmd", wait=FALSE)
see ?system in online help.
Thomas P.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/
Using system() is theer a way to make the R interpreter not wait for
the command to finish?
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.htm
Hi
what is wrong with matplot
sines <- outer(1:20, 1:4, function(x, y) sin(x / 20 * pi * y))
matplot(sines, pch = 1:4, type = "o", col = rainbow(ncol(sines)))
so you can use aditional parameters to exactly specify what type of
point and/or line and in what colour you will plot.
Or with plot/po
Le 07.09.2005 16:00, Stéphane Mattei a écrit :
>Thank you all for your answers.
>
>I eventually use the points command
>
>plot(MATRIX[,1])
>points(MATRIX[,2])
>points(MATRIX[,3])
>...
>
>with matplot I had numbers instead of points with type="p" and par(new=TRUE)
>makes complications
>with the ax
How can I pass parameters to an R script from the command line. And
how can I read them from within the script?
This is how I want to invoke the script:
R CMD BATCH r.in r.out
The script with read in the input values, process them and spit the
output to r.out.
_
Olshansky,Moshe wrote:
> What is the "classic" R function for solving a (possibly over
> determined) system of non-linear equations?
>
>
>
> Thank you!
>
>
>
> Moshe Olshansky
>
> e-mail: [EMAIL PROTECTED]
>
>
>
I'm not sure what your definition of '"classic"' is, but there are
se
Selon Gabor Grothendieck <[EMAIL PROTECTED]>:
> On 9/7/05, Stéphane Mattei <[EMAIL PROTECTED]> wrote:
> > Hello !
> >
> >
> > There is something quite simple I want to do with R but I found nowhere in
> > the help how to do
> it.
> > I just want to plot data which are in a matrix, every column be
Ronaldo,
Further to my previous posting on your Glycogen nested aov model.
Having read Douglas Bates' response and Reflected on his lmer analysis
output of your aov nested model example as given.The Glycogen treatment has
to be a Fixed Effect.If a 'treatment' isn't a Fixed Effect what is ? If
D
On 9/7/05, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Hello.
>
> This is my first post, so allow me to introduce myself.
>
> But first, I'd like to thank all the authors and contributors to the R
> software,
> as I think that it is truly a great and very useful package.
>
> I am the author
On Wed, 7 Sep 2005, Basile Chaix wrote:
> Dear all,
> I would have some questions on the coxph function for survival analysis,
> which I use with frailty terms.
>
> My model is:
> mdcox<-coxph(Surv(time,censor)~ gender + age + frailty(area, dist='gauss'),
> data)
> I have a very large proportion o
What is the "classic" R function for solving a (possibly over
determined) system of non-linear equations?
Thank you!
Moshe Olshansky
e-mail: [EMAIL PROTECTED]
The information transmitted is intended only for the person(s) or entity to
which it is addressed and may contain confident
some time ago I've written a function for the Hotelling test, maybe
you could give it a try:
hotel.test <- function(x, y = NULL, mu = 0) {
if(!is.numeric(x) || !is.matrix(x))
stop("'x' must be a numeric matrix")
n <- nrow(x)
p <- ncol(x)
xbar <- colMeans(x, na.rm = TRUE)
On Wed, 7 Sep 2005, Chris Buddenhagen wrote:
>
> 1) I have been really pleased with R as a means of doing and learning
> statistics. I work in a Spanish speaking country- and I wanted to pass on
> the benefits of R to my Spanish speaking colleagues. There are a couple of
> introductions to R in Spa
On 9/7/05, Stéphane Mattei <[EMAIL PROTECTED]> wrote:
> Hello !
>
>
> There is something quite simple I want to do with R but I found nowhere in
> the help how to do it.
> I just want to plot data which are in a matrix, every column being a data set
> and having the same
> x-axis (just an index
Dear all,
I would have some questions on the coxph function for survival analysis,
which I use with frailty terms.
My model is:
mdcox<-coxph(Surv(time,censor)~ gender + age + frailty(area, dist='gauss'),
data)
I have a very large proportion of censored observations.
- If I understand correctly,
Thank you all for your answers.
I eventually use the points command
plot(MATRIX[,1])
points(MATRIX[,2])
points(MATRIX[,3])
...
with matplot I had numbers instead of points with type="p" and par(new=TRUE)
makes complications
with the axis.
__
R-help@s
Hello R-users,
I've been looking for a function performing one and two sample Hotelling
test for testing equality of mean vectors. Has anyone implemented such a
function in R?
thanks a lot,
Bill
==
Bill Donner
Statistician
__
R-help@stat
Have a look at ?matplot
Stéphane Mattei wrote:
> Hello !
>
>
> There is something quite simple I want to do with R but I found nowhere in
> the help how to do it.
> I just want to plot data which are in a matrix, every column being a data set
> and having the same
> x-axis (just an index).
>
Le 07.09.2005 15:28, Stéphane Mattei a écrit :
>Hello !
>
>
>There is something quite simple I want to do with R but I found nowhere in the
>help how to do it.
>I just want to plot data which are in a matrix, every column being a data set
>and having the same
>x-axis (just an index).
>
>So for e
Chris Buddenhagen, Botany Department, Charles Darwin Research Station, Santa
Cruz,Galapagos. Mail: Charles Darwin Foundation, Casilla 17-01-3891 Avenida
6 de Diciembre N36-109 y Pasaje California Quito, ECUADOR
Dear all
1) I have been really pleased with R as a means of doing and learning
stat
On 9/7/2005 9:28 AM, Stéphane Mattei wrote:
> Hello !
>
>
> There is something quite simple I want to do with R but I found nowhere in
> the help how to do it.
> I just want to plot data which are in a matrix, every column being a data set
> and having the same
> x-axis (just an index).
>
> So
Hello !
There is something quite simple I want to do with R but I found nowhere in the
help how to do it.
I just want to plot data which are in a matrix, every column being a data set
and having the same
x-axis (just an index).
So for example if I have a 50 x 6 matrix I want 6 set of points on
Patrick Connolly wrote:
> I've never had this problem before and can't see what could be
> different from other times I've used keys with lattice.
>
>
>
> It appears that auto.key is being taken as TRUE when I specify a key
> list. The list I specify seems to be ignored.
>
> Where can I plac
funny optim message:
$MLE
$MLE$par
[1] -0.09554688 1.13100488 0.06651340
$MLE$value
[1] 48.93381
$MLE$counts
function gradient
100 100
$MLE$convergence
[1] 52
$MLE$message
[1] "ERROR: ABNORMAL_TERMINATION_IN_LNSRCH"
WHAT DOES THIS ERROR MESSAGE MEAN???
hope some one can help.
Not an R-response, but see this reference:
Dolan CV, van der Maas HLJ, Molenaar PCM
A framework for ML estimation of parameters of (mixtures of) common reaction
time distributions given optional truncation or censoring
BEHAVIOR RESEARCH METHODS INSTRUMENTS & COMPUTERS 34 (3): 304-323 AUG 2002
o
Hello,
I want to fit a distribution to a dataset. Important is not the "overall"
fitting but the fitting in the tail (e.g. all observations > x or the n
highest values). Standard ML-estimation sometimes doesn't work here very
well. We see that especially when we have truncated datasets the algori
hi,
when I use bclust in R, bclust(dat,centers=5,minsize=3,base.centers=4)
dat has 25 rows, there is an error as following:
Error in knn1(object$allcenters, x, factor(1:nrow(object$allcenters))) :
train and class have different lengths
when I debug this function , I found the err
Ronaldo ,
It looks as though you have specified you model incorrectly.
In the Rats example ,the Treatment is the only fixed effect,Rat and Liver
are random effects
In aov testing for sig of 'Means' of Random Effects is pointless and that is
why 'p' values are not given.Further more the interacti
Murray Jorgensen wrote:
> > thanks <- read.fwf("C:\\Files\\Reading\\thankyou.txt", c(43,37))
[CCing R-help again: I have looked at Murray Jorgensen's data in the
meantime]
tord <- order(thanks$V2)
sink("C:\\thanks.txt")
thanks[tord,]
sink()
Works for me with R-2.1.1.
--> Please upgrade your
Hello.
This is my first post, so allow me to introduce myself.
But first, I'd like to thank all the authors and contributors to the R software,
as I think that it is truly a great and very useful package.
I am the author of moodss, a GPL modular monitoring application
(http://moodss.sourceforge.
Uwe Ligges wrote:
> I guess there is just too much space or some special characters in your
> variables that cause problems when printing ...
> Hence you have to "debug" your data yourself.
>
> Uwe Ligges
However the problem persists when I don't try to print the fram "thanks"
to a file.
thank
Hi,
I read your answer and the message you pointed me
at, and you talked about the page 347 of the book MASS
3 in your posting as a place where the Fisher's method
was mentioned. The thing is that I don't have that
book, so I would like to ask you if you can give me
that information. If you don
Salang Pan wrote:
> hi，
>
> Is it possible to draw a string text in a rectangle according the width of
> this rectangle? that is, the fontsize of this string text can be adjusted
> according the width of the rectangle.
> How to set the cex parameter in text function?
>
> text (x, y = N
75 matches
Mail list logo
