Re: [R] if statement

2007-06-12 Thread Simon Blomberg
My solutions are usually too baroque, but does this do what you want? x <- rnorm(100) quants <- quantile(x, c(.3, .7)) Case <- rep(2, length(x)) # 2 lies in the middle of the distribution Case[x <= quants[1]] <- 0 Case[x >= quants[2]] <- 1 Case Cheers, Simon. On Mon, 2007-06-11 at 15:14 -0700,

Re: [R] PCA for Binary data

2007-06-12 Thread ssls sddd
Dear Prof Brian Ripley, Would you also recommend some packages for non-binary data to do variable and feature selection? Thanks a lot! Alex On 6/12/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > > On Tue, 12 Jun 2007, Spencer Graves wrote: > > > The problem with applying prcomp to bin

[R] Trouble making JRI.jar with Ubuntu and Java6

2007-06-12 Thread taivo
Hi, Forum newb here, looking for some help. Have been trying to install an R-Java interface to make R calls from Java. JRI's configure script runs fine, but when it comes to make, I get the "error: too few arguments to function 'R_ParseVector'" Java runs fine. R runs fine. But I can't get this .

[R] Panel data

2007-06-12 Thread Megh Dal
Dear all R users, I have a small doubt about panel data analysis. My basic understanding on Panel data is a type of data that is collected over time and subjects. Vector Autoregressive Model (VAR) model used on this type of data. Therefore can I say that, one of statistical tools used for analy

[R] if statement

2007-06-12 Thread Jiong Zhang, PhD
Hi all, I have a rather naive question. I have the height of 100 individuals in a table and I want to assign the tallest 30% as Case=1 and the bottom 30% as Case=0. How do I do that? thanks. jiong The email message (and any attachments) is for the sole use of the intended recipient(s) and may

Re: [R] Viewing a data object

2007-06-12 Thread Stephen Tucker
Hi Horace, I have also thought that it may be useful but I don't know of any Object Explorer available for R. However, (you may alread know this but) (1) you can view your list of objects in R with objects(), (2) view objects in a spreadsheet-like table (if they are matrices or data frames) wit

Re: [R] PCA for Binary data

2007-06-12 Thread Prof Brian Ripley
On Tue, 12 Jun 2007, Spencer Graves wrote: > The problem with applying prcomp to binary data is that it's not > clear what problem you are solving. > > The standard principal components and factor analysis models > assume that the observations are linear combinations of unobserved > "com

[R] passing (or obtaining) index or element name of list to FUN in lapply()

2007-06-12 Thread Stephen Tucker
Hello everyone, I wonder if there is a way to pass the index or name of a list to a user-specified function in lapply(). For instance, my desired effect is something like the output of > L <- list(jack=4098,sape=4139) > lapply(seq(along=L),function(i,x) if(i==1) "jack" else "sape",x=L) [[1]] [1]

Re: [R] Using dll with Visual Studio Compiler

2007-06-12 Thread Prof Brian Ripley
On Tue, 12 Jun 2007, Ian McCarthy wrote: > Hi, > > I have created a dll in Fortran and used the Visual Studio 2005 > Compiler. I've read that certain problems can arise based on the Since Visual Studio does not create Fortran, that's a bit vague. And let us assume you are on Win32, but you fai

Re: [R] Awk and Vilno

2007-06-12 Thread Christophe Pallier
On 6/13/07, Robert Wilkins <[EMAIL PROTECTED]> wrote: > > The point is : there are lots of data preparation scenarios where > large numbers of merges need to be done. This is an example where > Vilno and SAS are easier to use than the competition. I'm sure an Awk > programmer can come up with somet

[R] Design library installation problem

2007-06-12 Thread Ian Watson
Dear Listers I have tried to install Frank Harrell's two libaries: Hmisc and Design. I found that Hmisc was listed in the list of packages from the Install Packages command on the Packages menu, but Design was not. I installed Hmisc from this list, and when I issued the library(Hmisc) comma

Re: [R] PCA for Binary data

2007-06-12 Thread Simon Blomberg
You might try (detrended) correspondence analysis, which is designed for "count" data, if it makes sense to treat your binary data that way. I've used ade4 and also vegan, and they are both good packages for these types of ordinations. You could also look at non-metric multidimensional scaling. Th

