[R] How could plot a chart with some Chinese's words text?

2003-06-30 Thread Ken Lee
Dear all, I want to plot a chart title with some Chinese words at unix system, but I do not how to do. first, it is no problem when I use write.table like: a- write.table(a,file=z.txt,sep=\t) second, I try bitmap(file=z.png) plot(1:10,main=a) dev.off() I can't get a chart with the

Re: [R] How could plot a chart with some Chinese's words text?

2003-06-30 Thread Prof Brian Ripley
On Mon, 30 Jun 2003, Ken Lee wrote: Dear all, I want to plot a chart title with some Chinese words at unix system, but I do not how to do. first, it is no problem when I use write.table like: a-中文 write.table(a,file=z.txt,sep=\t) Not, that is not Chinese for me and I suspect

Re: [R] ./configure needs /sw/lib?

2003-06-30 Thread Prof Brian Ripley
On Sun, 29 Jun 2003, Tony Marlboro wrote: Hello, I am using R on Mac OS X. I have tried to install the package netCDF, but have run into a problem. The install.packages command fails during a call to the configure shell script in the package build directory, because it cannot find

Re: [R] Fitting inter-arrival time data

2003-06-30 Thread Prof Brian Ripley
On Sun, 29 Jun 2003, M. Edward Borasky wrote: I have a collection of data which includes inter-arrival times of requests to a server. What I've done so far with it is use sm.density to explore the distribution, which found two large peaks. However, the peaks are made up of Gaussians, and

Re: [R] Returning contour co-ordinates

