[R] lattice: cumsum and xyplot

2004-06-11 Thread Wolfram Fischer
I want to display cumulative summary functions with lattice. First I tried to get cumulated data: library(lattice) data(barley) d.cum - with( barley, by( yield, INDICES=list(site=site,year=year), FUN=cumsum ) ) I got a list of vectors. I tried to get a dataframe which I could use in

Re: [R] lattice: cumsum and xyplot

2004-06-11 Thread Wolski
Hi! To get a data.frame as.data.frame(do.call(rbind,d.cum)) Sincerely Eryk *** REPLY SEPARATOR *** On 6/11/2004 at 9:17 AM Wolfram Fischer wrote: I want to display cumulative summary functions with lattice. First I tried to get cumulated data: library(lattice)

[R] rownames of single row matrices

2004-06-11 Thread Robin Hankin
Hi I want to extract rows of a matrix, and preserve rownames even if only one row is selected. Toy example: R a - matrix(1:9,3,3) R rownames(a) - letters[1:3] R colnames(a) - LETTERS[1:3] R a A B C a 1 4 7 b 2 5 8 c 3 6 9 Extract the first two rows: R wanted - 1:2 R a[wanted,] A B C a 1 4 7 b

Re: [R] rownames of single row matrices

2004-06-11 Thread Eric Lecoutre
Hi Robin Have a look at: help([) The fact that dimensions are lost when extracting is a feature of the language. What you need is the drop option. a[1,,drop=FALSE] A B C a 1 4 7 Eric At 10:09 11/06/2004, Robin Hankin wrote: Hi I want to extract rows of a matrix, and preserve rownames even if

Re: [R] rownames of single row matrices

2004-06-11 Thread Jason Turner
On Fri, 2004-06-11 at 20:09, Robin Hankin wrote: Hi I want to extract rows of a matrix, and preserve rownames even if only one row is selected. The infamous drop=TRUE default. help([) goes into this. Your toy example, two ways: a - matrix(1:9,3,3) rownames(a) - letters[1:3]

Re: [R] Questions about Preserving registers

2004-06-11 Thread David Lennartsson
Duncan Murdoch wrote: On Thu, 10 Jun 2004 09:12:05 -0700 (PDT), Thomas Lumley [EMAIL PROTECTED] wrote : DLL attempted to change FPU control word from 9001f to 90003 I read the instruction on Duncan Murdochs website about preserving registers, but I still dont understand it. For example,

[R] setValidity() changes Extends?

2004-06-11 Thread Matthias Kohl
Hi, I'm using Version 1.9.0 (2004-04-12) on Windows NT/98/2000 and found the following difference between using setClass(..., valdity = ), respectively using setValidity() afterwards: setValidity() changes the Extends-part of a derived class, is this intended or a bug or am I missing

Re: [R] rownames of single row matrices

2004-06-11 Thread Robin Hankin
Go list! The examples on help([) were all I needed. The examples don't explicitly say that rownames are preserved; perhaps this is obvious. It'd be nice if matrix m on the help page had rownames. best wishes rksh -- Robin Hankin Uncertainty Analyst Southampton Oceanography Centre SO14 3ZH tel

Re: [R] running R UNIX in a mac computer

2004-06-11 Thread Jari Oksanen
On Fri, 2004-06-11 at 03:49, Tiago R Magalhaes wrote: Hi to you all My question is: there is a package written in UNIX for which there is no Mac version. I would like to know if it's possible to install the R UNIX version on the MacOSX and run that UNIX package on my Mac (through this

Re: [R] lattice: cumsum and xyplot

2004-06-11 Thread Wolfram Fischer
--- In reply to: --- Date:11.06.04 09:33 (+0200) From:Wolski [EMAIL PROTECTED] Subject: Re: [R] lattice: cumsum and xyplot Hi! To get a data.frame as.data.frame(do.call(rbind,d.cum)) Thanks for this interesting hint. To get the lost names of factors again, I tried now the following

[R] comparing regression slopes

2004-06-11 Thread Nathan Weisz
Dear List, I used rlm to calculate two regression models for two data sets (rlm due to two outlying values in one of the data sets). Now I want to compare the two regression slopes. I came across some R-code of Spencer Graves in reply to a similar problem: http://www.mail-archive.com/[EMAIL

Re: [R] comparing regression slopes

2004-06-11 Thread Prof Brian Ripley
It's not OK. Treat the results from rlm as having infinite df since the theory is asymptotic (use pnorm not pt), and don't expect anything useful in samples as small as 10 cases (30+, preferably 100+ would be OK). On Fri, 11 Jun 2004, Nathan Weisz wrote: Dear List, I used rlm to calculate

Re: [R] comparing regression slopes

2004-06-11 Thread Spencer Graves
Fortunately, you don't have to accept no for an answer: You can dream up something that you think is sensible, simulate scenarios you want to detect plus scenarios with no difference, and find out what your procedure produces, the distributions of your sample statistics, etc., e.g., as

[R] Continuous-lag Markov chains?

2004-06-11 Thread Monica Palaseanu-Lovejoy
Hi again, I am wondering if there is any R package which does continuous- lag Markov chains modelling? I discovered MCMCpack but i am not sure it does continuous-lags. Thanks for your answers, Monica Manchester University __ [EMAIL PROTECTED]

[R] set user-agent for download

2004-06-11 Thread Hans Kestler
Hi, does anybody know if it is possible to set a user-agent for download using a proxy? Thank you. Hans Kestler__ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] how to install yags in R

2004-06-11 Thread Tu Yu-Kang
Dear R users, I search the R archives and noted that the same problem has been posted but without solution. I know there IS instructions by the author of yags, but I just couldn't figured out. I know gee and geepack can also perform generalized estimating equation, but the reason why I need

Re: [R] set user-agent for download

2004-06-11 Thread Uwe Ligges
Hans Kestler wrote: Hi, does anybody know if it is possible to set a user-agent for download using a proxy? See ?download.file for details. Uwe Ligges Thank you. Hans Kestler __

Re: [R] Informal discussion group about R

2004-06-11 Thread Paul Lemmens
Hoi Harold, --On donderdag 10 juni 2004 7:55 -0700 Baize, Harold [EMAIL PROTECTED] wrote: I've started a tribe for discussing R and sharing scripts. Tribe.net is one of the popular I hope it will be of help to newbies, although I'm new to R myself. Here is the url:

Re: [R] Continuous-lag Markov chains?

2004-06-11 Thread Chris Jackson
Monica Palaseanu-Lovejoy wrote: Hi again, I am wondering if there is any R package which does continuous- lag Markov chains modelling? I discovered MCMCpack but i am not sure it does continuous-lags. MCMCpack fits Bayesian models using Markov Chain Monte Carlo simulation, which is a very

[R] [RSPerl] Loading problem

2004-06-11 Thread Emmanuel Engelhart
Hi I'm a new user with R, and I try to use it from Perl. So I have to use RSPerl. I tried to load it with library(RSPerl) : library(RSPerl) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library /usr/local/lib/R/library/RSPerl/libs/RSPerl.so:

[R] Sweave and multiple graphs

2004-06-11 Thread Gavin Simpson
Dear list, I am using Sweave to build a small report. I want to produce a series of figures, each figure containing a number of plots and then have them included in the Sweave file. An example would be to : postscript(file = ANCbwplot%03d.eps, onefile = FALSE, other options...) oldpar -

Re: [R] Sweave and multiple graphs

2004-06-11 Thread tobias . verbeke
[EMAIL PROTECTED] wrote on 11/06/2004 14:07:18: Dear list, I am using Sweave to build a small report. I want to produce a series of figures, each figure containing a number of plots and then have them included in the Sweave file. An example would be to : postscript(file =

[R] interaction plot with error bars based on TukeyHSD intervals

2004-06-11 Thread Manuel López-Ibáñez
Hello, maybe the previous message was not very clear. Thus, I will try to explain the problem in a better way. I would like to have an interaction plot with error bars based on TukeyHSD intervals. In an experiment with two factors A and B, each of them with two levels A={a1,a2} and B={b1,b2},

[R] Modifying Code in .rda based packages (e.g. lme4)

2004-06-11 Thread Dieter Menne
Dear List, assume I want to make a minor local change in a package that is supplied as .rda. For example, I want to get rid of the non-verbose-protected Iteration message in GLMM/lme4. Probably I have to load / change / save the package, but could someone help me to get the syntax right? Dieter

Re: [R] Modifying Code in .rda based packages (e.g. lme4)

2004-06-11 Thread Prof Brian Ripley
No standard R package is supplied as a .rda, including not lme4. You must be looking at a binary installation, and you would do best to reinstall from the sources. You could use R --vanilla load(/all.rda) fix(GLMM) save(ls(all=T), file=/all.rda, compress = TRUE) q() but we would not

Re: [R] Modifying Code in .rda based packages (e.g. lme4)

2004-06-11 Thread Roger D. Peng
You should probably download the source for the package from CRAN, modify the source, then reinstall. -roger Dieter Menne wrote: Dear List, assume I want to make a minor local change in a package that is supplied as .rda. For example, I want to get rid of the non-verbose-protected Iteration

[R] own family for glm()

2004-06-11 Thread I M S White
Is it possible in R to create a customized family to be used with glm()? I see no mention of this possibility in the documentation. (S-plus has something called make.family). == I.White ICAPB, University of Edinburgh Ashworth Laboratories, West Mains Road

[R] dll file missing?

2004-06-11 Thread Kemp S E (Comp)
Hi, I am trying to do a dyn.load(), but I get the following error... dyn.load(fileGT.dll) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library C:/R_Files/fileGT.dll: LoadLibrary failure: The specified module could not be found. It states it

Re: [R] own family for glm()

2004-06-11 Thread Prof Brian Ripley
On Fri, 11 Jun 2004, I M S White wrote: Is it possible in R to create a customized family to be used with glm()? I see no mention of this possibility in the documentation. (S-plus has something called make.family). It is possible. Package MASS has two examples, the simpler of which is

Re: [R] dll file missing?

2004-06-11 Thread Prof Brian Ripley
On Fri, 11 Jun 2004, Kemp S E (Comp) wrote: Hi, I am trying to do a dyn.load(), but I get the following error... dyn.load(fileGT.dll) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library C:/R_Files/fileGT.dll: LoadLibrary failure: The

[R] Einlesen von Daten unter R

2004-06-11 Thread Oscar Oehler
Sehr geehrte Damen und Herren, An der Fachhochschule Winterthur wurde für mich im Rahmen einer Diplomarbeit (Prof. Ruckstuhl) ein R-Programm zur Auswertung von IR-Spektren (line-shape-Analyse zur quantitativen Analyse von Gasmischungen) entwickelt. Die Daten werden einer Excel-Tabelle

Re: [R] Modifying Code in .rda based packages (e.g. lme4)

2004-06-11 Thread Douglas Bates
Dieter Menne [EMAIL PROTECTED] writes: assume I want to make a minor local change in a package that is supplied as .rda. For example, I want to get rid of the non-verbose-protected Iteration message in GLMM/lme4. Probably I have to load / change / save the package, but could someone help

[R] [R-pkgs] New package: RUnit

2004-06-11 Thread Klaus Juenemann
We would like to announce the availability on CRAN of a new package: RUnit It contains a unit testing framework strongly inspired by Javas popular JUint package. In addition it contains some functionality to investigate the degree to which some function is covered by a test suite. The main aims

[R] ROC for threshold value, biometrics

2004-06-11 Thread myint . tin-tin-htar
Hello, I am just a beginner of R 1.9.0. I try to construct a predictive score for the development of liver cancer in cirrhotic patients. So dependant variable is binanry (cancer yes or no). Independant variables are biological data. The aim is to find out a cut-off value which differentiate

[R] *** NEWBIE QUESTION *** QUANTILE FUNCTIONS

2004-06-11 Thread Pier Luca Lanzi
Dear all, sorry this will sound as naive as it can be. I need to know whether there is a closed analytical form (even approximated) to the quantile function for the Binomial? and for the Poisson? if not, what is the best citation to use when stating this? Thank you, Pier Luca

Re: [R] Sweave and multiple graphs

2004-06-11 Thread Gavin Simpson
[EMAIL PROTECTED] wrote: Use list.files(). It has a path argument (to specify the directory) and a pattern argument to put the regular expression. mygraphs - list.files(path=./mygraphs, pattern=^ANCbwplot.*\\.eps) See ?list.files and maybe ?regex HTH, Tobias Hi Tobias. Thanks for the reply. I just

Re: [R] *** NEWBIE QUESTION *** QUANTILE FUNCTIONS

2004-06-11 Thread Sundar Dorai-Raj
Pier Luca Lanzi wrote: Dear all, sorry this will sound as naive as it can be. I need to know whether there is a closed analytical form (even approximated) to the quantile function for the Binomial? and for the Poisson? if not, what is the best citation to use when stating this? Thank you,

Re: [R] Einlesen von Daten unter R

2004-06-11 Thread Martin Wegmann
hello, sorry, I cannot help you but if you post it again in English (because it is an English speaking mailing list) your are likely to receive a helpful reply. Cheers Martin On Friday 11 June 2004 15:52, Oscar Oehler wrote: Sehr geehrte Damen und Herren, An der Fachhochschule Winterthur

[R] par specification inside formula vs. outside function call

2004-06-11 Thread Peter Flom
The other day I asked a question about changing the size of axis labels in a dotchart. Chuck Cleland provided the answer, which was to specify par(cex.axis = .8) outside the call to dotchart I had tried this inside the formula, and it did nothing (nor did it produce an error). I was

Re: [R] Einlesen von Daten unter R

2004-06-11 Thread Gabor Grothendieck
My German is not very good but I assume your question is how to interface Labview and R. You might want to have a look at the two R COM interfaces that are available: http://cran.r-project.org/other-software.html http://mailman.csd.univie.ac.at/mailman/listinfo/rcom-l

[R] My useR! slides: Doing Customer Intelligence with R

2004-06-11 Thread Jim Porzak
FYI for useR! attendees others. The slides for my talk Doing Customer Intelligence with R are up on our R weblog: R.LoyaltyMatrix.com Also note any books ordered through the links on the blog will benefit the R Foundation - see blog for details. I will be using the blog as an informal log of

[R] question about Rcmd SHLIB

2004-06-11 Thread Laura Holt
Dear R People: I'm trying to use the Rcmd SHLIB to produce a dll. Rcmd SHLIB -o test2.dll test2.f make[1]: *** [libR.a] Error 255 make: *** [libR] Error 2 Where do I go to find out about the make errors, please? I suspect that I might be missing something. I have the tools for creating new

RE: [R] par specification inside formula vs. outside function cal l

2004-06-11 Thread Liaw, Andy
From: Peter Flom The other day I asked a question about changing the size of axis labels in a dotchart. Chuck Cleland provided the answer, which was to specify par(cex.axis = .8) outside the call to dotchart I had tried this inside the formula, and it did nothing (nor did it

Re: [R] question about Rcmd SHLIB

2004-06-11 Thread Prof Brian Ripley
See readme.packages. It would help to run make libR.a first. On Fri, 11 Jun 2004, Laura Holt wrote: Dear R People: I'm trying to use the Rcmd SHLIB to produce a dll. Rcmd SHLIB -o test2.dll test2.f make[1]: *** [libR.a] Error 255 make: *** [libR] Error 2 Where do I go to find out

[R] lme newbie question

2004-06-11 Thread Christoph Lehmann
Hi I try to implement a simple 2-factorial repeated-measure anova in the lme framework and would be grateful for a short feedback -my dependent var is a reaction-time (rt), -as dependent var I have -the age-group (0/1) the subject belongs to (so this is a between-subject factor), and

Re: [R] question about Rcmd SHLIB

2004-06-11 Thread Peter Dalgaard
Laura Holt [EMAIL PROTECTED] writes: Dear R People: I'm trying to use the Rcmd SHLIB to produce a dll. Rcmd SHLIB -o test2.dll test2.f make[1]: *** [libR.a] Error 255 make: *** [libR] Error 2 Where do I go to find out about the make errors, please? I suspect that I might be missing

[R] [StatDataML] compile error

2004-06-11 Thread Patrick Drechsler
Hi, sorry if this is OT here. I'm trying to compile StatDataML[1] to import Matlab data into R. The R part works fine but compiling the Matlab part exits with this error message (following the INSTALL instructions: autoconf - ./configure - make): ,[ [EMAIL

Re: [R] lme newbie question

2004-06-11 Thread JRG
On 11 Jun 04, at 20:12, Christoph Lehmann wrote: Hi I try to implement a simple 2-factorial repeated-measure anova in the lme framework and would be grateful for a short feedback -my dependent var is a reaction-time (rt), -as dependent var I have -the age-group (0/1) the subject

[R] space-time k function problem

2004-06-11 Thread hkang
I am intersted in time-space clustering of local industry. I made the point file, polygons and time table then run 'stkhat' function in Splancs, but it generates only time k function. I noticed that it has only one ks value. ... Can anybody help me?. thanks in advance. countiespt -

Re: [R] space-time k function problem

2004-06-11 Thread Roger Bivand
On Fri, 11 Jun 2004, hkang wrote: I am intersted in time-space clustering of local industry. I made the point file, polygons and time table then run 'stkhat' function in Splancs, but it generates only time k function. I noticed that it has only one ks value. ... Can anybody help

[R] probabilistic neural networks

2004-06-11 Thread Rajarshi Guha
Hi, I'm working on a classification problem and one of the methods I'd like to use are neural networks. I've been using nnet to build a classification network. However I would like to have the probabilities associated with the prediction. Are there any implementations of probabilistic neural

Re: [R] probabilistic neural networks

2004-06-11 Thread Prof Brian Ripley
That's not what a `probabilistic neural network' is. However, nnet already does what you ask: do read the references as the posting guide asks. On Fri, 11 Jun 2004, Rajarshi Guha wrote: I'm working on a classification problem and one of the methods I'd like to use are neural networks. I've

[R] Error when I try to build / plot a tree using rpart()

2004-06-11 Thread Jude Ryan
Hi, I am using the rpart package to build a classification tree. I did manage to build a tree with data on a previous project. However, when attampting to build a tree on a project I am working on, I seem to be getting the error shown below: nhg3.rp - rpart(profitresp ~., nhg3, method=class)

Re: [R] ROC for threshold value, biometrics

2004-06-11 Thread Frank E Harrell Jr
[EMAIL PROTECTED] wrote: Hello, I am just a beginner of R 1.9.0. I try to construct a predictive score for the development of liver cancer in cirrhotic patients. So dependant variable is binanry (cancer yes or no). Independant variables are biological data. The aim is to find out a cut-off value

[R] memory allocation and interrupts

2004-06-11 Thread Vadim Ogranovich
Hi, A recent discussion on the list about tryCatch and signals made me think about memory allocation and signals in C extension modules. What happens to the memory allocated by R_alloc and Calloc if the user pressed Ctr-C during the call? R-ext doesn't seem to discuss this. I'd guess that

[R] Re: R-help Digest, Vol 16, Issue 11

2004-06-11 Thread Alan Cobo-Lewis
[EMAIL PROTECTED] writes: I make a study in health econometrics and have a categorical dependent variable (take value 1-5). I would like to fit an ordered probit or ordered logit but i didn't find a command or package who make that. Does anyone know if it's exists ? Try polr() from the MASS

RE: [R] Error when I try to build / plot a tree using rpart()

2004-06-11 Thread Liaw, Andy
You didn't get a tree. The output of print() tells you that you only have the root node. You may need to adjust some of the parameters with rpart(..., control=rpart.control(...)). See ?rpart.control. HTH, Andy From: Jude Ryan Hi, I am using the rpart package to build a classification