Re: [R] Attempting to recode elements contained in a list

2005-08-24 Thread Prof Brian Ripley
On Wed, 24 Aug 2005, Greg Blevins wrote: > Hello R-Masters, > > I have a list 's' with three elements, as shown below. I want to recode > a.a, a.a2, and a.a3 to NA if the value in a.a is less than 3. I think s is a data frame, and you are looking at its printout. There is almost certainly no

Re: [R] Packages

2005-08-24 Thread Uwe Ligges
Roger D. Peng wrote: > What version of R are you using? 'modreg' no longer exists---its > functions have been moved to the 'stats' package. Let me add that the URL cited below is not an "official" one. Plaese take a look at some CRAN mirror close to you. Uwe Ligges > -roger > > Amir Safari

[R] Attempting to recode elements contained in a list

2005-08-24 Thread Greg Blevins
Hello R-Masters, I have a list 's' with three elements, as shown below. I want to recode a.a, a.a2, and a.a3 to NA if the value in a.a is less than 3. I reivewed my Modern Applied Statistic Book, the online help and did some searching of R-help on the internet. I explored unlist and as.list

Re: [R] How to do log normal regression?

2005-08-24 Thread Spencer Graves
1. Have you considered "nls" or "optim"? The documentation for both includes useful examples. 2. What do your Y values represent? The almost universal standard today is maximum likelihood estimation. If you tell us what the Y values are, someone might help you write a l

Re: [R] GLM/GAM and unobserved heterogeneity

2005-08-24 Thread Spencer Graves
Have you considered "lmer" in library(lme4)? See for example sec/ 4 pm "Two-level models for binary data" in vignette("MlmSoftRev") wiht library(mlmRev) in addition to www.r-project.org -> "Documentation: Newsletter" -> "R News Volume 5/1" -> "Fitting Linear Mixed Models in R" by Dou

[R] question about custom contrasts in ANOVA

2005-08-24 Thread Scot W McNary
Hi, I have a problem in which I have test score data on students from a number of schools. In each school I have a measure of whether or not they received special programming. I am interested in the interaction between school and attendance to the programming, but in a very select set of co

Re: [R] barplots - text direction