2003-06-30 Thread Ted Harding
On 29-Jun-03 Paul Murrell wrote: There is a partial implementation of this available via a tiny R package (http://www.stat.auckland.ac.nz/~paul/R/clines_1.0.tar.gz) Excellent! Thanks. [...] This will become available via a call to contour() in the future. Even better! (Really tiny: I like

Re: [R] Fitting inter-arrival time data

2003-06-30 Thread Adelchi Azzalini
On Monday 30 June 2003 01:23, M. Edward Borasky wrote: I have a collection of data which includes inter-arrival times of requests to a server. What I've done so far with it is use sm.density to explore the distribution, which found two large peaks. However, the peaks are made up of Gaussians,

[R] R as COM client (rather than server)

2003-06-30 Thread DED (David George Edwards)
Hello all Can R be used as a COM client (rather than as a COM server)? The following (Splus) code shows the sort of thing I'd like to do: pMIM - create.ole.object(mim31.Server) NoOutputLines - call.ole.method(pMIM, SendCmdLine, show w) for (i in 1:NoOutputLines)

[R] symbol size on a plot

2003-06-30 Thread Philippe Hupé
Hi, I would like to get from a plot the size of the symbols plotted. Imagine I have the following plot function : plot(1:2,1:2, pch=15, cex=4) I would like the get the values SIZE1 and SIZE2 so that if I plot the following rectangle : rect(1.5,1.5, 1.5+SIZE1, 1.5+SIZE2) then the size of this

[R] Rmpi: send a function to a slave from master and let it run?

2003-06-30 Thread lun li
Hello All, Can send different functions to different slaves from master and let each slave run itself function using Rmpi? From the Rmpi version 0.4.2, I know how to send a function to all slaves and let slave run the function. mpi.send.Robj may send a function to a slave, but how let the

Re: [R] repeatedly applying function with matrix-rows as argument

2003-06-30 Thread Prof Brian Ripley
apply(mymatrix, 1, function(x) myfun(x[1], x[2], x[3])) On Mon, 30 Jun 2003, Maarten Speekenbrink wrote: Suppose I have a function which takes three arguments. I would like to repeatedly apply the function, using a matrix N*3 in which each row supplies the three argements for the function. Is

[R] french map

2003-06-30 Thread Peyrard Nathalie
Hello, I would like to know if (as the usa map with Splus), it is possible with R to plot the french map and to add points (representing towns for instance) on the figure in the appropriate (x,y) system. Thank you. Nathalie Peyrard __ [EMAIL

[R] R CMD check

2003-06-30 Thread John Bjørnar Bremnes
when using R CMD check mypkg I get the error message ... * checking R files for library.dynam ... OK * checking generic/method consistency ... WARNING Error in .loadPackageQuietly(package, lib.loc) : Error in parse(file, n, text, prompt) : syntax error on line 95 Execution halted *

Re: [R] R as COM client (rather than server)

2003-06-30 Thread Duncan Temple Lang
There is also the Omegahat collection of packages for DCOM which supports R as a client, user-definable server objects and event handlers. Together, RDCOMClient and RDCOMServer (and the SWinRegistry and SWinTypeLibs packages) provide a general way to handle arbitrary data types in R and COM.

RE: [R] R as COM client (rather than server)

2003-06-30 Thread DED (David George Edwards)
Hi Duncan Thanks for the information, this looks very nice. Best regards David PS Some of the links on http://www.omegahat.org/RDCOMClient are broken. -Original Message- From: Duncan Temple Lang [mailto:[EMAIL PROTECTED] Sent: 30 June 2003 15:58 To: Erich Neuwirth Cc: DED (David

Re: [R] ./configure needs /sw/lib?

2003-06-30 Thread Prof Brian Ripley
On Mon, 30 Jun 2003, Tony Marlboro wrote: Thank you for your reply, Professor Ripley. Yes, to use R CMD INSTALL directly: use R CMD INSTALL --help to find out how. You need to supply --with-netCDF=/sw, I believe (but then the headers may not be found). R CMD INSTALL

Re: [R] french map

2003-06-30 Thread Barry Rowlingson
Peyrard Nathalie wrote: I would like to know if (as the usa map with Splus), it is possible with R to plot the french map and to add points (representing towns for instance) on the figure in the appropriate (x,y) system. You can use my Rmap library

[R] Your Message to Bonzi.com has an Attachment [T2003063000K3]

2003-06-30 Thread TechSupport
=== IMPORTANT! The message that you have sent to BONZI Software has a file attached to it. Due to the high level of virus messages we are currently receiving, your message has been discarded. Please reply to this message, make sure that there are no attachments to the message,

Re: [R] ./configure needs /sw/lib?

2003-06-30 Thread Tony Marlboro
Thank you for your reply, Professor Ripley. Yes, to use R CMD INSTALL directly: use R CMD INSTALL --help to find out how. You need to supply --with-netCDF=/sw, I believe (but then the headers may not be found). R CMD INSTALL --configure-args=--with-netCDF=/sw netCDF_1.5.tar.gz was the

Re: [R] french map

2003-06-30 Thread Spencer Graves
Have you tried www.r-project.org - search - R Site Search? There are several geographers who follow this list. From the above mentioned search, I found an exchange with subject Mapping in R discussing such a thing dated mercredi 21 novembre 2001 17:21, from: Christophe DECLERCQ, MD

Re: [R] R CMD check

2003-06-30 Thread Prof Brian Ripley
You have an error in the R files of your package. Before you even do R CMD check, do try loading the package in R. On Mon, 30 Jun 2003, John Bjørnar Bremnes wrote: when using R CMD check mypkg I get the error message ... * checking R files for library.dynam ... OK * checking generic/method

[R] spatial correlation test

2003-06-30 Thread Martin Wegmann
hello, I want to do a test for spatial correlation. I tried it with geary.test() but I don't understand the required input. x= a numeric vector the same length as the neighbours list in listw (my sampled data, I assume) listw= a listw object created for example by nb2listw (well when I check

[R] if else statements in R

2003-06-30 Thread Michael Rennie
Hi, there I am a grad student struggling to get my syntax right in a model I am building in R. I am trying to specify a variable,W, in a loop, such that on the first iteration, it is equal to Wo, my starting value, which I have already specified with Wo-9.2 On any other iteration other

[R] sample from a list of names for random groups

2003-06-30 Thread Daniele Medri
Dear R-users, I need to sample from a list of names (people) and create 2 random groups with unique item. Sure sample() function is a good point to start but I can't find example to solve my work. any tips? example? Thank -- Daniele Medri __

[R] if else statements in R

2003-06-30 Thread Michael Rennie
Hi, there I am a grad student struggling to get my syntax right in a model I am building in R. I am trying to specify a variable,W, in a loop, such that on the first iteration, it is equal to Wo, my starting value, which I have already specified with Wo-9.2 On any other iteration other

[R] Novice Questions

2003-06-30 Thread Bruce Moore
I'm writing a program to perform linear regressions to estimate the number of bank teller transactions per hour of various types based upon day of week, time of day, week of month and several prices. I've got about 25,000 records in my dataset, 85 columns of transaction counts (used 1 at a time),

Re: [R] if else statements in R

2003-06-30 Thread Timur Elzhov
On Mon, Jun 30, 2003 at 11:21:52AM -0400, Michael Rennie wrote: W- function(w) { if (comp[i,1]=1) W[i]-Wo else == is comparison operator, but = is assignment one. This code also doesn't work: x - 0 if (x = 1) Error: syntax error But this works: if (x == 1) + Chack also help(==), help(=)

Re: [R] NLME Covariates

2003-06-30 Thread Douglas Bates
Harold Doran [EMAIL PROTECTED] writes: In HLM, one can specify a covariate at one of the levels. For example, if the data structure are repeated observations nested within students nested within schools, school size might be a covariate that is used at level 3, but not at the other levels. In

Re: [R] if else statements in R

2003-06-30 Thread Andrew Perrin
This is untested, so try it first, but my approach would be: W- function(w) { W[i]-ifelse(comp[i,1]==1, Wo, (work[i-1,5]work[i-1,13]*/Ef)) } Best, Andy Perrin -- Andrew J Perrin - http://www.unc.edu/~aperrin Assistant Professor

RE: [R] french map

2003-06-30 Thread Christophe Declercq
Nathalie, You could have a look at the web site of the Institut Géographique National (http://www.ign.fr) where you can download a free SHAPEFILE with polygon data for the French départements. You can choose one of two projections: LAMBERT II étendu or LAMBERT 93

Re: [R] spatial correlation test

2003-06-30 Thread Stephane DRAY
At 18:48 30/06/2003, Martin Wegmann wrote: hello, I want to do a test for spatial correlation. I tried it with geary.test() but I don't understand the required input. x= a numeric vector the same length as the neighbours list in listw (my sampled data, I assume) listw= a listw object created for

Re: [R] within group variance of the coeficients in LME

2003-06-30 Thread Douglas Bates
J.R. Lockwood [EMAIL PROTECTED] writes: Dear listers, I can't find the variance or se of the coefficients in a multilevel model using lme. The component of an lme() object called apVar provides the estimated asymptotic covariance matrix of a particular

[R] Seemingly unrelated regressions

2003-06-30 Thread Richard A. Hinrichsen
Is there an implementation of SUR in R? These are used commonly in economics. Thanks -- Richard A. Hinrichsen __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] sample from a list of names for random groups

2003-06-30 Thread Spencer Graves
n1 - 2 n2 - 3 (S - sample(letters, n1+n2)) [1] f a d n l (S1 - sample(S, n1)) [1] l n (S2 - S[!is.element(S, S1)]) [1] f a d hope this helps. spencer graves Daniele Medri wrote: Dear R-users, I need to sample from a list of names (people) and create 2 random groups with unique item. Sure

Re: [R] if else statements in R

2003-06-30 Thread kjetil brinchmann halvorsen
On 30 Jun 2003 at 11:24, Michael Rennie wrote: Hola! Hi, there I am a grad student struggling to get my syntax right in a model I am building in R. I am trying to specify a variable,W, in a loop, such that on the first iteration, it is equal to Wo, my starting value, which I have

Re: [R] degrees of freedom in a LME model

2003-06-30 Thread Douglas Bates
Federico Calboli [EMAIL PROTECTED] writes: Dear All, I am analysing some data for a colleague (not my data, gotta be published so I cannot divulge). My response variable is the number of matings observed per day for some fruitlies. My factors are: Day: the observations were taken on

[R] make 1.7.1 crashes in Mac OS X

2003-06-30 Thread Michael Herron
I trying to compile R-.1.7.1 under OS X. Configure worked fine; then make crashes with: /sw/lib/libg2c.a(err.o) definition of common _f__formatted (size 4) /sw/lib/libg2c.a(err.o) definition of common _f__hiwater (size 4) /sw/lib/libg2c.a(err.o) definition of common _f__putn (size 4)

Re: [R] make 1.7.1 crashes in Mac OS X

2003-06-30 Thread Prof Brian Ripley
This is described in the MacOS X section of the R-admin manual. The file INSTALL does ask you to read it: please don't expect R-help readers to do so for you. On Mon, 30 Jun 2003, Michael Herron wrote: I trying to compile R-.1.7.1 under OS X. Configure worked fine; then make crashes with:

Re: [R] dropping factor levels in subset

2003-06-30 Thread Douglas Bates
Adaikalavan Ramasamy [EMAIL PROTECTED] writes: Another option is pruneLevels() in library nlme. x - factor( c( 0,1,2,1,2 ) ) x [1] 0 1 2 1 2 Levels: 0 1 2 pruneLevels( x[-1] ) [1] 1 2 1 2 Levels: 1 2 That function has been removed from the latest release of the nlme package because

Re: [R] spatial correlation test

2003-06-30 Thread Martin Wegmann
On Monday 30 June 2003 15:23, Barry Rowlingson wrote: Think you may be looking at the wrong sort of spatial correlation! For Geary tests you are comparing 'adjacent' objects, where adjacency is defined however you want - N-nearest neighbours, shared border between regions etc etc. sorry, I

[R] pairs and points

2003-06-30 Thread Christian Hennig
Hi, I am plotting a matrix/data frame with 10 variables by pairs. Is it possible to add another 10-dimensional point afterwards by something like points does for two-dimensional plots? Best, Christian -- *** Christian Hennig

Re: [R] degrees of freedom in a LME model

2003-06-30 Thread Federico Calboli
Dear Prof. Bates, Thank you for your reply. I did actually check the book at pages 89-92, but I have to say I found it a bit cryptic, if not downright confusing, for a genetist like me. Any day of the week. To use the example in the book, I cannot see why MACHINE is inner to WORKER. Or where

Re: [R] spatial correlation test

2003-06-30 Thread Roger Bivand
(I think Barry Rowlingson replied off-list, both his advice and Stéphane Dray's reply are relevant) On Mon, 30 Jun 2003, Martin Wegmann wrote: On Monday 30 June 2003 15:23, Barry Rowlingson wrote: Think you may be looking at the wrong sort of spatial correlation! For Geary tests you are

[R] date to age

2003-06-30 Thread William Briggs
I have files which have columns of data that look like this: DOB 9/27/1964 ... That is, dates in month/day/year format. When variables like DOB are read in, they are converted to factors. So, I wrote this to convert from date to age in years: age-function(x)

Re: [R] date to age

2003-06-30 Thread Prof Brian Ripley
On Mon, 30 Jun 2003, William Briggs wrote: I have files which have columns of data that look like this: DOB 9/27/1964 ... That is, dates in month/day/year format. When variables like DOB are read in, they are converted to factors. Why? You can avoid that, of course. So, I wrote

Re: [R] Seemingly unrelated regressions

2003-06-30 Thread SATO Masaei
Is there an implementation of SUR in R? These are used commonly in economics. Thanks -- Richard A. Hinrichsen You can use the function 'sur' in the library 'systemfit'. Masaei Sato __ [EMAIL PROTECTED] mailing list

[R] Deciphering an error message

2003-06-30 Thread Suzanne E. Blatt
Hello, I'm working in spatstat and having difficulty with the ppp objects. I can get a ppp object for one set of data, but with the same code applied to a second data set (all I am changing is the field identifier), I get the following error message: Error in switch(w$type, rectangle={:

Re: [R] Novice Questions

2003-06-30 Thread Spencer Graves
Hi, Bruce: I'm overwhelmed. Questioners seem to get quicker, more informative responses to questions that are short, well written and easily and quickly understood, and preferably include a toy example that someone else can run quickly to reproduce the problem and to evaluate alternative

[R] dev.copy()

2003-06-30 Thread David Brahm
In a batch script (i.e. there is no screen device), I would like to create both a PDF file and a PNG file from the same plot. I thought this would be the way to do it, but it isn't: R bitmap(copy.png, png16m, res=300) R pdf(copy.pdf) R plot(1:10, 1:10) # Plot into pdf file R

[R] Constructing loops in which i+1st element depends on ith

2003-06-30 Thread Michael Rennie
I feel greedy posting for help twice in one day- please forgive me, but the thesis can't wait. I have been trying to get an “if…else” statement working in a loop I am writing, in which I specify a variable value to 9.2 (Wo) on the first iteration, but for subsequent iterations, to adopt a

Re: [R] sample from a list of names for random groups

2003-06-30 Thread kjetil brinchmann halvorsen
On 30 Jun 2003 at 17:11, Daniele Medri wrote: You should really try to specify better what you want. Lets try names # character vector with your names m - length(names) n1 - 5 n2 - 10 muestra - sample(names, m, replace=FALSE) muestra1 - muestra[1:n1] muestra2 - muestra[(n1+1):(n1+n2)] Is this

[R] crossed random effects

2003-06-30 Thread Sarah Mclean
Hi, I have a data set on germination and plant growth with the following variables: dataset=fm mass (response) sub (fixed effect) moist (fixed effect) pop (fixed effect) mum (random effect nested within population) iheight (covariate) plot (random effect- whole plot factor for split-plot

[R] Creating a loop that works....

2003-06-30 Thread Michael Rennie
Hi there, First off, thanks to everyone who has helped me so far. Sorry to keep pestering you all. I'm including my code here, and I will comment down it where it is that I am having problems figuring out how to write this damn thing. temper - scan(temp2.dat, na.strings = ., list(Day=0,

RE: [R] Fitting inter-arrival time data

2003-06-30 Thread M. Edward Borasky
Unfortunately, the data are *non-negative*, not strictly positive. Zero is a valid and frequent inter-arrival time. It is, IIRC, the most likely value of a (negative) exponential distribution. -- M. Edward (Ed) Borasky mailto:[EMAIL PROTECTED] http://www.borasky-research.net Suppose that

RE: [R] Fitting inter-arrival time data

2003-06-30 Thread M. Edward Borasky
Thanks!! It does look like the easiest thing is direct ML; the code for a normal mixture is in the book, so all I have to do is modify that for a sum of a hyper-exponential, for which I have an approximate mean and CV, and a normal, for which I have an approximate mean and SD. I have two big

[R] crossed random effects

2003-06-30 Thread Sarah Mclean
Hi, if I have posted this twice, please ignore this. I'm not sure if I sent it to the correct e-mail address the first time. I have a data set on germination and plant growth with the following variables: dataset=fm mass (response) sub (fixed effect) moist (fixed effect) pop (fixed effect) mum