Re: [R] How do I obtain standard error of each estimated coefficients in polr

2007-06-12 Thread Spencer Graves
I'm confused: Have you considered the 'examples' in the 'polr' help file? The first example ends "summary(house.plr)". The print of this summary includes standard errors. If you want those numbers for subsequent computations, you can try str(summary(house.plr)) or names(summary

Re: [R] data from graph [SEC=UNCLASSIFIED]

2007-06-12 Thread Augusto.Sanabria
Jam, If you use "fpot" the GPD of package "evd" to calculate your return periods, you can generate the table you want by using the following commands: RP <- c(10,100,500,1000) #return periods rf <- your clean vector of rainfall data thres <- the threshold for your data nopy <- number of observ

Re: [R] PCA for Binary data

2007-06-12 Thread Spencer Graves
The problem with applying prcomp to binary data is that it's not clear what problem you are solving. The standard principal components and factor analysis models assume that the observations are linear combinations of unobserved "common" factors (shared variability), normally distr

[R] Stock Price Correlation to Index Price Levels

2007-06-12 Thread Josh Kalish
Thanks to all of the people who responded. What I was trying to do is to turn my matrix or frame containing index level returns and stock returns into a matrix of "betas". I don't really need to worry about risk-free interest rates. I just need to be able to come up with a number that shows the

Re: [R] pretty report

2007-06-12 Thread Robert A LaBudde
At 09:13 PM 6/12/2007, Don wrote: >At 5:01 PM -0400 6/12/07, Weiwei Shi wrote: > >Dear Listers: > > > >I have a couple of data frames to report and each corresponds to > >different condtions, e.g. conditions=c(10, 15, 20, 25). In this > >examples, four data frames need to be exported in a "pretty"

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-12 Thread Gabor Grothendieck
Here is yet another simplification. This one uses na.locf from the zoo package to shorten it further and also make it easier to understand. Below we have one line to read in the .ini file, one line to transform the characters [ and ] to = and =, the read.table line parses the result and the next

[R] specify constraints in maximum likelihood

2007-06-12 Thread adschai
Hi,I know only mle function but it seems that in mle one can only specify the bound of the unknowns forming the likelihood function. But I would like to specify something like, a = 2b or a <= 2b where 'a' and 'b' could be my parameters in the likelihood function. Any help would be really appreci

[R] Awk and Vilno

2007-06-12 Thread Robert Wilkins
In clinical trial data preparation and many other data situations, the statistical programmer needs to merge and re-merge multiple input files countless times. A syntax for merging files that is clear and concise is very important for the statistical programmer's productivity. Here is how Vilno do

Re: [R] pretty report

2007-06-12 Thread Don MacQueen
At 5:01 PM -0400 6/12/07, Weiwei Shi wrote: >Dear Listers: > >I have a couple of data frames to report and each corresponds to >different condtions, e.g. conditions=c(10, 15, 20, 25). In this >examples, four data frames need to be exported in a "pretty" report. > >I knew Perl has some module for ex

Re: [R] PCA for Binary data

2007-06-12 Thread Josh Gilbert
I don't understand, what's wrong with using prcomp in this situation? On Sunday 10 June 2007 12:50 pm, Ranga Chandra Gudivada wrote: > Hi, > > I was wondering whether there is any package implementing Principal > Component Analysis for Binary data > >

[R] Setting a minimum number of observations within an individual cluster

2007-06-12 Thread Norm.Good
Hi I'm trying to cluster a continuous dataset with a varying number of clusters and with a restriction that each cluster must have more than 'x' number of observations. I have tried the clara function, using silhouette to give me the neighbouring cluster mediod of each observation, then mergi

Re: [R] Generating artificial datasets with a specific correlati

2007-06-12 Thread Ted Harding
On 12-Jun-07 20:54:05, Ken Knoblauch wrote: > see mvrnorm in MASS and especially the empirical argument > > James Milks wright.edu> writes: > > >> I need to create artificial datasets with specific correlation >> coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples >> for a

Re: [R] Appropriate regression model for categorical variables

