thanx for the reply. i understood that the function found a maximum. i
was just a bit worried about the message. i assumed that it was an
ERROR message.
i see now that it is some sort of stopping rule. does this make sense?
/
allan
Douglas Bates wrote:
On 9/6/05, Clark Allan [EMAIL
Murray Jorgensen wrote:
[Using 2.0.1 under Windows XP]
There are a few pages on the internet that list equivalents of
thank you in many languages. I downloaded one from a Google search
and I thought that it would be interesting and a good R exercise to
sort the file into the order of the
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 = NULL,
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
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.
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
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 version
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
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
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
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
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.
/
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 place a
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 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 for
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
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 example if
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).
So for
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
__
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.
__
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 Spanish
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, 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 of
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 of moodss,
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
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 being a data
set
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
several options in
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 input values
The script with read in the input values, process them and spit the
output to r.out.
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 axis.
So,
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
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!
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
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
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 pointing out
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 Hotelling
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@stat.math.ethz.ch
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-devel only. For R
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
On Wed, 7 Sep 2005, Douglas Bates wrote:
On 9/7/05, Doran, Harold [EMAIL PROTECTED] wrote:
You could use eval=FALSE= 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
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
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
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, there is a
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
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 over more than
[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 time
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 some
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 lines) when long axis values
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 that are
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 are
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.
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 we
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
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
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
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 = NULL,
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 of
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
#
dftot
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,
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
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 get
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 PROTECTED]
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)
62 matches
Mail list logo