[R] pairs: altering pch options on upper and lower panel of pairwise scatter plots

2005-01-17 Thread Paul Hewson
Hello, I can't figure out how to use the upper.panel and lower.panel options in pairs to alter the label options for either panel independently of the other. I would like to be able to show the pairwise scatter plots for the data as they are (a vanilla pairs plot?) but separately to be able to

[R] 3d bar plot

2005-01-17 Thread R user
This graph - http://www.math.hope.edu/~tanis/dallas/images/disth36.gif is an example I found at http://www.math.hope.edu/~tanis/dallas/disth1.html created by Maple. Does anybody know how to create something similar in R? I have a feeling it could be possible using scatterplot3d (perhaps with

Re: [R] RWinEdt install problem

2005-01-17 Thread Chuck Cleland
I am using WinEdt 5.4 (Build: 20040513 (v. 5.4), on XP service pack 2), and the following approach works for me: 1) Install RWinEdt_1.6-2.zip using RGui menus. 2) Include the following line in the Rprofile file. library(RWinEdt) This will launch WinEdt with R customizations whenever RGui is

[R] help wanted using R in a classroom

2005-01-17 Thread Sam Parvaneh
Hi everyone! I'm using R 2.0.1 for Mac OS X in a classroom with 40 eMacs running Mac OS X version 10.3.6. These Macs are network based, meaning that the students log in to an XServe G4 where their user accounts and home directories are stored. The problem that I'm having each time a group of

RE: [R] 3d bar plot

2005-01-17 Thread Marwan Khawaja
You can check these packages, ?scatterplot3d ?scatter3d Best Marwan --- Marwan Khawaja http://staff.aub.edu.lb/~mk36/ --- -Original Message- From:

[R] bold face labelling/expression

2005-01-17 Thread hagric
Dear colleagues, I have a great problem in using expression for axes labels. I want the labels in bold face (i.e.: par(font.lab=2)). When typing boxplot(y ~ groups, names = , xlab = , ylab = , axes = F) axis(side=1, at=c(1,2), xlab=c(expression(H[2]*O),others)), I do not get a bold face

re:[R] Multiple plots in one screen

2005-01-17 Thread Niels Steen Krogh
see ?par look for par(new=T) (if you want plot-overlay) or just ?lines (if you want several lines in one plot) __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] randomForest: too many element specified?

2005-01-17 Thread luk
When I run randonForest with a 169453x5 matrix, I got the following message. Error in matrix(0, n, n) : matrix: too many elements specified Can you please advise me how to solve this problem? Thanks, Lu Uwe Ligges [EMAIL PROTECTED] wrote: Vera Hofer wrote: Dear colleagues, I have a

[R] Question about time series

2005-01-17 Thread Dennis Fisher
I have data in the following format: DATE [1] 01/13/2004 In order to find the difference between two data points, I presently use brute force to calculate the day of the year: strptime(DATE, format=%m/%d/%Y)$yday [1] 12 Although this works, it may not be robust over different years. I assume

RE: [R] randomForest: too many element specified?

2005-01-17 Thread Liaw, Andy
From: luk When I run randonForest with a 169453x5 matrix, I got the following message. Error in matrix(0, n, n) : matrix: too many elements specified Can you please advise me how to solve this problem? Thanks, Lu 1. When asking new questions, please don't reply to other

RE: [R] 3d bar plot

2005-01-17 Thread John Fox
Dear Marwan and Jonne, I don't think that there's a scatter3d package, so perhaps Marwan is referring to the scatter3d() function in the Rcmdr package. If so, that function won't make the kind of 3D graph that Jonne wants -- though the rgl package, on which scatter3d() is based, should be able to

[R] find source code

2005-01-17 Thread r . ghezzo
I am using R 2.0.2 on a WinXP I am trying to get the code of the Kruskal-Wallis test but kruskal.test function (x, ...) UseMethod(kruskal.test) environment: namespace:stats ls(3) [1] acf acf2AR add.scope .. [181] kruskal.test ks.test

Re: [R] find source code

2005-01-17 Thread Simon Wood
stats:::kruskal.test.default On Mon, 17 Jan 2005 [EMAIL PROTECTED] wrote: I am using R 2.0.2 on a WinXP I am trying to get the code of the Kruskal-Wallis test but kruskal.test function (x, ...) UseMethod(kruskal.test) environment: namespace:stats ls(3) [1] acf