2005-08-24 Thread Patrick Connolly
On Wed, 24-Aug-2005 at 01:54PM +1200, David Scott wrote: |> On Wed, 24 Aug 2005, Murray Pung wrote: |> |> > If the variable names are too long to allow room for each to be displayed on a barplot, how can the direction of the text be changed? |> > |> > a <- cbind(2.4,2.4,2.5,2.6,2.6,2.6,2.6,2.6,2

[R] Panel regression in R

2005-08-24 Thread Haibo Huang
I am currently trying to replicate the results I got from RATS for a panel regression. The codes in RATS looks like this: * Final equation for Office Cap Rate Spread * Regression, Panel Data preg(effects=time, method=random) CapRate # CapRate{1} RentCycle{1} VacancyChangeYTY InflationYTY RealGDPy

Re: [R] How to collect better estimations of a logistic model parameters, by using bootstrapping things ?

2005-08-24 Thread Prof Brian Ripley
I presume an `IC' is a confidence interval, but what is an `IP'? I think you need to think hard about the assumptions you want to make. The usual way to do logistic regression is via glm(), and the confint() function will give you confidence intervals based on profile likelihoods that are rath

[R] How to collect better estimations of a logistic model parameters, by using bootstrapping things ?

2005-08-24 Thread Laurent Valdes
Dear all, I know that when using R, people should have a sufficient level in statistics. As well, I'm not a genius, when dealing with logistic regressions. I would like to construct ICs, IPs, for a logistic regression, but the point is I have just 41 observations. I had a look at the Design

[R] Lacunarity

2005-08-24 Thread Rich Grenyer
Hello R-help, Does anyone know if the sliding-box lacunarity analysis of Allain & Cloitre [1] / Plotnick et al. [2] has been implemented by anyone in R? I've come up blank on searches of help files, mailing lists and from googling r-project.org. Thanks in advance, Rich [1] C. Allain and M. C

[R] summary: Call SAS from R

2005-08-24 Thread Shin, David
Thanks for all the insightful responses that help me solve my problem. The original question is: how to call SAS from R The solution is as below: If my sas code, test.sas, is in directory c:\ncme06\sas codes, then the following codes works > system('"c:\\program files\\sas institute\\sas\\v8\\sa

Re: [R] Packages

2005-08-24 Thread Roger D. Peng
What version of R are you using? 'modreg' no longer exists---its functions have been moved to the 'stats' package. -roger Amir Safari wrote: > > Hi Dear All, > I want to instal some packages, but I can not find them in both lists of > Package Index in internet ( > http://www.maths.lth.se/h

Re: [R] Remove NAs from Barplot

2005-08-24 Thread Doran, Harold
Thanks, Marc. It works. My apologies for not making the code reproducible. -Original Message- From: Marc Schwartz (via MN) [mailto:[EMAIL PROTECTED] Sent: Wednesday, August 24, 2005 1:29 PM To: Doran, Harold Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Remove NAs from Barplot On Wed,

[R] Packages

2005-08-24 Thread Amir Safari
Hi Dear All, I want to instal some packages, but I can not find them in both lists of Package Index in internet ( http://www.maths.lth.se/help/R/.R/doc/html/packages.html ) and of Instalation of Packages window in the softwar. For example "modreg". Are they included in larger packages ? Th

Re: [R] Call SAS from R

2005-08-24 Thread Nordlund, Dan
David, 1. I assume you are working on a MS Windows OS given your system command. 2. You need to determine what command will work to run your program from the command line in Windows. If your path is correct, then you can run SAS from the command line like "c:\program files\sas institute\

Re: [R] Call SAS from R

2005-08-24 Thread Prof Brian Ripley
Date: Wed, 24 Aug 2005 18:05:44 +0100 (BST) From: Prof Brian Ripley <[EMAIL PROTECTED]> To: David Forrest <[EMAIL PROTECTED]> Cc: "Shin, David" <[EMAIL PROTECTED]>, r-help@stat.math.ethz.ch, 'Don MacQueen' <[EMAIL PROTECTED]>, [EMAIL PROTECTED] Subject: Re: [R] Call SAS from R On Wed, 24 Aug

Re: [R] Remove NAs from Barplot

2005-08-24 Thread Marc Schwartz (via MN)
On Wed, 2005-08-24 at 13:06 -0400, Doran, Harold wrote: > Dear List: > > I'm creating a series of barplots using Sweave that must assume a > standard format. This is student achievement data and the x-axis must > include all grades 3 to 8. In some cases, the data for a grade (or more > than one gr

Re: [R] plotting GAM

2005-08-24 Thread Sean O'Riordain
Emil, I've no clue what gam() does nor how it works, but saying that is it possible to use lines(,lwd=???) in someway and plot on top of the original line? cheers, Sean On 24/08/05, Tkadlec Emil <[EMAIL PROTECTED]> wrote: > Sear Sean, > here is the more working example (I use the latest version

[R] Remove NAs from Barplot

2005-08-24 Thread Doran, Harold
Dear List: I'm creating a series of barplots using Sweave that must assume a standard format. This is student achievement data and the x-axis must include all grades 3 to 8. In some cases, the data for a grade (or more than one grade) are missing in the vector math.bar, but are never missing for t

Re: [R] Call SAS from R

2005-08-24 Thread David Forrest
On Wed, 24 Aug 2005, Shin, David wrote: ... > > system("c:\\program files\\sas institute\\v8\\sas.exe test") > Warning message: > c:\program not found Escape the spaces too. It is trying to run the program c:\program or c:\progra~1\sas Dave -- Dr. David Forrest [EMAIL PROTECTED]

Re: [R] Call SAS from R

2005-08-24 Thread Shin, David
Thanks for the responses. Here is the command I used and the error message I got. > system("c:\\program files\\sas institute\\v8\\sas.exe test") Warning message: c:\program not found if I change "program files" to "progra~1" then the output is: > system("c:\\progra~1\\sas institute\\v8\\sas.exe

Re: [R] Call SAS from R

2005-08-24 Thread Don MacQueen
How does it not work? Do you get error messages? Nothing at all? (it's difficult to help without more specific information, so this is a good time to suggest reading the posting guide: http://www.R-project.org/posting-guide.html) What happens if you try to run that exact command outside of R? Do

Re: [R] Call SAS from R

2005-08-24 Thread Shin, David
Thanks for some kind responses from the list. Maybe I should present my data here a little bit to explain why I want to use SAS. Below is part of my data set. First, I need to skip the first 9 lines which I can use skip=9 to do this. Second, I need to read to rows of data. the first row starts w

Re: [R] Call SAS from R

2005-08-24 Thread vincent
perhaps it would be easier to read directly your datafile from R ? read.table(...) is your friend, and quite easy to use. type ?read.table under R. hih __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read t

[R] Call SAS from R

2005-08-24 Thread Shin, David
Hi All, I am new to post question on this list. I apologize if this question is too easy or irrelevant. I am doing a simulation study and I need to read a data file that can be easily read by SAS. So, what I try to do is to execute SAS in R and then read the output of SAS to R. I try the followi

[R] studentship till January 2007

2005-08-24 Thread Ramon Diaz-Uriarte
Please pass along, and apologies for double posting. We have money to support a student till January 2007 (with a salary of about 1000 euro/month). Most of the work will focus on classification/prediction using microarray data. The work will involve both methodological research (mainly computat

Re: [R] How to get the list of the files when you read zip file by gzfile or unz

2005-08-24 Thread Prof Brian Ripley
On Wed, 24 Aug 2005, Shawn Lee wrote: > Does somebody know how to get the filename lists in the ziped when > "gzfile" or "unz" command is used ? There is no list when gzfile() is used: a gzipped file is not a zip archive. unz() opens a single specified file in a zip archive. AFAIK R does not di

[R] [R-pkgs] For2R - Interface from Fortran to R

2005-08-24 Thread Mike Prager
We announce the availability of the Fortran 95 module "For2R", which makes it easy to write complicated R data structures (i.e., lists of lists, matrices, dataframes, ...) from Fortran programs. By a series of subroutine calls, program outputs can be written to a file readable by R with a singl

[R] How to get the list of the files when you read zip file by gzfile or unz

2005-08-24 Thread Shawn Lee
Dear List, Does somebody know how to get the filename lists in the ziped when "gzfile" or "unz" command is used ? Thanks for your help. Shawn Lee -- The information contained in this communication and any atta...{{dropped}} __ R-help@stat.math.e

[R] The betareg package don't supporte large data file

2005-08-24 Thread justin bem
Hi I use betareg package to estimate part of food expenditure in household. In have a data frame of 10992 cases with 5 variables. I have this message "could not allow a vector of 43783 ko" In have use a laptop with - Intel PIV at 2.66GHZ - 512 MO of RAM I use set memory >memory.limit(size=1

[R] lm.ridge

2005-08-24 Thread daniel
Hello, I have posted this mail a few days ago but I did it wrong, I hope is right now: I have the following doubts related with lm.ridge, from MASS package. To show the problem using the Longley example, I have the following doubts: First: I think coefficients from lm(Employed~.,data=longley) sho

Re: [R] help with unknown function

2005-08-24 Thread Mathieu Basille
Hello. Sorry for the delay, (not so) long vacations... Seems to me that it has something to deal with the 'adehabitat' package which is very useful for radio-tracking data. ('adehabitat' is available on CRAN) BUT the 'getareahr' function doesn't exist. Could it be a mix between 'getverticeshr

Re: [R] correlation by subset

2005-08-24 Thread Krishna
Hi Mr. Cleland Could you please detail the function. Thanks and best regards snvk On 8/24/05, Chuck Cleland <[EMAIL PROTECTED]> wrote: > by(mydata, mydata$GROUPING, function(x) cor(x$x, x$y, use="pair")) > > ?by > > Krishna wrote: > > Hi all > > > > Having searched the available documentation

[R] (Fwd) Re: priority of operators in the FOR ( ) statement

2005-08-24 Thread Ravi . Vishnu
Dear List, Thanks to all those who responded to my questions. I now know how the R experts reason, which is good to know when faced with the a problem in R. I also became aware of R commands like syntax, browser, debug, traceback, xml etc. I certainly aim to try out these commands. Some have urg

Re: [R] correlation by subset

2005-08-24 Thread Chuck Cleland
by(mydata, mydata$GROUPING, function(x) cor(x$x, x$y, use="pair")) ?by Krishna wrote: > Hi all > > Having searched the available documentation on R, I request for help > in sorting out the underlying problem. > > I have a huge dataset containing 2 variables x and y, which is a daily > price ser

Re: [R] randomize

2005-08-24 Thread Uwe Ligges
Luis Tercero wrote: > Dear helpeRs, > > I am looking for an R equivalent of the S function "randomize" but > cannot find it anywhere (?randomize, help.search("randomize")). There > are some references to "Randomize()" in the list archives, but it is not > on my system (R 2.1.1 on Windows XP).

Re: [R] Wilcoxon Signed Rank Test

2005-08-24 Thread Torsten Hothorn
> Dear list, > > Please forgive me if this is a dumb question. I want to compare a set of > paired data (pre and late). There are ties in pre and late. I searched > on line, some document says use wilcox.exact in R instead of > wilcox.test. Anyone has any experience using wilcoxon signed rank test

Re: [R] randomize

2005-08-24 Thread Sean Davis
On 8/24/05 7:26 AM, "Luis Tercero" <[EMAIL PROTECTED]> wrote: > Dear helpeRs, > > I am looking for an R equivalent of the S function "randomize" but > cannot find it anywhere (?randomize, help.search("randomize")). There > are some references to "Randomize()" in the list archives, but it is not

[R] histogram method for S4 class.

2005-08-24 Thread ernesto
Hi, I'm trying to develop an histogram method for a class called "FLQuant" which is used by the package FLCore (http://flr-project.org). FLQuant is an extension to "array". There is an as.data.frame method that coerces flquant into a data.frame suitable for lattice plotting. The problem is that wh

[R] (Fwd) Re: priority of operators in the FOR ( ) statement

2005-08-24 Thread Petr Pikal
Hi On 23 Aug 2005 at 12:03, [EMAIL PROTECTED] wrote: > Dear All, > I spent an entire evening in debugging a small, fairly simple program > in R - without success. It was my Guru in Bayesian Analysis, Thomas > Fridtjof, who was able to diagonose the problem. He said that it took > a long

[R] correlation by subset

2005-08-24 Thread Krishna
Hi all Having searched the available documentation on R, I request for help in sorting out the underlying problem. I have a huge dataset containing 2 variables x and y, which is a daily price series. I would like to observe the quarterly correlations among these two variables. Is there anyway wh

[R] randomize

2005-08-24 Thread Luis Tercero
Dear helpeRs, I am looking for an R equivalent of the S function "randomize" but cannot find it anywhere (?randomize, help.search("randomize")). There are some references to "Randomize()" in the list archives, but it is not on my system (R 2.1.1 on Windows XP)... what package is it in? Thanks

[R] [Fwd: Re: reading "special" text files]

2005-08-24 Thread Jan Wiener
-- Dr. Jan Malte Wiener Max-Planck-Institute for Biological Cybernetics Spemannstr. 38, 72076 Tuebingen, Germany tel.: +49 7071 601 631 email: [EMAIL PROTECTED] url: www.kyb.tuebingen.mpg.de/~malte --- Begin Message --- well, the problem is that my txt file looks like this 1 1 2

Re: [R] reading "special" text files

2005-08-24 Thread Prof Brian Ripley
On Wed, 24 Aug 2005, Jan Wiener wrote: > hello, > > assume you have a txt-file (sep="\t" or "") looking as follows that you > want to read into a variable: Into a variable or a data frame? > 1 1 2 3 > 2 1 2 3 4 > 3 1 5

Re: [R] reading "special" text files

2005-08-24 Thread Uwe Ligges
Jan Wiener wrote: > hello, > > assume you have a txt-file (sep="\t" or "") looking as follows that you > want to read into a variable: > > 1 1 2 3 > 2 1 2 3 4 > 3 1 5 6 7 8 > > the result should look like

[R] reading "special" text files

2005-08-24 Thread Jan Wiener
hello, assume you have a txt-file (sep="\t" or "") looking as follows that you want to read into a variable: 1 1 2 3 2 1 2 3 4 3 1 5 6 7 8 the result should look like this (e.g., in a data.frame): 1

[R] plotting gam

2005-08-24 Thread Tkadlec Emil
Hi Sean, my working example is simple. I am modeling a response variable, say y, using two quantitative predictors, say x1 and x2: model<-gam(y~s(x1)+s(x2) Then, I plot the model using plot(model, residuals=T,se=T) If I set par(mfrow=c(1,2)) before that, I get two nice figures for both predictors

Re: [R] plotting GAM

2005-08-24 Thread Sean O'Riordain
Hi Emil, can you give us a working example of what you're trying to do? cheers! Sean ps. as per the posting-guide... :-) On 24/08/05, Tkadlec Emil <[EMAIL PROTECTED]> wrote: > Dear colleagues, > I would like to have GAM regression lines (package gam) thicker than the > default > setting does. I

[R] plotting GAM

2005-08-24 Thread Tkadlec Emil
Dear colleagues, I would like to have GAM regression lines (package gam) thicker than the default setting does. Is there any way to change the width of regression line when plotting gam.objects from the package GAM with more than one predictor? Changing lwd parameter in plot function controls

Re: [R] priority of operators in the FOR ( ) statement

2005-08-24 Thread Petr Pikal
__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] Problem with RCMD build ...

2005-08-24 Thread Martin Maechler
> "Carlos" == Carlos J Gil Bellosta <[EMAIL PROTECTED]> > on Tue, 23 Aug 2005 23:51:34 +0200 writes: Carlos> Dear Justin, Carlos> I also had similar problems recently... In fact, I have just created a Carlos> "package" using package.skeleton and if I try to build it with