[R] anova model refinement/clustering question

2003-10-23 Thread Murad Nayal
Hi, I am trying to refine models of a continuous response variable and a number of categorical predictor variables. I know of some model refinement tools available in R that help in the selection of model terms like dropterm and addterm from MASS etc. However, I would also like to try to refine

Re: [R] do.call() and aperm()

2003-10-23 Thread Martin Maechler
Tony == Tony Plate [EMAIL PROTECTED] on Tue, 21 Oct 2003 10:17:05 -0600 writes: Tony I've also been thinking about how to specify that 'along' should be Tony length(dim)+1. At the moment one can specify any number from 0 up to Tony length(dim)+1, but as you point out you have

[R] Can you create a MySQL database with RMySQL?

2003-10-23 Thread Barnet Wagman
Is it possible to create a database in MySQL via RMySQL? Also, is the format for the authorization field 'userName/[EMAIL PROTECTED]'? I saw an example like this somewhere in the documenation, but I haven't found the actual specification. Thanks, Barnet Wagman

Re: [R] R-help login page error

2003-10-23 Thread Martin Maechler
Héctor == Héctor Villafuerte D [EMAIL PROTECTED] on Tue, 21 Oct 2003 17:29:27 -0800 writes: Héctor Hi all, Héctor I'm trying to access my account at Héctor https://www.stat.math.ethz.ch/mailman/options/r-help Héctor but the following appears: Héctor *Error:

[R] multivariate ARMA analysis with R