2007-06-12 Thread Ted Harding
On 12-Jun-07 17:45:44, Tirthadeep wrote: > > Dear users, > In my psychometric test i have applied logistic regression > on my data. My data consists of 50 predictors (22 continuous > and 28 categorical) plus a binary response. > > Using glm(), stepAIC() i didn't get satisfactory result as > misc

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-12 Thread Gabor Grothendieck
In thinking about this a bit more here is an even shorter solution where Lines.raw is as before: # Lines <- readLines("myfile.ini") Lines <- readLines(textConnection(Lines.raw)) Lines2 <- chartr("[]", "==", Lines) DF <- read.table(textConnection(Lines2), as.is = TRUE, sep = "=", fill = TRUE) L <-

[R] Using dll with Visual Studio Compiler

2007-06-12 Thread Ian McCarthy
Hi, I have created a dll in Fortran and used the Visual Studio 2005 Compiler. I've read that certain problems can arise based on the compiler used and that these problems can sometimes cause R not to be able to access everything appropriately, but I've not found how to fix any of these problems

Re: [R] Can strptime handle milliseconds or AM/PM?

2007-06-12 Thread Dirk Eddelbuettel
On 12 June 2007 at 16:54, Earl F. Glynn wrote: | Any clues what I'm doing wrong with strptime? This seems to fail the same | way under Linux or Windows. | | For ?strptime would it make sense to explain %OS3 somewhere besides the | Examples? | | > # Why does %OS3 work here? | > format(Sys.ti

[R] data from graph

2007-06-12 Thread umarporn charusombat
hi i just learn how to use R in my research i used extreame value package to get the return level of rainfall. the output i got as the graph plot between return level and period. i wonder how i can get the value from the graph as a table format. please help me thanks jam [[alternative HTML

Re: [R] R Book Advice Needed

2007-06-12 Thread Cody_Hamilton
Alain, Can you tell us what you plan to use R for? Regards, -Cody [EMAIL PROTECTED] a écrit : > I am new to using R and would appreciate some advice on > which books to start with to get up to speed on using R. > > My Background: > 1-C# programmer. > 2-Programmed directly using IMSL (Now Visua

Re: [R] 2 Trellis graphics question

2007-06-12 Thread Deepayan Sarkar
On 6/12/07, Alan S Barnett <[EMAIL PROTECTED]> wrote: > 1) I have a data that includes some "bad" data. I want to make a > trellis plot where each panel contains > a) A scatter plot of the "good" data > b) A scatter plot of the "bad" data in a different color > c) A best fit line of all the data,

Re: [R] R Book Advice Needed

2007-06-12 Thread Alain Reymond
[EMAIL PROTECTED] a écrit : > I am new to using R and would appreciate some advice on > which books to start with to get up to speed on using R. > > My Background: > 1-C# programmer. > 2-Programmed directly using IMSL (Now Visual Numerics). > 3- Used in past SPSS and Statistica. > > I put together

Re: [R] pretty report

2007-06-12 Thread Gabor Grothendieck
Generating Excel or reports seems to be two different questions: 1. Reports. You can use xtable together with Sweave. See figure 1 in this link for an example: http://www.ci.tuwien.ac.at/~leisch/Sweave/Sweave-Rnews-2002-3.pdf 2. Regarding writing Excel files this can be done on Windows using th

[R] Can strptime handle milliseconds or AM/PM?

2007-06-12 Thread Earl F. Glynn
I'm trying to proess date/time fields from files that were given to me to analyze. Any clues what I'm doing wrong with strptime? This seems to fail the same way under Linux or Windows. For ?strptime would it make sense to explain %OS3 somewhere besides the Examples? > # Why does %OS3 work

[R] JGR and big list of packages.

2007-06-12 Thread Ronaldo Reis Junior
Hi, I have all CRAN packages installed on my Linux. Now I have problems with JGR. When I make a plot and close a device, the device dont work anymore, I nedd to use before javaGD() and after plot(). When I try do close JGR and save a session, It return a erro and dont close. Look: Exception i

[R] R Book Advice Needed

2007-06-12 Thread ngottlieb
I am new to using R and would appreciate some advice on which books to start with to get up to speed on using R. My Background: 1-C# programmer. 2-Programmed directly using IMSL (Now Visual Numerics). 3- Used in past SPSS and Statistica. I put together a list but would like to pick the "best of"

