[R] Plot ticks and tick labels: thickness, colour?

2003-08-10 Thread Dirk Eddelbuettel

I am displaying several series in one plot, and would like to make
them distinct without having to employ a legend. 

I managed to color tickmarks, but have been unsuccessful with either one of 

a) making tickmarks thicker (without increasing the axis at the same time).
   From reading ?axis:
   lty, lwd: line type, width for the axis line and the tick marks.
   in would appear that I cannot obtain the one _without_ the other; or
   
b) displaying the tick label in a different colour. Again, ?axis reads
 col: color for the axis line and the tick marks. [..]
   indicates that I can only set the tick mark, not the annotation.
   
Any ideas or suggestions?

Thanks in advance,  Dirk

-- 
Those are my principles, and if you don't like them... well, I have others.
-- Groucho Marx

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] Question about 'NA'

2003-08-10 Thread Christian Mora
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) ??, then
Data3 - omit.exclude(Data1$column2) and so on 

I tried the code above but with no results

Thanks for any help

CM

__

Yahoo! SiteBuilder - Free, easy-to-use web site design software

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] gregmisc

2003-08-10 Thread Stefano Iacus
In Carbon R 1.7.1 recommnded packages are already installed. VR itself 
is a bundle of package and not a package itself.
MASS, rpart, spatial and nnet are already present in Carbon R.

BTW, From 1.8.0 MacOS X users can take advange of the new (forthcoming) 
RAqua (http://www.economia.unimi.it/R/ for a preview)

stefano

On Venerdì, ago 8, 2003, at 04:10 Europe/Rome, Marc Schwartz wrote:

On Thu, 2003-08-07 at 18:42, atsuya fujito wrote:
Thank you, Marc;

I installed car, R commander and Bioconductor.
I am using R 1.7.1, on MAC OSX 10.2.6
[Long listing of packages snipped]

Atsuya,

Thanks for the listing. If that is the complete list, then you are
missing the VR bundle of packages, which includes MASS and is
typically part of the 'recommended' package bundle, at least under
Windows and Linux/Unix.
I did a search of the Mac FAQ (**which appears to be for 1.6.0**), the 
R
Admin manual and the r-help archive, since I do not have hands on
experience with Macs. I may have to defer to other Mac OSX users for
detailed guidance here.

I also just checked Jan de Leeuw's web site at UCLA
(http://gifi.stat.ucla.edu/pub/index.php) based upon a review of the
README at http://cran.r-project.org/bin/macosx/ReadMe.txt. The
indication is that the Mac port for V1.7.0 includes the base and
recommended packages, however V1.7.1 is the base package set only, if 
my
read is correct.

If that is indeed the case, the solution to your problem may be as
simple as using:
install.packages(VR)

being sure that you have appropriate R/W access permissions to the
library directory tree on your system when you use that R command.
If that solves the problem, then you should be ok and can load
'gregmisc' after VR is installed.
I am copying both Jan de Leeuw and Stefano Iacus on this communication
to be sure that I am offering appropriate advice here. Given the 
various
time zones involved here, hopefully someone can expediently confirm the
resolution to your problem.

I hope this helps.

Regards,

Marc Schwartz


__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] probability plot correlation coefficient

2003-08-10 Thread anne
As a newbie to R, I'm still rather at a loss for finding information 
(the commands names can be rather arcane)so I'm just posting my question:
I would like to estimate   the shape coefficient of diverse 
distributions (Weibull, gamma and Tukey-Lambda specifically, but other 
could be of interest)
- Does R have a PPCC utility to estimate such parameter?(maximum value 
of correlation coef)
- If yes how does one retrieve the numerical value from the graph? (see 
graphical example below)
- The retrieval of numerical values is also a problem for me from the 
probability plots



Thank for any help!
Anne
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] installing add-on packages on OS X

2003-08-10 Thread Philippe Glaziou
mpie [EMAIL PROTECTED] wrote:
 This is probably a stupid question, but since I'm a
 newbie on R, here it goes. I got R 1.6.2 running on OS
 10.1.5 and I'm trying to install an add-on package.
 
 I couldn't find info on installing this in a mac. Can
 anyone give me some pointers?


See R-admin. 

It may be a good idea to consider upgrading your R. It
compiles very well on MacOSX, provided that you install all
the missing soft that Apple could not be bothered to 
include in the OS. 

-- 
Philippe Glaziou
Pasteur Institute of Cambodia

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] 3D matrix

