Re: [R] R-etiquette

2005-01-10 Thread Anne
well, all these responses address part of my question! - what about citing people who helped with comments/ suggestions even code on the list? I do think they deserve aknowledgement! Thanks Anne - Original Message - From: Romain François [EMAIL PROTECTED] To: Gabor Grothendieck

[R] Votre message ne sera pas lu !

2005-01-10 Thread embauche
En raison de l'augmentation du pourriel (spam) sur cette adresse, nous ne pouvons plus consulter les messages reçus. Pour nous contacter merci de choisir le lien écrire au webmestre de notre site web ou d'utiliser le formulaire accessible via le lien ci-dessous

[R] Multiple comparisons following nlme

2005-01-10 Thread Stefan Uhlig
Dear Madam or Sir, I need to do multiple comparisons following nlme analysis (Compare the effects of different treatments on a response measured repeatedly over time; fixed = response ~ treat*time). On the web I found the notion that one might use the L argument from ANOVA. Do you have an

[R] Graphical table in R

2005-01-10 Thread Dan Bolser
I would like R to produce some tabulated data in a graphical output. When I say tabulated data, what I mean is a table with rows and columns. This would be useful when reading in a big file, performing some analysis on it, and then wanting to display the results as a table. Something like

[R] Invisible plot using RSvgDevice

2005-01-10 Thread Jan Hummel
Dear list members, I have a probably simple question concerning the RSvgDevice. After upgrading from R 1.9.0 to R 2.0.1 the computet svg files looking empty. Each time the RSvgDevice 0.5.3 were used. Scales and headers are printed but the plots are missing: rect x=433.10 y=246.13 width=0.93

[R] Partial wireframe plots

2005-01-10 Thread Rob Knell
Dear R-helpers Can anyone direct me to a method for plotting what you might call a partial wireframe plot? I have two explanatory variables in a dataset which give me a significant interaction term when I fit a model. The two variables are correlated with each other to a moderate degree, and

RE: [R] R-etiquette

2005-01-10 Thread Liaw, Andy
I'd thank them in the acknowledgement section. I think some (most?) journal will allow such a section, and most people use that to thank their research funding sources and/or collaborators who had not made the list of authors. Andy From: Anne well, all these responses address part of my

Re: [R] R-etiquette

2005-01-10 Thread Prof Brian Ripley
On Mon, 10 Jan 2005, Anne wrote: well, all these responses address part of my question! - what about citing people who helped with comments/ suggestions even code on the list? I do think they deserve aknowledgement! That depends on the journal and its conventions. Most commonly they are

Re: [R] Invisible plot using RSvgDevice

2005-01-10 Thread Prof Brian Ripley
Have you contacted the author? The conventions for R graphics devices changed in R 2.0.0, in particular the way the alpha channel is represented. I doubt if the same version of a graphics device works in both R 2.0.0 and R = 2.0.0, and this one appears to predate R 2.0.0. On Mon, 10 Jan

Re: [R] R-etiquette

2005-01-10 Thread Frank E Harrell Jr
Liaw, Andy wrote: I'd thank them in the acknowledgement section. I think some (most?) journal will allow such a section, and most people use that to thank their research funding sources and/or collaborators who had not made the list of authors. Andy You do need to get permission from people to

Re: [R] I have some problem about GLM function.

2005-01-10 Thread Prof Brian Ripley
On Sun, 9 Jan 2005, may sawasdee wrote: Dear R-Help I 'm using GLM function to Modelling. But when I used Gamma Family in GLM, then I can't run. It was error glm(DamageRatio~MinTEMP+MaxTEMP+DayRain+Group1+Group2+Group3+Year,family=Gamma()) Error in eval(expr, envir, enclos) : Non-positive values

Re: [R] Partial wireframe plots