Re: [R] Subset and logical operator error

2007-06-12 Thread Ken Knoblauch
Sébastien free.fr> writes: > > Can you please point to me my syntax mistake or indicate a method to get > this type of data.frame subset ? > > Thank you in advance > > ID value > 1 1 1.2 > 2 2 1.2 > 3 3 1.2 > 4 4 1.2 > 5 5 A > 6 6 A > 7 7 A > 8 8 A > subset

[R] ML, REML and several random effects

2007-06-12 Thread Arnaud Mosnier
Hello everyone, Hope that my question could have interest for more than one people here. I know that in order to compare mixed models with different fixed effect I need to use ML method. But ... what about comparing models with the same fixed effects and different random effects ? Thanks in adv

Re: [R] pretty report

2007-06-12 Thread Weiwei Shi
Just realized append=T might help for xtable, but I think it still not very "pretty" :) On 6/12/07, Weiwei Shi <[EMAIL PROTECTED]> wrote: > Dear Listers: > > I have a couple of data frames to report and each corresponds to > different condtions, e.g. conditions=c(10, 15, 20, 25). In this > example

[R] Subset and logical operator error

2007-06-12 Thread Sébastien
Hello, It looks to me as if the ! logical operator cannot be called when subsetting a data.frame. In the example below, the value column has two factor levels (but my typical datasets have more), and what I am trying to do is to exclude all lines for which the "value" is different from "A". I

[R] pretty report

2007-06-12 Thread Weiwei Shi
Dear Listers: I have a couple of data frames to report and each corresponds to different condtions, e.g. conditions=c(10, 15, 20, 25). In this examples, four data frames need to be exported in a "pretty" report. I knew Perl has some module for exporting data to Excel and after googling, I found R

Re: [R] Generating artificial datasets with a specific correlationcoefficient.

2007-06-12 Thread Ken Knoblauch
see mvrnorm in MASS and especially the empirical argument James Milks wright.edu> writes: > I need to create artificial datasets with specific correlation > coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples > for a lab I am teaching this summer. Is there a way to do tha

Re: [R] Overlaying lattice graphs

2007-06-12 Thread hadley wickham
On 6/12/07, Seb <[EMAIL PROTECTED]> wrote: > Hello > > I apologize in advance if this question has already be posted on the > list, although I could not find a relevant thread in the archives. > > I would like to overlay xyplots using different datasets for each plot. > I typically work on the foll

[R] Generating artificial datasets with a specific correlation coefficient.

2007-06-12 Thread James Milks
I need to create artificial datasets with specific correlation coefficients (i.e. a dataset that returns r = 0.30, etc.) as examples for a lab I am teaching this summer. Is there a way to do that in R? Thanks. Jim Milks Graduate Student Environmental Sciences Ph.D. Program 136 Biological Sc

Re: [R] Overlaying lattice graphs

2007-06-12 Thread Sébastien
Thanks for the information. These commands work perfectly fine and the ?panel.superpose help was most informative. If you don't mind, I will certainly come back to you as there will certainly be additionnal complexities in my datasets that I will be unable to handle (e.g. data in more than one d

Re: [R] how to find how many modes in 2 dimensions case

2007-06-12 Thread Ravi Varadhan
Patrick, Here is an example closer to what you are interested in - locating bumps in kernel density estimator. I am using the example from package KernSmooth, using the function bkde2D(). # Another example for locating maxima in kernel density estimation data(geyser, package="MASS") x <- cbind

Re: [R] Stacked barchart color

2007-06-12 Thread Dieter Menne
hadley wickham gmail.com> writes: > > On 6/12/07, Dieter Menne menne-biomed.de> wrote: > > Dear Latticer, > > > > I want to give individual colors to all elements in a simple stacked > > barchart. I know why the example below does not work (and it is a excellent > > default), but is there any w

Re: [R] how to find how many modes in 2 dimensions case

2007-06-12 Thread Ravi Varadhan
Hi Patrick, Here is a simple R code for locating ALL the local maximum of a bivariate function, which is known on a rectangular grid. I have illustrated it with a function called the Branin function, which is commonly used as a test function in the global optimization literature. It has 6 local

Re: [R] Build Windows pkgs from source - online

