Re: [R] plot.baysian error = only 0's may mix with negative subscripts

2004-10-29 Thread Uwe Ligges
Konstantinos Liolios wrote: I did make a typo and I meant plot.bayesian() from sma. Hopefully most people can easily guess what I meant. Forgive me, my brain is too small! I do not know which of the millions of function is in which of the 400 CRAN packages. It's is your homework to say which

Re: [R] array problem and for looping

2004-10-29 Thread Petr Pikal
Hi Beter not to give a same name your values (variables) as is function name. R is quite clever and sample - rnorm(10) sample - sample(sample,3) works as expected, but it is not a rule. Cheers Petr On 28 Oct 2004 at 17:54, Kunal Shetty wrote: Dear R- users and Helpers Is there some way

Re: [R][OT] Latex Integration question

2004-10-29 Thread Göran Broström
On Thu, Oct 28, 2004 at 05:13:10PM -0500, Laura Holt wrote: Hi R users! I have a Latex Question: When I use the following: \begin{eqnarray*} F(x) = -\lambda e^{-\lambda t} \mid_{0}^{x} \end{eqnarray*} the output looks ok. But I would like the limits of integration line to be

[R] Interfacing R_project from program languages

2004-10-29 Thread Colombo Giovanni (FIA)
Hello, I need an information: is it possible to call the statistical R_project funtions from a C or Fortran application? In the r_exts.pdf manual I found that is possible to interface R-project with compiled C/Fortran routines, but I didn't found anything about the opposite situation. Thanks

Re: [R] Interfacing R_project from program languages

2004-10-29 Thread Prof Brian Ripley
On Fri, 29 Oct 2004, Colombo Giovanni (FIA) wrote: I need an information: is it possible to call the statistical R_project funtions from a C or Fortran application? In the r_exts.pdf manual I found that is possible to interface R-project with compiled C/Fortran routines, but I didn't found

[R] Re: reproducible rodbc crash