RE: [R] find source code

2005-01-17 Thread Liaw, Andy
From: [EMAIL PROTECTED] I am using R 2.0.2 on a WinXP Does that version of R exist? I am trying to get the code of the Kruskal-Wallis test but kruskal.test function (x, ...) UseMethod(kruskal.test) environment: namespace:stats ls(3) [1] acf acf2AR

[R] [basic ?] Merging result of by processing with a data frame

2005-01-17 Thread Jean-Louis Abitbol
Dear All, I would like to merge a data frame such as: basetab subject dose cmax 1 1031 50 21.8 with the result of a by processing such as: tmax-by(pkga,subject,f.tmax) tmax subject: 1031 [1] 6 -- f.tmax being a function

Re: [R] Question about time series

2005-01-17 Thread Christian Schulz
Hi TODAY - as.POSIXlt('2004-01-17') LYTODAY - as.POSIXlt('2003-01-17') DAYS - TODAY-LYTODAY DAYS Time difference of 365 days DAYS[[1]] [1] 365 perhaps it helps, Christian Dennis Fisher wrote: I have data in the following format: DATE [1] 01/13/2004 In order to find the difference

Re: [R] find source code

2005-01-17 Thread Uwe Ligges
Simon Wood wrote: stats:::kruskal.test.default and how to get there: methods(kruskal.test) # note, you probably want the default method! getS3method(kruskal.test, default) Uwe On Mon, 17 Jan 2005 [EMAIL PROTECTED] wrote: I am using R 2.0.2 on a WinXP I am trying to get the code of the

Re: [R] 3d bar plot

2005-01-17 Thread Uwe Ligges
John Fox wrote: Dear Marwan and Jonne, I don't think that there's a scatter3d package, so perhaps Marwan is referring to the scatter3d() function in the Rcmdr package. If so, that function won't make the kind of 3D graph that Jonne wants -- though the rgl package, on which scatter3d() is based,

[R] kmeans with weighted variables

2005-01-17 Thread HYRIEN Matthieu
Hello, I'm having troubles finding how to weight some variables for a k-means clustering. Can you please let me know what arguments I should had to: z - kmeans(A, 3, 10) if I want to put different weights for 3 of the 6 variables I'm using to classify these individuals. Many

Re: [R] Setting the width and height of a Sweave figure

2005-01-17 Thread Bjørn-Helge Mevik
I haven't found any other solution than using fig=TRUE,height=7,width=14 theCode() @ (but of course that doesn't have any effect when theCode() is used interactively). -- Bjørn-Helge Mevik __ R-help@stat.math.ethz.ch mailing list

[R] Omitting constant in ols() from Design

2005-01-17 Thread Tobias Muhlhofer
Hi! I need to run ols regressions with Huber-White sandwich estimators and the correponding standard errors, without an intercept. What I'm trying to do is create an ols object and then use the robcov() function, on the order of: f - ols(depvar ~ ind1 + ind2, x=TRUE) robcov(f) However, when I

Re: [R] find source code

2005-01-17 Thread TEMPL Matthias
But how can I see the *documented* source code in Windows? I think i.e. in the /flexmix/library/man/all.rda file is the documented(?) code for all functions of Package flexmix, but the standard Windows XP cannot open .rda files correct. Is there a way to read the documented source code of a

RE: [R] find source code

2005-01-17 Thread Liaw, Andy
From: TEMPL Matthias But how can I see the *documented* source code in Windows? I think i.e. in the /flexmix/library/man/all.rda file is the documented(?) code for all functions of Package flexmix, but the standard Windows XP cannot open .rda files correct. Is there a way to read the

Re: [R] find source code

2005-01-17 Thread Peter Dalgaard
TEMPL Matthias [EMAIL PROTECTED] writes: But how can I see the *documented* source code in Windows? I think i.e. in the /flexmix/library/man/all.rda file is the documented(?) code for all functions of Package flexmix, but the standard Windows XP cannot open .rda files correct. No that's not

[R] How do I format something as 0.000?

