Re: [R] Installing R on Red Hat Linux, Tcl/Tk

2003-09-05 Thread Marc Schwartz
On Fri, 2003-09-05 at 21:37, Dirk Eddelbuettel wrote: > On Fri, 2003-09-05 at 07:13, James Wettenhall wrote: > > I've been trying to install R on Red Hat Linux 9 for some > > potential users of my R/TclTk application. I tried using the > > rpm for R 1.7.1 for Red Hat Linux 9. It told me that I n

Re: [R] Installing R on Red Hat Linux, Tcl/Tk

2003-09-05 Thread Dirk Eddelbuettel
On Fri, 2003-09-05 at 07:13, James Wettenhall wrote: > I've been trying to install R on Red Hat Linux 9 for some > potential users of my R/TclTk application. I tried using the > rpm for R 1.7.1 for Red Hat Linux 9. It told me that I needed: > libtcl8.3.so > so I looked for a binary release of Tc

[R] Using subset with describe in Hmisc

2003-09-05 Thread Greg Blevins
Hello, This is an embarassingly simple question, but I cannot get subset to work with describe. dataframe is attached and called Ph >describe(q1, subset=qs3a==1) qs3a is numeric. This runs, but no subset takes place. Thanks in advance. Greg Blevins The Market Solutions Group [[alter

RE: [R] R Documentation with Emphasis on Econometrics

2003-09-05 Thread len vir
It's been published in: Journal of Applied Econometrics 17 : 175-189 (2002) Bart Benthul <[EMAIL PROTECTED]> wrote: Short and basic: http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/2001/10-01.php -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Fr

Re: [R] fit data with skew t distribution

2003-09-05 Thread kjetil brinchmann halvorsen
On 5 Sep 2003 at 16:20, tzhou1 wrote: I don't know what a skew t distribution is, but help.search("skew t") gives dmst(sn)Multivariate skew-t distribution dst(sn) Skew-t Distribution dst2.plot(sn) Plot of Bivariate Skew-t Density Function msn.affine(sn)

Re: [R] Dotchart question

2003-09-05 Thread kjetil brinchmann halvorsen
On 5 Sep 2003 at 14:30, Peter Flom wrote: prof <- c('Accountant', 'Administrative assistant', 'Garment worker', 'Cook', 'Dentist', 'General practictioner', 'Graduate student', 'High level manager', 'Low level manager', 'Mechanical engineer', 'Mechanic', 'Minister/priest/rabbi', 'Nurse'

[R] Right formula for a mixed-design anova

2003-09-05 Thread Rafael Laboissiere
Hi, I am having some trouble to find the correct syntax in R for analysing my results. I am running a experiment in which we measure the force exerted by the subjects (the independent variable) under two different fixed factors. Each factor has two levels, so that each subject participated to sev

[R] fit data with skew t distribution

2003-09-05 Thread tzhou1
Hi, Is there a function in R that I can use to fit the data with skew t distribution? Speaking in detail, I first used the kernel density estimation to fit my data, then I drew the skew t using my specified location, scale, shape, and df to make it close to the kernel density. Now I want to get

Re: [R] stack overflow

2003-09-05 Thread Peter Dalgaard BSA
William Noble <[EMAIL PROTECTED]> writes: > Hello, > > I am trying to do an ANOVA on a microarray data set consisting of > 22690 elements. The ANOVA is fine, but when I try to put the data in > a frame in order to exporting it, I get a stack overflow. I have > found documentation on dynamic mem

Re: [R] Dotchart question

2003-09-05 Thread h . wickham
Hi Peter, Dotchart plots the points in the order that it receives them - to plot them in order of increasing mol, you need to sort mol first. To make sure you keep the right prof associated to the right mol, I'd use names to connect the two first (this also means you don't need to provide an expl

Re: [R] stack overflow

2003-09-05 Thread Thomas W Blackwell
Bill - Here's what I would do, starting after your display of anovaresults[[1]]. temp.1 <- unlist(lapply(anovaresults, function(x) { x["Pr(>F)"][1:3],] })) temp.2 <- matrix(temp.1, length(anovaresults), 3, byrow=T) dimnames(temp.2) <- list(names(anovaresults), dimnames(

[R] stack overflow

2003-09-05 Thread William Noble
Hello, I am trying to do an ANOVA on a microarray data set consisting of 22690 elements. The ANOVA is fine, but when I try to put the data in a frame in order to exporting it, I get a stack overflow. I have found documentation on dynamic memory in R, but not on how to increase the stack size.

Re: [R] Dotchart question

2003-09-05 Thread Tito de Morais Luis
This one is ugly, but works... molprof<-data.frame(mol,prof) molprof <- molprof[with(molprof,order(mol)), ] dotchart(molprof$mol, labels = as.character(molprof$prof), main = 'Dot chart', xlab = 'Meaning of life score') HTH Tito Le ven 05/09/2003 à 18:30, Peter Flom a écrit : > Sorry to keep as

Re: [R] laplace transform

2003-09-05 Thread Ravi Varadhan
Hi: I have written a Fortran program based on the Gaver-Stehfest algorithm, which uses only real numbers (as opposed to the more powerful methods using complex numbers). However, this can't be used in R since the function specifying the inverse of the Laplace transform must also be written in

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Ted Harding
On 05-Sep-03 Deepayan Sarkar wrote: > > Oops, forgot the attachment. > > On Friday 05 September 2003 12:45 pm, Deepayan Sarkar wrote: >> The prepanel function returns separate limits for x and y axes. This >> does not translate to splom, since each limit is used on both the x >> and y axes. >> Ho

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread Spencer Graves
One must estimate 2 coefficients for a 3-level factor. I therefore prefer to look for a plausible order between the 3 levels so I can code them as -1, 0, +1 and then estimate linear and quadratic coefficients. Then for interactions, I look first for interactions between the linear effect of

[R] Dotchart question

2003-09-05 Thread Peter Flom
Sorry to keep asking elementary questions..I appreciate the help. I am trying to create a dotchart with the rows sorted according to the values, rather than the labels. When I try prof <- c('Accountant', 'Administrative assistant', 'Garment worker', 'Cook', 'Dentist', 'General practictione

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Ted Harding
Hi Deepayan, Thanks again, especially for the careful explanation. Some comments (in snipped context) below: On 05-Sep-03 Deepayan Sarkar wrote: > On Friday 05 September 2003 10:04 am, Ted Harding wrote: >> Another query: >> >> I'm now trying to have the x- and y-axes all on the same scale >> (0:1

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread Rolf Turner
In response to a question from Francisco J. Bido, about how to create dummy variables, Doug Bates and others essentially said ``Don't.'' Which is good advice, but Recently I encountered a problem involving a linear model with a three level factor (levels low, medium, and high) crossed with l

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Deepayan Sarkar
Oops, forgot the attachment. On Friday 05 September 2003 12:45 pm, Deepayan Sarkar wrote: > The prepanel function returns separate limits for x and y axes. This does > not translate to splom, since each limit is used on both the x and y axes. > However, it is natural to add a new optional argume

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Deepayan Sarkar
On Friday 05 September 2003 10:04 am, Ted Harding wrote: > Another query: > > I'm now trying to have the x- and y-axes all on the same scale > (0:15) in every panel, whereas the default behaviour of splom > is to scale them according to the ranges of the individual > variables in each panel. > > So

RE: [R] R Documentation with Emphasis on Econometrics

2003-09-05 Thread Bart Benthul
Short and basic: http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/2001/10-01.php -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Francisco J. Bido Sent: Friday, September 05, 2003 10:52 AM To: [EMAIL PROTECTED] Subject: [R] R Documentation with Emphasi

[R] R Documentation with Emphasis on Econometrics

2003-09-05 Thread Francisco J. Bido
Hi Everyone, There's a lot of good documentation out there but I haven't come across one that makes econometrics it's primary focus. Do any of you know of any? Thanks, -Francisco __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread kjetil brinchmann halvorsen
On 5 Sep 2003 at 8:31, Francisco J. Bido wrote: Yes, model matrix is the answer, and if it has many arguments, it also has many reasonable defaults. When I am trying out a new function, I just accept the dafaults for a starter. > > x <- c(1:10) #data to be broken up into dummy variables

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread Francisco J. Bido
Thanks Douglas, I see what your saying. One of the reasons that I ask this question (besides being a complete R rookie) is to obtain good form and habits. R seems to be extremely capable and flexible (and of course also pretty dense) thank God for the mailing list. Your example and the feedb

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread Peter Dalgaard BSA
Douglas Bates <[EMAIL PROTECTED]> writes: > > x <- c(1:10) #data to be broken up into dummy variables > > v <- c(3,5,7) #breakpoints > > p = 1#drop this column to avoid dummy variable trap > > > > How can I get a matrix "y" that has the associated dummy variables for > >

Re: [R] eliminating a large subset of data from a frame

2003-09-05 Thread Uwe Ligges
Peter Whiting wrote: I have a data frame with 155,000 rows. One of the columns represents the user id (of which about 10,000 are unique). I am able to isolate 1000 of these user ids (stored in a list) that I want to eliminate from the data set, but I don't know of an efficient way to do this. Cer

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread Douglas Bates
"Francisco J. Bido" <[EMAIL PROTECTED]> writes: > Hi There, > > While looking through the mailing list archive, I did not come across > a simple minded example regarding the creation of dummy variables. > The Gauss language provides the command "y = dummydn(x,v,p)" for > creating dummy variables.

RE: [R] eliminating a large subset of data from a frame

2003-09-05 Thread Simon Fear
I can't see what would be wrong with newdf <- subset(df, uid %in% listofgoodusers) or newdf <- subset(df, !(uid %in% listofbadusers)) Is this what you want? Please note the code you supplied will not run at all, let alone slowly, so it is not easy to know exactly what you are trying to achieve >

[R] [Off] Ecological Modelling's Book Suggestion.

2003-09-05 Thread Ronaldo Reis Jr.
Hi, Anybody can suggest a good Ecological Modelling's Book. Thanks Ronaldo -- O céu sueco. -- palíndromo -- |> // | \\ [***] | ( õ õ ) [Ronaldo Reis Júnior] |> V [UFV/DBA-Entomologia] |/ \

Re: [R] Basic Dummy Variable Creation

2003-09-05 Thread John Fox
Dear Francisco, At 08:31 AM 9/5/2003 -0500, Francisco J. Bido wrote: Hi There, While looking through the mailing list archive, I did not come across a simple minded example regarding the creation of dummy variables. The Gauss language provides the command "y = dummydn(x,v,p)" for creating dum

[R] eliminating a large subset of data from a frame

2003-09-05 Thread Peter Whiting
I have a data frame with 155,000 rows. One of the columns represents the user id (of which about 10,000 are unique). I am able to isolate 1000 of these user ids (stored in a list) that I want to eliminate from the data set, but I don't know of an efficient way to do this. Certainly this would be s

Re: [R] all values from a data frame

2003-09-05 Thread partha_bagchi
If I understand you correctly, you want to stack the 15 columns on top of one another? I assuming all the data is numeric? In this case, convert the data.frame to a matrix and set the dim of the matrix to NULL. If df is the data.frame, df1 <- as.matrix(df) dim(df1) <- NULL <[EMAIL PROTECT

Re: [R] all values from a data frame

2003-09-05 Thread John Fox
Dear Arne, At 05:07 PM 9/5/2003 +0200, [EMAIL PROTECTED] wrote: Hello, I've a data frame with 15 colums and 6000 rows, and I need the data in a single vector of size 9 for ttest. Is there such a conversion function in R, or would I have to write my own loop over the colums? thanks for

RE: [R] all values from a data frame

2003-09-05 Thread Simon Fear
Use unlist. In any case, don't write an explicit "loop over the columns". See ?lapply instead. > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > Sent: 05 September 2003 16:08 > To: [EMAIL PROTECTED] > Subject: [R] all values from a data frame > > > Security War

RE: [R] Putting regression lines on SPLOM

2003-09-05 Thread Ted Harding
Another query: I'm now trying to have the x- and y-axes all on the same scale (0:15) in every panel, whereas the default behaviour of splom is to scale them according to the ranges of the individual variables in each panel. So I tried (emulating the responses to my earlier query): splom(log(1+

[R] all values from a data frame

2003-09-05 Thread Arne.Muller
Hello, I've a data frame with 15 colums and 6000 rows, and I need the data in a single vector of size 9 for ttest. Is there such a conversion function in R, or would I have to write my own loop over the colums? thanks for your help + kind regards Arne ___

[R] skins in R?

2003-09-05 Thread solares
Hello, it wanted to know if exists in R any package that carry out skins on the widget of tcl, rounded button style window XP. Thanks Ruben __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Two questions about building R add-on package

2003-09-05 Thread A.J. Rossini
Brian D Ripley <[EMAIL PROTECTED]> writes: > From <[EMAIL PROTECTED]> > > Windows Versions of R and R packages under Intel Linux$)A!1 by Jun > > Yan and A. J. Rossini and used their Makefile. But it kept giving error > We quite explicitly say you need to build R first. If Yun-Rossini do not

Re: [R] scatter.smooth warning "k-d tree" --> "loess bowels"

2003-09-05 Thread Martin Maechler
> "Peter" == Peter Flom <[EMAIL PROTECTED]> > on Thu, 04 Sep 2003 13:28:12 -0400 writes: Peter> Hello When I run Peter> scatter.smooth(jitter(weight), jitter(height2), span Peter>= .25, evaluation = 50, pch = '.') Peter> I get the type of graph I thoug

RE: [R] Comparison of SAS & R/Splus

2003-09-05 Thread rpugh
Sorry - couldn't resist chipping in. Firstly, this sort of conversation has been done over and over again on the S-News list, and I'd look in the archives for more info. My background: I was a SAS "statistical programmer" in the pharma industry before I joined Insightful (S-PLUS guys). I now wor

[R] multivariate spec.ar

2003-09-05 Thread Hugo Vilela Mendes
Im trying to perform an autoregressive spectral analysis to a multivariate time series in order to analyse cross spectrum. The multivariate case is not yet implemented; i would like to know if there are some package or add-on that allows to perform spec.ar analysis on multivariate time series

Re: [R] Overlaying graphs

2003-09-05 Thread Paul Meagher
From: "Damon Wischik" <[EMAIL PROTECTED]> > Paul Meagher wrote: > > 2. Does R have a suite of "best-fit" tools for finding the best > > fitting-probability distribution for any observed probability distribution? > > I think that the best-fitting probability distribution for an observed > probabili

RE: [R] Comparison of SAS & R/Splus

2003-09-05 Thread Thomas Lumley
On Fri, 5 Sep 2003, Paul, David A wrote: > > d. Because SAS is commercial software, a posteriori errors found in > clinical trials analyses (and due to software issues) can be > attributed by the NDA applicants to the SAS Institute. > Lawyers really like this. Of course, Splus is also commercial

Re: [R] Installing R on Red Hat Linux, Tcl/Tk

2003-09-05 Thread Marc Schwartz
On Fri, 2003-09-05 at 07:13, James Wettenhall wrote: > Hi, > > I've been trying to install R on Red Hat Linux 9 for some > potential users of my R/TclTk application. I tried using the > rpm for R 1.7.1 for Red Hat Linux 9. It told me that I needed: > libtcl8.3.so > so I looked for a binary rele

[R] Basic Dummy Variable Creation

2003-09-05 Thread Francisco J. Bido
Hi There, While looking through the mailing list archive, I did not come across a simple minded example regarding the creation of dummy variables. The Gauss language provides the command "y = dummydn(x,v,p)" for creating dummy variables. Here: x = Nx1 vector of data to be broken up into dummy

RE: [R] Comparison of SAS & R/Splus

2003-09-05 Thread Paul, David A
On Fri, 5 Sep 2003, Brian D. Ripley wrote: > In general I find such discussions irrelvant. > I bet those users make far, far more errors then any > of these packages do so. However, without having the discussions with my colleagues, nothing will ever change. The perception of SAS' "bestness

[R] Installing R on Red Hat Linux, Tcl/Tk

2003-09-05 Thread James Wettenhall
Hi, I've been trying to install R on Red Hat Linux 9 for some potential users of my R/TclTk application. I tried using the rpm for R 1.7.1 for Red Hat Linux 9. It told me that I needed: libtcl8.3.so so I looked for a binary release of Tcl 8.3.x on http://www.tcl.tk/software/tcltk/8.3.html but

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Ted Harding
Yes, that's it! Well spotted (and maybe new glasses for me ...). Thanks! Andy's solution now also works beautifully. Ted. On 05-Sep-03 Andrew C. Ward wrote: > I wonder if it's as simple as the two dots ("..") in the > splom line rather than three? > > Regards, > > Andrew C. Ward >> However, th

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Andrew C. Ward
I wonder if it's as simple as the two dots ("..") in the splom line rather than three? Regards, Andrew C. Ward CAPE Centre Department of Chemical Engineering The University of Queensland Brisbane Qld 4072 Australia [EMAIL PROTECTED] Quoting Ted Harding <[EMAIL PROTECTED]>: > Thanks, Deepayan

Re: [R] gnome gui question

2003-09-05 Thread Martin Maechler
> "Anne" == Anne York <[EMAIL PROTECTED]> > on Thu, 4 Sep 2003 10:41:25 -0700 (PDT) writes: Anne> I have compiled R 1.7.1 with gnome and installed the Anne> gtkDevice on a Linux RH 8 system. This is my first use Anne> of this device and I have two questions about it. A

Re: [R] Putting regression lines on SPLOM

2003-09-05 Thread Ted Harding
Thanks, Deepayan and Andy, for suggesting solutions. Deepayan's solution seems to work beautifully, producing exactly what is wanted, as it stands: On 05-Sep-03 Deepayan Sarkar wrote: > splom(log(1 + DF), panel = function(x, y, ...) { panel.xyplot(x, y, ...) ok

Re: [R] Comparison of SAS & R/Splus

2003-09-05 Thread Brian D Ripley
> On Thu, 4 Sep 2003, Paul, David A wrote: > >> I am one of only 5 or 6 people in my organization making the >> effort to include R/Splus as an analysis tool in everyday work - the >> rest of my colleagues use SAS exclusively. >> >> Today, one of them made the assertion that he believes the >> nu

Re: [R] Two questions about building R add-on package

2003-09-05 Thread Brian D Ripley
> I am a new user of R and am building an R add-on package. I followed the > "Writing R Extension" manual from cran website but still met some > problems that I cannot solve. I build it under redhat linux 9.0 R1.7.0. > Say foo is the name of the package. Do read the instructions! > First, after