2004-10-29 Thread Jack Tanner
Naturally, the problem has already been fixed. update.packages() has an upgrade to RODBC 1.1-2, which does not crash. Huge thanks to Brian Ripley. Jack Tanner wrote: On Windows XP, R 2.0.0. library(RODBC) mydsn - odbcDriverConnect(DSN=mydsn) odbcClose(mydsn) sqlQuery(mydsn, select * from

[R] winDialog (equivalent) on Unix anyone ?

2004-10-29 Thread Steffen Neumann
Hi, We have a few R scripts, that have a minimalistic GUI Interface using winDialog() calls. Is there any cross-plattform equivalent we can port to, such that the code will run on both Windows and Unix ? I thought about either tcltk or gtk Stuff. I was astonished that I couldn't find anything

RE: [R] winDialog (equivalent) on Unix anyone ?

2004-10-29 Thread Henrik Bengtsson
Check out James Wettenhall's 'R TclTk Examples'; http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/ Henrik Bengtsson -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Steffen Neumann Sent: Friday, October 29, 2004 10:12 AM To: [EMAIL PROTECTED]

[R] missing values in logistic regression

2004-10-29 Thread Avril Coghlan
Dear R help list, I am trying to do a logistic regression where I have a categorical response variable Y and two numerical predictors X1 and X2. There are quite a lot of missing values for predictor X2. eg., Y X1 X2 red 0.6 0.2* red 0.5 0.2* red 0.5 NA red 0.5 NA

[R] Subsetting with more than one criteria

2004-10-29 Thread CG Pettersson
Hello all! I have been using R for two years, but still have trivial problems. For the moment, I run R1.9.1 on a W2k computer. I work with a dataset from a variety evaluation project. Every year 25 varieties have been tested. The tested cultivars have been the same each year in all trials, but

Re: [R] missing values in logistic regression

2004-10-29 Thread Peter Dalgaard
Avril Coghlan [EMAIL PROTECTED] writes: Dear R help list, I am trying to do a logistic regression where I have a categorical response variable Y and two numerical predictors X1 and X2. There are quite a lot of missing values for predictor X2. eg., Y X1 X2 red 0.6 0.2*

Re: [R] Subsetting with more than one criteria

2004-10-29 Thread Prof Brian Ripley
On Fri, 29 Oct 2004, CG Pettersson wrote: I have been using R for two years, but still have trivial problems. For the moment, I run R1.9.1 on a W2k computer. Time to read a book! This sort of thing is discussed at length in Chapter 2 of MASS4 (see the FAQ), for example. It even discusses why

Re: [R] missing values in logistic regression

2004-10-29 Thread Prof Brian Ripley
On 29 Oct 2004, Avril Coghlan wrote: Dear R help list, I am trying to do a logistic regression where I have a categorical response variable Y and two numerical predictors X1 and X2. There are quite a lot of missing values for predictor X2. eg., Y X1 X2 red 0.6 0.2*

[R] extraction of p-value

2004-10-29 Thread Luis Rideau Cruz
R-help, I wish to get the p-values of anova model as numeric and not as part of the standard output (print,summary,,,) I have tried several commands but to not vail Thank you __ [EMAIL PROTECTED] mailing list

[R] Installation problems with R.classes bundle

2004-10-29 Thread Michael Griffiths
Firstly hi to everyone on the list, I am new to this list ans also R so please forgive the simplicity of my questions over the next few months. I have version R 1.9.1 and want to perform z-scoring for a benchmarking procedure. I have tried to install Henrik Bengtsson's R.classes bundle

[R] How to obtain confidence intervals from a coxph object?

2004-10-29 Thread Vegard Andersen
Hi! I am performing an analysis which looks like this: cox.out - coxph(Survobject ~ x1 + x2 + x3) summary(cox.out) This works fine, but I have not been able to extract any of the output from summary(cox.out), for instance like: summary(cox.out)$coef or summary(cox.out)$lower So my question is

Re: [R] as.list.matrix

2004-10-29 Thread Kjetil Brinchmann Halvorsen
Prof Brian Ripley wrote: On Fri, 29 Oct 2004, Peter Alspach wrote: Kjetil Isn't a data.frame as special type of list, and thus one could use as.data.frame? Yes, or split(mat, col(mat)) works. So simple it hardly needs an as.list method. Thanks! I thought about split(), bit did'nt

Re: [R] Installation problems with R.classes bundle

2004-10-29 Thread Jan Goebel
On Fri, 29 Oct 2004, Michael Griffiths wrote: I have version R 1.9.1 and want to perform z-scoring for a benchmarking ---snip library(R.oo) Error in firstlib(which.lib.loc, package) : couldn't find function lazyLoad In addition: Warning message: package R.oo was built under

Re: [R] Installation problems with R.classes bundle

2004-10-29 Thread Uwe Ligges
Michael Griffiths wrote: Firstly hi to everyone on the list, I am new to this list ans also R so please forgive the simplicity of my questions over the next few months. I have version R 1.9.1 and want to perform z-scoring for a benchmarking procedure. I have tried to install Henrik Bengtsson's

Re: [R] meaning of frailty estimates

2004-10-29 Thread Emanuela Rossi
Sorry, but I have some difficulties to properly interpret these kind of hazard ratios without reference group. I'm studying the risk of introduction on a particular virus in poultry farms. My frailty variables is the specie (1, 2, 3, 4 are four types of poultry species which, I suppose, have

Re: [R] extraction of p-value

2004-10-29 Thread Prof Brian Ripley
On Fri, 29 Oct 2004, Luis Rideau Cruz wrote: I wish to get the p-values of anova model as numeric and not as part of the standard output (print,summary,,,) It is part of the standard output of summary. I have tried several commands but to not vail Please do as the posting guide asks and do

Re: [R] Installation problems with R.classes bundle

2004-10-29 Thread Prof Brian Ripley
On Fri, 29 Oct 2004, Michael Griffiths wrote: Firstly hi to everyone on the list, I am new to this list ans also R so please forgive the simplicity of my questions over the next few months. I have version R 1.9.1 and want to perform z-scoring for a benchmarking procedure. I have tried to

Re: [R] How to obtain confidence intervals from a coxph object?

2004-10-29 Thread Prof Brian Ripley
On Fri, 29 Oct 2004, Vegard Andersen wrote: I am performing an analysis which looks like this: cox.out - coxph(Survobject ~ x1 + x2 + x3) summary(cox.out) This works fine, but I have not been able to extract any of the output from summary(cox.out), for instance like:

Re: [R] extraction of p-value

2004-10-29 Thread Petr Pikal
Hi str(summary(fit)) will give you a whole picture how summary output list is structured. summary(fit)[n] where n is an appropriate number will extract parts of a list. Just try and choose the one you need. Maybe n=4? Cheers Petr On 29 Oct 2004 at 11:50, Luis Rideau Cruz wrote: R-help,

[R] glmmPQL and REML

2004-10-29 Thread sgrabasa
Hi, I am trying to use glmmPQL package for Generalized linear mixed models. This package works by repeated calls to lme. lme uses by default REML method for estimation. Then, does glmmmPQL use REML too? In contrast, how can I change it? I have tried it, writing : method=REML, but the

RE: [R] Installation problems with R.classes bundle

2004-10-29 Thread Henrik Bengtsson
Hi, I've added information on http://www.maths.lth.se/help/R/R.classes/ how to install Windows binaries of the latest R.classes on either R v2.0.0 or R v1.9.1. Quote: R v2.0.0: install.packages(R.classes, contriburl=http://www.maths.lth.se/help/R;) R v1.9.1: install.packages(R.classes,

RE: [R] missing values in logistic regression

2004-10-29 Thread Ted Harding
On 29-Oct-04 Avril Coghlan wrote: Dear R help list, I am trying to do a logistic regression where I have a categorical response variable Y and two numerical predictors X1 and X2. There are quite a lot of missing values for predictor X2. eg., Y X1 X2 red 0.6 0.2* red

[R] Datainput in R

2004-10-29 Thread christine mies
Sehr geehrte Damen und Herren, im Rahmen einer Seminararbeit an der Universität Augsburg soll ich Daten aus einer Excel-Datei in R einlesen und damit eine Faktorenanalyse, MDS und Procrustesanalyse durchführen. Nach einigen fehlgeschlagenen Versuchen richten wir nun folgende Fragen an Sie: Wie

Re: [R] glmmPQL and REML

2004-10-29 Thread Prof Brian Ripley
On Fri, 29 Oct 2004 [EMAIL PROTECTED] wrote: I am trying to use glmmPQL package for Generalized linear mixed models. This package works by repeated calls to lme. lme uses by default REML method for estimation. Then, does glmmmPQL use REML too? In contrast, how can I change it? I have

RE: [R] Installation problems with R.classes bundle

2004-10-29 Thread Henrik Bengtsson
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley Sent: Friday, October 29, 2004 2:28 PM To: Michael Griffiths Cc: [EMAIL PROTECTED] Subject: Re: [R] Installation problems with R.classes bundle On Fri, 29 Oct 2004, Michael

Re: [R] Datainput in R

2004-10-29 Thread Uwe Ligges
christine mies wrote: Sehr geehrte Damen und Herren, im Rahmen einer Seminararbeit an der Universitt Augsburg soll ich Daten aus einer Excel-Datei in R einlesen und damit eine Faktorenanalyse, MDS und Procrustesanalyse durchfhren. Nach einigen fehlgeschlagenen Versuchen richten wir nun folgende

[R] R list

2004-10-29 Thread Leonardo L Miceli
Hi Is there any function to get the name of the components of a given list object? ok. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

RE: [R] transitivity

2004-10-29 Thread Huntsinger, Reid
For transitivity, if P is your incidence matrix of 0s and 1s, then P%*%P has i,j entry equal to the number of k with Pik = 1 and Pkj = 1. So you want to check that when P%*%P has a nonzero (i,j) entry, then Pij = 1 too. Acyclicity: more generally the nth matrix power of P counts the number of

RE: [R] R list

2004-10-29 Thread Robert Sams
names() is what you want, if i understand your question correctly. cheers, robert -Original Message- From: Leonardo L Miceli [mailto:[EMAIL PROTECTED] Sent: Friday, October 29, 2004 3:22 PM To: [EMAIL PROTECTED] Subject: [R] R list Hi Is there any function to get the name of the

RE: [R] R list

2004-10-29 Thread Liaw, Andy
Yes, names(). Andy From: Leonardo L Miceli Hi Is there any function to get the name of the components of a given list object? ok. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

Re: [R] R list

2004-10-29 Thread Giovanni Petris
What about names() ? Giovanni Date: Fri, 29 Oct 2004 11:21:35 -0300 From: Leonardo L Miceli [EMAIL PROTECTED] Sender: [EMAIL PROTECTED] Cc: Precedence: list Hi Is there any function to get the name of the components of a given list object? ok. [[alternative HTML

[R] question about R on Linux Cluster

2004-10-29 Thread apollo wong
Hi, I am in trying to get a new Linux Cluster. I am thinking about putting R on it. Does anyone have any experience with running R on a distributed Linux Cluster? Or does R support a cluster environment at all. Any pitfall I should really watch out for? Thanks in advance Apollo

[R] Error in PDF output in R 2.0.0

2004-10-29 Thread james . holtman
The following script works fine in R 1.9.1. It was creating a PDF file with the graphs in it. In R 2.0.0, I got the error message below. I tried the same script just outputting to Windows and postscript and the output was OK. The error message only showed up when trying to create a PDF

RE: [R] lag variable addition to data frame question

2004-10-29 Thread Austin, Matt
?diff --Matt -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of s viswanath Sent: Friday, October 29, 2004 7:22 AM To: [EMAIL PROTECTED] Subject: [R] lag variable addition to data frame question Hi, I was wondering if there is a more efficient way of

Re: [R] question about R on Linux Cluster

2004-10-29 Thread Roger D. Peng
As far as I know, there are R interfaces to MPI and PVM. There are also the `snow' and `snowFT' packages on CRAN. -roger apollo wong wrote: Hi, I am in trying to get a new Linux Cluster. I am thinking about putting R on it. Does anyone have any experience with running R on a distributed Linux

RE: [R] lag variable addition to data frame question

2004-10-29 Thread Pfaff, Bernhard
Hello Sri, how about embed() in stats? Citing from the doc's detail section: Each row of the resulting matrix consists of sequences 'x[t]', 'x[t-1]', ..., 'x[t-dimension+1]', where 't' is the original index of 'x'. If 'x' is a matrix, i.e., 'x' contains more than one variable,

Re: [R] question about R on Linux Cluster

2004-10-29 Thread Seth Falcon
I've used an openMosix cluster to load balance simultaneously running R scripts. For the most part this has worked as expected, but sometimes the R processes don't migrate. Recently, I converted one of my analyses to utilize the SNOW package (using PVM on top of openMosix). I've been quite

Re: [R] Error in PDF output in R 2.0.0

2004-10-29 Thread Deepayan Sarkar
On Fri, 29 Oct 2004 [EMAIL PROTECTED] wrote: The following script works fine in R 1.9.1. It was creating a PDF file with the graphs in it. In R 2.0.0, I got the error message below. I tried the same script just outputting to Windows and postscript and the output was OK. The error message

Re: [R] ifelse() question

2004-10-29 Thread F Z
Thanks for you reply Peter. I tried using as.character and then converting to factors but it did not work since it generates missing values for all the dat[,4]==POR. See: dat[1:5,4] [1] BOV POR BOV POR BOV Levels: BOV CAP CER OVI POR test-dat test[,4]-as.character(test[,4])

[R] Poisson regression

2004-10-29 Thread Federico Gherardini
Dear all, First of all sorry if this is a dumb question... I'm trying to fit a linear model between the logarithm of two numerical variables (log(y) ~ log(x)). A log-log plot shows that the variance of log(y) is decreasing with the mean of log(x), in other words the points are quite dispersed

[R] Giving column names to a matrix

2004-10-29 Thread jmoreira
Heloo, I have the following problem: orig.data - NULL Inside a loop I have instructions like: orig.data - rbind(orig.data, ...) After that I do: colnames(orig.data)-c('Data','InicioViagem', ...) Everything works fine. For example, the first line of the matrix is: orig.data[1,]

Re: [R] ifelse() question

2004-10-29 Thread Adaikalavan Ramasamy
Francisco, a more reproducible example would have helped but see if the following helps in your understanding # Create dataset df - data.frame( type=1:5, species=LETTERS[1:5] ) df type species 11 A 22 B 33 C 44 D 55 E df[ ,2] == E

Re: [R] Giving column names to a matrix

2004-10-29 Thread Sundar Dorai-Raj
[EMAIL PROTECTED] wrote: Heloo, I have the following problem: orig.data - NULL Inside a loop I have instructions like: orig.data - rbind(orig.data, ...) After that I do: colnames(orig.data)-c('Data','InicioViagem', ...) Everything works fine. For example, the first line of the matrix is:

[R] Warning message if the plot statement breaks into 2 lines

2004-10-29 Thread partha_bagchi
Here is a curious observation. In version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major2 minor0.0 year 2004 month10 day 04 language R Try the following: Plot(c(1:100), type = l + ) Warning message:

[R] pattern search

2004-10-29 Thread Sean Liang
hi, I like to find a pattern within a giver sequence. There might be multiple occurences of the pattern. I like to know the number of occurences and the positions if possible. grep can tell me if a pattern exists but can't give me the information I need. Does anyone know any function that I can

[R] fitting linear mixed model for incomplete block design

2004-10-29 Thread Caimiao Wei
Dear R developers and users: I have the following data, x is the response vaiable, nsample(individual) nested within trt, and subsample nested within nsample, I want to fit trt as fixed effect, and block, nsample(trt) as random effects using lme, is the following coding correct? dat$vgrp

Re: [R] Warning message if the plot statement breaks into 2 lines

2004-10-29 Thread Mike Meyer
It has all to do with a string going over two lines and hence the line terminator is part of the string. If you close the string on the first line then the command is just fine. At least the newer version of the parser allows for embedded newline characters. --Mike On Fri, 29 Oct 2004

Re: [R] Warning message if the plot statement breaks into 2 lines

2004-10-29 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Here is a curious observation. In version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major2 minor0.0 year 2004 month10 day 04 language R Try the following: Plot(c(1:100), type = l +

Re: [R] pattern search

2004-10-29 Thread Uwe Ligges
Sean Liang wrote: hi, I like to find a pattern within a giver sequence. There might be multiple occurences of the pattern. I like to know the number of occurences and the positions if possible. grep can tell me if a pattern exists but can't give me the information I need. Does anyone know any

Re: [R] pattern search

2004-10-29 Thread Sean Liang
Uwe: Thanks for getting back to me. My situation is a bit more complicated. I have a vector of sequences each of which might contain any number of a given pattern (e.g. pat=c(ATCGTTTGCTAC, GGCTAATGCATTGC); grep (TGC, pat) [1] 1 2 grep only tells me the position of first occurrence in each

RE: [R] question about R on Linux Cluster

2004-10-29 Thread Huntsinger, Reid
We use R on a Linux/openMosix cluster and it works very well. Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of apollo wong Sent: Friday, October 29, 2004 10:48 AM To: [EMAIL PROTECTED] Subject: [R] question about R on Linux Cluster Hi,

[R] [rmetasim] Need help deciphering this error msg... targeted to those who use rmetasim...

2004-10-29 Thread Brian Seok
Hello, I am trying to do some simulation using the rmetasim package and I've run to this problem. --beginning of error msg-- Error in [-(`*tmp*`, slice[l, ], slice[l, ], value = c(0.20002980232, : number of items to replace is not a multiple of replacement length --end of error

[R][rmetasim] Need help deciphering this error msg... targeted to those who use rmetasim...

2004-10-29 Thread Brian Seok
Hello, I am trying to do some simulation using the rmetasim package and I've run to this problem. --beginning of error msg-- Error in [-(`*tmp*`, slice[l, ], slice[l, ], value = c(0.20002980232, : number of items to replace is not a multiple of replacement length --end of error

Re: [R] plot.baysian error = only 0's may mix with negative subscripts

2004-10-29 Thread Konstantinos Liolios
Thanx a lot my friend, you have been helpfull enough. Dinos On Fri, 2004-10-29 at 01:36, Uwe Ligges wrote: Konstantinos Liolios wrote: I did make a typo and I meant plot.bayesian() from sma. Hopefully most people can easily guess what I meant. Forgive me, my brain is too small! I

Re: [R] Warning message if the plot statement breaks into 2 lines

2004-10-29 Thread Prof Brian Ripley
This is a consequence of a change described in NEWS to allow multi-line *strings*. On Fri, 29 Oct 2004 [EMAIL PROTECTED] wrote: Here is a curious observation. In version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major2

Re: [R] fitting linear mixed model for incomplete block design

2004-10-29 Thread Douglas Bates
Caimiao Wei wrote: Dear R developers and users: I have the following data, x is the response variable, sample(individual) nested within trt, and subsample nested within nsample, I want to fit trt as fixed effect, and block, nsample(trt) as random effects using lme, is the following coding

Re: [R] missing values in logistic regression

2004-10-29 Thread Frank E Harrell Jr
(Ted Harding) wrote: On 29-Oct-04 Avril Coghlan wrote: Dear R help list, I am trying to do a logistic regression where I have a categorical response variable Y and two numerical predictors X1 and X2. There are quite a lot of missing values for predictor X2. eg., Y X1 X2 red 0.6 0.2*

[R] pairs(x)

2004-10-29 Thread John-J . Smith
Hi, I am using pairs(x) , but would like to change the point font and color for the bottom half of x. I am using windows and opened a graphics page using x11(), then tried: pairs(x) points(x[110:114,1],x[110:114,2],col=6,pch=8) points(x[115:119,1],x[115:119,2],col=4,pch=17) Nothing happened.

[R] Mac OS X binary package naming convention

2004-10-29 Thread Paul Roebuck
Did the convention change during the R 2.0.0 Cocoa update? When I issue the build command: $ R CMD build --binary mypkg I expected to get 'mypkg_0.9-0.tgz' but instead got 'mypkg_0.9-0_R_powerpc-apple-darwin6.8.tar.gz'. Or should I rename it manually? platformpowerpc-apple-darwin6.8

[R] test

2004-10-29 Thread Rajdeep Das
[[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] why should you set the mode in a vector?

2004-10-29 Thread Joel Bremson
Hi all, If I write v = vector(mode=numeric,length=10) I'm still allowed to assign non-numerics to v. Furthermore, R figures out what kind of vector I've got anyway when I use the mode() function. So what is it that assigning a mode does? Joel __

Re: [R] ifelse() question

2004-10-29 Thread F Z
Many thanks to Peter, Sundar and Adaikalavan for their useful help. As Adaikalavan correctly pointed out, the key was to pass as.character() as an argument within the ifelse() function. Here is the soution using a subset of my data: dat[1:10,] id long lat species type size 1

Re: [R] why should you set the mode in a vector?

2004-10-29 Thread Tony Plate
It's useful when you need to be certain of the mode of a vector. One such situation is when you are about to call a C-language function using the .C() interface. As you point out, some assignments (even just to vector elements) can change the mode of the entire vector. This is why it's

Re: [R] pattern search

2004-10-29 Thread Gabor Grothendieck
Sean Liang SLiang at wyeth.com writes: I have a vector of sequences each of which might contain any number of a given pattern (e.g. pat=c(ATCGTTTGCTAC, GGCTAATGCATTGC); grep (TGC, pat) [1] 1 2 grep only tells me the position of first occurrence in each element whereas the second

[R] Can R import data from firebird relational database

2004-10-29 Thread Unung Istopo Hartanto
Dear Users, Can R import data or process data from firebird relational database (firebird.sourceforge.net). Cause i still using it on linux and i'll process my data using R. Thanks all, regards, Unung __ [EMAIL PROTECTED] mailing list