Re: [R] POSIXct bug?

2005-02-28 Thread Martin Maechler
Spencer == Spencer Graves [EMAIL PROTECTED] on Sun, 27 Feb 2005 07:15:44 -0800 writes: Spencer Hi, Gabor: Of course: time zone vs. GMT. SpencerNext question: Might a simple example that Spencer illustrates this be added to the help file for Spencer as.POSIXct, and

[R] How to plot

2005-02-28 Thread WeiQiang . Li
Hello, I would like to know if R provides the similar function to produce the 4-in-1 graph (Normal Plot of Residuals, I Chart of Residuals, Histogram of Residuals Residuals vs Fits) as MiniTab. If do not have, could anyone please tell me how to produce above-mentioned individual

Re: [R] How to plot

2005-02-28 Thread Dimitris Rizopoulos
you could use directly the plot function on an `lm' object, e.g., x - runif(100, -3, 3) y - 1 + 2*x + rnorm(100) m - lm(y~x) par(mfrow=c(2,2)) plot(m) This doesn't produce the I Chart you wanted but it gives a very good picture about any strange behaviour of the residuals. I hope it helps. Best,

RE: [R] averaging within columns

2005-02-28 Thread Rau, Roland
Hi, I'm hoping there is something like mean(dataname$wtime[name]) which will just create a column with length equal to the number of different names (levels) and an average wtime for each. So far though, I haven't had much luck figuring that one out. Does the following code do

Re: [R] How to plot

2005-02-28 Thread Petr Pikal
On 28 Feb 2005 at 17:35, [EMAIL PROTECTED] wrote: Hello, I would like to know if R provides the similar function to produce the 4-in-1 graph (Normal Plot of Residuals, I Chart of Residuals, Histogram of Residuals Residuals vs Fits) as MiniTab. If do not have, could

[R] Unable to install packages

2005-02-28 Thread MARTIN CALMARZA AGUSTIN
Dear Sirs , Today I've downloaded the last release of the R-Program (I an a novel user). I am interested on performing time series analysys (regression,ARIMA,ARCH) to some data-sets, therefore I would like to use some of the packages of R-Contributors (like the tseries one).When I try to

Re: [R] Unable to install packages

2005-02-28 Thread Peter Dalgaard
MARTIN CALMARZA AGUSTIN [EMAIL PROTECTED] writes: Dear Sirs , Today I've downloaded the last release of the R-Program (I an a novel user). I am interested on performing time series analysys (regression,ARIMA,ARCH) to some data-sets, therefore I would like to use some of the packages

Re: [R] Problems with downloading

2005-02-28 Thread Uwe Ligges
Maria Befring Hovda wrote: Hello I'm using your program R in a course I'm taking at the University of Oslo, and thereby I need to download it to my PC. Unfortunately I do have some problems, I do not know whish files to download and how I do it. Can you please send me a list of what to download

Re: [R] Problems with downloading

2005-02-28 Thread Adaikalavan Ramasamy
I am assuming you are talking about a Windows machine. If so, reading http://cran.r-project.org/bin/windows/base/rw-FAQ.html#Installation-and- Usage will point you the following windows installer http://cran.r-project.org/bin/windows/base/rw2001.exe However, you may need to ask your IT person

[R] Four parameter Lognormal Dis.

2005-02-28 Thread F Monadjemi
Dear Sir, I want to generate a random observations from four parameter lognormal distribution. Would you please tell me how I can do that? Thanks Farinaz __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Problems with downloading

2005-02-28 Thread Uwe Ligges
Adaikalavan Ramasamy wrote: I am assuming you are talking about a Windows machine. If so, reading http://cran.r-project.org/bin/windows/base/rw-FAQ.html#Installation-and- Usage will point you the following windows installer http://cran.r-project.org/bin/windows/base/rw2001.exe However, you may

[R] persistance of factor levels in a data frame

2005-02-28 Thread Lefebure Tristan
Hi, Just something I don't understand: data - data.frame(V1=c(1:12),F1=c(rep(a,4),rep(b,4),rep(c,4))) data_ac - data[which(data$F1 !=b), ] levels(data_ac$F1) Why the level b is always present ? thanks Tristan, R 2.0.1 for Linux Fedora 3 --

Re: [R] persistance of factor levels in a data frame

2005-02-28 Thread Dimitris Rizopoulos
look at ?[.data.frame and also check this: dat - data.frame(V1=c(1:12), F1=rep(letters[1:3], each=4)) dat.ac - dat[dat$F1 !=b, ] ### dat.ac$F1 dat.ac$F1[, drop=TRUE] ### dat.ac$F1 - dat.ac$F1[, drop=TRUE] levels(dat.ac$F1) I hope it helps. best, Dimitris Dimitris

Re: [R] persistance of factor levels in a data frame

2005-02-28 Thread Peter Dalgaard
Lefebure Tristan [EMAIL PROTECTED] writes: Hi, Just something I don't understand: data - data.frame(V1=c(1:12),F1=c(rep(a,4),rep(b,4),rep(c,4))) data_ac - data[which(data$F1 !=b), ] levels(data_ac$F1) Why the level b is always present ? Because it is a property of the definition,

Re: [R] persistance of factor levels in a data frame

2005-02-28 Thread Douglas Bates
Lefebure Tristan wrote: Hi, Just something I don't understand: data - data.frame(V1=c(1:12),F1=c(rep(a,4),rep(b,4),rep(c,4))) data_ac - data[which(data$F1 !=b), ] levels(data_ac$F1) Why the level b is always present ? thanks Tristan, R 2.0.1 for Linux Fedora 3 You must explicitly drop

RE: [R] 3d scatterplots of more than 1 data set

2005-02-28 Thread John Fox
Dear Saurav, The scatter3d() function in the Rcmdr package will plot by groups, using different colours for the groups. The function can be used directly or called via the Rcmdr menus. scatter3d() appears to do what you want, but in any event is a straightforward function, and you should be able

Re: [R] persistance of factor levels in a data frame

2005-02-28 Thread Petr Pikal
On 28 Feb 2005 at 14:07, Lefebure Tristan wrote: Hi, Just something I don't understand: data - data.frame(V1=c(1:12),F1=c(rep(a,4),rep(b,4),rep(c,4))) data_ac - data[which(data$F1 !=b), ] levels(data_ac$F1) Why the level b is always present ? H Tristan from ?[.factor Extract or

Re: [R] persistance of factor levels in a data frame

2005-02-28 Thread Marc Schwartz
On Mon, 2005-02-28 at 14:07 +0100, Lefebure Tristan wrote: Hi, Just something I don't understand: data - data.frame(V1=c(1:12),F1=c(rep(a,4),rep(b,4),rep(c,4))) data_ac - data[which(data$F1 !=b), ] levels(data_ac$F1) Why the level b is always present ? thanks Tristan, R 2.0.1

[R] [R-pkgs] New package: ROCR (Visualizing classifier performance)

2005-02-28 Thread Tobias Sing
Dear R users, we are glad to announce the release of our new R package ROCR, for visualizing the performance of scoring classifiers (available on CRAN). We hope that the package might be useful for those of you working on classification problems. For details, see the package description below,

[R] Legal issues

2005-02-28 Thread Kenneth Cabrera
Hi R Users: I have some legal questions about R development. R is under GNU Licence, version 2. How the core team deal with their employer (most of the cases it Universities, I don't know if public or private) about their contribution to R development? Do you have to make all the work outside the

Re: [R] prediction, gam, mgcv

2005-02-28 Thread Simon Wood
My model is of the form: mod-gam(y~s(x0)+s(x1)+s(x2),family=poisson). But I want to get estimate and standard error of the difference of two fitted values. The following code shows you how to do this (a) for differences on the scale of the linear predictor, and (b) for differences on the

[R] A problem about outer()

2005-02-28 Thread Feng Chen
Dear all, I have something about function outer() that I can't understand. Just see the following example. The two NaNs are due to 0/0, but I can't figure out the cause of the last two errors. I wonder if some one can explain this for me.

[R] Ask for your help

2005-02-28 Thread sczhang
Dear Sir(or madam), I got following error information when I run the haplo.glm program.Could you tell me which wrong was happened in my program?Many thanks. Shanchun fit.gaus - haplo.glm(y ~ SEX+geno, family = gaussian, + data=my.data, locus.label=label, control = +

Re: [R] A problem about outer()

2005-02-28 Thread Adaikalavan Ramasamy
You might want to read (or re-read) the posting guide about giving a simple example. See comments below. On Mon, 2005-02-28 at 23:03 +0800, Feng Chen wrote: Dear all, I have something about function outer() that I can't understand. Just see the following example. The two NaNs are due to

[R] [R-pkgs] phpSerialize 0.8

2005-02-28 Thread Dieter Menne
New on CRAN: phpSerialize Version 0.8 Dieter Menne, [EMAIL PROTECTED] Serializes R objects for PHP import into an associative array. Main use is for building web pages with R-support. Has mainly been tested with lm,lme, nlme and their summaries. A web example is provide, showing -- How to

[R] Using mutiply imputed data in NLME

2005-02-28 Thread Shige Song
Dear All, I am doing a growth modeling using NLME. I have three levels in my data: observation, individual, household. About half of my total sample have missing values in my household-level covariates. Under this situation, the best way to go is probably to multiply impute the data (for, say, 5

Re: [R] tkRplot under macos x

2005-02-28 Thread Luke Tierney
It only works under X11; if that is what you want this posting should tell you what needs to be done: https://stat.ethz.ch/pipermail/r-sig-mac/2004-December/001465.html Best, luke On Mon, 28 Feb 2005, Sean Davis wrote: I have not successfully gotten tkRplot to install on macos 10.3.8, R

Re: [R] and [ESS] Starting ESS

2005-02-28 Thread B Suresh Krishna
Hi, Yes, I defined the HOME environment variable as instructed in the installation instructions. Also, I can happily switch buffers in Emacs (C-x C-b or using C-x C-o) to the R-window and back to the .r file I am editing without any problem... If it matters, I am using Windows XP on a Dell

[R] memory problem with mac os X

2005-02-28 Thread Edouard Henrion
Dear list, I am using R.2.0.1 on a G5 biprocessor 2.5GHz with 2Go RAM (Mac OS X 10.3.8). I'm trying to calculate an object of type dist. I am getting the following memory error : *** malloc: vm_allocate(size=1295929344) failed (error code=3) *** malloc[25960]: error: Can't allocate region

[R] draw random samples from empirical distribution

2005-02-28 Thread bogdan romocea
Dear useRs, I have an empirical distribution (not normal etc) and I want to draw random samples from it. One solution I can think of is to compute let's say 100 quantiles, then use runif() to draw a random number Q between 1 and 100, and finally run runif() again to pull a random value from the

[R] formatting output

2005-02-28 Thread Peyuco Porras Porras .
Dear R-users A basic question that I wasn't able to solve: Is it possible to get the results of the function 'quantile' expressed as data.frame? What I'm doing is to apply the following code to get the quantiles in a particular dataset:

[R] R Reference Card (especially useful for Newbies)

2005-02-28 Thread Berton Gunter
Newbies (and others!) may find the R Reference Card made available by Tom Short and Rpad at http://www.rpad.org/Rpad/Rpad-refcard.pdf or on the Contributed link on CRAN (where some other reference cards are also linked) useful. It categorizes and organizes a bunch of R's (S's) basic, most used

[R] fitting gamma using glm

2005-02-28 Thread Daniel Bogod
The data in Table 1, which is Table T.1 in Cox Snell's Applied Statistics, gives the intervals in service-hours between failures of the air-conditioning eequipment in 10 Boeing 720 jet aircraft. The following possible models are under consideration: (a) separate gamma distributions fitted to all

[R] number formatting

2005-02-28 Thread Laura Holt
Dear R People: I have used the command round(x,3) to produce values with 3 places to the right of the decimal. Is there any command to remove the leading zero before the decimal point, please: that is, if I have 0.375, how do I produce just .375, please? Thanks in advance R 2.0.1 for Windows

RE: [R] number formatting

2005-02-28 Thread McGehee, Robert
x - .001 x [1] 0.001 sub(+0, , x) [1] .001 Of course, that means you'll be storing this number as a character, but if you're looking to format how the number is printed, that's probably what you want. -Original Message- From: Laura Holt [mailto:[EMAIL PROTECTED] Sent: Monday,

[R] problems with Rd format

2005-02-28 Thread Renaud Lancelot
Dear R-helpers, I have an Rd file with a section: \details{ The model is:\cr \eqn{y | \lambda ~ Binomial(n, \lambda)},\cr [snip] } I would like the equation line to be indented, with a \tab character for example. How can I do that ? Moreover, the panel of greek letters available in HTML

[R] negative multinomial regression models

2005-02-28 Thread Jhilbe
I am trying to find out the log-likelihood function of the negative multinomial. I would like to program a regression model for this distribution for a book I am writing. Any help I can obtain will be most appreciated. Joseph Hilbe [[alternative HTML version deleted]]

Re: [R] formatting output

2005-02-28 Thread Marc Schwartz
On Mon, 2005-02-28 at 14:05 -0400, Peyuco Porras Porras . wrote: Dear R-users A basic question that I wasn't able to solve: Is it possible to get the results of the function 'quantile' expressed as data.frame? What I'm doing is to apply the following code to get the

[R] Re: R-help Digest, Vol 24, Issue 28

2005-02-28 Thread John Maindonald
You've omitted a comma. races2000 is a data frame, which for purposes of extracting rows behaves like a 2-dimenional object. The following works fine: hills2000 - races2000[races2000$type == 'hill', ] Additionally, you might like to ponder type - races2000[names(races2000)==type]

Re: [R] problems with Rd format

2005-02-28 Thread Prof Brian Ripley
On Mon, 28 Feb 2005, Renaud Lancelot wrote: Dear R-helpers, I have an Rd file with a section: \details{ The model is:\cr \eqn{y | \lambda ~ Binomial(n, \lambda)},\cr [snip] } I would like the equation line to be indented, with a \tab character for example. How can I do that ? Well, use markup

[R] Journal Quality R Graphs?

2005-02-28 Thread Werner Wernersen
Hi! I have browsed the help archives but did not find anything on the subject: How to make publication quality graphs with R best? Is there some document about that topic out there? The problem is that the graphs look nice on the screen but when printed in black and white every color apart

[R] (no subject)

2005-02-28 Thread Brett Stansfield
Dear R Can you tell me how to change the working directory of R It's just that I have some text files that I wish to save separately from the R filing structure eg. into C:/my documents and need to change the working directory of R so that it reads these files . This means if I ever upgrade the

RE: [R] (no subject)

2005-02-28 Thread Berton Gunter
help.search(working directory) The R folks have worked very hard to provide good documentation. Please make use of the Help system and other resources before posting, as the Posting Guide asks (below). -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the

RE: [R] Journal Quality R Graphs?

2005-02-28 Thread Berton Gunter
Please read the posting guide. As a minimum, you'll need to specify your OS and R version plus probably other particulars about what form the output must be in, etc. I am surprised that you were unable to find anything in the archives, as this is a frequent discussion topic. -- Bert Gunter

Re: [R] change working directory [was: (no subject)]

2005-02-28 Thread Spencer Graves
'help.search(working directory)' identifies getwd, the help file for which describes setwd. Similarly, www.r-project.org - search - R site search for working directory produced 564 hits. I checked the first 4, and they all mentioned either getwd or setwd. hope this helps.

Re: [R] formatting output

2005-02-28 Thread Gabor Grothendieck
Marc Schwartz MSchwartz at MedAnalytics.com writes: : : On Mon, 2005-02-28 at 14:05 -0400, Peyuco Porras Porras . wrote: : Dear R-users : : A basic question that I wasn't able to solve: Is it possible to get : the results of the function 'quantile' expressed as data.frame?

Re: [R] formatting output

2005-02-28 Thread Marc Schwartz
On Tue, 2005-03-01 at 01:21 +, Gabor Grothendieck wrote: snip Here is a slight variation of Marc's solution that avoids explicit setting of the row names: R my.tmp - split(warpbreaks$breaks, + list(warpbreaks$wool, warpbreaks$tension)) R my.tmp - lapply(my.tmp, quantile) R

[R] packages masking other objects

2005-02-28 Thread James Smith
hello all, I am trying to use the function getCovariateFormula(nlme) in conjunction with the library lme4. When I load both packages I get the following message and the getCovariateFormula function no longer works: library(nlme) library(lme4) Attaching package 'lme4': The following

[R] packages masking other objects

2005-02-28 Thread James Smith
hello all, I am trying to use the function getCovariateFormula(nlme) in conjunction with the library lme4. When I load both packages I get the following message and the getCovariateFormula function no longer works: library(nlme) library(lme4) Attaching package 'lme4':

Re: [R] packages masking other objects

2005-02-28 Thread Gabor Grothendieck
James Smith james.smith at cdu.edu.au writes: : : hello all, : : I am trying to use the function getCovariateFormula(nlme) in conjunction with the library lme4. When I : load both packages I get the following message and the getCovariateFormula function no longer works: : : library(nlme) :

[R] Problems Building Ron AIX 5.2.0.0 (Solved)

2005-02-28 Thread paul . boutros
Happily I got this to work, largely by trial-and-error. In hopes that this will help somebody else, my config.site ended up being: OBJECT_MODE=64 R_PAPERSIZE=letter CC=/usr/local/bin/gcc MAIN_LDFLAGS=-Wl,-brtl SHLIB_LDFLAGS=-Wl,-G Which is virtually identical to that recommended in R-admin:

Re: [R] (no subject)

2005-02-28 Thread Adaikalavan Ramasamy
Reading the posting guide http://www.R-project.org/posting-guide.html will tell you to specify an informative subject line in your postings. You can change working directory by using setwd(), see help(setwd). For example setwd(c:\My Documents) might work. I am not sure if you need to truncate the

[R] variable importance in random forest

2005-02-28 Thread Murad Nayal
Hello, In Breiman papers on random forests 4 variable importance measures are described. as far as I can tell only two are available in the random forest R package. reduction in accuracy when the variable is permuted, and the mean decrease in the gini index due to the variable (no permutation).

Re: [R] Journal Quality R Graphs?

2005-02-28 Thread Adaikalavan Ramasamy
Searching for graph publication on http://maths.newcastle.edu.au/~rking/R/ gave me the following hit : http://tolstoy.newcastle.edu.au/R/help/04/03/0202.html which suggests postscript(). Have you tried printing other documents in black and white on the same printer or tried different printers for

Re: [R] (no subject)

2005-02-28 Thread Adaikalavan Ramasamy
Sorry, all backslashes have to be doubled in R as mentioned in FAQ 2.14 http://cran.r-project.org/bin/windows/base/rw-FAQ.html#R-can_0027t-find- my-file Regards, Adai On Tue, 2005-03-01 at 04:33 +, Adaikalavan Ramasamy wrote: Reading the posting guide

[R] Creating matrices for Mantel test

2005-02-28 Thread scott . meers
This seems like a real simple question - but I am not finding the answer in the help files and manuals. I am not a real sophisticateduser of R so maybe that is why I cant seem to get this. Anyway - I ne= ed to create matrices for doing a Mantel test on my data which was

Re: [R] problems with Rd format

2005-02-28 Thread Renaud Lancelot
Prof Brian Ripley a crit : On Mon, 28 Feb 2005, Renaud Lancelot wrote: Dear R-helpers, I have an Rd file with a section: \details{ The model is:\cr \eqn{y | \lambda ~ Binomial(n, \lambda)},\cr [snip] } I would like the equation line to be indented, with a \tab character for example. How can I

Re: [R] Sweave and \input or \include LaTeX commands

2005-02-28 Thread Gregor GORJANC
Hi! What is wrong if there would be the same command? Recall my example from previous posts and at the end of this mail. If I have a file a.Rnw and this one inputs file a1.Rnw. My idea was that Sweave would check \input{a1} or \include{a1} statements. If file a1 would have extension .Rnw it

Re: [R] packages masking other objects

2005-02-28 Thread Uwe Ligges
James Smith wrote: hello all, I am trying to use the function getCovariateFormula(nlme) in conjunction with the library lme4. When I load both packages I get the following message and the getCovariateFormula function no longer works: library(nlme) library(lme4) Attaching package 'lme4':