2005-01-17 Thread Christian Hennig
Hi, I would like to use the format function to get numbers all with three digits to the right of the decimal point, even in cases where there is no significant digit left. For example, I would like to get c(0.3456789,0.053) as 0.346 0.000. It seems that it is not possible to force format to

Re: [R] Omitting constant in ols() from Design

2005-01-17 Thread Achim Zeileis
On Mon, 17 Jan 2005 14:57:33 + Tobias Muhlhofer wrote: Hi! I need to run ols regressions with Huber-White sandwich estimators and the correponding standard errors, without an intercept. What I'm trying to do is create an ols object and then use the robcov() function, on the order of:

Re: [R] How do I format something as 0.000?

2005-01-17 Thread Achim Zeileis
On Mon, 17 Jan 2005 16:49:35 +0100 (MET) Christian Hennig wrote: Hi, I would like to use the format function to get numbers all with three digits to the right of the decimal point, even in cases where there is no significant digit left. For example, I would like to get

Re: [R] Omitting constant in ols() from Design

2005-01-17 Thread Frank E Harrell Jr
Tobias Muhlhofer wrote: Hi! I need to run ols regressions with Huber-White sandwich estimators and the correponding standard errors, without an intercept. What I'm trying to do is create an ols object and then use the robcov() function, on the order of: f - ols(depvar ~ ind1 + ind2, x=TRUE)

[R] merge data frames taking mean/mode of multiple macthes

2005-01-17 Thread avneet singh
Hello :) I have two data frames, one has properties taken on a piece by piece basis and the other has performance on a lot by lot basis. I wish to combine these two data frames but the problem is that each lot has multiple pieces and hence i need to take a mean of the properties of multiple

Re: [R] RWinEdt install problem

2005-01-17 Thread Murray Eisenberg
The problem with RWinEdt can be resolved by NOT using User Profiles in WinEdt: I did a clean re-install of WinEdt 5.4, build 20050114. Then, executing library(RWinEdt) in the RGui did start up a properly configured R-WinEdt session of WinEdt. No error messages; the expected R buttons were

Re: [R] find source code

