RE: [R] Newbie

2003-08-14 Thread Liaw, Andy
Try typing R --gui=gnome at a terminal prompt and see if that works for you. I don't know if the binary for RedHat has the GNOME support compiled in. Andy -Original Message- From: Marcos Llobera [mailto:[EMAIL PROTECTED] Sent: Thursday, August 07, 2003 9:42 AM To: [EMAIL

Re: [R] A question on orthogonal basis vectors

2003-08-14 Thread kjetil brinchmann halvorsen
On 13 Aug 2003 at 10:50, John Fox wrote: Dear Fred, If I understand correctly what you want, the answer is not unique. Think about the 3D case where you start with one vector. (I assume, by the way, that you mean orthonormal and that you mean unique up to a reflection.) There are

Re: [R] Rcgi

2003-08-14 Thread Chris Evans
On 30 Jul 2003 at 19:01, Chris Evans wrote: I am keen to look at Rcgi as I want to put up some simple bits of R to do prescribed tasks on HTML form input. Rweb is overkill and worryingly flexible for what I want and it sounds as if Rcgi is more what I need. However, I can't get any of the

Re: [R] problems with lda

2003-08-14 Thread Prof Brian Ripley
On Wed, 6 Aug 2003, Nicolaas Busscher wrote: enclosed a simple R script (and a data file) , with calls lda similar to the example with the iris data in the documentation. it is not working and i dont understand the error message. can anybody help me? Nothing arrived attached: please include

Re: [R] read.spss doesn't work anymore [more info]

2003-08-14 Thread Hedderik van Rijn
Some more information about my read.spss issues: on a cleanly installed Debian Woody machine: ~/TMP % R R : Copyright 2003, The R Development Core Team Version 1.7.1 (2003-06-16) [...] library(foreign); x - read.spss(dataDef.sav) Error in read.spss(dataDef.sav) : Calloc could not allocate

[R] spdep error message

2003-08-14 Thread rnaidoo
Hello, I have been using the package spdep to run spatial regressions on a data set with about 2500 observations. It has performed well up until now, but the following code resulted in an error: load(Panel.90s.ok.R) attach(Panel.90s.ok) neighs-dnearneigh(cbind(x,y),0,5)

Re: [R] graph for selected lines in stars()

2003-08-14 Thread Peter Dalgaard BSA
Tito de Morais Luis [EMAIL PROTECTED] writes: Dear listers, The following command (derived from the example in the ?stars help page) works : data(mtcars) stars(mtcars[, 1:7]) But the following gives an error: stars(mtcars[1, 1:7]) Error in s.y[i, ] : incorrect number of dimensions

[R] Mapping function dependencies

2003-08-14 Thread David Khabie-Zeitoune
Hi Does anyone know of any automatic way of mapping out the dependencies of a function -- i.e. automatically listing which other functions that function calls? I tried using all.names or all.vars, but I can't get it to work on a function: e.g. test = function(x, y) return(sin(x + y)*cos(x + y))

[R] How to get the pseudo left inverse of a singular square matrix?

