[R] netcdf read/write package

2003-06-28 Thread Viet Nguyen
Hi everyone,

I'm looking for an R package that is capable of _write_ as well as read 
in netcdf format.  Could someone offer advice on this?

I have some data in netcdf format and I only need to make simple 
modifications to it.
I thought doing that via R could be handy and would save me reading up 
on C++/Fortran netcdf lib.

Thanks and regards,
viet nguyen
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help


[R] question

2003-06-28 Thread Jan de Leeuw
Is it correct that

lsfit(matrix(0,10,1),1:10,intercept=FALSE)

returns zero residuals ?
===
Jan de Leeuw; Professor and Chair, UCLA Department of Statistics;
Editor: Journal of Multivariate Analysis, Journal of Statistical  
Software
US mail: 9432 Boelter Hall, Box 951554, Los Angeles, CA 90095-1554
phone (310)-825-9550;  fax (310)-206-5658;  email: [EMAIL PROTECTED]
homepage: http://gifi.stat.ucla.edu
  
 
-
  No matter where you go, there you are. --- Buckaroo Banzai
   http://gifi.stat.ucla.edu/sounds/nomatter.au

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


Re: [R] histogram cells

2003-06-28 Thread Prof Brian Ripley
?hist says

 breaks: one of:

...
 *  a single number giving the number of cells for the
histogram,

...

  In the last three cases the number is a suggestion only. 


On Fri, 27 Jun 2003, Tommy E. Cathey wrote:

 Why does the following code generate a Histogram of 7 cells
 instead of 5?

  data(islands)
  str(hist(sqrt(islands), br = 5, col=lightblue, border=pink))
 List of 7


-- 
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] question

2003-06-28 Thread Prof Brian Ripley
On Sat, 28 Jun 2003, Jan de Leeuw wrote:

 Is it correct that
 
 lsfit(matrix(0,10,1),1:10,intercept=FALSE)
 
 returns zero residuals ?

No.  I would use lm.fit for such problems, and that seems to behave 
(probably rather coincidentally).

-- 
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] sprintf(%c,...)

2003-06-28 Thread Uwe Ligges
Gabor Grothendieck wrote:
How does one get the character that corresponds to a number?  Also
the inverse?
The %c format code is not supported by the R sprintf function and I
could not find another way to do it.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
The ASCII character table:

ascii - sapply(1:255, function(i)
  parse(text=paste(\\\, structure(i, class = octmode), \,
sep = ))[[1]])
Uwe Ligges

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


Re: [R] R-help Digest, Vol 4, Issue 27 ( -Reply)

2003-06-28 Thread Peter Dalgaard BSA
Leo Wang-Kit Cheung [EMAIL PROTECTED] writes:

 Hi,
 
 I am out of town and will get back to you on the 13th of July.
 
 Leo
 
  [EMAIL PROTECTED] 06/27/03 00:32 
 
 Send R-help mailing list submissions to
   [EMAIL PROTECTED]
 
 To subscribe or unsubscribe via the World Wide Web, visit
   https://www.stat.math.ethz.ch/mailman/listinfo/r-help
 or, via email, send a message with subject or body 'help' to
   [EMAIL PROTECTED]
 
 You can reach the person managing the list at
   [EMAIL PROTECTED]
 
 When replying, please edit your Subject line so it is more specific
 than Re: Contents of R-help digest...
 
 
 Today's Topics:
 
1. create help files ([EMAIL PROTECTED])

...and a full week of digested messages gets quoted back to the list.
Let's hope that not every single digest-subscriber does likewise when
he/she goes on holiday!

-- 
   O__   Peter Dalgaard Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics 2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark  Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907

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


[R] R and RMySQL - segmentation fault

2003-06-28 Thread Xavier Fernández i Marín
Hello, 

I want to import data from a database in Mysql to a R data frame.
MySQL is running ok on my computer.

I can start library RMySQL, and connect to the database, but when I try to 
make an statement, R says segmentation fault and stops:

 rs - dbExecStatement (con, select * from labels)
Segmentation fault
[EMAIL PROTECTED]

I'm running R 1.7.1 on a Linux Mandrake 9.0.

Thank you,



Xavier Fernández i Marín
Universitat Pompeu Fabra - Barcelona

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


[R] Re: Application

2003-06-28 Thread NJStar Postmaster
Hi,

Thanks for your email, but your email contains attachement of 
executable program which may contain virus.
 
Please re-send the file in a zip file.

NJStar Mail Processor.

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


[R] Truncating y axis

2003-06-28 Thread Shravan Vasishth
Hi folks,

I want to have a y-axis using matplot where the y axis is truncated as
shown schematically below:

(msecs)

3000 |   x
 ~ 
 ~  
800  |x
 |
 | x 
 |
 |||-
 123  
  Position

The ~ is supposed to stand for a curved line showing a break. The issue is
that I want to plot the points shown by the x's, but if I plot a
continuous y axis then the confidence intervals about the points at
positions 2 and 3 are not clearly visible...

Can someone suggest how this might be done? I did study the manuals and
contributed manuals first, and searched the archives and FAQs, just
couldn't find any discussion of such a situation.

Thanks much in advance,

-- 
Shravan Vasishth   Phone: +49 (681) 302 4504
Computational Linguistics, Universität des Saarlandes, Postfach 15 11 50
D-66041 Saarbrücken, Germany http://www.coli.uni-sb.de/~vasishth

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


Re: [R] Truncating y axis

2003-06-28 Thread Frank E Harrell Jr
On Sat, 28 Jun 2003 14:03:27 +0200 (MEST)
Shravan Vasishth [EMAIL PROTECTED] wrote:

 Hi folks,
 
 I want to have a y-axis using matplot where the y axis is truncated as
 shown schematically below:
 
 (msecs)
 
 3000 |   x
  ~ 
  ~  
 800  |x
  |
  | x 
  |
  |||-
  123  
   Position
 
 The ~ is supposed to stand for a curved line showing a break. The issue is
 that I want to plot the points shown by the x's, but if I plot a
 continuous y axis then the confidence intervals about the points at
 positions 2 and 3 are not clearly visible...
 
 Can someone suggest how this might be done? I did study the manuals and
 contributed manuals first, and searched the archives and FAQs, just
 couldn't find any discussion of such a situation.
 
 Thanks much in advance,
 
 -- 
 Shravan Vasishth   Phone: +49 (681) 302 4504

See W. Cleveland The Elements of Graphing Data (Hobart Press) for reasons not to do 
this.Consider using a nonlinear scale or making the graph taller.

---
Frank E Harrell Jr  Prof. of Biostatistics  Statistics
Div. of Biostatistics  Epidem. Dept. of Health Evaluation Sciences
U. Virginia School of Medicine  http://hesweb1.med.virginia.edu/biostat

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


Re: [R] Truncating y axis

2003-06-28 Thread Wolfgang Viechtbauer
 I want to have a y-axis using matplot where the y axis is truncated as
 shown schematically below:

One (not so elegant) way is this:

# original data with one very large y value
x - c(1, 2, 3, 4, 5)
y - c(1000, 120, 110, 108, 104)
plot(x,y)

# truncated y-axis plot
y - c(140, 120, 110, 108, 104)
plot(x,y, yaxt=n)
axis(side=2, at=c(100, 105, 110, 115, 120, 140), labels=c(100, 105, 110, 115, 120, 
1000))
rect(0, 130, 1, 131, col=white, border=white)
par(xpd=T)
lines(x=c(0.7,1), y=c(130, 130))
lines(x=c(0.7,1), y=c(131, 131))

--
Wolfgang Viechtbauer

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