2005-01-10 Thread Prof Brian Ripley
On Mon, 10 Jan 2005, Rob Knell wrote: Can anyone direct me to a method for plotting what you might call a partial wireframe plot? I have two explanatory variables in a dataset which give me a significant interaction term when I fit a model. The two variables are correlated with each other to a

Re: [R] R-etiquette

2005-01-10 Thread Duncan Murdoch
On Mon, 10 Jan 2005 10:03:33 +0100, Anne [EMAIL PROTECTED] wrote : well, all these responses address part of my question! - what about citing people who helped with comments/ suggestions even code on the list? I do think they deserve aknowledgement! Mention them in an acknowledgments section at

[R] contrasts involving random terms in lme

2005-01-10 Thread Joris DeWolf
Hello, Has somebody an idea of how to fit contrast involving random terms and obtain their standard errors with lme? I am referring to Welham, S.; Cullis, B.; Gogel, B.; Gilmour, A. and Thompson, R. (2004). Prediction in linear mixed models. Australian and New Zealand Journal of Statistics 46,

[R] doing many commands within R

2005-01-10 Thread Cserhti Mtys
Dear all, I'm new to this list, so let me greet everyone. My problem is that I have several thousand data files which I want to perform a lot of R commands on, which are found in a seperate .R script. Now, what I did was within the R prompt, I used to read in a list of the data files.

Re: [R] doing many commands within R

2005-01-10 Thread Uwe Ligges
Cserhti Mtys wrote: Dear all, I'm new to this list, so let me greet everyone. My problem is that I have several thousand data files which I want to perform a lot of R commands on, which are found in a seperate .R script. Now, what I did was within the R prompt, I used to read in a list of the

[R] mouse position with package tcltk

