[R] The 'test.terms' argument in 'regTermTest' in package 'survey'

2004-08-19 Thread Almirall, Daniel
This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X, say, how do the two test specifications

[R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread Jacob Wegelin
Alternate title: How can I persuade my students that R is for them? Alternate title: Can R replace SAS, SPSS or Stata for clinicians? I am teaching introductory statistics to twelve physicians and two veterinarians who have enrolled in a Mentored Clinical Research Training Program. My course

Re: [R] glmmPQL in R and S-PLUS 6 - differing results

2004-08-19 Thread Prof Brian Ripley
glmmPQL calls lme, and it is that which differs between S-PLUS and R. These are optimization problems with multiple local maxima, and like any complex statistical fitting problem you should not expect all programs to give the same answer. The short answer is to believe all of them, and to

RE: [R] Do you know if you can map a large minimum spanning tree in R?

2004-08-19 Thread Roger Bivand
Thanks Mike. My data has longitude and latitude coords and I used distAB {clim.pact} then mst {ape} to calculate my minimum spanning tree. The nodes are telecoms sites from all over Australia. My goal is to determine the minimum cost of linking them via cabling, and I'm starting by

RE: [R] Do you know if you can map a large minimum spanning tree in R?

2004-08-19 Thread Roger Bivand
Thanks Mike. My data has longitude and latitude coords and I used distAB {clim.pact} then mst {ape} to calculate my minimum spanning tree. The nodes are telecoms sites from all over Australia. My goal is to determine the minimum cost of linking them via cabling, and I'm starting by

Re: [R] OS X specific question: help.start() won't launch

2004-08-19 Thread Ingmar Visser
Hi Dave, This has come up a couple of times recently. I don't remember by heart but searching the archives should help you find the answer. It's best to try R-SIG-MAC list archive. Bye, ingmar On 8/19/04 4:34 AM, David L. Van Brunt, Ph.D. [EMAIL PROTECTED] wrote: It's been a while since I used

Re: [R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread Ernesto Jardim
On Thu, 2004-08-19 at 07:45, Jacob Wegelin wrote: Alternate title: How can I persuade my students that R is for them? Alternate title: Can R replace SAS, SPSS or Stata for clinicians? I am teaching introductory statistics to twelve physicians and two veterinarians who have enrolled in a

[R] how to extract data from excel sheet?

2004-08-19 Thread Sivan Aldor
i want to upload different sheets from one excel file. does anyone know of a way to do so? (i use read.csv to load the excel but i would like to load form the sheets themselves.) thanks sivan __ [EMAIL PROTECTED] mailing list

Re: [R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread A.J. Rossini
One thought -- the first course I took, I taught physicians XLispStat. They appreciated it since it allowed them to do logistic regression without spending large $$$. Now, that isn't quite true any more. So might depend on how times have changed. Jacob Wegelin [EMAIL PROTECTED] writes:

Re: [R] how to extract data from excel sheet?

2004-08-19 Thread Uwe Ligges
Sivan Aldor wrote: i want to upload different sheets from one excel file. does anyone know of a way to do so? Yes, the R Data Import/Export manual knows. Uwe Ligges (i use read.csv to load the excel but i would like to load form the sheets themselves.) thanks sivan

[R] sample selection problem, inverse mills ratio (Heckman, Lewbel, ...)

2004-08-19 Thread Wildi Marc, wia
-Ursprüngliche Nachricht- Von: Wildi Marc, wia Gesendet: Mittwoch, 18. August 2004 10:11 An: [EMAIL PROTECTED] Betreff: Hi Does anybody know from an R-package devoted to sample selection problems

[R] Erlang distribution

2004-08-19 Thread Nils Aschenbruck
Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: 7bit Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted sender hosts) X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on

Re: [R] how to extract data from excel sheet?

2004-08-19 Thread Kevin Wang
Hi, Take a look at read.xls() from gregmisc package. HTH, Kevin On Thu, 19 Aug 2004, Sivan Aldor wrote: i want to upload different sheets from one excel file. does anyone know of a way to do so? (i use read.csv to load the excel but i would like to load form the sheets themselves.) thanks

RE: [R] OS X specific question: help.start() won't launch

2004-08-19 Thread Brunsdon C (Comp)
I've got the same problem with help.start(). It seems to be a problem with the 'open' command when executed through R's 'system' function - ie system('open something.html') gives an error message, whilst opening a new xterm (or console) and entering open something.html works fine (here

RE: [R] Erlang distribution

2004-08-19 Thread Trenkler, Dietrich
-Original Message- From: Nils Aschenbruck Sent: Thursday, August 19, 2004 11:51 AM To: [EMAIL PROTECTED] Subject: [R] Erlang distribution Hello, is there a packet that supports the Erlang distribution? I want to use this distribution for tests against empirical

Re: [R] how to extract data from excel sheet?

2004-08-19 Thread Chuck Cleland
Here is one way to read the mysheet1 worksheet from the mydata.xls file: library(RODBC) z - odbcConnectExcel(c:\\mydata.xls) myframe - sqlFetch(z, mysheet1) close(z) Your data will be in myframe. Thanks to Michael Lapsley and Brian Ripley for the RODBC package. See the R

Re: [R] distance to cluster center

2004-08-19 Thread Christian Hennig
On Wed, 18 Aug 2004, Moises Hassan wrote: Is there a way to get the distance between each point and the center of the cluster it was assigned to in cluster methods such as agnes, pam, and clara. I don't think that it's directly available, but for pam and clara you can easily get it from your

[R] Use Font Frutiger in PDF-File

2004-08-19 Thread Helga Neidlinger
Hi, I'd like to use the font Frutiger in my pdf files generated by R. but I only manage this if I create a postscript file with postscript() and the option family = c( LT_14361.afm, LT_14363.afm, LT_14373.afm, LT_14375.afm ) and convert this postscript file to pdf. If I specify the option

Re: [R] how to extract data from excel sheet?

2004-08-19 Thread Petr Pikal
On 19 Aug 2004 at 12:05, Sivan Aldor wrote: i want to upload different sheets from one excel file. does anyone know of a way to do so? (i use read.csv to load the excel but i would like to load form the sheets themselves.) thanks sivan Hi Open Excel, go to the sheet, select rectangular

[R] Clipping of display in Lattice graphics

2004-08-19 Thread Matthew Walker
I'm baffled as to how the Lattice package achieves clipping. Would someone mind explaining this to me? Firstly, my attempt using just the grid package: x-seq(0,3,by=0.3)/2.8 y-seq(0,1,by=0.1) grid.newpage() grid.rect(gp=gpar(fill=pink)) vp-viewport(width=0.8, height=0.8) pushViewport(vp)

Re: [R] Use Font Frutiger in PDF-File

2004-08-19 Thread Prof Brian Ripley
Have you read the help file? It says (R 1.9.1) family: the font family to be used, one of 'AvantGarde', 'Bookman', 'Courier', 'Helvetica', 'Helvetica-Narrow', 'NewCenturySchoolbook', 'Palatino' or 'Times'. Now, it would seem if you can get away with specifying a

[R] sorting character vectors

2004-08-19 Thread andreas . krause
The following is not what I expected in sorting characters (single letters and the same letters with preceding spaces). Can someone enlighten me as to why the following might be a correct result for sorting? ; x - c(LETTERS[1:3], paste( , LETTERS[1:3], sep=)) ; x [1] A B C A B C ; sort(x)

[R] pdf font encoding

2004-08-19 Thread Petr Pikal
Dear all. This is probably more PDF then R issue but anyway, maybe somebody can help me. I try to persuade pdf device to use eastern european (Czech) encoding but I was not successful yet. Working on W 2000, R 1.9.0, lattice plots. pdf(redprov.pdf, width = 5, height = 5, bg=white,

Re: [R] sorting character vectors

2004-08-19 Thread Prof Brian Ripley
It is documented to depend on your locale. I get sort(x) [1] A B C A B C in the C locale. The help page does say so: The sort order for character vectors will depend on the collating sequence of the locale in use: see 'Comparison'. The default collation sequences for

[R] Fortran 90 in R?

2004-08-19 Thread Ben Yip
Hi there, I am using Fortran 77 g77 compiler to write Fortran subprogram to R. I would like to change and write Fortran 90 code instead (no more specification of matrix!). Does gcc support Fortran 90? And if Yes, does it works in R? Ben __ [EMAIL

Re: [R] sorting character vectors

2004-08-19 Thread Roger D. Peng
Sorting depends on the locale. For example, I get x - c(LETTERS[1:3], paste( , LETTERS[1:3], sep=)) x [1] A B C A B C sort(x) [1] A B C A B C On Linux (Fedora Core 1) I set LANG=C, which is necessary for some other (non-R) things to work. -roger [EMAIL PROTECTED] wrote: The

RE: [R] sorting character vectors

2004-08-19 Thread Liaw, Andy
I do get (R-1.9.1 on WinXPPro): x - c(LETTERS[1:3], paste( , LETTERS[1:3], sep=)) sort(x) [1] A B C A B C The `right' sequence is locale dependent, and I wonder if that's the discrepancy here. Andy From: [EMAIL PROTECTED] The following is not what I expected in sorting characters

RE: [R] A question about external time-dependent covariates in co x model

2004-08-19 Thread Hanke, Alex
Dear Rui, From my understanding of time-dependent covariates (not an expert but have been working on a similar problem), it would appear that the coding of the status column is not correct. Unless you have observed an event at each interval you should only have status=1 for the last interval. In

Re: [R] labeled break statements in R?

2004-08-19 Thread Gabor Grothendieck
Gabor Grothendieck ggrothendieck at myway.com writes: Roger Levy rog at stanford.edu writes: : Are there labeled break statements in R? i.e., something along the : lines of : : TOPLOOP: for(i in 1:m) { : for(j in 1:n) { : ... : if(condition) { :break TOPLOOP : }

Re: [R] Fortran 90 in R?

2004-08-19 Thread Prof Brian Ripley
On Thu, 19 Aug 2004, Ben Yip wrote: I am using Fortran 77 g77 compiler to write Fortran subprogram to R. I would like to change and write Fortran 90 code instead (no more specification of matrix!). Does gcc support Fortran 90? Isn't that a question for a gcc-help list? Take a look at

Re: [R] sorting character vectors

2004-08-19 Thread andreas . krause
Thank you very much to Brian Ripley, Roger Peng, and Andy Liaw. Everyone pointed out the same solution. Setting LC_COLLATE=C did it. This default setting is indeed odd to me. Thanks again! Andreas Krause __ [EMAIL PROTECTED] mailing list

RE: [R] OS X specific question: help.start() won't launch

2004-08-19 Thread Prof Brian Ripley
You need to find out what put /usr/local/lib/libxml2.2.dylib there and correct it. Look at PR#7063 in the R-bugs archive. Stefano Iacus said there Please remove your copy of libxml2 in /usr/local/lib, you no longer need it on Panther and the one included in 10.3.4 is incompatible with the

[R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread White, Charles E WRAIR-Wash DC
I'm not in your target audience but I consult with and teach to your target audience, so I'll chime in anyway. 1. Look at the user friendliness of the alternatives you offered to R. I haven't used SPSS recently or Stata at all but I know SAS is NOT that user friendly. Note that my institute

[R] Suggestion for posting guide

2004-08-19 Thread Gabor Grothendieck
I have a suggestion for the posting guide. One problem with some posts is that they do not provide an example that can be reproduced. I think that many people just do not know how to easily specify some data and some technical assistance should be provided in the posting guide. If the problem

Re: [R] Fortran 90 in R?

2004-08-19 Thread Christian T. Steigies
On Thu, Aug 19, 2004 at 01:50:47PM +0100, Prof Brian Ripley wrote: On Thu, 19 Aug 2004, Ben Yip wrote: I am using Fortran 77 g77 compiler to write Fortran subprogram to R. I would like to change and write Fortran 90 code instead (no more specification of matrix!). Does gcc support

Re: [R] sorting character vectors

2004-08-19 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: Thank you very much to Brian Ripley, Roger Peng, and Andy Liaw. Everyone pointed out the same solution. Setting LC_COLLATE=C did it. This default setting is indeed odd to me. It's odd to everyone, except perhaps to the POSIX locale mafia... Once upon a time all

[R] How to randomize a set of integers in R

2004-08-19 Thread Susanne Osterkamp
Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: 7bit Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted sender hosts) X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on

Re: [R] Suggestion for posting guide

2004-08-19 Thread Rolf Turner
Technical point: Gabor Grothendieck wrote ``data(iris) will attach data set iris'' This is not true --- data(iris) put ``iris'' into the Global environment; if one wants it attached, one has to do the follow up by attach(iris). cheers,

[R] Ques about the 'test.terms' argument in 'regTermTest'

2004-08-19 Thread Almirall, Daniel
Hi, Z is a three level factor variables and ZB and ZC are its corresponding binary dummy variables. In a 'glm' of Y on Z and X, say, how do the two test specifications test.terms = c(ZB:X,ZC:X) # and test.terms = ~ ZB:X + ZC:X in 'regTermTest' differ? I thought that both

Re: [R] How to randomize a set of integers in R

2004-08-19 Thread Olaf Mersmann
Hi Susanne, * Susanne Osterkamp [EMAIL PROTECTED] [040819 15:52]: Hello, I want to validate my data for working with CART. I need to now how to create a randomized column of some 6800 intgers within R. Do you want to sample 6800 random integers? In that case try ?runif (for uniform

Re: [R] A question about external time-dependent covariates in co x model

2004-08-19 Thread Göran Broström
On Thu, Aug 19, 2004 at 09:36:22AM -0300, Hanke, Alex wrote: Dear Rui, From my understanding of time-dependent covariates (not an expert but have been working on a similar problem), it would appear that the coding of the status column is not correct. Unless you have observed an event at each

[R] RMySQL, Solaris and getopt????

2004-08-19 Thread Roggie Boone
Hi, I've been struggling with RMySQL installation for several days. I'm running Solaris 2.9 and MySQL 4.0. RMySQL (v0.4.3) seems to install OK, but when I load it in R I get the message: --- library(RMySQL) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load

Re: [R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread Jonathan Baron
A few comments. 1. What students use. I've gotten a few undergraduate psychology students to use R to analyze their data for independent study projects. They do fine. I prepare the data for them and show them how to do t tests and correlations, and how to define new variables from old ones.

Re: [R] Clipping of display in Lattice graphics

2004-08-19 Thread Deepayan Sarkar
On Thursday 19 August 2004 06:38, Matthew Walker wrote: I'm baffled as to how the Lattice package achieves clipping. Would someone mind explaining this to me? Firstly, my attempt using just the grid package: x-seq(0,3,by=0.3)/2.8 y-seq(0,1,by=0.1) grid.newpage()

[R] convert strings to object names

2004-08-19 Thread Gudrun Jonasdottir
Dear R-Help list, I have a problem with convertions of strings. I want to use the function paste() to create an object name and then use that character string to call on that object. So, for example: dat99 - matrix(rbind(1,1,2),3,3) no - 99 dat - paste(dat,no,sep=) dat [1] dat99 What should I

RE: [R] How to randomize a set of integers in R

2004-08-19 Thread Liaw, Andy
Use sample() in either case: ## randomly shuffle the elements of the vector `x': rx - sample(x, replace=FALSE) ## randomly draw 6800 integers between 1 and 6800: x - sample(6800, 6800, replace=TRUE) Andy From: Olaf Mersmann Hi Susanne, * Susanne Osterkamp [EMAIL PROTECTED] [040819 15:52]:

RE: [R] convert strings to object names

2004-08-19 Thread Liaw, Andy
See ?get. Andy From: Gudrun Jonasdottir Dear R-Help list, I have a problem with convertions of strings. I want to use the function paste() to create an object name and then use that character string to call on that object. So, for example: dat99 - matrix(rbind(1,1,2),3,3) no -

Re: [R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread Marc Schwartz
On Thu, 2004-08-19 at 01:45, Jacob Wegelin wrote: Alternate title: How can I persuade my students that R is for them? Alternate title: Can R replace SAS, SPSS or Stata for clinicians? I am teaching introductory statistics to twelve physicians and two veterinarians who have enrolled in a

Re: [R] convert strings to object names

2004-08-19 Thread Dimitris Rizopoulos
Hi Gudrun, try to use eval(parse(text=dat)) Best, Dimitris Dimitris Rizopoulos Doctoral Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/396887 Fax: +32/16/337015 Web:

Re: [R] convert strings to object names

2004-08-19 Thread Uwe Ligges
Gudrun Jonasdottir wrote: Dear R-Help list, I have a problem with convertions of strings. I want to use the function paste() to create an object name and then use that character string to call on that object. So, for example: dat99 - matrix(rbind(1,1,2),3,3) no - 99 dat - paste(dat,no,sep=) dat

Re: [R] convert strings to object names

2004-08-19 Thread Adaikalavan Ramasamy
assign(dat, get(paste(dat, no, sep=))) or simply dat - get(paste(dat, no, sep=)) On Thu, 2004-08-19 at 15:15, Gudrun Jonasdottir wrote: Dear R-Help list, I have a problem with convertions of strings. I want to use the function paste() to create an object name and then use that

Re: [R] convert strings to object names

2004-08-19 Thread Thomas Lumley
And this is all in the FAQ. -thomas On Thu, 19 Aug 2004, Adaikalavan Ramasamy wrote: assign(dat, get(paste(dat, no, sep=))) or simply dat - get(paste(dat, no, sep=)) On Thu, 2004-08-19 at 15:15, Gudrun Jonasdottir wrote: Dear R-Help list, I have a problem with

Re: [R] A question about external time-dependent covariates in co x model

2004-08-19 Thread Thomas Lumley
On Thu, 19 Aug 2004, Göran Broström wrote: On Thu, Aug 19, 2004 at 09:36:22AM -0300, Hanke, Alex wrote: Dear Rui, From my understanding of time-dependent covariates (not an expert but have been working on a similar problem), it would appear that the coding of the status column is not

Re: [R] The 'test.terms' argument in 'regTermTest' in package 'survey'

2004-08-19 Thread Thomas Lumley
On Thu, 19 Aug 2004, Almirall, Daniel wrote: This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X,

RE: [R] A question about external time-dependent covariates in co x model

2004-08-19 Thread Rui Song
Dear Alex, Thanks a lot for your help! The end is not proportional to the covariate in my data set. I recoded the status as you suggested, but the problem still exists. I tried to fit the model with partial data, it has no problem when I use the first 25 rows (from the data in the form start,

[R] GEEs for time series data

2004-08-19 Thread Duncan Lee
I want to run a GEE for a time series of counts. The data are daily respiratory mortality counts and so there aren't any 'clusters' in the longitudinal sense. Neither the gee or geese packages work. The gee one wont run at all and the geese one produces NaNs or just runs indefinitely depending on

[R] rgdal under windows?

2004-08-19 Thread Benjamin Lloyd-Hughes
Hi, Has anyone had any joy getting the rgdal package to compile under windows? I've been trying with MinGW using Duncan Murdoch's Rtools but to no avail. After unpacking the tarball I put it here: D:\rw1091\src\library\rgdalls DESCRIPTION Rconfigure inst man INDEX

[R] precision problems in testing with Intel compilers

2004-08-19 Thread Samuelson, Frank*
I compiled the 1.9.1 src.rpm with the standard gnu tools and it works. I tried compiling the 1.9.1 src.rpm with the Intel 8 C and FORTRAN compilers and it bombs out during the testing phase: comparing 'd-p-q-r-tests.Rout' to './d-p-q-r-tests.Rout.save' ...267c267 df = 0.5[1] Mean

Re: [R] How to convert a vector into a list

2004-08-19 Thread S Peri
Dear group, Apologies for asking the most chomped FAQ. I have a file with list of gene names(genes.txt): EGF EGFR PTPN6 TIEG2 MAPK1 I have another object in R, I do not know the data type for that object that looks like this (lidnames): RABGGTA MAPK3 TIE CYP2C19 lidnames[1:10]

[R] NLME: Holding constant the across group correlational structure of the fixed effects in nlme

2004-08-19 Thread [EMAIL PROTECTED]
Hello all. I was wondering if there is a way to hold constant the fixed effects correlation structure across multiple groups? For example, I have two groups and I fit a three parameter logistic growth curve where the fixed effects are free to vary across the groups. I'll paste in the code as a

Re: [R] convert strings to object names

2004-08-19 Thread Gardar Johannesson
Use the get() function: dat99 - matrix(c(1,1,2),3,3) no - 99 dat.name - paste(dat,no,sep=) get(dat.name) [,1] [,2] [,3] [1,]111 [2,]111 [3,]222 At 04:15 PM 8/19/2004 +0200, Gudrun Jonasdottir wrote: Dear R-Help list, I have a problem with convertions

Re: [R] Is R good for not-professional-statistician, , un-mathematical clinical researchers?

2004-08-19 Thread Frank E Harrell Jr
We are finding more and more clinical researchers interested in learning R and are starting to teach R to clinicians. We have put some teaching material on our site: http://biostat.mc.vanderbilt.edu -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] How to convert a vector into a list

2004-08-19 Thread Uwe Ligges
S Peri wrote: Dear group, Apologies for asking the most chomped FAQ. I have a file with list of gene names(genes.txt): EGF EGFR PTPN6 TIEG2 MAPK1 I have another object in R, I do not know the data type for that object that looks like this (lidnames): RABGGTA MAPK3 TIE CYP2C19

Re: [R] paired t-test vs pairwise t-test

2004-08-19 Thread Rolf Turner
You wrote: What's the difference between t.test(x, y) and pairwise.t.test()? Is it just that the former takes two vectors, whereas the latter takes a vector and a factor? No. The pairwise.t.test() function (according to the help file) does a multiplicity of t-tests, on more

Re: [R] pdf font encoding

2004-08-19 Thread M.Kondrin
Petr Pikal wrote: Dear all. This is probably more PDF then R issue but anyway, maybe somebody can help me. I try to persuade pdf device to use eastern european (Czech) encoding but I was not successful yet. Working on W 2000, R 1.9.0, lattice plots. pdf(redprov.pdf, width = 5, height = 5,

Re: [R] IDE or an Editor for R

2004-08-19 Thread S Peri
Hi, Is there any IDE or any editor or any pice of code for .vimrc is available? Please let me know. Thanks SP ___ Win 1 of 4,000 free domain names from Yahoo! Enter now. __ [EMAIL PROTECTED] mailing

Re: [R] paired t-test vs pairwise t-test

2004-08-19 Thread Jack Tanner
From: Rolf Turner [EMAIL PROTECTED] No. The pairwise.t.test() function (according to the help file) does a multiplicity of t-tests, on more than two samples, adjusting the p-value to compensate for the multiplicity by various methods. Thank you, Rolf, that's helpful. So pairwise.t.test() is used

Re: [R] paired t-test vs pairwise t-test

2004-08-19 Thread Duncan Murdoch
On Thu, 19 Aug 2004 13:42:21 -0300 (ADT), Rolf Turner [EMAIL PROTECTED] wrote : You wrote: What's the difference between t.test(x, y) and pairwise.t.test()? Is it just that the former takes two vectors, whereas the latter takes a vector and a factor? No. The pairwise.t.test()

[R] Clustering and the test for proportional hazards

2004-08-19 Thread Hanke, Alex
Dear List, Is the test for proportional hazards valid when the model contains a cluster variable? The output looks strange with the cluster variable. My intervals are based on calendar time and the clustering variable is related to the season the event occurs in.

RE: [R] paired t-test vs pairwise t-test

2004-08-19 Thread Liaw, Andy
From: Duncan Murdoch On Thu, 19 Aug 2004 13:42:21 -0300 (ADT), Rolf Turner [EMAIL PROTECTED] wrote : You wrote: What's the difference between t.test(x, y) and pairwise.t.test()? Is it just that the former takes two vectors, whereas the latter takes a vector and a factor?

[R] save() load()

2004-08-19 Thread Harris A. Jaffee
Is there a way to list the contents of a save file, i.e. the variable names, without load'ing from it first? __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

[R] Question on TukeyHSD

2004-08-19 Thread Hem Nalini Morzaria-Luna
Hi, I am running a ANOVA on a factorial design, and using TukeyHSD for post hoc comparisons. I have 2 factors with three levels each: Factor B Factor A1 2 3 1 2 3 When I look at the Tukey output (on the interaction of the factors) the

RE: [R] More precision problems in testing with Intel compilers

2004-08-19 Thread Samuelson, Frank*
The Intel compiled version also fails the below test: ### Very big and very small umach - unlist(.Machine)[paste(double.x, c(min,max), sep='')] xmin - umach[1] xmax - umach[2] tx - unique(outer(-1:1,c(.1,1e-3,1e-7)))# 7 values (out of 9) tx - unique(sort(c(outer(umach,1+tx#

Re: [R] paired t-test vs pairwise t-test

2004-08-19 Thread Peter Dalgaard
Jack Tanner [EMAIL PROTECTED] writes: From: Rolf Turner [EMAIL PROTECTED] No. The pairwise.t.test() function (according to the help file) does a multiplicity of t-tests, on more than two samples, adjusting the p-value to compensate for the multiplicity by various methods. Thank you,

[R] Package rmutil

2004-08-19 Thread Savano
Users, Where I can to download rmutil package? thanks __ [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

Re: [R] How to convert a vector into a list

2004-08-19 Thread Adaikalavan Ramasamy
Peri, Please stop posting to both R-help and BioConductor simultaneously. Please read the posting guide and decide which list is more appropriate. http://www.R-project.org/posting-guide.html On Thu, 2004-08-19 at 17:20, S Peri wrote: Dear group, Apologies for asking the most chomped

Re: [R] save() load()

2004-08-19 Thread Prof Brian Ripley
On Thu, 19 Aug 2004, Harris A. Jaffee wrote: Is there a way to list the contents of a save file, i.e. the variable names, without load'ing from it first? I believe not. But that's not too bad: x - system.file(package=MASS, data, hills.rda) print(load(x , new.env())) [1] hills -- Brian D.

Re: [R] Question on TukeyHSD

2004-08-19 Thread Prof Brian Ripley
On Thu, 19 Aug 2004, Hem Nalini Morzaria-Luna wrote: Hi, I am running a ANOVA on a factorial design, and using TukeyHSD for post hoc comparisons. I have 2 factors with three levels each: Factor B Factor A1 2 3 1 2 3 When I look at the

Re: [R] IDE or an Editor for R

2004-08-19 Thread Adaikalavan Ramasamy
I am not sure what does an IDE consist besides syntax highlighting and ability to parse lines from script into console. Nor am I familiar with VIM but IMHO emacs with ESS is the best for R. On Thu, 2004-08-19 at 19:13, S Peri wrote: Hi, Is there any IDE or any editor or any pice of code for

Re: [R] Package rmutil

2004-08-19 Thread Jeff Gentry
Where I can to download rmutil package? You can use the reposTools package from Bioconductor: install.packages2(rmutil) There's also the website: http://popgen0146uns50.unimaas.nl/~jlindsey/rcode.html __ [EMAIL PROTECTED] mailing list

Re: [R] Package rmutil

2004-08-19 Thread Prof Brian Ripley
On Thu, 19 Aug 2004, Savano wrote: Where I can to download rmutil package? See the FAQ 5.1.5. PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics,

RE: [R] paired t-test vs pairwise t-test

2004-08-19 Thread Marc Schwartz
On Thu, 2004-08-19 at 14:42, Liaw, Andy wrote: From: Duncan Murdoch On Thu, 19 Aug 2004 13:42:21 -0300 (ADT), Rolf Turner [EMAIL PROTECTED] wrote : You wrote: What's the difference between t.test(x, y) and pairwise.t.test()? Is it just that the former takes two

[R] Coagulation data by Box/Hunter/Hunter

2004-08-19 Thread T. Murlidharan Nair
Hi!! Has anyone used the coagulation data for statistical analysis ? I managed to get the data from the web but unsure of the way its supposed to read. I am new to R so trying to gets myself familiarized with the statistical tools using available data. I am appending the data I got of the web.

[R] ROracle and vector elements

2004-08-19 Thread Laura Holt
Hi there! Is ROracle available for Windows, please? I found a download site, but it's really for UNIX/Linux. Here is a thought question, please: Why do the vector elements start at location 1 rather than zero, as C does? Thanks in advance! R Version 1.9.1 Windows Sincerely, Laura Holt mailto:

[R] A question about memory size

2004-08-19 Thread Stephanie Tsung
Hi I have a question about how to increase my memory size, could someone answer it for me?? I am using Bioconductor in R to calculate gene expression values with mas5, dchip, and mas4. I have only 18 samples, all from Affymetrix U133A Plus 2 arrays, which has ~54,000 genes. My machine

[R] List dimention labels to plots of components

2004-08-19 Thread White, Charles E WRAIR-Wash DC
It is frustrating to see the labels I want in the dimensions of a list but not be able to extract those labels into titles for plots generated from component objects. If someone could set me straight, I would appreciate it. For your amusement, I have provided an example of the Byzantine code I

Re: [R] A question about memory size

2004-08-19 Thread Prof Brian Ripley
Please read the rw-FAQ, as the posting guide asks. The answer is there. On Thu, 19 Aug 2004, Stephanie Tsung wrote: I have a question about how to increase my memory size, could someone answer it for me?? I am using Bioconductor in R to calculate gene expression values with mas5, dchip,

Re: [R] ROracle and vector elements

2004-08-19 Thread Peter Dalgaard
Laura Holt [EMAIL PROTECTED] writes: Hi there! Is ROracle available for Windows, please? I found a download site, but it's really for UNIX/Linux. As I understand it, a Windows version can easily (?) be built, but to do so you need Oracle libraries and header files, which you can't have

[R] How do you read in a table with numeric and dates fields without it defaulting to reading in the dates as factors?

2004-08-19 Thread Briggs, Meredith M
__ [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

Re: [R] Is R good for not-professional-statistician, un-mathematical clinical researchers?

2004-08-19 Thread Andrew Robinson
Jake, this is a great question. I have experience teaching R to forestry and natural resources graduate and undergraduate students. These students are amongst the least numerically comfortable students I've taught. I have found that the numerically comfortable among them leap into the

Re: [R] Clipping of display in Lattice graphics

2004-08-19 Thread Paul Murrell
Hi Deepayan Sarkar wrote: On Thursday 19 August 2004 06:38, Matthew Walker wrote: I'm baffled as to how the Lattice package achieves clipping. Would someone mind explaining this to me? Firstly, my attempt using just the grid package: x-seq(0,3,by=0.3)/2.8 y-seq(0,1,by=0.1) grid.newpage()

[R] column names in data.frame

2004-08-19 Thread Louize Hill
Dear R-help, Please can someone explain how to put a column name on an output data.frame. ##Starting with a data.frame with 3 columns (d$Year, d$NoIndiv, d$wtd_tl) yr_ind - split (d$NoIndiv, d$Year) yr_tl - split (d$wtd_tl, d$Year) ann_ind - sapply (yr_ind, sum) ann_tl - sapply (yr_tl, sum)

Re: [R] glmmPQL in R and S-PLUS 6 - differing results

2004-08-19 Thread Andrew Robinson
Professor Ripley, These are optimization problems with multiple local maxima, and like any complex statistical fitting problem you should not expect all programs to give the same answer. Examining the log-likelihoods from the anova statements, while fit1 and m1 are similar, the maximum

RE: [R] column names in data.frame

2004-08-19 Thread Steve Upton
Is names(d2) - c(year,av_t1) what you're looking for? steve -Original Message- From: [EMAIL PROTECTED] [mailto:r-help- [EMAIL PROTECTED] On Behalf Of Louize Hill Sent: Thursday, August 19, 2004 7:26 PM To: [EMAIL PROTECTED] Subject: [R] column names in data.frame Dear R-help,

[R] probability histogram question

2004-08-19 Thread Joseph LeBouton
Hello, all; I get an unexpected result when trying to plot a probability histogram with R1.9.1 on windows xp: #with the following code: x - runif(100,0,1) hist(x) hist(x, freq=F) h - hist(x, freq=F) summary(h) #Length Class Mode #breaks 11 -none- numeric #counts 10

Re: [R] column names in data.frame

2004-08-19 Thread Louize Hill
apparently not that gives me the following error: Error in names-.default(`*tmp*`, value = c(year, av_t1)) : names attribute [2] must be the same length as the vector [1] presumably because the 2nd column (av_tl) already has a column name, and therefore seems longer? - Original

Re: [R] probability histogram question

2004-08-19 Thread Alec Stephenson
Alec Stephenson Department of Statistics Macquarie University NSW 2109, Australia Joseph LeBouton [EMAIL PROTECTED] 08/20/04 09:56am Hello, all; I get an unexpected result when trying to plot a probability histogram with R1.9.1 on windows xp:

Re: [R] probability histogram question

2004-08-19 Thread Joseph LeBouton
Alec, Thanks for your reply. I guess what I'm getting at is that I to plot the histogram such that the HEIGHT of each bar represents the proportion of that class in the sample. From your reply I gather that the AREA of each bar is currently representing the proportion. My current work-around

RE: [R] column names in data.frame

2004-08-19 Thread Steve Upton
Hi Louize, Sorry to lead you astray, but I misinterpreted what you were doing (I could still be leading you astray :-(. You should be able to make a names attribute for d2 by names(d2) - c(year, col2namewhatyouwant, col3namewhatyouwant, ...) to the length of av_t1. When you created the data frame

Re: [R] column names in data.frame

2004-08-19 Thread Alec Stephenson
I suspect you are interpreting the column names as a column, whereas in fact you only have a single column data frame, so using names(d2) - whatever will work. Thanks, Alec Alec Stephenson Department of Statistics Macquarie University NSW 2109,

  1   2   >