2003-08-14 Thread Feng Zhang
Dear R-listers, I have a dxr matrix Z, where d r. And the product Z*Z' is a singular square matrix. The problem is how to get the left inverse U of this singular matrix Z*Z', such that U*(Z*Z') = I? Is there any to figure it out using matrix decomposition method? Thanks a lot for your help.

Re: [R] nls confidence intervals

2003-08-14 Thread Duncan Murdoch
On Thu, 14 Aug 2003 09:08:26 -0700, Spencer Graves [EMAIL PROTECTED] wrote : This seems to identify a possible bug in R 1.7.1 under Windows 2000: tstDf - data.frame(y = 1:11, x=1:11) fit - nls(y~a/x, data=tstDf, start=list(a=1)) predict(fit, se.fit=TRUE) [1] 7.0601879 3.5300939 2.3533960

Re: [R] nls confidence intervals

2003-08-14 Thread Torsten Hothorn
On Thu, 14 Aug 2003, Spencer Graves wrote: This seems to identify a possible bug in R 1.7.1 under Windows 2000: tstDf - data.frame(y = 1:11, x=1:11) fit - nls(y~a/x, data=tstDf, start=list(a=1)) predict(fit, se.fit=TRUE) [1] 7.0601879 3.5300939 2.3533960 1.7650470 1.4120376

Re: [R] How to get the pseudo left inverse of a singular squarematrix?

2003-08-14 Thread Feng Zhang
Thank, Jerome The question is if this generalized inverse can make their product to be identity matrix? - Original Message - From: Jerome Asselin [EMAIL PROTECTED] To: Feng Zhang [EMAIL PROTECTED]; R-Help [EMAIL PROTECTED] Sent: Thursday, August 14, 2003 11:52 AM Subject: Re: [R] How to

[R] X11 not configured on linux RedHat 9

2003-08-14 Thread Justin Fay
I'm using Linux version 2.4.20-8smp (gcc version 3.2.2, RedHat 9). I installed R-1.7.1-1.src.rpm following INSTALL. However, no X11 device is found when I run R. X11() Error in X11() : X11 is not available In the configuration, I found checking for X... no How do I configure so that X11

RE: [R] Integer precision etc.

2003-08-14 Thread Ted Harding
Thanks to James Holtman for the confirmation of the IEEE definition, and to Marc Schwartz and Roger Koenker for pointing out .Machine which I had not been aware of! For the latter, the information I wanted is in .Machine$double.digits [1] 53 so that the largest integer exactly represented is

Re: [R] Regexpr with .

2003-08-14 Thread Stephen C. Upton
Trevor, The . is a regex meta-character that matches any character. In order to look specifically for a ., the you must escape it with a \, and that \ must also be escaped, thus, regexpr(\\., Female.Alabama) [1] 7 attr(,match.length) [1] 1 HTH steve Thompson, Trevor wrote: I'm trying to

Re: [R] vectorization question

2003-08-14 Thread Alberto Murta
Thank you very much. I just would expect that 'as.matrix' would have the same behaviour as 'data.matrix' when all columns in a data frame are numeric. Regards Alberto On Thursday 14 August 2003 16:41, Liaw, Andy wrote: If you look at the structure, you'll see: x$V4 - 0 str(x)

[R] stars graphs

2003-08-14 Thread Tito de Morais Luis
Hi listers, A few days ago I posted a question about the use of the stars function on selected lines of a frame. Thanks to two helpers, a closer look at the scale argument allowed to partially solve the problem. Yet I still have a problem with stars. Allow me to explain what I intend to do

RE: [R] placing labels in polygon center ?

2003-08-14 Thread Liaw, Andy
This is rather simple-minded: rot - function(x, k=1) { k - k %% length(x) x[c((k+1):length(x), 1:k)] } Andy -Original Message- From: Richard A. O'Keefe [mailto:[EMAIL PROTECTED] Sent: Thursday, August 14, 2003 1:28 AM To: [EMAIL PROTECTED] Subject: Re: [R] placing labels in

Re: [R] A question on orthogonal basis vectors

2003-08-14 Thread John Fox
Dear Fred, If I understand correctly what you want, the answer is not unique. Think about the 3D case where you start with one vector. (I assume, by the way, that you mean orthonormal and that you mean unique up to a reflection.) There are infinitely many pairs of orthonormal basis vectors for

Re: [R] vectorization question

2003-08-14 Thread Tony Plate
From ?data.frame: Details: A data frame is a list of variables of the same length with unique row names, given class `data.frame'. Your example constructs an object that does not conform to the definition of a data frame (the new column is not the same length as the old columns).

[R] Removing variables

2003-08-14 Thread Aurora Torrente
Hi all, Is there any way to remove at the same time several variables that share a suffix, for example (similar to * or ? in DOS) ? Thanks in advance, Aurora __ [EMAIL PROTECTED] mailing list

[R] leave-one-out

2003-08-14 Thread array chip
Hi, is there a package for performing leave-one-out cross validation in R? __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

RE: [R] Mapping function dependencies

2003-08-14 Thread Martin Maechler
TL == Thomas Lumley [EMAIL PROTECTED] on Wed, 13 Aug 2003 07:16:43 -0700 (PDT) writes: TL On Wed, 13 Aug 2003, David Khabie-Zeitoune wrote: With regards to my question below, I've found a way to do it: all.names(body(test)) TL This also gives all the variable

[R] please advise re: data mining in Germany

2003-08-14 Thread Lisa Solomon
Our CEO, Dr. Dan Steinberg, is planning to visit Germany in September. He would like the opportunity to introduce statisticians (and statistically minded people) to data mining, data mining applications and to forefront data mining tools. Our algorithms are probably familiar to many

Re: [R] Regexpr with .

2003-08-14 Thread Jeff Gentry
I'm trying to use the regexpr function to locate the decimal in a character string. Regardless of the position of the decimal, the function returns 1. You need to escape it. gsub(\\.,,,Female.Alabama) [1] Female,Alabama __ [EMAIL PROTECTED]

Re: [R] How to copy and paste a R plot onto Word (or PowerPoint)

2003-08-14 Thread James MacDonald
Either ?savePlot and import using Import -- Picture -- File in Word/PowerPoint or right click, copy as..., and paste in Word/Power point Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109

[R] Binary file reading problem

2003-08-14 Thread Jonathan Williams
I am trying to read data from a binary file. The file has a header of mixed data modes:- Bytes Length Mode 1-4 4 Integer 5-6 2 integer 7-8 2 integer .. 15-18 4 integer 19-20

Re: [R] Newbie

2003-08-14 Thread Tobias Verbeke
On Thu, 7 Aug 2003 14:41:51 +0100 Marcos Llobera [EMAIL PROTECTED] wrote: I just installed R's RPM on a Linux box using red hat 8.0. I would like to use some graphical interface (e.g. gnome) but I have not managed to do it. type R at the prompt in a gnome terminal and you can use R.

[R] error message in fitdistr

2003-08-14 Thread vincent . stoliaroff
Hi R lovers Here is a numerical vector test test [1] 206 53 124 112 92 77 118 75 48 176 90 74 107 126 99 84 114 147 99 114 99 84 99 99 99 99 99 104 1 159 100 53 [33] 132 82 85 106 136 99 110 82 99 99 89 107 99 68 130 99 99 110 99 95 153 93 136 51 103 95 99

Re: [R] Importing Data

2003-08-14 Thread TyagiAnupam
These functions are in the library foreign. You may not have loaded the library. Load it first, before using the functions. I tried the following on R1.7.1 on Windows; the functions are available. library(foreign) help(read.xport) help(read.table) In a message dated 8/6/03 11:48:56 AM

Re: [R] Standard error of standard deviation: bootstrap or theoreticalresults?

2003-08-14 Thread Prof Brian Ripley
On Wed, 6 Aug 2003, Thomas W Blackwell wrote: Huan - The difference between the empirical (bootstrap') result and the theoretical results shows evidence for autocorrelation in the time series data. I don't think that's where the correlation is (and for positive autocorrelation I would

Re: [R] Sorting a dataframe

2003-08-14 Thread Spencer Graves
Why won't order solve your problem? ?order in R 1.7.1 for Windows includes an example sorting lexicographically 3 variables in ascending order. spencer graves Paul, David A wrote: Undoubtedly a simple question: I've looked at order() and sort() in the help pages for R1.7.1. It doesn't

Re: [R] Plot ticks and tick labels: thickness, colour?

2003-08-14 Thread Prof Brian Ripley
On Wed, 6 Aug 2003, Thomas W Blackwell wrote: You can certainly make tickmarks thinner than the axis line by multiple calls to axis() with different values for lwd. MAYBE you can overwrite an earlier call by setting col.axis=white (and no tickmarks) but I've never tried this. mtext()

[R] Thanks: batch package install advice

2003-08-14 Thread Al Piszcz
Thanks for the help, example follows. pkg-c(gregmisc,e1071,car,foreign,tree) install.packages(pkg,destdir=/tmp) __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Filtering Time Series / apply STL method

2003-08-14 Thread Prof Brian Ripley
No graph was shown: R-help strips most attachments. On Wed, 6 Aug 2003, Jan Verbesselt wrote: Dear R Helpers, Which technique can I apply to get the noise out of the following graph? The data displayed shows the seasonal variation of NDVI (vegetation photosynthesis) extracted from the

[R] bwplot colours

2003-08-14 Thread Federico Calboli
Dear All, is any way I can change the colours of the box and umbrella of a bwplot without having to go to the length of: box.rectangle-trellis.par.get(box.rectangle) box.rectangle$col-black trellis.par.set(box.rectangle, box.rectangle) etc... but straight from the call: bwplot(y ~ x | z,

RE: [R] capturing output from Win 98 shell

2003-08-14 Thread Simon Fear
Prof Brian Ripley writes: If you insist on using a long-obselete OS, please help us continue ot support it by supplying patches and workarounds. I think we (like Microsoft) are going to have to think seriously about withdrawing support from non-NT-based versions of Windows fairly soon.

[R] A little problem

2003-08-14 Thread Unternährer Thomas, uth
Hi, Does anybody knows a easy way (without for-loops, maybe with something like match) to solve this problem: x - rep(1,3) y - c(0,1,0,1,0,1,1,1,0,0,0,0,1,0) if (x is a part of y){ find out where it is and do something } Thanks a lot Thomas __

Re: [R] Importing Data

2003-08-14 Thread Uwe Ligges
Hadassa Brunschwig wrote: Im trying to import data from an excel sheet or a sas file to R...im not succeeding. Apparently the function read.xport for reading a SAS file doesnt exist. What do i have to type in EXACTLY to read from an excel sheet(i guess i would be using read.table?)? -

Re: [R] think i found it, problems with lda ,

2003-08-14 Thread Marc Feldesman
Your problem is in the use of the underscore character. Replace it with something else. In your version of R, it is synonymous with the - operator. On Wed, 6 Aug 2003 18:04:56 +0200 amazing electrons exploded from the fingers of Nicolaas Busscher [EMAIL PROTECTED]: i think i found my

Re: [R] Marginal (type II) SS for powers of continuous variablesin a linear model?

2003-08-14 Thread Spencer Graves
I'm confused. Consider the following example: Df - data.frame(x=1:9, y=rep(c(-1,1), length=9)) anova(lm(y~x, Df)) Analysis of Variance Table Response: y DfSum Sq Mean Sq F value Pr(F) x 1 2.861e-34 2.861e-34 2.253e-34 1 Residuals 78.88891.2698

Re: [R] Comprehesive Package/Library list?

2003-08-14 Thread A.J. Rossini
Al Piszcz [EMAIL PROTECTED] writes: Is there a list of all contributed R libraries available through CRAN? Ideally it would include a one or two line description. I am looking for a packages() command similar to library() but that would access the CRAN repository and provide a listing of

[R] new version of Rcmdr package

2003-08-14 Thread John Fox
Dear list members, I've uploaded a new version of the Rcmdr package to CRAN. There are many additions and (I hope) improvements, as indicated in the relevant portion of the CHANGES file, reproduced below. As usual, comments, suggestions, and bug reports are appreciated. John

Re: [R] Alignment when rotating text and symbols

2003-08-14 Thread p.b.pynsent
Dear Prof Ripley, Many thanks for this response. I do not understand how your example addresses my problem. Basically what concerns me is the relationship between pos=1 and pos=3 upon rotation, if I take your example and modify it slightly to be at srt=180 I would expect to get an inverted

Re: [R] high memory allocation

2003-08-14 Thread Andrew C. Ward
Dear Georg, You may prefer to try clara() which uses less memory than the other cluster routines (and stands for clustering large applications). The documentation for clara() says that all variables must be numeric, which could be a problem if you have nominal or ordinal variables. Regards,

Re: [R] Regexpr with .

2003-08-14 Thread Chuck Cleland
Thompson, Trevor wrote: I'm trying to use the regexpr function to locate the decimal in a character string. Regardless of the position of the decimal, the function returns 1. For example, regexpr(., Female.Alabama) You probably want backslashes to indicate that . should not be treated as a

Re: [R] gregmisc

2003-08-14 Thread atsuya fujito
Thank you, Marc; I installed car, R commander and Bioconductor. I am using R 1.7.1, on MAC OSX 10.2.6 library() Packages in library '/usr/local/lib/R/library': AnnBuilder Bioconductor annotation data package builder Biobase Biobase: Base functions for Bioconductor

Re: [R] L10N and i18n of R

2003-08-14 Thread M.Kondrin
Shigeru Mase wrote: Nakama's patches can also handle other languages than Japanese. He kindly built Korean and Russian version of R (1.7.1) as rpm binaries (although he can understand neither Korean nor Russian). It is interesting, but with Russian R works fine as is (can say nothing about

[R] gnls - Step halving....

2003-08-14 Thread Christian Mora
Hi all, I'm working with a dataset from 10 treatments, each treatment with 30 subjects, each subject measured 5 times. The plot of the dataset suggests that a 3-parameter logistic could be a reasonable function to describe the data. When I try to fit the model using gnls I got the message 'Step

Re: [R] problems with lda , data included, can somebody test withthe new version

2003-08-14 Thread Ko-Kang Kevin Wang
On Wed, 6 Aug 2003, Nicolaas Busscher wrote: enclosed a simple R script (and a data file, and the output) , with calls lda similar to the example with the iris data in the documentation. it is not working and i dont understand the error message. can anybody help me? i am using R 1.5.1

Re: [R] size of a pixmap image

2003-08-14 Thread Roger Bivand
On 7 Aug 2003, Tito de Morais Luis wrote: Dear listers, I loaded a pixmap image with: brobo - read.pnm(brobo.pnm) I can have the characteristics of the image with: brobo Pixmap image Type : pixmapRGB Size : 609x682 Resolution: 1x1 Bounding box : 0 0

Re: [R] Plot ticks and tick labels: thickness, colour?

2003-08-14 Thread Thomas W Blackwell
Dirk - You can certainly make tickmarks thinner than the axis line by multiple calls to axis() with different values for lwd. MAYBE you can overwrite an earlier call by setting col.axis=white (and no tickmarks) but I've never tried this. mtext() allows building custom tick labels. In

[R] subscripts in lists

2003-08-14 Thread Chris Knight
I am tying myself in knots over subscripts when applied to lists I have a list along the lines of: lis-list(c(a,b,next,want1,c),c(d, next, want2, a)) From which I want to extract the values following next in each member of the list, i.e. something along the lines of answer-c( want1, want2).

Re: [R] gls function

2003-08-14 Thread Douglas Bates
Prodromos Zanis [EMAIL PROTECTED] writes: I use the gls function but in contrast to the lm function in which when I type summary(lm(...))$coef I receive all the coefficients (estimate, Std. Error, t-value and pvalue), with gls when I type summary(gls(...))$coef I only receive the estimate of

[R] Newbie

2003-08-14 Thread Marcos Llobera
I just installed R's RPM on a Linux box using red hat 8.0. I would like to use some graphical interface (e.g. gnome) but I have not managed to do it. Marcos [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

Re: [R] R like SPSS

2003-08-14 Thread Jonathan Baron
On 08/06/03 13:02, Gerhard Prade wrote: Hello all, (Sorry, my bad english) i am searching a alternativ for spss that i can use with linux. I studied sociology in germany and work sometimes for some hours in the marketresearch. There they use statistics like anova and crosstabs. I see that r

[R] code speed help? -- example and results provided

2003-08-14 Thread Debruicker, Paul Andrew
I have the following piece of code that combines lists comprised of components of varying length into a list with components of constant length. I have found 2 ways to do it, and the faster of the two is posted below along with sample results. Do you have any suggestions on how to decrease

[R] Levene test of homogeneity of variance

2003-08-14 Thread Haynes, Maurice (NIH/NICHD)
Has the Levene test of homogeneity of variance been implemented in any library in R? Thanks, Maurice Haynes National Institute of Child Health and Human Development Child and Family Research Section 6705 Rockledge Drive Bethesda, MD 20892 Voice: 301-496-8180 Fax: 301-496-2766 E-Mail: [EMAIL

Re: [R] problem with download.packages

2003-08-14 Thread Prof Brian Ripley
On Wed, 6 Aug 2003, David Parkhurst wrote: R Version 1.7.0, under windows XP pro The help page for update.packages says that { `download.packages' takes a list of package names and a ^^^ destination directory, downloads the newest versions of

Re: [R] problems with lda , data included,can somebody test with the new version

2003-08-14 Thread Philippe Glaziou
Nicolaas Busscher [EMAIL PROTECTED] wrote: enclosed a simple R script (and a data file, and the output) , with calls lda similar to the example with the iris data in the documentation. it is not working and i dont understand the error message. I tested your data and commands without any

[R] problem with download.packages

2003-08-14 Thread David Parkhurst
R Version 1.7.0, under windows XP pro The help page for update.packages says that { `download.packages' takes a list of package names and a destination directory, downloads the newest versions of the package sources and saves them in `destdir'. If the list of available packages is not given as

RE: [R] na.action in randomForest --- Summary

2003-08-14 Thread Liaw, Andy
To make it clear: Version 3.3 and older of Breiman's code do not handle NAs at all: you need to exclude them before running random forest. One can easily do so by using na.action=na.omit (which, despite the default, seemed to be what's used). The airquality data used in the examples in the

Re: [R] ginv vs. solve

2003-08-14 Thread Peter Dalgaard BSA
Angel [EMAIL PROTECTED] writes: Why do x-b%*%ginv(A) and x-solve(A,b) give different results?. It seems that I am missing some basic feature of matrix indexing. e.g.: A-matrix(c(0,-4,4,0),nrow=2,ncol=2) b-c(-16,0) x-b%*%ginv(A);x x-solve(A,b);x [ginv() is from MASS, please remember

[R] Startup for RExcel

2003-08-14 Thread David Khabie-Zeitoune
Hi When I start an instance of an R stats server using the R(D)COM RExcel Addin by Thomas Baier, it appears to start the server in a temporary directory. Is there any way I can change the directory in which the server starts in so that I can, for example, take advantage of Rprofile or Renviron

Re: [R] Weighted SUR, 2SLS regressions

2003-08-14 Thread TyagiAnupam
An (the?) early reference for SUR is: * Arnold Zellner (1962). An efficient method for estimating seemingly unrelated regressions, and tests for aggregation bias, Journal of the American Statistical Association, Vol 57, pp348--68. * For a discusssion, generalization and related references see,

RE: [R] subscripts in lists

2003-08-14 Thread Liaw, Andy
sapply(lis, function(x) x[which(x == next) + 1]) [1] want1 want2 HTH, Andy From: Chris Knight I am tying myself in knots over subscripts when applied to lists I have a list along the lines of: lis-list(c(a,b,next,want1,c),c(d, next, want2, a)) From which I want to extract the

Re: [R] X11 not configured on linux RedHat 9

2003-08-14 Thread Marc Schwartz
On Wed, 2003-08-13 at 12:53, Justin Fay wrote: I'm using Linux version 2.4.20-8smp (gcc version 3.2.2, RedHat 9). I installed R-1.7.1-1.src.rpm following INSTALL. However, no X11 device is found when I run R. X11() Error in X11() : X11 is not available In the configuration, I found

Re: [R] Manipulating text

2003-08-14 Thread Peter Dalgaard BSA
Uwe Ligges [EMAIL PROTECTED] writes: Dennis Fisher wrote: In labeling axes, I want to combine symbols and text/superscripts. Examples include: m2 (m, followed by a superscripted 2) * µg (micrograms) How can I accomplish this in R? See

Re: [R] Princomp function in R

2003-08-14 Thread Ko-Kang Kevin Wang
It is in package mva. On Tue, 12 Aug 2003, dai wrote: Date: Tue, 12 Aug 2003 13:51:10 -0700 From: dai [EMAIL PROTECTED] To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Subject: [R] Princomp function in R Hi, I want to use Princomp function in R, but com up an error as Error: couldn't

[R] tsdiag and tsStructure for np,ns,nt and nl determination

2003-08-14 Thread Jan Verbesselt
Hi R-Helpers, I'm dealing with the STL procedure and trying to apply the tsdiag and StructTS onto the ts object to analyse the different parameters which need to be set. How can I use the tsStructure tsdiag to create a seasonal, trend and cycle subseries plot so that I can select analyse the

Re: [R] Problems with addition in big POSIX dates

2003-08-14 Thread Martin Maechler
Whit == Whit Armstrong [EMAIL PROTECTED] on Wed, 13 Aug 2003 15:23:58 -0400 writes: Whit Have you noticed any problems with big dates (=1/1/2040) in R? Whit Here is the bit of code that I'm having trouble with: test.date - strptime(1/1/2040,format=%m/%d/%Y)

Re: [R] error message in sem

2003-08-14 Thread John Fox
Dear Remko, Can you supply the input that produced this error message? John At 04:50 PM 8/11/2003 -0700, Remko Duursma wrote: Dear R-helpers, i get a mysterious error message when using sem, it reads: Error in tapply(grad.P[arrows.2.free], ram[ram[, 1] == 2 ram[, 4] != : arguments

Re: [R] Changing default browser in options()

2003-08-14 Thread Tobias Verbeke
so I have to type every single R session: options(browser='dillo') Is there anyway I can change this globally? If you put the line options(browser='dillo') in the .Rprofile of your home directory, it will fire up dillo the next time. HTH, Tobias

[R] Startup for RExcel

2003-08-14 Thread David Khabie-Zeitoune
Hi When I start an instance of an R stats server using the R(D)COM RExcel Addin by Thomas Baier, it appears to start the server in a temporary directory. Is there any way I can change the directory in which the server starts in so that I can, for example, take advantage of Rprofile or Renviron

[R] ^-operation exponentiation problem

2003-08-14 Thread Alexander Singer
Hi, There seems to be a problem with the power-operation of the form c(..., -a, ...)^x.yyy Four code examples for the same calculation, which perform differently: -3^3.2 Result: [1] -33.63474 c(-3,6)^3.2 Result: [1] NaN 309.0893 But: c(-3.62,3)^3 Result: [1] -47.43793 27.0 I can

[R] default directory RGui for windows NT

2003-08-14 Thread Roy Werkman
Hello, Can anyone tell me how to change the default directory in RGui for windows NT? Thanx, Roy __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Re: www.acd.ucar.edu mailing list memberships reminder

2003-08-14 Thread Philippe Glaziou
Daniel Bayard [EMAIL PROTECTED] wrote: how do you add x and y error bars on a plot. Have queried the search engine? help.search(error bar) -- Philippe Glaziou Pasteur Institute of Cambodia __ [EMAIL PROTECTED] mailing list

Re: [R] placing labels in polygon center ?

2003-08-14 Thread Barry Rowlingson
Richard A. O'Keefe wrote: I wrote: I found myself wishing for a function to rotate a vector. Is there one? I know about ?lag, but help.search(rotate) didn't find anything to the point. Here I was regarding a vector as a _sequence_. The (one-step) rotation of c(u,v,w,x,y,z)

Re: [R] What does m$... mean?

2003-08-14 Thread Prof Brian Ripley
On Tue, 12 Aug 2003, Spencer Graves wrote: a$b = a[['b']] = attribute b of list a. (Not quite always. First, it is `component' not `attribute' and second $ and [[ ]] do behave differently, e.g. for data frames in 1.7.x.) A basic object in R is a list, and the $ operator provides one means

[R] Comprehesive Package/Library list?

2003-08-14 Thread Al Piszcz
Is there a list of all contributed R libraries available through CRAN? Ideally it would include a one or two line description. I am looking for a packages() command similar to library() but that would access the CRAN repository and provide a listing of the current libraries, and version.

Re: [R] function lme

2003-08-14 Thread Jason Turner
Hadassa Brunschwig wrote: Thanks to everyone who replied to my first problem. Here is the second one: Is the function lme (Mixed Model) also in a foreign library yes. When you encounter this sort of problem, try help.search(lme) This will tell you which installed library has the function.

Re: [R] default directory RGui for windows NT

2003-08-14 Thread Ko-Kang Kevin Wang
On Mon, 11 Aug 2003, Roy Werkman wrote: Can anyone tell me how to change the default directory in RGui for windows NT? I take it you mean the working directory? At least two ways: 1) Right click on the Rgui shortcut, click on Properties. Enter the path into the Starting Location (or

RE: [R] error message in fitdistr

2003-08-14 Thread Ted Harding
On 05-Aug-03 [EMAIL PROTECTED] wrote: Here is a numerical vector test test [1] 206 53 124 112 92 77 118 75 48 176 90 74 107 126 99 84 ... 89 94 69 74 99 97 91 92 Assuming it follows a lognormal distribution I'd like to determine the mean and the sd thanks to maximum

Re: [R] Marginal (type II) SS for powers of continuous variables ina linear model?

2003-08-14 Thread Bjørn-Helge Mevik
Prof Brian Ripley [EMAIL PROTECTED] writes: drop1 is the part of R that does type II sum of squares, and it works in your example. So does Anova in the current car: I'm sorry, I should have included an example to clarify what I meant (or point out my misunderstandings :-). I'll do that

Re: [R] How to get the pseudo left inverse of a singular squarematrix?

2003-08-14 Thread Gabor Grothendieck
As already pointed out ZZ' is not invertible. That is, it is not one-to-one and onto. What we can do is restrict the domain and range of ZZ' to the range of ZZ' or equivalently to the range of Z. In operational terms this means if y = ZZ'x then we only consider y's expressable as y = Zu for

[R] economic woes?

2003-08-14 Thread badams4
hiya! check this: it's very special.only the banks know about it.. I personally couldnt have got out of the mess I was in without this site --- http://btrack.iwon.com/r.pl?redir=http://[EMAIL PROTECTED]/index.asp?RefID=198478 --- dont want any more?

Re: [R] code speed help? -- example and results provided

2003-08-14 Thread Jerome Asselin
It doesn't seem that you need data frame to do this since x1 and x2 are always numeric. I think your final data frame object has N rows. This should be much faster. n - 5 N - 10 life - c(3,2,7,8,10) ind - N-life matrix(rep(rep(0:1,n),c(rbind(life,ind))),ncol=N,byrow=T) BTW, it would have been

RE: [R] evaluating and walking in names

2003-08-14 Thread Ted Harding
On 06-Aug-03 Cezar Augusto de Freitas Anselmo wrote: Hi, all. Suppose I have an object with names (like a data.frame) and I want to walk in a loop with your names. How can I do this? The idea is like this: my.data-data.frame(matrix(runif(6),ncol=2)) names(my.data) [1] X1 X2 for(i in

Re: [R] The desire for prettier and cheaper images ...

2003-08-14 Thread Martin Maechler
ChRi == Ritter, Christian C MCIL-CTGAS [EMAIL PROTECTED] on Wed, 13 Aug 2003 11:49:03 +0200 writes: ChRi This might be easy (but I didn't find an answer in the archives). ChRi I'm trying to make nice looking images (using image()). ChRi To make them look nice (not jagged), it

[R] (no subject)

2003-08-14 Thread Nolwenn Le Meur
Hi everybody, Hope your are not all on holyday because I've got a problem that is going to drive me crazy... I would like to remove some rows from a dataframe. The rows correspond to some specific indexes which I can get by looking at the name in the first column of my dataset. But I manage to

[R] Problems with addition in big POSIX dates

2003-08-14 Thread Whit Armstrong
Have you noticed any problems with big dates (=1/1/2040) in R? Here is the bit of code that I'm having trouble with: test.date - strptime(1/1/2040,format=%m/%d/%Y) unlist(test.date) sec min hour mday mon year wday yday isdst 0 0 0 1 0 140 0 0 0

Re: [R] Replacing underscore character in Windows GUI

2003-08-14 Thread kjetil brinchmann halvorsen
On 12 Aug 2003 at 17:38, nels.tomlinson.1 wrote: I thought one of the reasons for removing the underscore as synonym to - was that at some later point it could be reintroduced in the grammar, for example as a valid part of variable names. If this is the case, this suggestion is a bad idea.

[R] pb of importation data: truncation...

2003-08-14 Thread Pinaud David
Dear all, I'm importing some data (with read.table, file .txt) but I have some problem: R truncates my numeric data at 2 digits during the importation... It's the first time that I see that. Thanks for an answer. David __ [EMAIL PROTECTED] mailing list

Re: [R] newbie question: fitting power law

2003-08-14 Thread Douglas Bates
Thomas Hoffmann [EMAIL PROTECTED] writes: Dear R-helpers, dear Andrew Andrew C. Ward schrieb: Dear Thomas, I wonder if you have any NAs in xy$x or xy$y. The Dataframe looks like the following: xy x y 1 3.0 121.50951 2 2.0 30.61258 3 4.0 323.14837 4 8.2

RE: [R] capturing output from Win 98 shell

2003-08-14 Thread Adaikalavan Ramasamy
Why should it take months to get these installed ? Both R and emacs are simple and free to install. If high speed internet access is a problem, you can always download these from some Internet Café or a friend, burn it into CD and install onto the machines. If you have red tapes on installing

RE: [R] capturing output from Win 98 shell

2003-08-14 Thread Simon Fear
Peter, thanks for quick reply, but I can't drop it: system(dir /b, intern=T, show.output.on.console=T) Error in system(dir /b, intern = T, show.output.on.console = T) : dir not found From: Peter Dalgaard BSA [mailto:[EMAIL PROTECTED] snip Does it help if you drop the COMMAND.COM

[R] Regexpr with .

2003-08-14 Thread Thompson, Trevor
I'm trying to use the regexpr function to locate the decimal in a character string. Regardless of the position of the decimal, the function returns 1. For example, regexpr(., Female.Alabama) [1] 1 attr(,match.length) [1] 1 In trying to figure out what was going on here, I tried the below

[R] Integer precision etc.

2003-08-14 Thread Ted Harding
Hi Folks, With a bit of experimentation I have determined (I think) that on my R implementation the largest positive integer that is exactly represented is (2^53 - 1), based on (((2^53)-1)+1) - ((2^53)-1) [1] 1 ((2^53)+1) - (2^53) [1] 0 System: platform i686-pc-linux-gnu arch i686 os

[R] converting GLM from SAS to R

2003-08-14 Thread Monder, Harvey
Dear List members, I am in the process of translating a number of SAS programs into R scripts. The translation is surprisingly straight-forward. However, I want to convert some PROC GLM code to glm in R. Unfortunately my understanding of the procedure has become rust through the

Re: [R] re: two dimentional hierarchical clustering algorithm

2003-08-14 Thread liping
Dear Dr. Liaw Andy: I have a few more questions about your heatmap function. actually heatmap is what I am looking for. heatmap(x, Rowv, Colv, distfun = dist, hclustfun = hclust, add.expr, scale=c(row, column, none), na.rm = TRUE, ...) my data is a XNEW, dim(XNEW) [1] 554 335

  1   2   3   >