Justin Fay wrote:
I'm using Linux version 2.4.20-8smp (gcc version 3.2.2, RedHat 9). I
installed R-1.7.1-1.src.rpm following INSTALL. However, no X11 device is
found when I run R.
X11()
Error in X11() : X11 is not available
In the configuration, I found checking for X... no
How do I
Hello.
Is there an easy-to-use contour plot function analogous to scatterplot3d
that can draw handle a dataset of arbitrary (x,y,z) triplets? That is,
say x, y, and z are each measured quanties, and exhibit neither order
nor regularity.
I looked at the lattice package function contourplot but
# I need to generate some data. I'm modeling some time series that follow
a
# negative exponential decay (mostly). I have 20 samples that can easily
be fit with cubic splines.
# What I want to do is generate many thousands of similar samples using
the parameters from the splines
# For instance
On Wed, 13 Aug 2003, anne wrote:
As a newbie to R, I need to learn my way around (no previous experience
of S).What books, doc for R are recommended? I'm interested primarly in
non linear regression and process modelling (and have downloaded the R
documentation from the site).
So you have
Aurora Torrente wrote:
Is there any way to remove at the same time several variables that share
a suffix, for example (similar to * or ? in DOS) ?
To remove objects sharing a suffix use something like:
remove(list=objects(pattern=SUFFIX$))
To remove variables sharing a suffix from a
I wrote:
I found myself wishing for a function to rotate a vector.
Is there one? I know about ?lag, but help.search(rotate)
didn't find anything to the point.
Here I was regarding a vector as a _sequence_.
The (one-step) rotation of c(u,v,w,x,y,z) is c(v,w,x,y,z,u).
This
This might be easy (but I didn't find an answer in the archives).
I'm trying to make nice looking images (using image()).
To make them look nice (not jagged), it usually takes me at least 100x100 points. This
can be slow for frequent redraws.
Is there a smarter (less point intensive) way?
Hello everyone
(Im not sure if i already posted this problem)
I have grouped data and i used a Trellis plot to show the curve of a fitted
model for each group. Additionally i have a vector whose values are the
highest points the curve can reach (the curve rises exponetially). Each
value of the
Dear Dassy,
I'm not entirely sure of the specifics of what you want to
achieve, but for similar sounding things I typically make
use of the subscripts argument to panel functions. Check out
the following simple example.
x - rnorm(500)
y - 0.9*x + sqrt(1-0.9*0.9)*rnorm(500)
g -
On Tue, 2003-08-12 at 11:38, Paul, David A wrote:
Undoubtedly a simple question:
I've looked at order() and sort() in the help pages for
R1.7.1. It doesn't appear that these functions are immediately
suited to doing the same thing as
PROC SORT DATA = BLAH;
BY X Y Z;
RUN;
in
Alexander Singer wrote:
Hi,
There seems to be a problem with the power-operation of the form c(...,
-a, ...)^x.yyy
Four code examples for the same calculation, which perform differently:
-3^3.2
Result: [1] -33.63474
c(-3,6)^3.2
Result: [1] NaN 309.0893
But:
c(-3.62,3)^3
Result: [1] -47.43793
Using some simulated data:
A - rnorm(7, mean=1); B - rnorm(9, mean=2); C - rnorm(13, mean=2.5)
y - c(A, B, C)
f - factor(rep(1:3, c(7, 9, 13)))
TukeyHSD(aov(y~f))
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = y ~ f)
$f
difflwr
On 13 Aug 03, at 11:26, Harold Doran wrote:
I believe it is in the Rcmdr package, which requires the car library to be loaded.
You can also perform an ANOVA using the absolute value of the deviations from each
respective group mean, which is what Levene's Test does.
The function just
Dear all,
is there any function to calculate the center of a polygon mass in R?
Actually I need to find the best location within polygons to place labels.
Thanks for any hint
Jens Oehlschlägel
--
COMPUTERBILD 15/03: Premium-e-mail-Dienste im Test\ ...{{dropped}}
a$b = a[['b']] = attribute b of list a.
A basic object in R is a list, and the $ operator provides one means
of accessing named attributes of a list.
Beginning with R 1.7, objects can also have slots, which are
accessed as [EMAIL PROTECTED]. I have yet to understand why slots were
I would like to do a one way anova and then a tukeyHSD. I have three vectors
A,B and C. In a previous help message, I was told to do the following for
the anova:
y = c(A,B,C)
group = factor(rep(a:3,c(7,9,13))) #provided there a 7 elements in A,9 in B
and 13 in C
and then
Hi,
Sorry, I was to fast in sending the mail to the list.
I could find the error myself.
^operator binds stronger than -operator, so -3^3.2 calculates as
-(3)^3.2 and not as (-3)^3.2 which results correctely in NaN.
Alex
__
[EMAIL PROTECTED] mailing
With regards to my question below, I've found a way to do it:
all.names(body(test))
Seems to do the trick. Still, if there are any tools already developed
for mapping the dependencies between a set of functions I be very
interested to hear before I reinvent the wheel.
Thanks,
Dave
Hey, R-listers,
I have a question about determining the orthogonal
basis vectors.
In the d-dimensinonal space, if I already know
the first r orthogonal basis vectors, should I be
able to determine the remaining d-r orthognal basis
vectors automatically?
Or the answer is not unique?
Thanks for
I have grouped data and i used a Trellis plot to show the curve of a
fitted model for each group. Additionally i have a vector whose values
are the highest points the curve can reach (the curve rises
exponetially). Each value of the vector is the corresponding highest
point for each
Anna H. Pryor [EMAIL PROTECTED] writes:
Does anyone know of some methods already programmed up ( ie freeware
We prefer the term Open Source to freeware.
: ) ) in R for the post hoc methods in anova. Particularly
Scheffe's method or Tukey's?
?TukeyHSD
See also package multcomp.
On Wed, 2003-08-13 at 07:55, [EMAIL PROTECTED] wrote:
Hi Folks,
With a bit of experimentation I have determined (I think)
that on my R implementation the largest positive integer
that is exactly represented is (2^53 - 1), based on
(((2^53)-1)+1) - ((2^53)-1)
[1] 1
((2^53)+1) - (2^53)
Please read the INSTALL file that comes with the package.
Actually reading the documentation can be wonderfully liberating, and is a
skill that can work wonders in life. Please get some practice in.
You might even work out which list supports Omegahat packages.
On Tue, 12 Aug 2003, Jin, Ying
Then the best book is Venables Ripley (2002) for R.
A.S.
Alessandro Semeria
Models and Simulations Laboratory
Montecatini Environmental Research Center (Edison Group),
Via Ciro Menotti 48,
48023 Marina di Ravenna (RA), Italy
Tel. +39 544 536811
Fax. +39 544 538663
On Wed, 13 Aug 2003, Paul Sorenson wrote:
How can I have SQL queries (ODBB/Windows) return dates (cf text) for
date columns?
I can successfully query and ODBC database with R. I can subsequently
convert the text representations of dates in columns using as.POSIXxx
functions but I figured
The following code works, to gsub single quotes to double quotes:
line - gsub(', '', line)
(that's a single quote within doubles then a double within singles if
your
viewer's font is not good).
But The R Language Manual tells me that
Quotes and other special characters within strings
are
Hi all,
I want to conduct a cluster analysis with quantitative variables.
More precisely, it concerns binary and non-ordered categorical
variables. For such data, various
similarity measures have been proposed, such as the Jaccard index or the
simple matching index.
So, is there a package
Usman Shehu wrote:
Dear R- Users,
Please, what is wrong ? I am trying to run one sample Kolmogorov-Smirov test but I always get an error as shown below.
ks.test(A, normal, 0,1)
For sure, you mean
ks.test(A, pnorm, 0,1)
because pnorm() is the distribution function of the normal (Gauss)
Thompson, Trevor wrote:
It looks like R is treating every character in the string as if it were
decimal. I didn't see anything in the help file about . being some kind
of special character. Any idea why R is treating a decimal this way in
these functions? Any suggestions how to get around
Dear all,
Thank for your answers and your help. My problem was that, during
exportation with write.table, numeric variable precision corresponds
to the printed presentation and was then too low for me. But you can set
this option options(digits=XX) before exporting.
Thank again
David
Try with write.matrix (MASS pkg)...
good luck
A.S.
Alessandro Semeria
Models and Simulations Laboratory
Montecatini Environmental Research Center (Edison Group),
Via Ciro Menotti 48,
48023 Marina di Ravenna (RA), Italy
Tel. +39 544 536811
Fax. +39 544 538663
Try
regexpr(\\., Female.Alabama) and gsub(\\., ,, Female.Alabama)
X-Sybari-Trust: 9293cd92 d90ef28b 235e1558 093d
From: Thompson, Trevor [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Date: Wed, 13 Aug 2003 10:46:45 -0400
MIME-Version: 1.0
X-Virus-Scanned: by amavisd-milter (http://amavis.org/)
Greetings,
I am trying to use nlme to model a large data set
of pharmacokinetic concentrations with sparse data per cluster.
The base model is a one compartment model with first order absorption.
The program runs but the PNLS step never converges. The parameter
estimates are reasonable.
This is becoming a relevant ESS topic, hence I divert it to the
ESS-help mailing list, please continue there
{archives at https://www.stat.math.ethz.ch/pipermail/ess-help/ )
tony == A J Rossini [EMAIL PROTECTED]
on Tue, 12 Aug 2003 20:07:35 -0700 writes:
tony kjetil brinchmann
In addition to Gelman et al, there is also Jeff Gill's BAYESIAN METHODS FOR
THE SOCIAL AND BEHAVIORAL SCIENCES with accompanying R/SPLUS code. You
have to be careful though, at least in the first printing, it isn't always
clear when the code is for R versus SPlus. The code is available on
his
Hello everyone
I have data grouped by subjects and i used the trellis plot to show the curve for each
subject. I have another vector a which gives me the highest points the curves for
each subjects can achieve ( the curve is exponential ). Now i want to draw a
horizontal line, whose values are
This is covered in the FAQ.
The most comprehensive accounts of non-linear regression I know of
are in Chambers Hastie (1992, for S) and Venables Ripley (2002 and
earlier editions, for R or S).
On Wed, 13 Aug 2003, anne wrote:
As a newbie to R, I need to learn my way around (no previous
On Wed, 13 Aug 2003, Richard A. O'Keefe wrote:
While reading through some of the R source code, I have come
across forms such as
m$...
m$... - e
and I wondered what they meant.
?$ mentions x$name, but not $...
... is a name, a special one in the context of function calls.
All
Petr Pikal [EMAIL PROTECTED] wrote (embedded in much XML):
Not sure about efficency and it is not very general solution but
you maybe can use embed() function
It's an interesting suggestion, but I don't see *how*.
Here's the function I want, only I'd like it to be something built
On Tue, 12 Aug 2003, Uwe Ligges wrote:
Simon Fear wrote:
Uwe, you suggest
So you almost got it, just read ?shell carefully enough:
filenames - shell(dir D:\\tmp\\*.sasb7dat /b, intern = TRUE)
Unfortunately, this does not work and is why I wrote to Rhelp (OK I used
Dibakar Ray [EMAIL PROTECTED] writes:
I am very new to R. I was trying to load some publicly available Expression
data in to R.
I used the following commands
mydata-read.table(dataALLAMLtrain.txt, header=TRUE, sep
=\t,row.names=NULL)
It reads data without any error
Now if I use
Thank you. Single backslash version, first thing I tried (I thought)
works
just fine when I copy and paste, ergo I must have got confused by some
stupid
typo of mine. Sorry to waste everyone's time over this. (Still, I am
probably not the only confused user when it comes to RE handling - I
hope
On Tue, 12 Aug 2003, Simon Fear wrote:
The following code works, to gsub single quotes to double quotes:
line - gsub(', '', line)
(that's a single quote within doubles then a double within singles if
your
viewer's font is not good).
But The R Language Manual tells me that
Quotes
On Tue, 12 Aug 2003, Gerhard Prade wrote:
i think i am to silly. I have installed R 1.7.1 (2003-06-16). Installed
some packages like xtables ore xml. I tried out this to installing
packages. Then i tried to make a crosstable like i know it from spss.
They say in this list that it would be
Hello all,
i think i am to silly. I have installed R 1.7.1 (2003-06-16). Installed
some packages like xtables ore xml. I tried out this to installing
packages. Then i tried to make a crosstable like i know it from spss.
They say in this list that it would be going.
I made a table in asci-format,
Hi.
(B
(BI've been very interested in your seao seao-gui package for R.
(B
(BHowever manuals of seao seao-gui package for R have doesn't include
(Bexamples, so I can't understand how
(Bto use them.
(B
(BCould you tell me availability of tutorials of them in English.
(B
(B Regards.
(B
nels.tomlinson.1 [EMAIL PROTECTED] writes:
I'd like to propose that now that the underscore-as-assignment-operator
is to be removed from R (good thing, too), that the Windows GUI should
replace the underscore ``_'' with the proper assignment operator ``-''
when you type in the underscore
What about ...
unlist(lis)[which(unlist(lis)==next)+1]
[1] want1 want2
... to avoid the loop in sapply?
Richard A. O'Keefe wrote:
Chris Knight [EMAIL PROTECTED] has
lis-list(c(a,b,next,want1,c),c(d, next, want2, a))
and wants c(want1,want2)
Step 1:
inx -
Dear All,
I am trying to install R using debian's apt-get system. When I do 'apt-get install
r-base', it appears to install everything correctly. However, when I start up R I get
the error message:
[EMAIL PROTECTED]:/knoppix$ R
cannot find system RenvironError in options(...) : invalid editor
Prof Brian D Ripley [EMAIL PROTECTED] writes:
On Tue, 12 Aug 2003, [iso-8859-1] Bjørn-Helge Mevik wrote:
Also, is this example (lm(y~x+I(x^2), Df)) really balanced? I think
No, and I did not use summary,aov on it!
And I didn't say you did!
This gives the SSs R(x | A, B, A:B, x^2), R(x^2
On Tue, 12 Aug 2003, Ciril Rozman wrote:
Sirs,
I have recently been interested in your Maxstat.
For questions related to add-on packages please at least cc to the
maintainer since some might not read r-help regulary.
I have computed
with my own programme the ranks (by using the
Have you considered adding noise to predict(spline.model)$y? If this
won't solve your problem, then I think I don't understand what you want
to do.
spencer graves
[EMAIL PROTECTED] wrote:
# I need to generate some data. I'm modeling some time series that follow
a
# negative exponential decay
Hi,
I'm trying to use the command glm.nb in library(MASS) to test for a significant
difference in the aggregation parameter theta between the three levels of a factor.
Any help gratefully received!
Martin.
Martin Hoyle,
School of Life and Environmental Sciences,
University of Nottingham,
I have a time series of data from an electroencephalogram (EEG).
I wish to filter the data to get rid of 50Hz mains 'hum'. I have
'designed' a combination bandpass and notch filter using a web-
site. The site returns the filter in ANSI C source code. It is:-
/* Digital filter designed by
Hello All
I am currently looking at spatial data - Chorophyll A concentration in sea
water over a wide geographic area. These data are used to determine the
location of ocean fronts and hence where tuna are located. A front is
identified by a steep gradient in the change in chloroA
Johnathan Williams wrote:
Would anyone be so kind as to write a routine to time
mouse button presses in R to the nearest millisecond?
If R had a timer of this kind and a few basic screen
handling routines (to write characters or graphics of
different sizes and colours at precise times)
Jim Lemon [EMAIL PROTECTED] writes:
Johnathan Williams wrote:
Would anyone be so kind as to write a routine to time
mouse button presses in R to the nearest millisecond?
If R had a timer of this kind and a few basic screen
handling routines (to write characters or graphics of
Anova != anova.
drop1 is the part of R that does type II sum of squares, and it works in
your example. So does Anova in the current car:
drop1(lm(y~x+I(x^2), Df)) # add test=F if you like
Single term deletions
Model:
y ~ x + I(x^2)
Df Sum of SqRSSAIC
none 8.3117
1. Have you considered exporting to a txt or csv file and then using
read.table?
2. Have you considered hist?
3. For other questions, I suggest you explore www.r-project.org -
search - R site search before emailing this list.
hope this helps.
spencer graves
Aimin Yan wrote:
Hello,
How to
Barry Rowlingson [EMAIL PROTECTED] wrote:
Do you want:
c(x[-1],x[1]) for a one-step 'rotation'?
That's just the kind of thing I did, except that it's ugly.
I've browsed src/main/subscript.c and got rather lost, but
it looks very much as though x[-1] starts by
Dear Aimin,
Your question about Excel comes up rather often on the list
and answers can be found in the archives in the R Data
Import/Export PDF file installed with R. In short, you can
read CSV files through read.csv or you can get read
directly from Excel spreadsheets through the excellent
Thank you to all who helped me with my previous question regarding random
numbers. I have read up on the subject and found a whole new interesting world.
I found some code wich I translated into java for my project and now have a few
pseudorandom number generators. I was wondering if R had any
Dennis Fisher wrote:
In labeling axes, I want to combine symbols and text/superscripts.
Examples include:
m2 (m, followed by a superscripted 2)
µg (micrograms)
How can I accomplish this in R?
See ?plotmath, e.g.:
plot(1:10, xlab=expression(m[2]), ylab=expression(mu*g))
Uwe Ligges
In labeling axes, I want to combine symbols and text/superscripts.
Examples include:
m2 (m, followed by a superscripted 2)
µg (micrograms)
How can I accomplish this in R?
Dennis Fisher MD
P (The P Less Than Company)
Phone: 1-866-PLessThan (1-866-753-8864)
Fax: 1-415-564-2220
www.PLessThan.com
Christian Mora wrote:
Hi all,
Ive got a database with 10 columns (different
variables) for 100 subjects, each column with
different # of NA's. I'd like to know if it is
possible to use a function to exclude the NA's using
only a specific column, lets say:
Data2 - omit.exclude(Data1$column1) ??,
Dear all,
I've worked with optim before but never encountered this error message:
Nelder-Mead direct search function minimizer
0.23 0Error: subscript out of bounds
The error seems to depend on the initial parameter values. However,
strangely (I think), I recieve this error message when I have
In
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major1
minor7.1
year 2003
month06
day 16
language R
On Thu, 2003-08-07 at 11:07, Al Piszcz wrote:
I am performing this sequence
barplot
title
legend
abline
When abline renders the lines they appare to be in the layer
above the bars in the graph. Is there a way to make them
render first or 'behind' the bars?
Thanks.
Are you using
### First, is there a way to access a univariate time series as a matrix
instead of a vector?
# For example:
data(UKLungDeaths)
# If I do
apply(mdeaths,1,cumsum)
# Gives an error as mdeaths is not a matrix but a vector, although when I
look at it :
mdeaths
# the ts object has a matrix like
Hi,
I'm working on a set of demo files for a package, and I'm having a problem because I
have two demo files (one called logit, the other called blogit) and when I type
demo(logit), it gives me demo(blogit) instead. (And calls it the the demo for logit.)
I've figured out that this is
dat - data.frame(x=c(1,NA,2,2),y=c(3,2,NA,1))
dat
dat[rownames(na.omit(dat[,y,drop=F])),]
dat - data.frame(x=c(1,NA,2,2),y=c(3,2,NA,1))
dat
x y
1 1 3
2 NA 2
3 2 NA
4 2 1
dat[rownames(na.omit(dat[,y,drop=F])),]
x y
1 1 3
2 NA 2
4 2 1
HTH,
Jerome
On August 7, 2003 01:05 pm,
Curso GIS y Teledeteccion
Para solicitar información mandar mail a [EMAIL PROTECTED]
CENTRO DE TRANSFERENCIA TECNOLÓGICA EN SISTEMAS DE INFORMACIÓN GEOGRÁFICA Y
TELEDETECCIÓN
---TUTORIAL GRATUITO
Entendiendo la proyección de los mapas - SISTEMA Gauss-Krüger
---CLASES GRATUITAS: Primer módulo
On Wed, 2003-08-06 at 16:38, Marc Schwartz wrote:
On Wed, 2003-08-06 at 16:21, Frank E Harrell Jr wrote:
In
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major1
minor
Dear List,
I have a problem building an R package using R 1.7.1 that I cannot resolve
myself. R CMD check tells me that functions I documented properly in Rd-files
are not documented and gives a TeX error at the end, altough the Rd files
check is OK. I have three fucntions, one gets through
enclosed a simple R script (and a data file, and the output) , with
calls lda similar to the example with the iris data in the
documentation. it is not working and i dont understand the error
message. can anybody help me? i am using R 1.5.1 (2002.06.17) on
debian woody stable.
I would like to
Hi
Can anyone help with the technique of obtaining leverages from a
conditional logistic regression model? The code lm.influence does not seem
to work for this data.
Thanks
Jane Murray
__
[EMAIL PROTECTED] mailing list
Those are not predicted values, they are fitted values. Try predicting on
the same set of variables as you printed.
Please do try to give a small reproducible example so we can see what you
actually did.
On Thu, 7 Aug 2003, orkun wrote:
I carried out a logistic regression and found
Yao, Minghua wrote:
All,
Anybody can tell how to export a R plot onto Word (or Power Point)?
Many thanks in advance.
-MY
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Since you are presumably on Windows (and
Dear R-Helpers,
I try to fit my x and y vector-data with a power law using a the
following command:
test - nls(y ~ A*x^B, xy, start=list(A=0.5,B=0.8))
and I get the error message:
Error in numericDeriv(form[[3]], names(ind), env) :
Missing value or an Infinity produced when evaluating
As previous replies have suggested -- you could use an array
E.g.
X = array(0, dim=c(100,12,12))
Then X[i, , ] returns the i-th 12 by 12 matrix slice in this array,
e.g.
X[1, ,]
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
[1,]00000000
201 - 280 of 280 matches
Mail list logo