2007-06-12 Thread Roger Bivand
On Tue, 12 Jun 2007, Benilton Carvalho wrote: > Hi, > > First of all, I apologize for sending out this message, as I'm sure > the answer is on the archives, but I just can't find it. > > Not long ago, there was a discussion about building Windows packages > from the source code and someone posted

Re: [R] Stacked barchart color

2007-06-12 Thread Deepayan Sarkar
On 6/12/07, Dieter Menne <[EMAIL PROTECTED]> wrote: > Dear Latticer, > > I want to give individual colors to all elements in a simple stacked > barchart. I know why the example below does not work (and it is a excellent > default), but is there any workaround for this? > > Dieter > > > # This only

Re: [R] Stacked barchart color

2007-06-12 Thread hadley wickham
On 6/12/07, Dieter Menne <[EMAIL PROTECTED]> wrote: > Dear Latticer, > > I want to give individual colors to all elements in a simple stacked > barchart. I know why the example below does not work (and it is a excellent > default), but is there any workaround for this? > > Dieter > > > # This only

Re: [R] LASSO coefficients for a specific s

2007-06-12 Thread Bruno C\.
Hy, no need to find the best s value. CV does it for you: cvres<-cv.lars(X,Y,K=10,type='lasso') sAtBest<-cvres$fraction[which.min(cvres$cv)] fits <- predict.lars(object, type="coefficients", s=sAtBest, mode="fraction") ... Ciao Bruno > Hello, > > I have a question about the lars package. I am us

Re: [R] format.dates, chron and Hmisc