2003-08-10 Thread Prof Brian Ripley
Do you mean a 3-dimensional array?  If so see the help for array().
In S/R terminology, a matrix is 2D, by definition.

On Tue, 5 Aug 2003, Marcel Vieira wrote:

 Is it possible to define a 3 dimension matrix in R?
 (without using list - if possible)

-- 
Brian D. Ripley,  [EMAIL PROTECTED]
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel:  +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax:  +44 1865 272595

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] robust regression

2003-08-10 Thread Joerg Schaber
Hi,

the package is called quantreg and yes, it solved my problem. LAD 
regression was actually the method I was looking for, thanks a lot.
However, some of my problems are rather large and even if I use the 
method 'fn' and 'pfn' recommended for large problems in 'rq' I get an 
error:

 res - rq(o ~ y + s 
-1,tau=0.5,method=fn,contrasts=list(s=(contr.sum)))
Error: cannot allocate vector of size 1950212 Kb

Is there are way to avoid that? The least square procedures are able to 
handle these large problems.

best wishes,

joerg

Patrick Burns wrote:

Two possible alternatives are:

1) Use least absolute deviation regression, which you can
get from Roger Koenker's package that I think is called
quantile.
2) The LAD regression is essentially equivalent to median
polish (see, for instance, Understanding Robust and Exploratory
Data Analysis by Hoaglin, Mosteller and Tukey).  To gain some
more efficiency for nearly Gaussian data, you could replace the
median by a location estimator that is more efficient, such as a
trimmed mean.
Good luck,

Patrick Burns

Burns Statistics
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S User)
Joerg Schaber wrote:

Hi,

trying to do a robust regression of a two-way linear model, I keep 
getting the following error:

 lqs(obs ~ y + s -1,method=lms, contrasts=list(s=(contr.sum)))
Error: lqs failed: all the samples were singular
Robust regression with M-estimators works (also regular least square 
fits, of course):
rlm.formula(formula = obs ~ y + s - 1, method = M, contrasts = 
list(s = (contr.sum)))

I tried an exact sampling (psamp=exact), but I keep getting syntax 
errors. Any idea how I can make the first one work?

Thanks,

joerg




--
--
Jörg Schaber
Instituto Cavanilles de Biodiversidad y Biologia Evolutiva
Universidad de Valencia   Tel.: ++34 96 354 3666
A.C. 22085Fax.: ++34 96 354 3670
46071 Valencia, Españaemail : [EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] L10N and i18n of R

2003-08-10 Thread Shigeru Mase
Dear R-users of non-English speaking countires.

Maybe a good news for those who want to use R in their
local languages. Recently two Japanese, E. Nakama and
M. Okada, succeeded in making R speak Japanese. At
present, R can handle Japanese (as well as other languages,
I guess) character strings if one use consoles which
can understand Japanese. Also P. Murrell kindly
provided a facility to plot hundreads of Japaense
characaters as graphical symbols (Hershey vector fonts).

Nakama and Okada's L10N (Localization) and i18n
(internationalization) patches make it possible to
use Japanese object names and to display Japanese
characters on graphical devices (although still limited).
If you are interested, please visit URL:

http://www.okada.jp.org/RWiki/index.php?%5B%5Bi18n_of_R%5D%5D

which is the only English page of RjpWiki, a Wiki-based
collaboration site supported by Japanese R users.

Nakama's patches can also handle other languages
than Japanese. He kindly built Korean and Russian
version of R (1.7.1) as rpm binaries (although
he can understand neither Korean nor Russian).

We are glad if these patches can give hints to R
users who want to use local languages. Please note
their patches are by no means complete and may
potentially give your systems troubles. They and
related binaries are offered without no warranty.
Feedbacks and reports are welcome, but mere inquiries
and complains not. Since the version up of R is so
frequent, we cannot promise that these patches will
be updated in future.

Yours

Shigeru Mase
Dept. Math. and Comp. Sciences,
Tokyo Institute of Technology, Tokyo, Japan

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


Re: [R] packages on R on a Mac

2003-08-10 Thread Thomas Lumley
On Tue, 22 Jul 2003, Roger Coupal wrote:

 Hello, i am new to R and have it on my Mac with OSX 10.2. I downloaded
 the systemfit package and tried to run it and it didn't work. I think i
 need to install the package, or update R to let it know that I have
 that in the library. (i simply placed the systemfit folder in the R
 library.) When i run update.packages(systemfit) or
 install.packages(systemfit) i get the following error:

 Error: couldn't find function install.packages