2005-01-17 Thread TEMPL Matthias
-Ursprüngliche Nachricht- Von: Liaw, Andy [mailto:[EMAIL PROTECTED] Gesendet: Montag, 17. Jänner 2005 16:33 An: TEMPL Matthias; r-help@stat.math.ethz.ch Betreff: RE: [R] find source code From: TEMPL Matthias But how can I see the *documented* source code in Windows? I

Re: [R] merge data frames taking mean/mode of multiple macthes

2005-01-17 Thread Christoph Buser
Dear Avneet the function aggregate (see also ?aggregate) could be useful for your problem. See the short example I've written below. dat1 - data.frame(lot = factor(1:10),y1 = rnorm(10)) str(dat1) dat2 - data.frame(nr = 1:100, lot = factor(rep(1:10, each = 10)),y2 = rnorm(100)) str(dat2)

Re: [R] 3d bar plot

2005-01-17 Thread Duncan Murdoch
On Mon, 17 Jan 2005 08:49:06 -0500, John Fox [EMAIL PROTECTED] wrote : Dear Marwan and Jonne, I don't think that there's a scatter3d package, so perhaps Marwan is referring to the scatter3d() function in the Rcmdr package. If so, that function won't make the kind of 3D graph that Jonne wants --

[R] debian /etc/apt/sources.list for CRAN?

2005-01-17 Thread Michael Friendly
When I set up my debian linux (unstable) system, the only R debian packages were on CRAN, so I followed directions and added the last line to my /etc/apt/sources.list euclid: # m /etc/apt/sources.list deb http://debian.yorku.ca/debian/ unstable main non-free contrib deb

[R] Skewness test

2005-01-17 Thread Christian Hennig
Hi, is there a test for the H0 skewness=0 (or with skewness as test statistic and normality as H0) implemented in R? Thank you, Christian *** Christian Hennig Fachbereich Mathematik-SPST/ZMS, Universitaet Hamburg [EMAIL

Re: [R] find source code

2005-01-17 Thread r . ghezzo
Thanks to all who answered my query, I forgot completely to call methods() first to check the true whole name of the function. Heberto Ghezzo Quoting Uwe Ligges [EMAIL PROTECTED]: Simon Wood wrote: stats:::kruskal.test.default and how to get there: methods(kruskal.test) # note, you

[R] C-c C-c and tcltk cause R to get a segmentation fault

2005-01-17 Thread Jean Eid
when tcltk is called and a C-c C-c is issued a message of Illegal instruction is issued and the R process is killed. To reproduce please do the following require(tcltk) C-c C-c It will say Illegal instruction I use a debian 2.4 kernel with R 2.0.1 tcltck version 2.0.1. I use tcltk to have gui

Re: [R] debian /etc/apt/sources.list for CRAN?

2005-01-17 Thread Dirk Eddelbuettel
On Mon, Jan 17, 2005 at 01:39:23PM -0500, Michael Friendly wrote: Should I just delete the cran.r-project.org line above, I think so. Plus, if you already point to Debian unstable, CRAN has nothing you wouldn't have otherwise. or are there debian R packages elsewhere on CRAN that I should

Re: [R] Skewness test

2005-01-17 Thread Dirk Eddelbuettel
On Mon, Jan 17, 2005 at 07:59:32PM +0100, Christian Hennig wrote: Hi, is there a test for the H0 skewness=0 (or with skewness as test statistic and normality as H0) implemented in R? Not that I know of, but tseries has the standard Jarque-Bera test that combines 3rd and 4th moment tests into

Re: [R] bwplot: how not to draw outliers

2005-01-17 Thread Sundar Dorai-Raj
RenE J.V. Bertin wrote: Hello, and (somewhat belated) best wishes for 2005. Can one order not to draw outliers in bwplot, or at least exclude them from the vertical axis scaling? If so, how (or what doc do I need to consult)? The options that have this effect in boxplot() do not appear to have

RE: [R] subsetting like in SAS

2005-01-17 Thread Denis Chabot
I want to thank Petr Pikal, Robert Balshaw and Na Li for suggesting the use of unique or !duplicated on a subset of my data where unwanted variables have been removed. This worked perfectly. Denis Chabot On 13 Jan 2005 at 11:52, Denis Chabot wrote: Hi, Being in the process of translating some

[R] problem installing RSPython

2005-01-17 Thread Robert Sams
Hi, I'm trying to install RSPython v0.5-4 on a debian machine (woody, testing) but am having the following problem. $R CMD INSTALL -c --library=/usr/lib/R/library RSPython_0.5-4.tar.gz 2err But then... $python import RS Error in .PythonInit() : Error in Python call: values Error in

Re: [R] bwplot: how not to draw outliers

2005-01-17 Thread Deepayan Sarkar
On Monday 17 January 2005 14:51, Sundar Dorai-Raj wrote: RenE J.V. Bertin wrote: Hello, and (somewhat belated) best wishes for 2005. Can one order not to draw outliers in bwplot, or at least exclude them from the vertical axis scaling? If so, how (or what doc do I need to consult)? The

Re: [R] bwplot: how not to draw outliers

2005-01-17 Thread Frank E Harrell Jr
Sundar Dorai-Raj wrote: RenE J.V. Bertin wrote: Hello, and (somewhat belated) best wishes for 2005. Can one order not to draw outliers in bwplot, or at least exclude them from the vertical axis scaling? If so, how (or what doc do I need to consult)? The options that have this effect in

Re: [R] How do I format something as 0.000?

2005-01-17 Thread Marc Schwartz
On Mon, 2005-01-17 at 14:19 -0600, Dirk Eddelbuettel wrote: On Mon, Jan 17, 2005 at 04:49:35PM +0100, Christian Hennig wrote: Hi, I would like to use the format function to get numbers all with three digits to the right of the decimal point, even in cases where there is no significant

[R] a question of mixed effect in R

2005-01-17 Thread chang10000.tw
Dear all, I have a question about mixed effect model in R. The data set has 5 variables, X(response),subject, times, repeat, indicator The model is X_hijk=a_h+Z_h*b_i+r(ij)+e_hijk , where h=0,1(indicator), i=1,...,n(subject), j=1,...,n_i(times within subject; nested effect),k=1,2,3(repeat).

Re: [R] debian /etc/apt/sources.list for CRAN?

2005-01-17 Thread Christian T. Steigies
On Mon, Jan 17, 2005 at 02:14:04PM -0600, Dirk Eddelbuettel wrote: On Mon, Jan 17, 2005 at 01:39:23PM -0500, Michael Friendly wrote: Should I just delete the cran.r-project.org line above, I think so. Plus, if you already point to Debian unstable, CRAN has nothing you wouldn't have

[R] rpart

2005-01-17 Thread Weiwei Shi
Hi, there: I am working on a classification problem by using rpart. when my response variable y is binary, the trees grow very fast, but if I add one more case to y, that is making y has 3 cases, the tree growing cannot be finished. the command looks like: x-rpart(r0$V142~.,data=r0[,1:141],

Re: [R] debian /etc/apt/sources.list for CRAN?

2005-01-17 Thread Julien Damon
I think the debian packages are now maintened in the official debian repository. http://packages.debian.org/cgi-bin/search_packages.pl?keywords=r-cransearchon=namessubword=1version=unstablerelease=all

Re: [R] debian /etc/apt/sources.list for CRAN?

2005-01-17 Thread Dirk Eddelbuettel
On Mon, Jan 17, 2005 at 09:00:53PM +0100, Julien Damon wrote: I think the debian packages are now maintened in the official debian repository. Yes, as they have been since December 1997. Dirk -- Better to have an approximate answer to the right question than a precise answer to the wrong

[R] Weighted least squares

2005-01-17 Thread Ming Hsu
Hi, I would like to run a weighted least squares with the the weighting matrix W. I ran the following two regressions, (W^-1)Y = Xb + e Y = WXb+ We In both cases, E[bhat] = b. I used the following commands in R lm1 - lm(Y/W ~ X) lm2 - lm(Y ~ W:X, weights = W) where Y - rnorm(10,1) X - Y +

[R] lme confusion

2005-01-17 Thread Ming Hsu
Hi, this is my first time using the nlme package, and I ran into the following puzzling problem. I estimated a mixed effects model using lme, once using groupedData, once explicitly stating the equations. I had the following outputs. All the coefficients were similar, but they're always

[R] Why can't run Gamma Distribution by GLM?

2005-01-17 Thread may sawasdee
y-cy-c(15,16,17,0,20,50,0,11,60,12) x-c(1,6,3,11,13,20,12,14,15,3) gamma-glm(y~x,family=Gamma()) Error in eval(expr, envir, enclos) : Non-positive values not allowed for the gamma family but if i change data of y is y1-c(15,16,17,12,20,50,51,11,60,12) ## doesn't have zero then I can

[R] X11 installation problem Linux X86_64

2005-01-17 Thread Robert Nelson
Hello all, This is my first time posting, so forgive my ignorance. I've searched high and low but haven't come across the specific answer to my questions. configure runs fine, then make gives this error: ... gcc -shared -L/usr/local/lib -o R_X11.so dataentry.lo devX11.lo rotated.lo rbitmap.lo

RE: [R] 3d bar plot

2005-01-17 Thread Lorenz . Gygax
This graph - http://www.math.hope.edu/~tanis/dallas/images/disth36.gif is an example I found at http://www.math.hope.edu/~tanis/dallas/disth1.html created by Maple. Does anybody know how to create something similar in R? I have a feeling it could be possible using scatterplot3d

RE: [R] 3d bar plot

2005-01-17 Thread Marwan Khawaja
Dear John, Yes, I meant the scatter3d function in the Rcmdr package -- I will 'behave' next time :-) Jonne asked for 'something similar in R' -- hence the suggestion to also use the package scatterplot3d. Best Marwan --- Marwan

Re: [R] rpart

2005-01-17 Thread Prof Brian Ripley
On Mon, 17 Jan 2005, Weiwei Shi wrote: I am working on a classification problem by using rpart. when my response variable y is binary, the trees grow very fast, but if I add one more case to y, that is making y has 3 cases, Do you mean 3 classes?: you have many more than 3 cases below. the tree

Re: [R] Weighted least squares

2005-01-17 Thread Prof Brian Ripley
On Mon, 17 Jan 2005, Ming Hsu wrote: I would like to run a weighted least squares with the the weighting matrix W. This is generalized not weighted least squares if W really is a matrix and not a vector of case-by-case weights. I ran the following two regressions, (W^-1)Y = Xb + e Y = WXb+ We