2007-06-12 Thread Deepayan Sarkar
On 6/10/07, R.H. Koning <[EMAIL PROTECTED]> wrote: > Hello, I have some problems in using chron, Hmisc, and lattice. First, > using both chron and Hmisc, I get an error message when describing data: > > df$Date <- chron(df$Date,format=c("d/m/y")) > > ll <- latex(describe(df),file="..//text//df.tex

Re: [R] dyn.load( ) problem

2007-06-12 Thread Vladimir Eremeev
Ian McCarthy wrote: > > I am trying to access a dll with dyn.load, but I get an error message box > titled "R Console: Rgui.exe - Unable to Locate Component". The error > message itself states "this application has failed to start because > libifcoremdd.dll was not found. Re-installing the appli

Re: [R] Linux equivalent to windows menus and script window

2007-06-12 Thread Stefan Grosse
There is, two candidates are ESS-emacs and JGR as a gui. Have a look at: http://www.r-project.org/GUI , on Tuesday 12 June 2007 19:42:33 John Sorkin wrote: JS > R 2.5 under windows XP JS > R 2.5 under Linux FC 6 JS > JS > I currently use R under windows and when I do so, I get to work in a nice

Re: [R] Appropriate regression model for categorical variables

2007-06-12 Thread Robert A LaBudde
At 01:45 PM 6/12/2007, Tirtha wrote: >Dear users, >In my psychometric test i have applied logistic regression on my data. My >data consists of 50 predictors (22 continuous and 28 categorical) plus a >binary response. > >Using glm(), stepAIC() i didn't get satisfactory result as misclassification >r

[R] Stacked barchart color

2007-06-12 Thread Dieter Menne
Dear Latticer, I want to give individual colors to all elements in a simple stacked barchart. I know why the example below does not work (and it is a excellent default), but is there any workaround for this? Dieter # This only colors red and green, but I want blue and gray for Peatland. barcha

Re: [R] getting R2.4 (Win)

2007-06-12 Thread Chuck Cleland
Jeremy Miles wrote: > Hi, > > I would like to get hold of the R version 2.4.0 windows installer, it > doesn't seem to be available on CRAN (except the source, which needs > compiling). Does anyone know if it's still available anywhere? > > Thanks, > > jeremy Hi Jeremy: The windows binary for

[R] Viewing a data object

2007-06-12 Thread Horace Tso
Dear list, First apologize that this is trivial and just betrays my slothfulness at the keyboard. I'm sick of having to type a long name just to get a glimpse of something. For example, if my data frame is named 'AuroraStochasticRunsJune1.df" and I want to see what the middle looks like, I hav

[R] Build Windows pkgs from source - online

2007-06-12 Thread Benilton Carvalho
Hi, First of all, I apologize for sending out this message, as I'm sure the answer is on the archives, but I just can't find it. Not long ago, there was a discussion about building Windows packages from the source code and someone posted a link to a website to which we could submit the sou

[R] Linux equivalent to windows menus and script window

2007-06-12 Thread John Sorkin
R 2.5 under windows XP R 2.5 under Linux FC 6 I currently use R under windows and when I do so, I get to work in a nice windows environment complete with a very useful menu bar that includes a FILE, EDIT, MISC, PACKAGE, WINDOWS, and HELP menus. I have begun to use R under LINUX. When I run und

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-12 Thread Gabor Grothendieck
Here is some code. It replaces [ and ] with = sign and reads the result into a data frame, DF. DF2 is similar except the section is now in V3. DF3 is like like DF2 except sections are carried forward and finally we remove the rows which only had sections. Lines.raw <- "[Section1] var1=value1 var2

[R] Appropriate regression model for categorical variables

2007-06-12 Thread Tirthadeep
Dear users, In my psychometric test i have applied logistic regression on my data. My data consists of 50 predictors (22 continuous and 28 categorical) plus a binary response. Using glm(), stepAIC() i didn't get satisfactory result as misclassification rate is too high. I think categorical varia

[R] Cause of error message in cov function?

2007-06-12 Thread Matthew Keller
Hi all, I have written a script in R that simulates genetically informative data - it is posted on my website and available to the public. This is my first time to write a script for use by others and am learning that it isn't as easy as it seems. To the issue. My script runs fine on my machine a

[R] dyn.load( ) problem

2007-06-12 Thread Ian McCarthy
I am trying to access a dll with dyn.load, but I get an error message box titled "R Console: Rgui.exe - Unable to Locate Component". The error message itself states "this application has failed to start because libifcoremdd.dll was not found. Re-installing the application may fix this problem."

[R] getting R2.4 (Win)

2007-06-12 Thread Jeremy Miles
Hi, I would like to get hold of the R version 2.4.0 windows installer, it doesn't seem to be available on CRAN (except the source, which needs compiling). Does anyone know if it's still available anywhere? Thanks, jeremy -- Jeremy Miles Learning statistics blog: www.jeremymiles.co.uk/learni

[R] Stock Price Correlation to Index Price Levels

2007-06-12 Thread Kalish, Josh
Hi, This is probably trivial to most people out there, but I'm struggling with this. I have a data set which contains the closing prices (properly adjusted for dividends and splits) for several hundred securities and the closing prices for a general stock market index (S&P 500). I have no prob

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-12 Thread ngottlieb
Earl: .Ini files are, for lack of a better description, ancient. There are old windows functions such as GetProfileString. However you will have to make reference to load these from the windows Kernel.dll. Probably not worth the effort to code really old things as .ini files. >From what I see of

[R] Pareto Distribution

2007-06-12 Thread livia
I would like to fit a Pareto Distribution and I am using the following codes. First, I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the "fitted" turns out to be the same value. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit="c"

[R] Read Windows-like .INI files into R data structure?

2007-06-12 Thread Earl F. Glynn
I need to process some datasets where the configuration information was stored in .INI-like files, i.e., text files with sections like this: [Section1] var1=value1 var2=value2 [Section2] A=value3 B=value4 ... >From Google and other searches I haven't found any package, or function within a pac

Re: [R] [OT]Web-Based Data Brushing

2007-06-12 Thread hadley wickham
On 6/12/07, Roy Mendelssohn <[EMAIL PROTECTED]> wrote: > I apologize for the off-topic post, but my Google search did not turn > up much and I thought people on this list my have knowledge of this. > I am looking for examples of data brushing (i.e. dynmaic linked > plots) either on a web site, or

[R] nlme model

2007-06-12 Thread Daniel O'Shea
I am having trouble figuring out the right form for the nlme arguments. I do have examples in Modern and Applied Statistics with S and from other sources, but I still can't figure it out. I am trying to estimate species richness (sr) in streams across minnesota. My predictor variables are de

[R] [OT]Web-Based Data Brushing

2007-06-12 Thread Roy Mendelssohn
I apologize for the off-topic post, but my Google search did not turn up much and I thought people on this list my have knowledge of this. I am looking for examples of data brushing (i.e. dynmaic linked plots) either on a web site, or in a web-based application, such as an AJAX app. Eve

[R] 2 Trellis graphics question

2007-06-12 Thread Alan S Barnett
1) I have a data that includes some "bad" data. I want to make a trellis plot where each panel contains a) A scatter plot of the "good" data b) A scatter plot of the "bad" data in a different color c) A best fit line of all the data, and d) A best fit line of the "good" data. I tried using xyplo