2003-10-23 Thread Vincent Spiesser
Hi, I try to make a multivariate ARMA analysis. Function arima (from ts packages only accept univariate argument. Does an R function exist for multivariate analysis ? Thanks Vincent Spiesser __ [EMAIL PROTECTED] mailing list

Re: [R] multivariate ARMA analysis with R

2003-10-23 Thread Prof Brian Ripley
Look in package bundle DSE for multiple ARMA processes. (I would say multivariate and multiple time series are different things, but guess multiple is what you want.) On Thu, 23 Oct 2003, Vincent Spiesser wrote: Hi, I try to make a multivariate ARMA analysis. Function arima (from ts

AW: [R] Can you create a MySQL database with RMySQL?

2003-10-23 Thread Christian Schulz
library(DBI) library(RMySQL) drv - dbDriver(MySQL) con - dbConnect(drv,test,root,) dbListTables(con) [1] dmjanuar dmjulidmmai model0ndata newtable [7] uebergang weka mysqlQuickSQL(con, Create table newTable( attr1 int, attr2 int)) dbListTables(con) [1] dmjanuar dmjulidmmai

Re: [R] nlm, hessian, and derivatives in obj function?

2003-10-23 Thread Arne Henningsen
Hi, I don't know much about non-linear models, but there is another possibility to fit these models: 1) get some starting values for the parameters 2) take the derivatives of the model with respect to the parameters at the point of the starting values of the parameters 3) perform a linear

[R] GIS re-mapping / polygon overlap

2003-10-23 Thread Jens Oehlschlägel
In Germany the Unemployment Agency uses a sectioning of the german map that is different from the usual Administrative Boundaries. Some demographic data are available in Administrative Boundaries only, some in Unemployment Boundaries only. I would like to generate estimates in one boundary

[R] OOP like handling of lists?

2003-10-23 Thread Thomas Petzoldt
Hello, I am writing a package with a collection of several models. In order to allow users to play interactively with the models (in contrast to hacking lengthy scripts), I want to put all what is needed to run a particular model into a single list object for each model. Then there will be a

Re: [R] GIS re-mapping / polygon overlap

2003-10-23 Thread Martin Wegmann
On Thursday 23 October 2003 13:36, Jens Oehlschlägel wrote: In Germany the Unemployment Agency uses a sectioning of the german map that is different from the usual Administrative Boundaries. Some demographic data are available in Administrative Boundaries only, some in Unemployment Boundaries

[R] R command for input/output with parallel port

2003-10-23 Thread solares
Hello, my question is,Exist a command for to make input from parallel port , i have a sensor connected to parallel port, and i need have the data in real time,i need obtain the data and put in .txt file. Thank Ruben __ [EMAIL PROTECTED] mailing list

Re: [R] OOP like handling of lists?

2003-10-23 Thread Patrick Burns
I think what you want to write is: setParms(lvmodel) - list(k1=0.5) which means that you need to define an assignment function. If you are using S3 style OOP (which you appear to be), then do something along the lines of: 'setParms-' - function(x, value) UseMethod('setParms-')

[R] Writing and running a R program

2003-10-23 Thread Subramanian Karthikeyan
Is there a way I can combine multiple lines of R commands (see below) into a little code snippet or a program in a text file, and run it in R to do my analysis? sink(mysink.txt) for (..) { code for creating a dataframe from supplied data code for doing anova from selected data } Thanks very

[R] repeating colors in graph

2003-10-23 Thread Anna Pryor
I have a plot with 13 lines and R repeats colors. Is there any way to make R not repeat colors, or is this just a fundamental limitation on R's part to say have only 8 colors that it can use? Anna __ [EMAIL PROTECTED] mailing list

Re: [R] Writing and running a R program

2003-10-23 Thread Arne Henningsen
source(mycode.R) On Thursday 23 October 2003 15:29, Subramanian Karthikeyan wrote: Is there a way I can combine multiple lines of R commands (see below) into a little code snippet or a program in a text file, and run it in R to do my analysis? sink(mysink.txt) for (..) { code for

Re: [R] repeating colors in graph

2003-10-23 Thread Uwe Ligges
Anna Pryor wrote: I have a plot with 13 lines and R repeats colors. Is there any way to make R not repeat colors, or is this just a fundamental limitation on R's part to say have only 8 colors that it can use? See ?colors and the other help topics linked in its See Also Section. Uwe Ligges

Re: [R] Writing and running a R program

2003-10-23 Thread Uwe Ligges
Subramanian Karthikeyan wrote: Is there a way I can combine multiple lines of R commands (see below) into a little code snippet or a program in a text file, and run it in R to do my analysis? sink(mysink.txt) for (..) { code for creating a dataframe from supplied data code for doing anova

[R] Stochastic dynamic programming

2003-10-23 Thread Martin Biuw
Hello all, Does anyone know of any already written functions for carrying out stochastic dynamic programming in R? Thanks! Martin -- Martin Biuw Sea Mammal Research Unit Gatty Marine Laboratory, University of St Andrews St Andrews, Fife KY16 8PA Scotland Ph: +44-(0)1334-462637 Fax:

[R] Quantreg Package

2003-10-23 Thread Radaelli Paolo - Dottorati di Ricerca
I've just installed R 1.0.8 (for Windows) and I tried to install the package Quantreg directly from Cran but it's not in the list of downlodable packages. I tried also downloading the zip file and then install it but there is an error. How can I do it? Thank you [[alternative HTML version

Re: [R] Writing and running a R program

2003-10-23 Thread Philipp Pagel
On Thu, Oct 23, 2003 at 09:29:46AM -0400, Subramanian Karthikeyan wrote: Is there a way I can combine multiple lines of R commands (see below) into a little code snippet or a program in a text file, and run it in R to do my analysis? It's quite easy: Write your R code with your favorite text

Re: [R] Quantreg Package

2003-10-23 Thread Roger Koenker
I'm not sure why quantreg isn't in the windows binaries on CRAN, perhaps Uwe can explain this...if there is a problem I'd be happy to try to help, but I don't have access to a windows machine myself. url:www.econ.uiuc.edu/~roger/my.htmlRoger Koenker email [EMAIL PROTECTED]

Re: [R] Quantreg Package

2003-10-23 Thread Uwe Ligges
Radaelli Paolo - Dottorati di Ricerca wrote: I've just installed R 1.0.8 (for Windows) and I tried to install the package Quantreg directly from Cran but it's not in the list of downlodable packages. I tried also downloading the zip file and then install it but there is an error. How can I do

[R] estimating probit models

2003-10-23 Thread Millo Giovanni
Dear all, I am looking for a convenient way to model a binary choice variable in R. Would you suggest glm() with family=binomial(link=probit))? Is there something more focused on that kind of analysis? Or am I plainly wrong? Cheers and thanx for your answers Giovanni Giovanni Millo RD Dept.

Re: [R] nlm, hessian, and derivatives in obj function?

2003-10-23 Thread M.Kondrin
Arne Henningsen wrote: Hi, I don't know much about non-linear models, but there is another possibility to fit these models: 1) get some starting values for the parameters 2) take the derivatives of the model with respect to the parameters at the point of the starting values of the parameters

[R] Re: generic algorithm

2003-10-23 Thread zhu wang
Sorry, it is my bad. I mean genetic algorithm. It seems genoud {rgenoud} is the one I am looking for. On Thu, 2003-10-23 at 09:51, zhu wang wrote: Dear all, Is there any generic algorithm code for optimization implemented in R? I searched without success. Thanks, -- Zhu Wang

Re: [R] generic algorithm

2003-10-23 Thread alessandro . semeria
Are you sure? look at nlm on base pkg and at pkgs nls nlme ... all pkgs to minimize something optimize or better, what signify optimize for you? Best A.S. Alessandro Semeria Models and Simulations Laboratory Montecatini Environmental Research Center (Edison Group),

Re: [R] Quantreg Package

2003-10-23 Thread Uwe Ligges
Radaelli Paolo - Dottorati di Ricerca wrote: I saw the read-me but I didn't undersstand wich is the problem. I only know that in a previous version of R I installed on my pc it was all ok. Yes, but on the recent version it is *not* OK.

Re: [R] OOP like handling of lists?

2003-10-23 Thread M.Kondrin
Thomas Petzoldt wrote: Hello, I am writing a package with a collection of several models. In order to allow users to play interactively with the models (in contrast to hacking lengthy scripts), I want to put all what is needed to run a particular model into a single list object for each model.

Re: [R] Re: generic algorithm

2003-10-23 Thread alessandro . semeria
Look at pkg gafit Bye A.S. Alessandro Semeria Models and Simulations Laboratory Montecatini Environmental Research Center (Edison Group), Via Ciro Menotti 48, 48023 Marina di Ravenna (RA), Italy Tel. +39 544 536811 Fax. +39 544 538663 E-mail: [EMAIL PROTECTED]

Re: [R] Weighted Clustering

2003-10-23 Thread vishal goyal
Hi, I have read some of the papers on weighted clustering but those weigths are with respect to variables (say we are clustering height measured in cms and weights in tons etc) rather than each point itself and hence they just multiply the distances with the weights. In my case, it is something

Re: [R] Weighted Clustering

2003-10-23 Thread Christian Hennig
Hi, sorry, I do not have any references, it was only an intuitive idea (but there might be references somewhere out in space). The idea is that you can indeed weight the points. For this you have to multiply the contribution of each point to the target function by the weight, and you have to

[R] Variance-covariance matrix for beta hat and b hat from lme

2003-10-23 Thread Haifeng \(Kevin\) Xie
Dear all, Given a LME model (following the notation of Pinheiro and Bates 2000) y_i = X_i*beta + Z_i*b_i + e_i, is it possible to extract the variance-covariance matrix for the estimated beta_i hat and b_i hat from the lme fitted object? The reason for needing this is because I want to have

[R] getAttrib segfaults in Windows

2003-10-23 Thread Eric Blair
Dear R people: I've written a rather extensive simulation in R, much of which I've ported to C. I did all of this in Linux, where it works great. On a Windows machine, it compiles without errors or warnings, but when it runs, it segfaults. Here are the exact lines: void init(SEXP rho){

Re: [R] repeating colors in graph

2003-10-23 Thread Barry Rowlingson
Thomas Lumley wrote: See the help pages for color() and palette(), and perhaps the RColorBrewer package. Splus lets you define colour palettes in its Gui, using a nice little notation where you specify something like black 8 white 6 red to get a 17-colour palette with 8 colours between balck

Re: [R] repeating colors in graph

2003-10-23 Thread Thomas Lumley
On Thu, 23 Oct 2003, Barry Rowlingson wrote: Thomas Lumley wrote: See the help pages for color() and palette(), and perhaps the RColorBrewer package. Splus lets you define colour palettes in its Gui, using a nice little notation where you specify something like black 8 white 6 red to get

Re: [R] generic algorithm

2003-10-23 Thread Hisaji Ono
Hi. (B (B Here is GA. I've never used this. (B (B http://www.dst.unive.it/~claudio/R/index.html#ga (B (B But this can deal with GP(genetic programming)? (B (B Cheers. (B (B__ (B[EMAIL PROTECTED] mailing list

Re: [R] repeating colors in graph

2003-10-23 Thread Barry Rowlingson
Thomas Lumley wrote: It isn't, but it works fairly well with well chosen colors (for example, linearly interpolating the RColorBrewer sequential or diverging schemes gives quite nice color ramps). It's unlikely to work for producing palettes for 13 different lines. Indeedy. Perhaps the original

Re: [R] Variance-covariance matrix for beta hat and b hat from lme

2003-10-23 Thread Douglas Bates
Haifeng \(Kevin\) Xie [EMAIL PROTECTED] writes: Given a LME model (following the notation of Pinheiro and Bates 2000) y_i = X_i*beta + Z_i*b_i + e_i, is it possible to extract the variance-covariance matrix for the estimated beta_i hat and b_i hat from the lme fitted object? Not easily.

[R] installation problems Rgraphviz

2003-10-23 Thread Yectli A Huerta
i having issues installing Rgraphviz. i installed the prereqs from the att bell labs. the latest. i'm using R 1.8. below are my error messages. any help will be appreciated. thanks, -- Yectli * Installing *source* package 'Rgraphviz' ... ** libs `Rgraphviz.so' is up to date. ** R ** inst **

RE: [R] generic algorithm

2003-10-23 Thread Liaw, Andy
Could it be that you meant genetic (rather than generic) algorithm? If so, from the package descriptions, the gafit and seao packages on CRAN might be relevant. HTH, Andy -Original Message- From: zhu wang [mailto:[EMAIL PROTECTED] Sent: Thursday, October 23, 2003 10:51 AM To:

[R] Problem w/ SJava package

2003-10-23 Thread Humphreys, Darryl
Hello, Seems like it is almost working properly. I was able to run examples calling R from Java successfully for the java code below. However, some commands, such as objects when not commented out, give me the following error? If you could point me to the correct place to find a solution, I would

[R] what's going on here with substitute() ?

2003-10-23 Thread Tony Plate
I was trying to create a function with a value computed at creation time, using substitute(), but I got results I don't understand: this.is.R Error: Object this.is.R not found substitute(this.is.R - function() X, list(X=!is.null(options(CRAN)[[1]]))) this.is.R - function() TRUE # the above

Re: [R] multivariate ARMA analysis with R

2003-10-23 Thread Jason Turner
Vincent Spiesser wrote: Hi, I try to make a multivariate ARMA analysis. Function arima (from ts packages only accept univariate argument. Does an R function exist for multivariate analysis ? Thanks Vincent Spiesser The dse bundle on CRAN has functions for VARX multivariate time series models.

Re: [R-sig-Geo] RE: [R] adjacency matrix

2003-10-23 Thread danlinyu
Hi- As far as I know, some of the error messages results from the fact that you either did not install the package or did not load the package. If you did not install the package, you can install from R's GUI (my experiences are based on the Windows Version), and load them by type

Re: [R] adjacency matrix

2003-10-23 Thread Christophe Declercq
- Original Message - From: Orlando Zacarias [EMAIL PROTECTED] To: [EMAIL PROTECTED] Cc: [EMAIL PROTECTED]; [EMAIL PROTECTED] Sent: Thursday, October 23, 2003 6:56 PM Subject: RE: [R] adjacency matrix Hi once again, Now i have been trying to run the R changed example from

[R] What is the definition of curve continuity and smoothness

2003-10-23 Thread Feng Zhang
Dear R-listers, I am now using principal curves for data analysis. The definition of it required the accurate concept of curve continuity (C0, C1, ..., CK) and curve smoothness. So if anyone can introduce the exact definition of continuous for curves and what is the most popular textbooks for

Re: [R] repeating colors in graph 2

2003-10-23 Thread Thomas Lumley
On Thu, 23 Oct 2003, Anna Pryor wrote: I've tried looking at ?colors and ?palette and if I'm understanding it correctly, I'm supposed to type in (for example) palette(rainbow(13)) before I type in my plot (of 13 lines) if I want 13 different colors. But this does not work. Other things that

Re: [R] repeating colors in graph 2

2003-10-23 Thread James MacDonald
I believe you want plot(blah blah, col=palette(rainbow(13))) HTH, 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 734-647-5623 Anna Pryor [EMAIL PROTECTED] 10/23/03 02:19PM I've tried

Re: [R] what's going on here with substitute() ?

2003-10-23 Thread Peter Dalgaard
Tony Plate [EMAIL PROTECTED] writes: Why is the body of the function X when I substituted a different expression for X? Also, given that the body of the function is X, how does the function evaluate to TRUE since X is not defined anywhere (except in a list that should have been discarded.)

pre-compiled win binaries (was: Re: [R] Quantreg Package)

2003-10-23 Thread Achim Zeileis
It happened for the second time in a week that a Windows binary for a CRAN package was not availabe (Hmisc and quantreg) although the package maintainer would have been willing to try to fix the problems that prevented automatic pre-compilation if he would have realized that there is a problem

Re: pre-compiled win binaries (was: Re: [R] Quantreg Package)

2003-10-23 Thread Roger Koenker
You are right, of course, it would be nice to have notification, but I'm also sympathetic to Uwe's situation, and not everything that could be automated, could be _easily_ automated within the constraints imposed by the rest of world. The lesson I've drawn from this is that complaints will

[R] List of lm objects

2003-10-23 Thread CENDOYA, Gabriela
Hi R-Helpers: I’m trying to fit the same linear model to a bunch of variables in a data frame, so I was trying to adapt the codes John Fox, Spencer Graves and Peter Dalgaard proposed and discused yesterday on this e-mail list: for (y in df[, 3:5]) { mod = lm(y ~ Trt*Dose, data = x, contrasts =

Re: [R] repeating colors in graph 2

2003-10-23 Thread Thomas Lumley
On Thu, 23 Oct 2003, James MacDonald wrote: I believe you want plot(blah blah, col=palette(rainbow(13))) No. palette() sets the default palette, so eg palette(rainbow(13)) plot(1:13,col=1:13) will plot in the new palette. -thomas __

Re: [R] List of lm objects

2003-10-23 Thread Spencer Graves
How about the following: AList - vector(mode=list, length=3) for(i in 1:3){ y - df[,i]) mod - lm(y~Tr*Dose, data=x, contrasts=list(Tr=contr.sum, Dose=contr.sum)) AList[[i]] - Anova(mod, type=III) } hope this helps. spencer graves CENDOYA, Gabriela wrote: Hi R-Helpers: Im trying to fit the same

Re: [R] GIS re-mapping / polygon overlap

2003-10-23 Thread Ross Ihaka
Jens Oehlschlägel wrote: AND - is there a general function in R available to map demographic data between two boundary systems OR - is there an efficient R function that can determine the amount of area overlap between two general polygons (beyond gpclib) I am willing to share

Re: [R] List of lm objects

2003-10-23 Thread Peter Dalgaard
Peter Dalgaard [EMAIL PROTECTED] writes: CENDOYA, Gabriela [EMAIL PROTECTED] writes: for (myname in names(myframe)){ mycall - substitute(lm(myvar~etc.etc.),list(myvar=as.name(myname))) myfit - eval(mycall) print(summary(myfit)) } ## by Peter Dalgaard But instead of

Re: [R] GIS re-mapping / polygon overlap

2003-10-23 Thread Hisaji Ono
Hi. How about GEOS available at. http://geos.refractions.net/ GEOS supports hole for clipping. I could build it in MiGW. However I've only executed demo codes. - Original Message - From: Ross Ihaka [EMAIL PROTECTED] To: Jens Oehlschlägel [EMAIL PROTECTED] Cc: [EMAIL PROTECTED]

[R] Generating Data Sets -- Each with a Different Name

2003-10-23 Thread Ken Kelley
Hello all. I was wondering if anyone had insight into how I might generate a large number of data sets/vectors, where each of the data sets/vectors have a different name? For example, suppose I wanted to generate N data sets/vectors, each with a different name such as: Random.Data.1 -

Re: [R] Generating Data Sets -- Each with a Different Name

2003-10-23 Thread Spencer Graves
Have you considered using paste(..., sep=) to generate either the object names or file names you want? If you want objects in R, have you considered assign? If you want files, then have you considered write? hope this helps. spencer graves Ken Kelley wrote: Hello all.

Re: [R] Generating Data Sets -- Each with a Different Name

2003-10-23 Thread apjaworski
You need paste and assign functions. Paste will create a sequence of names and assign will give values to the variables with these names. This modification of your loop should work for(i in 1:N) { assign(paste(Random.Data., i, sep=''), rnorm(10, 0, 1)) } Cheers, Andy

Re: [R] GIS re-mapping / polygon overlap

2003-10-23 Thread Ross Ihaka
Hisaji Ono wrote: Hi. How about GEOS available at. http://geos.refractions.net/ GEOS supports hole for clipping. I could build it in MiGW. However I've only executed demo codes. Very nice reference. I will be taking a close look through it. The Greiner-Hormann algorithm is still pretty

[R] linear discriminant analysis

2003-10-23 Thread array chip
Hi all, When I use lda for discriminant analysis, should I normalized my data (variables) to mean 0, variance 1 before running lda if my variables might not be exactly on the same scale? I have this question because in principle component analysis, this is indeed an issue where we can choose

[R] IPFS Watch This St ock Trade

2003-10-23 Thread Hot Pennies
Hot Pennies.net Special Situation Alert Watch This Stock Trade!!! STOCK DETAILS Symbol: IPFS Exchange: OTCBB Shares Outstanding: 13,379,000 Current Price: $1.10 Industry: Wireless Manufacturing and Technology COMPANY DETAILS Integrated Performance

[R] Scan() question

2003-10-23 Thread Janice Warner
I am just beginning to use R 1.7.1 on Windows XP and can not do a simple scan of a .txt file with integer data separated by tabs. First I figured out that I need to open a connection to the file (which isn't mentioned in scan()). But now when I scan, I get Read 0 items. What can I do?

[R] Searching a Table

2003-10-23 Thread Nicole Baggette
Is there any way to search a column of a table? I read in a table and I want to create a subset based on the criteria that a certain column (filled with words) has a word which starts with a certain letter. The only method I have seen that searches this way is apropos(), but that doesn't seem

[R] R -- for help

2003-10-23 Thread jbwu
hi, guys: I have a question on how to read data from file in which the length of each line is unfixed. For instance, the data in sample.dat are: 1 2 3 4 5 6 7 8 I want to read this kind of data into a matrix as: 1 2 3 4 5 6 7 0 8 0 0 0 Thank you Jianbing Wu

Re: [R] Searching a Table

2003-10-23 Thread Jason Turner
Nicole Baggette wrote: Is there any way to search a column of a table? I read in a table and I want to create a subset based on the criteria that a certain column (filled with words) has a word which starts with a certain letter. The only method I have seen that searches this way is

Re: [R] Scan() question

2003-10-23 Thread Jason Turner
Janice Warner wrote: I am just beginning to use R 1.7.1 on Windows XP and can not do a simple scan of a .txt file with integer data separated by tabs. First I figured out that I need to open a connection to the file (which isn't mentioned in scan()). It isn't mentioned because it isn't