2005-01-10 Thread solares
HI, i need a know what is wrong with this code. If you run then not known, X0, W or XPlace. This code not accept parameter but in tcl yes. This code make a window panned. Thanks Ruben In R: library(tcltk) XPlace-function(fract){ tkcmd(place,nb1,relwidth=fract)

Re: [R] Graphical table in R

2005-01-10 Thread Marc Schwartz
On Mon, 2005-01-10 at 11:05 +, Dan Bolser wrote: I would like R to produce some tabulated data in a graphical output. When I say tabulated data, what I mean is a table with rows and columns. This would be useful when reading in a big file, performing some analysis on it, and then wanting

[R] Installation of XML library can't find libxml2.dll

2005-01-10 Thread michael watson \(IAH-C\)
Sorry to ask a (probably) dumb question, but I am trying to install XML package on Windows XP, R 2.0.1, and I get the error: This application has failed to start because libxml2.dll was not found. Re-installing the application may fix this problem library(XML) Error in dyn.load(x,

[R] Installing RCurl on Linux

2005-01-10 Thread michael watson \(IAH-C\)
Hi I'm not having much luck today! Running Suse Linux 8.2, R 2.0.0, trying to install RCurl_0.5-1.tar.gz. Perhaps there's an issue with my C compiler, but I got some errors: IAHC-LINUX03:/usr/users/mwatson # R CMD INSTALL RCurl_0.5-1.tar.gz * Installing *source* package 'RCurl' ... checking

Re: [R] Graphical table in R

2005-01-10 Thread Dan Bolser
On Mon, 10 Jan 2005, Marc Schwartz wrote: On Mon, 2005-01-10 at 11:05 +, Dan Bolser wrote: I would like R to produce some tabulated data in a graphical output. When I say tabulated data, what I mean is a table with rows and columns. This would be useful when reading in a big file,

[R] help diagnosing ftp problem

2005-01-10 Thread Liaw, Andy
Dear R-help, I've been compiling R from source on our Linux boxes for quite a while. One thing that bugs me is that I always get an error when make check-all where it choked up running tests/internet.R. That wasn't a big deal, as I can run install.packages()/update.packages() fine, and that's

Re: [R] how to do by-processing using weighted.mean

2005-01-10 Thread Thomas Lumley
On Sun, 9 Jan 2005, Denis Chabot wrote: Hi, I'd like to compute the weighted mean of some variables, all using the same weight variable, for each combination of 3 factor variables. I found how I could use summarize (from Hmisc) to do normal means for combinations of 3 factors, but I cannot find

Re: [R] plot.default and open ended limits

2005-01-10 Thread Thomas Lumley
On Sun, 9 Jan 2005, Dan Bolser wrote: Thats a good point. help.search(xlim) / ?xlim don't turn up any pages, and only plot.default seems to document [xy]lim. Aside from greping the source / docs how can I find out which plot functions use / abuse the [xy}lim parameter? That's actually good news.

Re: [R] doing many commands within R

2005-01-10 Thread Tim F Liao
Sounds like you've got a file referencing problem. You may want to try and input the file into R using a simple read.table command or something similar without yet running your long script to see if your file can be called up. Tim Original message Date: Mon, 10 Jan 2005 13:55:28 +0100

[R] Mixing portrait/landscape in a postscript file

2005-01-10 Thread Yves Brostaux
Dear list, I'm stuck with a little graphical problem. I'm generating several lattice plots which are printed in a single postcript device opened by trellis.device(postscript, theme=canonical.theme(postscript, color=F), file=an_phase2_graph.ps, paper=a4, pointsize = 10, onefile=TRUE,

RE: [R] Graphical table in R

2005-01-10 Thread Liaw, Andy
I recall there's a function in the gregmisc bundle that can write text to an R plot (was it textplot?). That might be of use here, I guess. Andy From: Dan Bolser On Mon, 10 Jan 2005, Marc Schwartz wrote: On Mon, 2005-01-10 at 11:05 +, Dan Bolser wrote: I would like R to produce

Re: [R] R-etiquette

2005-01-10 Thread Anne
In this case, it is not a publication for a journal (where I would follow the prescriptions of my academic institution, ie ETH) but an internal report for a local gvt agency. Looking at some such reports I was horrified to find out that credits were not made (in a case I do know that R was used,

Re: [R] Graphical table in R

2005-01-10 Thread Marc Schwartz
On Mon, 2005-01-10 at 14:52 +, Dan Bolser wrote: On Mon, 10 Jan 2005, Marc Schwartz wrote: On Mon, 2005-01-10 at 11:05 +, Dan Bolser wrote: I would like R to produce some tabulated data in a graphical output. When I say tabulated data, what I mean is a table with rows and columns.

Re: [R] help diagnosing ftp problem

2005-01-10 Thread Uwe Ligges
Liaw, Andy wrote: Dear R-help, I've been compiling R from source on our Linux boxes for quite a while. One thing that bugs me is that I always get an error when make check-all where it choked up running tests/internet.R. That wasn't a big deal, as I can run install.packages()/update.packages()

Re: [R] Installation of XML library can't find libxml2.dll

2005-01-10 Thread Prof Brian Ripley
Save yourself the problem and get a precompile package from http://www.stats.ox.ac.uk/pub/RWin On Mon, 10 Jan 2005, michael watson (IAH-C) wrote: Sorry to ask a (probably) dumb question, but I am trying to install XML package on Windows XP, R 2.0.1, and I get the error: This application has failed

Re: [R] Mixing portrait/landscape in a postscript file

2005-01-10 Thread Prof Brian Ripley
On Mon, 10 Jan 2005, Yves Brostaux wrote: Dear list, I'm stuck with a little graphical problem. I'm generating several lattice plots which are printed in a single postcript device opened by trellis.device(postscript, theme=canonical.theme(postscript, color=F), file=an_phase2_graph.ps, paper=a4,

[R] SSOAP and Rcurl with proxy server

2005-01-10 Thread michael watson \(IAH-C\)
Hi I'm trying to use a bioconductor package (KEGGSOAP) which relies on Rcurl and SSOAP. As an example, a function exists: list.organisms function () { orgs - matrix(unlist(.SOAP(KEGGserver, list_organisms, , action = KEGGaction, xmlns = KEGGxmlns), use.names = FALSE),

Re: [R] Graphical table in R

2005-01-10 Thread Anne
Well I saw the title graphical table... and had just a cursory glance However I do use this function which is good if not too many categories are represented. Anne - Original Message - From: Marc Schwartz [EMAIL PROTECTED] To: Anne [EMAIL PROTECTED] Sent: Monday, January 10, 2005 5:54

RE: [R] Installation of XML library can't find libxml2.dll

2005-01-10 Thread dr mike
My reading of this is that you may have installed the 'Windows binary' package from the Omegahat site - if so, it would appear that this is for versions of R 2.0 . Further to the pointer to the reply by Prof Ripley, hisa link provides the correct Windows binary for R = 2.0. Regards Mike

[R] Query: simple autocorrelation ina time series

2005-01-10 Thread Crabb, David
I hope you can help with what might be an easy query. I am doing some simple simulations where I am generating series of data from a normal distribution using rnorm. I am then treating these as a time series. I want to know how I can incorporate correlation in the series (autocorrelation) i.e.

Re: [R] Invisible plot using RSvgDevice

2005-01-10 Thread Csardi Gabor
The new 0.6 package is submitted to CRAN. You can grab it from the author's web site: http://darkridge.com/~jake/RSvg if you want it now. It works with R =2.0.0 correctly. Cheers, Gabor On Mon, Jan 10, 2005 at 12:13:44PM +, Prof Brian Ripley wrote: Have you contacted the author? The

RE: [R] Query: simple autocorrelation ina time series

2005-01-10 Thread Doran, Harold
You need mvrnorm() to do this. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Crabb, David Sent: Monday, January 10, 2005 1:04 PM To: r-help@stat.math.ethz.ch Subject: [R] Query: simple autocorrelation ina time series I hope you can help with what

Re: [R] dyn.load Excel add-in on Windows XP and XLCall32.dll

2005-01-10 Thread Prof Brian Ripley
Use pedump -i to find the depedencies (see README.packages). This is a Windows issue, and you will need to seek advice on a Windows forum. That message is coming (as it says) from the Windows routine LoadLibrary. On Mon, 10 Jan 2005 [EMAIL PROTECTED] wrote: I am on Win XP SP2 using R 2.0.1

Re: [R] Graphical table in R

2005-01-10 Thread Dan Bolser
Cheers. This is really me just being lazy (as usual). The latex function in Hmisc allows me to make a .ps file then grab a screen shot of that ps and make a .png file. I would just like to use plot so I can wrap it in a png command and not have to use the 'screen shot' in between. I can use

Re: [R] Query: simple autocorrelation ina time series

2005-01-10 Thread Rolf Turner
?arima.sim __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Mixing portrait/landscape in a postscript file

2005-01-10 Thread Ted Harding
On 10-Jan-05 Prof Brian Ripley wrote: On Mon, 10 Jan 2005, Yves Brostaux wrote: [...] trellis.device(postscript, theme=canonical.theme(postscript, color=F), file=an_phase2_graph.ps, paper=a4, pointsize = 10, onefile=TRUE, horizontal=TRUE) Everything works fine,but some of these plots

Re: [R] Graphical table in R

2005-01-10 Thread Peter Dalgaard
Dan Bolser [EMAIL PROTECTED] writes: Cheers. This is really me just being lazy (as usual). The latex function in Hmisc allows me to make a .ps file then grab a screen shot of that ps and make a .png file. I would just like to use plot so I can wrap it in a png command and not have to use

[R] mle() and with()

2005-01-10 Thread Ben Bolker
I'm trying to figure out the best way of fitting the same negative log-likelihood function to more than one set of data, using mle() from the stats4 package. Here's what I would have thought would work: -- library(stats4) ## simulate values r = rnorm(1000,mean=2) ## very basic

Re: [R] Graphical table in R

2005-01-10 Thread Thomas Lumley
On Mon, 10 Jan 2005, Peter Dalgaard wrote: A screen shot of a ps file? That sounds ... weird. If you can view it, presumably you have Ghostscript and that can do png files. However, will textual output do? plot(0,type=n,axes=FALSE, xlab=, ylab=) con - textConnection(txt,w) sink(con);

Re: [R] Mixing portrait/landscape in a postscript file

2005-01-10 Thread Prof Brian Ripley
On Mon, 10 Jan 2005 [EMAIL PROTECTED] wrote: On 10-Jan-05 Prof Brian Ripley wrote: On Mon, 10 Jan 2005, Yves Brostaux wrote: [...] trellis.device(postscript, theme=canonical.theme(postscript, color=F), file=an_phase2_graph.ps, paper=a4, pointsize = 10, onefile=TRUE, horizontal=TRUE) Everything

Re: [R] Graphical table in R

2005-01-10 Thread apjaworski
The outer two lines could also be done more easily (?) by: library(gplots) #part of gregmisc bundle textplot(txt) Check ?textplot for additional options. Andy __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail:

[R] defining lower part of distribution with covariate

2005-01-10 Thread Troels Ring
Dear friends, forgive me a simple question, possibly related to quantreg but I failed to get it done and hope for basic instruction. I have two sets of observed Glasgow coma scores at admission to ICU after operation, and accompanying time of anesthesia (in hours). Thio is cheap and perhaps old

[R] Stadard errors and boxplots with 632plus error estimator, errorest

2005-01-10 Thread Antoine Beuchat
Dear R-users, I'd like to estimate standard errors (for lda) and make a boxplot with the 632plus and boot error estimators included in package ipred (method: errorest). The boot estimator returns only a standard deviation but not the whole error data. Thank you in advance, regards, Antoine

[R] Data Set

2005-01-10 Thread Talita Leite
Hi everybody, I'm studying descriptive statistics using R and I have to make an important work about that. I need some help to choose a good data set to apply those statistics. Does anybody know a good data set I could work with? Thanx, Talita Perciano Costa Leite Graduanda em Ciência da

[R] Help requested .....

2005-01-10 Thread Sridhar Murthy
Happy New Year !!! I am working on a AIX machine with R Version 1.9.1 on AIX 5.2 ML-3. The perl.rte fileset is at 5.8.0.0 The following line is a part of a cgi script: open(T, | /usr/local/bin/R --silent --no-save /var/tmp/lll.input /var/tmp/messages.out); The contents of file

[R] How to obtain nls parameter estimates

2005-01-10 Thread Jagarlamudi, Choudary
How can I receive parameter estimates for a given curve that has been fit with the NLS function. The problem is that I have a n parameter curve and I want the optimal fit. The NLS procedure gives me final function values and not the individual parameter estimates that were used to define

Re: [R] Data Set

2005-01-10 Thread Spencer Graves
Much of the documentation including manuals, help files and other contributed descriptions include same data sets. What kinds of applications and techniques most interest you? That with the posting guide (http://www.R-project.org/posting-guide.html, especially the search at

RE: [R] How to obtain nls parameter estimates

2005-01-10 Thread Liaw, Andy
Use coefficients(nlsobject), or if you want standard errors, t-statistics, and p-values, summary(nlsobject)$paramters. Andy From: Jagarlamudi, Choudary How can I receive parameter estimates for a given curve that has been fit with the NLS function. The problem is that I have a n

[R] Data Set

2005-01-10 Thread Talita Leite
Hi, I'll try to be more specific asking my question. I want to apply some functions like mean(), median(), var(), sd(), mad(), quantile(), kurtosis(), skewness() and make some graphics like boxplot, barplot, histogram, stars... In order to do that I need a simple data set, simple but

RE: [R] Data Set

2005-01-10 Thread Liaw, Andy
There are a few packages on CRAN that are collections of data sets, some from intro textbooks. You might find some of them suitable. There are also datasets that come with R. Type data() at the R prompt to see a list. Andy From: Talita Leite Hi, I'll try to be more specific asking my

RE: [R] Windows package development: bad html links to functions in non-standard packages

2005-01-10 Thread Andrew Ward
Thank you to Professor Ripley for clearing this up for me. Embarrassingly, I read this FAQ several times, completely missing the point each time. Regards, Andrew C. Ward, Senior Analyst (Quantitative), Tel: +61 7 3864 0439 Queensland Studies Authority, Fax: +61 7 3229 3318 295

[R] Calculate Mean of Column Vectors?

2005-01-10 Thread Thomas Hopper
Hello, I've got an array defined as y - rnorm(3000), dim(y) - c(3, 1000). I'd like to produce a 1000-element vector z that is the mean of the corresponding elements of y (like z[1,1] - mean(y[1,1], y[2,1], y[3,1])), but being new to R, I'm not sure how to do this for all elements at once (or,

Re: [R] Calculate Mean of Column Vectors?

2005-01-10 Thread Marc Schwartz
On Mon, 2005-01-10 at 22:45 -0500, Thomas Hopper wrote: Hello, I've got an array defined as y - rnorm(3000), dim(y) - c(3, 1000). I'd like to produce a 1000-element vector z that is the mean of the corresponding elements of y (like z[1,1] - mean(y[1,1], y[2,1], y[3,1])), but being new

Re: [R] Calculate Mean of Column Vectors?

2005-01-10 Thread Tim F Liao
Tom, I often use a loop, such as: m - rep(0,1000) for (i in 1:1000) + m[i]-mean(y[,i]) m[1:20] [1] -0.04914724 -0.28253861 -0.31112690 -0.18034371 0.18839167 0.66448244 [7] 0.19769017 -1.28363405 -0.05167451 -0.95492534 -1.23285174 0.10288562 [13] -0.73792584 -0.19297468 -0.59059036

RE: [R] Calculate Mean of Column Vectors?

2005-01-10 Thread Mulholland, Tom
apply(y,2,mean) Is this what you are after. If it is I would suggest that you look at the examples not just for this but for what I call the apply family sapply, tapply, mapply. Once you get the hang of these they are really helpful. Tom. -Original Message- From: Thomas Hopper

Re: [R] Calculate Mean of Column Vectors?

2005-01-10 Thread Simon Blomberg
z - apply(y, 2, mean) Cheers, Simon. Hello, I've got an array defined as y - rnorm(3000), dim(y) - c(3, 1000). I'd like to produce a 1000-element vector z that is the mean of the corresponding elements of y (like z[1,1] - mean(y[1,1], y[2,1], y[3,1])), but being new to R, I'm not sure how

Re: [R] Calculate Mean of Column Vectors?

2005-01-10 Thread Spencer Graves
Have you considered the following: apply(y, 1, mean) Alternatively: y.means - rep(NA, 3) for(i in 1:3) y.means[i] - mean(y[i,]) hope this helps. spencer graves Thomas Hopper wrote: Hello, I've got an array defined as y - rnorm(3000), dim(y) - c(3, 1000). I'd like

[R] problem with fitted probability in glm

2005-01-10 Thread flavio . zanini
Good morning, Sorry to disturb you. I have a set of variable (not correlated) and with a stepwise selection procedure (step) on my multiple glm I obtain a result but with a warning message saying: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x[, jj, drop = FALSE], y, wt, offset

Re: [R] suse 9.1 x86_64 rpms?

2005-01-10 Thread Kristian Eric Markon
Thanks for the help. You were all right--I was missing the -devel versions of a few of the packages. Everything seems to be working well in 64-bit so far. Peter Dalgaard wrote: Kristian Eric Markon [EMAIL PROTECTED] writes: I noticed on the mailing list archives and through Google searches