[R] LASSO coefficients for a specific s

2007-06-12 Thread Seungho Huh
Hello, I have a question about the lars package. I am using this package to get the coefficients at a specific LASSO parameter s. data(diabetes) attach(diabetes) object <- lars(x,y,type="lasso") cvres<-cv.lars(x,y,K=10,fraction = seq(from = 0, to = 1, length = 100)) fits <- predict.lars(obje

Re: [R] Building packages with subroutine in fortran 90 under windows xp

2007-06-12 Thread Prof Brian Ripley
The tools do not include an F9x compiler: you need to edit MkRules to specify the path to one (after installing it) (and you did last March). If I do that, I can install and check packages such as oc and wnominate from CRAN that make use of F9x. On Tue, 12 Jun 2007, Cinzia Viroli wrote: > Hel

Re: [R] how to ignore error messages?

2007-06-12 Thread Jianping Jin
Dear Xiaohui, Thanks a lot for your help! It worked for me. Best regards! JP- --On Monday, June 11, 2007 9:25 AM -0700 Xiaohui <[EMAIL PROTECTED]> wrote: > see ?try > > Jianping Jin wrote: >> Dear group: >> >> I wrote a code to iterate a non-linear fit with a set of data. The >> entire proce

[R] Building packages with subroutine in fortran 90 under windows xp

2007-06-12 Thread Cinzia Viroli
Hello, I work under windows xp and I am trying to build a R package with a subroutine written in fortran 90. I have installed all the updated tools and I am working with R-2.4.0 or R-2.5.0. When I check a package with a subroutine in fortran 77 (and extension f) everything is ok. When I try

[R] distribution graph

2007-06-12 Thread Wentzel-Larsen, Tore
The following gives two functions for producing distribution graphs: distribution-graph produces a single graph, and multiple.distribution.graph produces a number of graphs side by side. Regards, Tore Wentzel-Larsen statistician Centre for Clinical research Armauer Hansen house Haukeland Univ

Re: [R] Package update announcements

2007-06-12 Thread Martin Henry H. Stevens
Hi Martin, Maybe we should appeal to each others' vanity ("Just think how many hundreds---even thousands---of people will use your package once they know about it") Hmmm. Maybe that would have the opposite of the intended effect. Cheers, Hank On Jun 12, 2007, at 3:16 AM, Martin Maechler w

Re: [R] Problems with Vista, R 2.5.0 and function save

2007-06-12 Thread Prof Brian Ripley
So you are running R somewhere in your file system where you do not have permission to write. Did you create a shortcut with a working directory set as the rw-FAQ advised you to? On Tue, 12 Jun 2007, "Maja Schröter" wrote: Hi everyone, I want to make use of the save function but it did not

[R] hedonic

2007-06-12 Thread christoph . bischoff
I would like to compute a direct hedonic prize-index, calld Time-Dummy-Variable-Method. Is this method implmented in R? and further indirect methods? Best regards [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list ht

Re: [R] Lines in dotchart & dotplot ?

2007-06-12 Thread John Kane
Thanks again. Both layouts look very usefull and certainly a lot better than I was getting on my own. --- Deepayan Sarkar <[EMAIL PROTECTED]> wrote: > On 6/11/07, John Kane <[EMAIL PROTECTED]> wrote: > > > > --- [EMAIL PROTECTED] wrote: > > > > > On 6/9/07, John Kane <[EMAIL PROTECTED]> wrote

[R] Problems with Vista, R 2.5.0 and function save

2007-06-12 Thread Maja Schröter
Hi everyone, I want to make use of the save function but it did not work. I'm using vista and R 2.5.0, winzip is installed too. Here's the code (from example ?save): > x <- runif(20) > y <- list(a = 1, b = TRUE, c = "oops") > save(x, y, file = "xy.Rdata") Fehler in gzfile(file, "wb")

[R] bitmap function in R 2.4.1 vs R 2.2.1

2007-06-12 Thread luc.jouneau
Hello, I work under windows with two versions of R : 2.4.1 and 2.2.1 Ghostscript tool is also installed (gswin32c Version 8.54). I have a small sample program: > bitmap("a.bmp") > m=(1) > plot(m) > dev.off() With R 2.2.1, it works well, but with R 2.4.1 it fails when gswin32c is launched with fo

Re: [R] Error using mgcv package

2007-06-12 Thread spime
Dear Mr. Oksanen, First of all thanks for your reply. I have solved this problem in this way. My data consists of some categorical(CAT..) predictors and also some numerical variables(NUM..) have only {0,1} 0r {0,1,2,3} values. For applying GAM i just didnot consider their splines. I had came to t

Re: [R] [ANN] Static and dynamic graphics course, July 2007, Salt Lake City

2007-06-12 Thread Berwin A Turlach
G'day Hadley, On Wed, 30 May 2007 10:57:54 +0200 "hadley wickham" <[EMAIL PROTECTED]> wrote: > We're pleased to announce a one day course covering static and dynamic > graphics using R, ggplot and GGobi. The course will be held just > before the JSM, on Saturday, 28 July 2007, in Salt Lake City.

[R] Post-hoc tests for interactions of between- and within-subject factors

2007-06-12 Thread Julius Verrel
Is there a standardized way in R to perform post-hoc comparisons for main adn interaction effects of between- and within-subject factors? For instance, I have a data set of performance of adults of different age groups (20-30, 60-70,70-80) performing a WM task (n-back, with n=1,2,3,4) in two diffe

Re: [R] Data transformation for chi-square test.

2007-06-12 Thread Christoph Buser
Dear Charlie dat <- data.frame(id = 1:10, member = c(4,2,3,5,6,2,4,6,3,4), cost = c(320,150,420,330,540,310,169,647,347,567)) dat[,"costF"] <- cut(dat[,"cost"], breaks = seq(100, 700, by=100)) table(dat[,"costF"], dat[,"member"]) This should create the table you like. Best r

Re: [R] Error using mgcv package

2007-06-12 Thread Jari Oksanen
spime gmail.com> writes: > > > Hi all, > > I need some solution in the following problem. The following error appears > when i use "mgcv" package for implementing GAM. But the same formula works > fine in "gam" package. > > > model.gam <- gam(formula = RES ~ > > CAT01+s(NUM01,5)+CAT02+CAT03+s

Re: [R] Data transformation for chi-square test.

2007-06-12 Thread Christophe Pallier
If you want to test whether ' member' has an effect on 'cost' (a continuous numerical variable), I do not recommend using a chi.square test, but rather a simple linear regression or a one-way analysis of variance. Chi.square are for categorical variables and unless you have a good reason, there is

Re: [R] Textpad help

2007-06-12 Thread Daniel Nordlund
> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > On Behalf Of Mary Royerr > Sent: Monday, June 11, 2007 8:14 PM > To: r-help@stat.math.ethz.ch > Subject: Re: [R] Textpad help > > Thx for the response John and Jim. > > I had heard before that textpad is a good edi

Re: [R] selecting characters from a line of text

2007-06-12 Thread Vladimir Eremeev
Tim Holland wrote: > > Is there a way in R to select certain characters from a line of text? I > have some data that is presently in a large number of text files, and I > would like to be able to select elements of each text file (elements are > always on the same line, in the same position) an

[R] Data transformation for chi-square test.

2007-06-12 Thread Charlie Chi
Dear all R users : I am a IT student with few statistical background and new R user for only have two month exprience. I have a data named medcost, import by read.table() as follow for example (real dataset has 500 cases), the heander id means case id, member means members in a family and cost

Re: [R] Package update announcements

2007-06-12 Thread Martin Maechler
> "MHHS" == Martin Henry H Stevens <[EMAIL PROTECTED]> > on Mon, 11 Jun 2007 17:23:46 -0400 writes: MHHS> Hi Folks, I was wondering what everyone thought about MHHS> adding a sentence to each package update announcement MHHS> that described what the package did. R extension