You don't say whether you are using the Carbon GUI version or the Darwin
version, but if install.packages() doesn't exist then I assume the former.

The Carbon version needs precompiled binary packages (In principle you
could compile them yourself, but that is not straightforward for this
version).  The precompiled binary packages live at
http://cran.r-project.org/bin/macos/contrib/

As systemfit has no compiled code then I believe a binary package created
on any system will work, so you could get a Windows zip file and unzip in
in the library directory.  Note that I have only tried this with a package
created under Darwin.

This system is clearly unsatisfactory, and that's why our Mac guru,
Stefano Iacus, is now working on a GUI for the Darwin version rather than
compiling packages.  In the current development version of RAqua
install.packages() works very nicely.


-thomas

Thomas Lumley   Assoc. Professor, Biostatistics
[EMAIL PROTECTED]   University of Washington, Seattle

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] read.spss doesn't work anymore

2003-08-10 Thread Hedderik van Rijn
A couple of months ago, probably using an older version of R, R used to 
run the following code just fine:

 library(foreign)
 data.exp1 - as.data.frame(read.spss(dataDef.sav))
Issuing the same commands now (after starting R using --vanilla), gives 
me the following behavior:

 library(foreign)
 x - read.spss(dataDef.sav)
Error in read.spss(dataDef.sav) : Calloc could not allocate 
(-2147483648 of 1) memory
 x - read.spss(dataDef.sav)
 x - read.spss(dataDef.sav)
Segmentation fault

The first two assignments return instantaneously, after the second x 
contains all original SPSS variable/column names but no data, after the 
last read.spss - it takes a while before R returns with the 
Segmentation fault. During that time, the harddisk seems to be working 
quite hard.

Does anyone know what the problem might be?

 - Hedderik.

 version
 _
platform i386-pc-linux-gnu
arch i386
os   linux-gnu
system   i386, linux-gnu
status
major1
minor7.1
year 2003
month06
day  16
language R
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


RE: [R] index all subsets of k of m items?

2003-08-10 Thread Liaw, Andy
help.search(subset) on my R installation gives:

nchoosek(vsn)List all subsets of size k from n objects

and vsn is one of the Bioconductor packages.  From that, the
index.subsets() you want shouldn't be too hard to construct...

HTH,
Andy

 -Original Message-
 From: Spencer Graves [mailto:[EMAIL PROTECTED] 
 Sent: Friday, August 01, 2003 6:56 AM
 To: R-help
 Subject: [R] index all subsets of k of m items?
 
 
 How can I efficiently index all choose(m, k) subsets 
 of m items taken 
 k at a time?  For example, with (m, k) = (3, 2), the subsets 
 are (1, 2), 
 (1, 3), and (2, 3).  I'd like a function something like 
 index.subsets(subset, k, m) that would return 1, 2 or 3 for these 3 
 subsets.  Examples:
 
 index.subsets(c(1,2), 2, 3) - 1
 index.subsets(c(1,3), 2, 3) - 2
 index.subsets(c(2,3), 2, 3) - 3
 
 index.subsets(c(1,2,3), 3, 5) - 1
 index.subsets(c(1,2,5), 3, 5) - 3
 index.subsets(c(3,4,5), 3, 5) - 10
 
 Thanks.
 Spencer Graves
 
 __
 [EMAIL PROTECTED] mailing list 
 https://www.stat.math.ethz.ch/mailman/listinfo /r-help
 

--
Notice:  This e-mail message, together with any attachments, contains
information of Merck  Co., Inc. (Whitehouse Station, New Jersey, USA), and/or
its affiliates (which may be known outside the United States as Merck Frosst,
Merck Sharp  Dohme or MSD) that may be confidential, proprietary copyrighted
and/or legally privileged, and is intended solely for the use of the
individual or entity named on this message.  If you are not the intended
recipient, and have received this message in error, please immediately return
this by e-mail and then delete it.

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] high memory allocation

2003-08-10 Thread gowuban
Hello,

I have trouble with my cluster analysis using package cluster. diana and agnes 
both seem to try to allocate memory directly, so I can not use virtual memory of my 
Windows2000 operation system.
I do have 320 MB of memory. But they claim about 600 MB. Do I have a chance to do the 
analysis with my amount of memory. 
Thanks for all comments, I did not find a way yet.

Regards
Georg

[[alternative HTML version deleted]]

__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help